The MWIR light detection has widespread applications in the fields of health monitoring, environmental protection, defense and national security, as well as space exploration and other fields. Existing technologies with high sensitivity are mainly based on semiconductor photo-detectors. In the past half-century, many semiconductor material systems have been intensively studied, and significant progress has been made   . Currently, the leading photodetector approaches are mainly based on HgCdTe (MCT)   III-V type-II superlattices, and InSb-based material systems. Among them, MCT is the premier material of interest for MWIR Focal Plane Arrays (FPA) applications. The challenges have been reducing the cost, size, weight, and power requirements, while improving the detection range and resolution. So far, these MWIR FPAs still require cryogenic cooling to achieve high detectivity, which is bulky and expensive. It is known that the major fundamental hurdle of the MWIR photodetectors at high operation temperature is the high Auger recombination rate. Among these detectors, the MCT technology has demonstrated its state-of-the-art performance, which is closely approaching the Auger theoretical limit.
It is well known that Auger coefficient in IV-VI semiconductors    is about an order of magnitude lower than those in Sb-based type-II QWs,    which are in turn significantly suppressed relative to other III-V and II-VI semiconductors such as MCT   for the same waveleng. Such low Auger recombination should result in superior device performance such as high detectivity for detectors at a high operating temperature. High specific detectivity (D*) of PbSe MWIR photoconductive (PC) detector at room temperature was reported, 2.8 × 1010 cm・Hz1/2/W and 4.2 × 1010 cm・Hz1/2/W at ~3.8 µm, with and without antireflective coating    . Northrop Grumman has demonstrated monolithic integrated PbSe PC detector FPA with Si read-out circuitry (ROIC) with thermoelectric (TE) cooling  . These exciting results have stimulated renewed interests in using Pb-salt cameras for TE-cooled or uncooled applications  .
It is worth nothing that photovoltaic (PV) detectors offer advantages over their PC detector counterparts, such as low power consumption, high pixel density, and small pixel size. In addition, a PV detector FPA operating at low bias with no or low flicker noise eases the design of ROIC which further reduces cost. Although the D* of the IV-VI Pb-salt semiconductor junction at a low temperature has been investigated via experimental and theoretical research, the performance limit of PbSe at high temperature has not been well investigated. In this paper, we investigate the performance limit of a PbSe n+-p junction detector operating at room temperature and TE-cooled temperature.
2. The Specific Detectivity D*
Specific detectivity, or D*, for a photodetector is an important figure of merit and is determined by  :
where A is the active area of the detector and NEP is the noise equivalent power (the power needed to generate the signal current which is equal to the noise current)
Here  is the radiant power incident on a detector.
In a junction device, the fluctuation of diffusion rates in the neutral region and the generation-recombination (g-r) fluctuation in both depletion region and quasi-neutral region are indistinguishable. All of them give rise to the shot noise which has the form of 
where is the diode current, is the dark reverse-bias saturation current of the device and is the electrical bandwidth.
β is the ideality factor which determines how the actual diode deviates from the ideal diode. In most of our later calculations, we assume the ideal diode case where .
For the diode under the illumination of photon density flux , which generates the photon current , the total noise contribution is:
where R0 is the zero-voltage resistance which is defined by
Under the zero-voltage case, the expression of shot noise in (3) has the same form as the Johnson-noise associated with the R0 of the junction, thus it is often called Johnson noise.
Substituting (1) and (2) to (5), the expression for the detectivity of the device is
In Equation (7), η is the quantum efficiency, λ is the wavelength of the incident radiation, q is the charge of the carrier, h is Planck’s constant, c is the speed of light, k is Boltzmann’s constant, R0 is the diode incremental resistance at 0 V, A is the detector sensitive area, T is the temperature of the detector in Kelvin, and is the photon flux incident on the detector from its surroundings. It is clear that both R0A and η are key factors on the specific detectivity D*. Therefore, the resistance-area (R0A) product at zero bias needs to be calculated to determine D* based on Equation (6).
3. R0A Analysis and Discussions
The total dark current density flowing through the junction can be found with the following equation:
In Equation (7), JDF is the diffusion current density, JGR is generation-recombination current density that is often dominated in the depletion layer, JT is tunneling current density and JL is leakage current density. JL may be due to bulk as well as surface defects of the material. When suitable diode technology and construction are used, the contribution of the JL component is negligible.
The maximum R0A is obtained from the one-sided abrupt junction  , giving the limitation on the specific detectivity D*. The minority carriers on the lightly doped side of the n+-p junction will dominant the diffusion current. The R0A product determined by the diffusion current in the case of radiative recombination and Auger recombination is
where ni is the intrinsic carrier concentration, NA(ND) is the concentration of donors (acceptors), τh(τe) is the minority carrier lifetime of holes (electrons) and μh(μe) is the minority carrier mobility.
The resistance-area product (R0A) of PbSe photo-voltaic detectors under conditions of zero bias is calculated for one-side abrupt junction. Highly doped n-side homojunction (n+-p) is discussed in this paper. At room temperature, the effective R0A product on the PbSe n+-p junction can be expressed by:
where is the R0A product due to radiative recombination, is the R0A product due to Auger recombination, and is the R0A product due to Shockley-Read-Hall recombination. Different R0A products due to each recombination mechanism will be discussed in the following section.
3.1. Radiative Recombination Contribution
The R0A product determined by the diffusion current in the case where radiative recombination has the relationship on the lifetime . The radiative lifetime for PbSe follows as  :
The radiative lifetime is inversely proportional to the total concentration of free carriers under all circumstances. The radiative lifetime decreases when the doping concentration increases.
GR is the capture probability for radiative recombination. The band gap of a semiconductor Eg is expressed in electron-volts. nr is the index of refraction, and is the density-of-states electron and hole effective masses. ni is intrinsic carrier concentration. n0 and p0 represent the equilibrium carrier concentrations. According to Equation (14), the (R0A)RA product can be found after the radiative recombination lifetime is calculated.
Figure 1 describes the radiative life time on the dopant concentrations for the PbSe junction at 300 K. When carrier concentration is between 1017 and 1018, the radiative recombination life time of PbSe is between 10−7 s and 10−8 s at 300 K.
3.2. Auger Recombination Contribution
The R0A product determined by the diffusion current in the case of Auger recombination has the relationship on the lifetime . The PbSe Auger’s recombination coefficient CA is given by  :
where and are the longitudinal and transverse effective mass, is high frequency dielectric constant, and is vacuum permittivity. During the Auger recombination process, the carrier lifetime is defined as:
Figure 2 shows the Auger life time of different dopant concentrations for the PbSe junction at room temperature. When carrier concentration is between 1017
Figure 1. The dependence of the radiative life time on the dopant concentrations for the PbSe junction at 300 K.
Figure 2. The dependence of the Auger life time on the dopant concentrations for the PbSe junction at 300 K.
and 1018, the Auger recombination life time of PbSe is between 10−7 s and 10−8 s at 300 K. According to Equation (17), the (R0A)Auger product can be derived when the Auger recombination lifetime is calculated.
3.3. Tunneling Current Contribution
The (R0A)T determined by tunneling is given by
Tunneling simulation of PbSe depends on the crystal orientation due to the difference in effective masses. The four effective masses to consider are: (conduction band, longitudinal), (conduction band, transverse), (transverse band, longitudinal) and (transverse band, transverse). is the electron resting mass  .
For PbSe  directions, the effective mass on different orientations can be determined by the Table 1.
In order to compare Auger, radiative and tunneling recombination, R0A production on those three recombination mechanisms versus different temperatures is analyzed. In the following calculation , a typical PbSe hole concentration is used. To satisfy n+-p junction, ND is chosen to be 1 × 1018 cm−3, which is five times higher than NA. Throughout this paper, we use ND ≥ 5NA for our n+-p junction calculation. Looking at Figure 3(a) for the given carrier concentration , at a higher temperature, (R0A)T increases significantly for the PbSe  junction so that the tunneling (R0A)T can be ignored. This is mainly because the energy bang gap of PbSe increases with temperature. When
Table 1. PbSe  effective mass in different directions.
Figure 3. The dependence of the tunneling R0A product on the dopant concentrations for the PbSe  junction at 160 K and 240 K. For Figure 3(a) the n side carrier concentration is 1 × 1018 cm−3, and in Figure 3(b) n side carrier concentration is 5 × 1018 cm−3.
the temperature is around 160 K, (R0A)T is comparable to the Auger R0A and radiative R0A. For even lower temperatures, (R0A)T becomes dominant. At temperatures above 180K K (R0A)T is at least five times higher than the other mechanisms and thus (R0A)T can be neglected relative to the total R0A product.
To keep the same NA and increase ND doping concentrations, Figure 3(b) shows that (R0A)T cannot be neglected even at room temperature. At about 275 K, (R0A)T is five time lower than the other mechanisms and at higher than 275 K tunneling needs to be considered in the total R0A product. Comparing Figure 3(a) and Figure 3(b), it is apparent that carrier concentration has a large influence on R0A, and therefore, the total device performance.
3.4. SRH Recombination Contribution
The carrier generation-recombination mechanisms in detector devices are distinguished as Auger, radiative, tunneling, and Shockley-Read-Hall’s (SRH) generation-recombination mechanisms. Auger, radiative, and tunneling mechanisms are determined by energy band structures. However, SRH’s mechanism is determined by the material quality. In Equation (4), the depletion generation and recombination R0A term is largely dominated by the lifetime associated with SRH centers. However, the value of this lifetime is extremely uncertain as has been noted in previous studies   . These studies often assume longer SRH lifetimes in their investigations, typically on the order of 10−8s25. Thus, using this assumed SRH lifetime, higher temperature operations produce high values for (R0A)GR which may be neglected in our case. However, SRH recombination may need to be considered with a shorter SRH lifetime.
3.5. Discussion of the Overall R0A Limitation
Utilizing the R0A expressions shown in the above sections, calculations of R0A for Auger, radiative, and tunneling were described in Figure 4. For a fixed operating temperature, R0A dependence on carrier concentration reveals the contribution of each mechanism, revealing an optimized condition for the maximum R0A limit. Three temperatures of interest are used. 300 K, 240 K, and 160 K represent temperatures of an uncooled, cut-off wavelength of 5 μm, and that the tunneling effect needs to be considered.
Table 2 shows simulation parameters at 300 K. Between 50 K and 300 K, the mobility varies as 
And the expression of energy gap is given as 
At 300 K, the total effective R0A product is around 1.3 × 10−4 ohm・m2. At 240 K, the total effective R0A product is around 1.6 × 10−4 ohm・m2. When temperature is higher, the total effective R0A product becomes smaller. According to Equation (10), Figure 4 gives the results of the total effective R0A product on the
Figure 4. The PbSe R0A product determined by the diffusion current in the case of Auger recombination and radiative recombination at different high temperature. The black curve is the total effective PbSe R0A product on the diffusion current in the case of Auger recombination and radiative recombination.
Table 2. The parameters are applied in the simulations for 300 K  .
diffusion current in the case of Auger recombination and radiative recombination at a different temperature. At 160 K, tunneling R0A value is comparable to the Auger and radiative R0A value.
Figure 4(a) shows the R0A limitation at 160 K when the n+ doping concentration is . In this case, is comparable with and . When the tunneling effect on the total R0A cannot be ignored, there exists an optimized carrier concentration for the largest R0A. In this case, the largest R0A occurs at . When the n side carrier concentration was increased from to , tunneling could not be ignored at 300 K, as shown in Figure 4(b). The contribution from tunneling recombination results in a lower total R0A in Figure 4(b). In Figure 4(c) and Figure 4(d), where tunneling can be ignored, it seems that higher PbSe doping concentration (NA) gives a larger total R0A. This is because a larger NA leads to a smaller minority carrier concentration, which in turn reduces the saturation current. However, for a n+-p junction n-doping concentration has to be increased, as p-doping concentration increases, which will increase the tunneling effect. In our simulations, the maximum R0A limit at 240 K (Figure 2(c)) and room temperature (Figure 2(d)) occurs when carrier concentration on the p side is at a maximum, while satisfying the condition for the n+-p junction where and . The recombination mechanisms in n+ side are not taken into consideration in this paper. Recombination in n+ side poses a further limitation of n+ concentration and thus p-doping concentration for optimized D*. Assuming a wide bandgap N+-doping layer with negligible Auger and Radiative recombination rate, and negligible tunneling effect, then higher p-type doping concentration could be used to further increase the R0A.
4. Thermal Noise and Background Photon Shot Noise
As can be seen from Equation (6), the two main mechanisms contributing to the homojunction device are the thermal noise and the background photon shot noise. By using calculated R0A when and , the temperature dependence current noise of these two sources is described in Figure 5.
From Figure 5, we can easily see that at room temperature the dominant source of noise in the device is the Johnson noise. For temperatures around 191 K to 244 K, the contribution of both noise sources is in the same order of magnitude. Below 200 K, photon noise becomes the dominating mechanism, with the background photon noise current five times greater than the Johnson noise. Therefore, if we ignore one noise source when it is five times lower than the other source, the detectivity expressed in (7) can be simplified to high temperature Johnson noise limited D* and background limited infrared photodetector (BLIP) D*.
Figure 5. The comparison between the thermal noise and photon background noise on PbSe PV detectors related to temperature.
5. Theoretical Limitation on the Specific Detectivity D* at High Temperature
To find the theoretical limit for PbSe photovoltaic detector D* at 300 K, we apply the R0A product from the above conditions with Equation (23) because photon shot noise can be negligible at room temperature.
At 240 K, to calculate the theoretical limit for PbSe photovoltaic detector D* following by:
Theoretically, the quantum efficiency should be calculated by the number of electron-hole pairs generated per incident photon. In order to get the D* limitation, the detector quantum efficiency is simplified to a constant for a different cutoff wavelength. The D* of the PbSe junction is simulated for a quantum efficiency of 75% and 100%, λ is the wavelength of the incident radiation, q is the charge of the carrier, h is Planck’s constant, c is the speed of light, k is Boltzmann’s constant, R0 is the diode incremental resistance at 0 V, A is the detector sensitive area, T is the temperature of the detector in Kelvin. With cutoff wavelength λ = 5 μm at T = 240 K, the photon flux equals to 1.23 × 1016 cm−2 sec−1. Figure 6 shows D* for the PbSe junction as it depends on the absorbing light wavelength at 240 K and 300 K. At 300 K with 100% quantum efficiency, the theoretical limit of D* for the photovoltaic PbSe n+-p junction detector is 2.8 × 1010 HZ1/2/W at 4.5 μm. At 240 K, the D* theoretical limit is 3.7 × 1010 HZ1/2/W at the cut-off wavelength of 5 μm.
Due to reabsorption in the crystal, the effect of radiative recombination may be reduced and carrier radiative life time becomes larger  . Photon reabsorption is an effect by which a photon created by the recombination of an electron-hole pair is reabsorbed to produce a new pair and the process is repeated over and over again until a photon can either escape from the sample, or is absorbed by a process which does not create electron-hole pairs, such as free carrier absorption  . Considering this, Auger recombination becomes the dominating mechanism at high temperature. The D* of the PbSe junction in this case is shown in Figure 7.
At 300 K, the limitation of D* now increases from 2.8 × 1010 HZ1/2/W to 3.9 × 1010 HZ1/2/W at 4.5 μm when radiative recombination is ignored.
Figure 6. The specific detectivity D* on the PbSe junction versus the absorbing light wavelength at 240 K and 300 K.
In this n+-p homojunction junction model, the n+ PbSe was not considered for its impact on diffusion R0A product. As discussed in previous sections, high carrier concentration will lead to high Auger and radiative recombination and thus reduced R0A product. To achieve the calculated D* in this paper, a heterojunction structure with a wide bandgap N+ doped layer with much reduced Auger and radiative recombination is needed.
6. Optimization on Carrier Concentration
According to Figure 4(a), when tunneling can be ignored, R0A value increases with increasing NA. However, to keep the device n+-p junction ND needs to be increased with NA, which increases the tunneling effect. Therefore, there exists an optimized NA which provides the maximum R0A at a given temperature. In the following calculation, ND/NA is set to be 5 to satisfy the condition for a n+-p junction. NA is then optimized at 300 K and 240 K, as show in Figure 8. The optimized R0A for PbSe homojunction at 300 K and 240 K are 2.5 × 1018 cm−3 and 1.5 × 1018, respectively. For optimized carrier concentration the detector reaches BLIP limited D* at 210 K, 10 degrees higher than the “typical” carrier concentration.
At 240 K, the optimized and are used to calculate maximum R0A product in Figure 9(a). By using the maximum R0A value, the limitation D* is increased to 5.6 × 1010 HZ1/2/W. At 300 K, the optimized and are applied to get the limitation D* which can be increased to 3 × 1010 HZ1/2/W in Figure 9(b). Therefore, at different temperatures the carrier concentration can be optimized to achieve the highest D*.
Figure 7. The specific detectivity D* on the PbSe junction with Auger limitation versus the absorbing light wavelength at 300 K.
Figure 8. The n side doping concentration dependent PbSe R0A products at 240 K and 300 K.
Figure 9. D* versus wavelength at 240 K and 300 K after carrier concentration optimization.
Performance limitation of the PbSe homojunction at high temperatures is theoretically studied in this paper. For such PV detectors thermal noise is dominating at temperatures higher than 240 K. BLIP limit could be achieved at temperatures around 210 K. The calculated peak D* are 3 × 1010 HZ1/2/W and 5.6 × 1010 HZ1/2/W at 300 K and 240 K respectively. Achieving D* of more than 1010 HZ1/2/W in a MWIR detector at uncooled temperatures allows for its use in practical applications of high sensitivity.
The authors thank helpful discussions with Dr. Sung-Shik Yoo from Northrop Grumman. This work is partially supported by the DARPA WIRED program through Northrop Grumman, US army research office (ARO) under Grant No. W911NF-14-1-0312. The views, opinions and/or findings expressed are those of the author and should not be interpreted as representing the official views or policies of the Department of Defense or the U.S. Government.