Suppression of compensating native defect formation during semiconductor processing via excess carriers

In many semiconductors, compensating defects set doping limits, decrease carrier mobility, and reduce minority carrier lifetime thus limiting their utility in devices. Native defects are often responsible. Suppressing the concentrations of compensating defects during processing close to thermal equilibrium is difficult because formation enthalpies are lowered as the Fermi level moves towards the majority band edge. Excess carriers, introduced for example by photogeneration, modify the formation enthalpy of semiconductor defects and thus can be harnessed during crystal growth or annealing to suppress defect populations. Herein we develop a rigorous and general model for defect formation in the presence of steady-state excess carrier concentrations by combining the standard quasi-chemical formalism with a detailed-balance description that is applicable for any defect state in the bandgap. Considering the quasi-Fermi levels as chemical potentials, we demonstrate that increasing the minority carrier concentration increases the formation enthalpy for typical compensating centers, thus suppressing their formation. This effect is illustrated for the specific example of GaSb. While our treatment is generalized for excess carrier injection or generation in semiconductors by any means, we provide a set of guidelines for applying the concept in photoassisted physical vapor deposition.

underpinning the improvements were not well developed. Understanding the relationships and interactions between defects, photogenerated carriers and the Fermi level is critical for assessing the potential for defect suppression via light stimulated processing techniques and can ultimately direct their implementation.
Herein, we develop a complete framework for determining the interplay of steady-state excess carrier generation, trapping and recombination, and native defect formation to allow prediction of native defect populations during semiconductor processing. Our theoretical framework self-consistently determines the free carrier concentrations, quasi-Fermi levels, recombination rates, defect occupations, and defect concentrations. In order to illustrate its utility, we apply this model to the concrete example of compensating native acceptor suppression during light-assisted epitaxy or annealing of nominally un-doped and intentionally n-type doped GaSb. This material system exhibits dominant native defects of either positive or negative charge depending on the prevailing chemical potentials during processing and is therefore ideal for assessing the influence of each. From these results, we have constructed a set of guidelines that can be applied to the light-assisted processing of semiconductors. Our findings have implications for physical vapor deposition techniques such as evaporation, hydride vapor phase epitaxy (HVPE) and molecular beam epitaxy (MBE) and also especially for understanding differences between furnace and lamp-heated annealing treatments such as rapid thermal processing (RTP).

Results
Theoretical Framework. At thermal and chemical equilibrium, the concentrations of point defects within a material phase describe a unique thermodynamic state and are found in general by minimizing the appropriate free energy. It is important to remember that this thermodynamic equilibrium is microscopically a state of dynamic equilibrium. Individual point defects are created and annihilated at certain rates resulting in steady-state concentrations with negligible fluctuations for large system sizes. When the dominant entropy is configurational in the dilute limit, the system can be treated within the quasi-chemical formalism 6,7 . The concentrations of N species of point defects can be computed if the enthalpies of formation or quasi-chemical reaction constants are measured or calculated. A system of N + 2 equations, consisting of N quasi-chemical defect creation reactions, the electron/hole pair generation reaction (or law of mass-action = np n i 2 ), and the charge balance equation can be solved at a given temperature by self-consistently determining the Fermi level. In particular, the dependence of the defect formation energy on the Fermi level is introduced through energy transfer between the charge transition level of the defect and the appropriate reservoir of hole or electron carriers (i.e. the valence or conduction band). The chemical potential that characterizes the reservoir is simply the Fermi level when the bands of the semiconductor are in thermal equilibrium.
Under conditions of steady-state photocarrier generation, we must further consider the variation of defects' formation energies when the electron and hole QFLs split. Since two quasi-Fermi levels instead of one Fermi level must be solved for, the problem is cast as a system of N + 3 equations. The extra needed equation comes from the requirement of steady state for either electrons or holes: requirement of steady state for either electrons or holes: dn/dt = dp/dt = G - Σ(U) = 0. It is well established (e.g. in analysis of open circuit voltage of solar cells) that the QFLs represent the electrochemical potentials of the electron and hole populations in the conduction and valence bands of the crystal once those populations have thermalized amongst themselves and with the lattice but prior to recombination 8 Defects may exist in multiple charge states and may generally capture both electrons and holes via the electron-phonon interaction, as described by their capture cross-sections 9 The occupancy of defect charge states depends on both the free electron and hole concentrations. Consequently, the formation free energies for all defects in the system are a function of both QFLs, weighted by the strength or frequency of interactions with the valence and conduction bands.
Because charged states arise through real forward and reverse exchange of carriers between the localized defect and the conduction and valence bands, the most accurate method of the weighting the dependence of the formation energy on the QFLs that govern this exchange is as a ratio of the rates of these transactions. To illustrate this point, consider the formation of a defect of charge state q = -1. This state can be produced by either the capture of an electron or the emission of a hole. The electron and hole QFLs would then be weighted based on these two exchange rates, respectively.
A consequence of this behavior is that the defect formation energy will be tied most strongly to the QFL that governs the dominant exchange rate. For instance, negatively charged defects are typically more likely to trap holes, and thus their concentration will be affected more by the energy of the hole QFL than the electron QFL. However, it is important to bear in mind that the defect's interaction with both bands and corresponding dependence on both QFLs is in stark contrast to a shallow dopant, which exchanges carriers exclusively with only one band. The full details of the theoretical approach are included in the Computational Methods section.

Defect Calculations.
To understand the potential effects of illumination on defect concentrations during light-assisted semiconductor processing, such as epitaxy or rapid thermal annealing, we have applied this framework to the model system of n-type GaSb. GaSb is an interesting but simple system in which to study the effects of light-based processing because of its propensity to form a single species of compensating defect over the full range of chemical potentials. Furthermore, excessive native defect concentrations and apparent doping limits impact its performance in practical applications.
Epitaxial growth is typically carried out in the temperature range 400-550 °C. Under Ga-rich growth conditions, negatively charged Ga Sb antisites form with relatively low energy (< 1.8 eV) and are the single most dominant defect in intentionally n-type doped single crystals [10][11][12] Positively charged Sb Ga antisites and Ga i interstitials also form at non-negligible concentrations in n-type GaSb grown under Sb-rich conditions 11 . Overall, the formation energies of negatively charged Ga Sb antisites are higher in Sb-rich vs Ga-rich conditions, suggesting that maintaining an Sb-rich environment during GaSb crystal growth or annealing is advantageous for controlling compensating native defect concentrations. In practice, achieving Sb-rich growth conditions by MBE can be somewhat difficult. Since Sb has a relatively low vapor pressure, the thermally evaporated molecules must be thermally cracked from its initial Sb 4 state, and deposited atoms tend to clump on the growth surface, leaving exposed Ga sites. Typical growths are carried out at modest V/III ratios (~2-6) close to those needed to achieve a Sb-stabilized surface. This realistically places the growth in an intermediate regime 13 Unintentionally doped bulk crystals and epilayers are usually p-type due to high levels of Ga Sb antisites. The presence of charged defect species also degrade the carrier mobility. GaSb has a relatively low bandgap (0.72 eV at room temperature), and the theoretical electron mobility can be as high as 44,000 cm -1 V s 14 However, much lower mobilities are typically observed in part because of scattering by ionized centers 15 Suppressing defect formation is therefore critical to realizing the theoretical mobilities of GaSb and enabling its widespread use.
Under all conditions, native defect concentrations rise exponentially with growth temperature, and the total concentrations decrease exponentially with defect formation free energy. The magnitude of any light-induced changes to the defect formation free energy will be bounded (approximately, in degenerate injection conditions) by the bandgap, similar to limits of open circuit voltage in photovoltaics. Likewise, the magnitude of any atomic chemical potential-induced changes to the defect formation free energy are limited by the formation free energy of the compound, which is typically a few eV per formula unit and in general correlates with bandgap. Variation of the photogeneration rate and the chemical potentials are therefore theoretically expected to lead to similar changes in the native defect concentrations. Practically, though, the effects of illumination on native defect populations may be more muted because of the large photon or other particle fluxes required to maintain large QFL splitting, especially in the presence of recombination-active centers and at high temperatures.
To evaluate the relative influence of photogenerated carriers and chemical potentials on the native defect concentrations, calculations were performed as a function of temperature for GaSb crystals under both Ga-rich and Sb-rich environments. The trap species, their ionization levels, and their formation enthalpies used in these calculations are listed in Table 1 11 . Two donor doping concentrations were considered: undoped and moderately doped (N d = 10 17 cm -3 ). Additionally, three photon fluxes were applied: 0, 10 23 and 5 × 10 24 cm -3 s -1 . Assuming a carrier lifetime on the order of 1 ns, these rates are expected to produce excess carrier concentrations in the range of 10 14 to 5 × 10 15 cm -3 , which is well below the extrinsic n-type doping concentration and are achievable with common lamp and laser light sources.    Figure 1 shows the calculated QFLs as functions of temperature and excess carrier generation rates for both doping cases. As expected, QFL splitting increases with increasing photocarrier generation rate. The splitting also decreases with increasing growth temperature as the intrinsic carrier concentrations begin to dominate the photo-injected electron and hole concentrations and recombination becomes faster.
The resulting general trends in the concentrations of typical native defects in the presence of excess carriers are determined primarily by a combination of their formation energies, their charge states and the degree of QFL splitting. Our first main finding is that excess carriers always reduce the concentration of charged native defects of both polarities (regardless of extrinsic doping) and that the degree of suppression increases with QFL splitting (Fig. 1). The behaviors of − Ga sb 2 and + Sb Ga 2 antisites exemplify the main trends observed in all of the defects listed in Table 1 and are displayed in Fig. 2. Under Ga-rich growth conditions, Ga Sb antisites are the dominant native defect. Their concentration rises with temperature in undoped GaSb grown under dark conditions as the electron QFL moves closer to the conduction band edge (see Fig. 1) and eventually exceeds 10 16 cm -3 . These defect concentrations calculated for dark growth conditions are comparable with the excess free hole concentrations measured in nominally un-doped GaSb crystals, which was linked to the presence of native acceptor defects 10,16 The higher free electron concentration in the doped case pushes the QFLs closer to the conduction band edge, driving up the − Ga sb 2 concentration at all temperatures to approach that of the donor concentration. Compared to these baseline equilibrium defect concentrations in the dark (generation rate = 0 cm -3 s -1 ), the suppression under steady-state excitation is most significant at the lowest temperatures due to the larger QFL splitting. In fact, the suppression mechanism is so dominant that the − Ga sb 2 concentrations are approximately equal between the undoped and doped cases at the highest generation rate. Under Sb-rich growth conditions, the formation energies of Sb Ga antisites are lower than Ga Sb antisites, making them the dominant defect.
It is interesting and instructive to compare the behavior of defects with different charge types. For the n-type extrinsic doping case considered here, the presence of photogenerated carriers induce a strong reduction in the compensating − Ga sb 2 antisite concentration. The suppression, which is predicted to be approximately 4 orders of magnitude at 300 °C, is most significant in doped GaSb given the higher − Ga sb 2 concentration induced by the high Fermi-level energy under dark growth conditions. The suppresion in the undoped case of approximately 3 orders of magnitude is also quite substantial and is only tempered by the lower concentration found under dark conditions. Thus in all cases, photogenerated carriers can be used to substantially reduce the formation of the otherwise compensating native defects. + Sb Ga 2 antisites, which have the same charge state as the ionized extrinsic donors, exhibit a more muted change in their concentration in the presence of photogenerated carriers. Additionally, the degree of reduction decreases rather than increases as the extrinsic shallow donor concentration is increased. The disparity in the response to the presence of photogenerated carriers between these two defects is largely due to their different dependences on the electron and hole QFLs. The − Ga sb 2 antisite formation energy is primarily sensitive to the hole QFL, which is significantly affected by the presence of a non-equilibrium concentration of holes at both donor doping levels. On the other hand, the formation energy of + Sb Ga 2 antisites is strongly tied to the electron QFL. Its energy is altered by the photocarrier generation in undoped GaSb due to the ability to achieve conditions of high injection with relatively low light intensity. In contrast, achieving high injection conditions requires much greater generation rates in more heavily-doped semiconductors. Photocarrier-induced suppression of + Sb Ga 2 antisites is therefore minimal for the generation rates considered here. Inverse conclusions can be drawn for p-type doped material.
Interestingly, there is a divergence in the QFLs between the Ga-rich and Sb-rich conditions at elevated temperatures. This divergence arises from a difference in the dominant native defect type and concentration, which because of their different charge states and cross sections alter the free electron and hole concentrations. The much higher concentrations of − Ga sb 2 and − Ga sb 1 defects (~10 15 -10 16 cm -3 ) generated under Ga-rich conditions effectively act as additional acceptors and push the QFLs toward the valence band. This result is consistent with the p-type behavior observed in undoped and lightly doped GaSb crystals. The degree of QFL splitting at the same doping level and temperature is also smaller in an Ga-rich vs Sb-rich environment. We attribute this to excessive carrier trapping at the higher concentration of defects generated under those conditions, which effectively regulates the excess carrier population. This result highlights the need to properly account for carrier trapping at defects when calculating the QFL splitting.
In the preceding sections, we have discussed the effects of excess carriers on native defect concentrations and thus dopant compensation and recombination rates. Additional improvements to optoelectronic performance are also realized in light-assisted processing by the higher mobilities made possible by lower compensation and ionized center concentrations. To demonstrate these effects, we calculated the total concentration of all ionized centers (native defects plus donor and acceptor impurities) for undoped and n-type doped GaSb processed under Ga-rich conditions and generation rates of 0 and 5 × 10 24 cm -2 s -1 . The calculation assumes perfect quenching of the total native defect populations from the growth or processing temperature, and the occupations of the charge states of each defect were recalculated at room temperatureFig. 3a 17 displays these concentrations as a function of growth temperature. The native defect concentration entirely dictates the number of ionized centers in undoped GaSb, and thus the suppression of defect formation by excess carriers leads to a corresponding decrease in the center concentration. A much smaller relative change is observed for n-type doped GaSb due to the additional presence of ionized extrinsic dopants. The extrinsic dopant concentration is of course independent of the presence of photogenerated carriers, and the doping concentration considered here (10 17 cm -3 ) is high enough to overwhelm the defect concentration. For ionized impurity scattering, the mobility scales inversely with the concentration of each defect multiplied by the square of their charge state. Figure 3b displays the expected mobility ratios resulting from a reduction in the defect concentration under a photocarrier generation rate of 5 × 10 24 cm -3 s -1 (compared to the dark case). The upper limit of mobility enhancement is set by the undoped case; for this generation rate and the parameters in Table 1 the ionized impurity mechanism can be suppressed by a factor of ~250 for a growth temperature around 300 °C. This ratio does not include any ionized centers due to shallow dopants, which will ultimately be present at some non-zero concentration in real samples. With increased doping, the mobility enhancement effect decreases, and the processing temperature at which the maximum occurs increases.

Discussion
When contemplating the implications and potential for processing semiconductors under a light bias, one should first consider that these calculations apply only to the case of thermodynamic equilibrium. Realization of local equilibrium conditions depends on several factors, including the processing temperature, diffusion coefficients for atoms at the crystal surface and in the bulk, and chemical potentials. Each material will behave differently. However, local equilibrium can be achieved for conventional growth techniques such as MBE and HVPE, as well as for thermal annealing, with the correct selection of temperatures and time scales (i.e. growth rate, soak time, etc.).
Defect reduction and a general improvement in crystalline quality have been experimentally observed in variety of semiconductors (i.e. CdTe, ZnSe, ZnS, GaN) grown under light stimulation by both molecular beam epitaxy and metal organic chemical vapor deposition 3,5,[18][19][20][21][22] The influence of excess photogenerated carrier populations on defect formation was first explored separately by Marfaing 23 and Ichimura, et al. 24 through simplified microscopic models of specific growth processes. More recently, a theoretical model proposed by Bryan, et al., explicitly considered the role of QFL splitting on defect populations 5 The concepts therein represent an advance in defining the intriguing mechanism coupling excess carrier populations to defect concentrations treating QFLs as chemical potentials. However, a complete, rigorous framework was not developed, nor were other mechanisms eliminated as possible explanations for their experimental observations. Our model offers a comprehensive and rigorous description of this specific mechanism of suppressing charged defect formation by altering the electron and hole chemical potentials via injection or generation of excess carriers, which couple to the defect formation enthalpy through capture and emission rates. It defines the conjugate variables to the QFLs (which are the chemical potentials for the free carriers in the bands) as the normalized capture and emission rates between the defect's energy level and appropriate bands (see Eq. 11a,b). This allows us to write the thermodynamic potential Δ G = Δ H -TΔ S (with Δ H from Eq. 3) which must be minimized in order to find the steady-state concentrations of defects in the light. The quasichemical formalism utilized herein is equivalent to assuming the dilute limit where Δ S is dominated by configurational entropy 25 .
It is important to emphasize that light or other methods of introducing excess carriers, such as radiation or electron beams 26 , may modify doping and compensation via many mechanisms. These include exciting or decomposing precursor molecules, modifying desorption rates 27 , reducing surface band bending 28 , or creating or annihilating defect complexes. We note the relevance of related bodies of work on photo-induced defects in semiconductors and insulators, for example, color-center formation in II-VI crystals, the Staebler-Wronski effect in amorphous Si, and radiation-enhanced diffusion in GaAs [29][30][31][32][33][34][35] Carefully-controlled and designed experiments measurements are required to differentiate between the various physical mechanisms. Focused experimental measurements should be made to determine whether enhancements in crystal quality and doping efficiency can be definitively assigned to the mechanism proposed here. Experiments comparing growth to annealing may be particularly helpful in differentiating bulk thermodynamic mechanisms such as the present one from kinetic ones of the growth interface.
The observations presented above naturally lead to a set of common guidelines for controlling native defect concentrations using light-assisted epitaxy or other processing methods during which excess carriers are generated. (1) A greater reduction in the concentration of native defects can be achieved with excess carriers at lower temperatures where the photogenerated carrier population is large enough to support substantial QFL splitting. It is subsequently suppressed when the intrinsic carrier concentration approaches the excess carrier density. From a thermodynamic standpoint, lower temperatures can be beneficial for controlling native defect populations and should be used when it is possible to also suppress defect formation due to kinetic limitations. (2) Increased photogeneration rates result in a greater reduction in native defect concentrations. Practically, there will be an upper limit to the light intensity and density of photogenerated carriers that can be introduced without initiating other processes including: melting, desorption, and crystal damage. There will therefore be an optimal light intensity for each material based on its melting temperature, extrinsic doping level and the vapor pressure of its constituent atoms.
(3) The largest changes in concentration are realized for defects that efficiently exchange carriers primarily with the minority carrier reservoir. In typical cases, these will be compensating native defects having shallower charge transition levels and the opposite charge polarity to the extrinsic dopant (and thus larger capture cross section for the minority carrier). In fact, techniques that target the injection of minority carriers (i.e. electrical injection) could also be very effective at reducing compensating defects. Very small effects will be expected in heavily-doped cases for defects of the same charge type as the dominant dopant. (4) In intrinsic or unintentionally-doped cases, the formation of defects of both charge polarities will be more strongly suppressed since both QFLs will change significantly with the introduction of excess carriers. This is also true for recombination at defects interacting equally with both carrier types. (5) The magnitude of compensating defect reduction using this technique will be greater for higher doping concentrations. This is because of the greater driving force for compensating defect formation under dark conditions in the first place. Ultimately, the defect concentration should reach the same minimum for all doping levels at high enough generation rates. (6) From a broader perspective, light-based processing techniques will exhibit the most dramatic effects in semiconductors having compensating defects with shallower charge transition levels, smaller formation enthalpies for the neutral defect to maximize the impact on the formation free energy for charged defects, and large bandgaps, which allow greater modulation of the QFLs. These observations help to explain why light-stimulated epitaxy has produced the most dramatic changes in the free carrier concentrations of doped II-VI and wide bandgap III-V semiconductors. An obvious benefit is that light-assisted processing techniques can be easily implemented. Our model now provides a way to evaluate the outcome of specific processing conditions and optimize them to achieve the greatest reduction in defect concentrations.
While the example of GaSb has been used herein, we reiterate that the model generally treats the behavior of the charge states and capture and emission behavior of the defects, which we assumed to be governed primarily by the charge states. Thus the findings are generally applicable to any semiconductor and any semiconductor deposition or processing technique. This includes the effects of photogenerated carrier populations on compensating defects formed by extrinsic dopants themselves. The model may also explain differences observed in defect concentrations resulting from thermal equilibrium furnace annealing and light-heated techniques such as RTP.

Computational Method
Here we develop the model framework for the simplified case of one extrinsic shallow dopant and one dominant compensating acceptor-like native defect; extension to multiple species is straightforward. We consider a binary compound semiconductor with formula unit AB doped n-type with a fixed density N D of extrinsic shallow donors and compensated by a fixed density of shallow extrinsic acceptors N A . We assume a thermal equilibrium concentration of a single dominant native defect on the B sublattice (X B ) having a dominant charge state, q: where g is the number of equivalent configurations, N B is the sublattice site density, and ∆H f a is the defect formation energy. ∆H f a can be calculated as 36 : where E VB is the valence band maximumand E D is the total energy of the supercell comprised of n A and n B atoms and a single charged defect. The chemical potentials of A and B are μ A and μ B , and μ F is the Fermi level, which represents the chemical potential for both free electrons and holes in thermal equilibrium. The second term in Eq. 2 accounts for the exchange of carriers between the defect and host semiconductor bands and strictly is written for 0 K, where the Fermi-Dirac distribution is a step function. The atomic chemical potentials are restricted to values less than their bulk values, μ A(bulk) and μ B(bulk) and the chemical potential of the AB compound is assumed to be the sum of μ A and μ B such that Eq. 2 can be re-written as: The defect's charge transition energies (denoted q n /q m , e.g. 0/+ for a donor) can be determined by equating the formation enthalpies for the charge states ∆ = ∆ H H f n f m . Under steady-state excess carrier generation (e.g. photon irradiation), the QFLs μ Fn and μ Fp independently describe the concentrations of band electrons and holes. Carriers are exchanged between a defect state and the extended bands through capture and emission processes. It is well known that the occupation of a defect state between the bands' QFLs is ambiguously described by the equilibrium Fermi-Dirac functions associated with the QFLs and that detailed balance rate calculations are required to determine the steady-state occupancy [37][38][39] .
Following Simmons and Taylor 39 , the electron and hole capture rates are: and the rates of emission are: gives the number of band free carriers, E t is the trap's charge transition energy, E C and E V are the energies of the conduction and valance band edges, k B is the Boltzmann constant, T is temperature, and N t is the concentration of traps. The steady-state occupation of a defect by an electron is given by: n n n n C C p p It is not immediately obvious how the formation free energy for a charged defect depends on quasi Fermi energies in the presence of excess carriers since the defect can exchange carriers with both bands. However, interpreting ∆H f a as a statistical average over a large number of identical defect formation events allows an unequivocal result to be derived. As a gedanken experiment, we consider the formation of a population of native defects in the q = -1 charge state. This charge state can be achieved starting from the charge-neutral ideal crystal by the formation of a neutral defect combined with 1) the emission of a hole to the valence band reservoir having carrier chemical potential μ Fp , or 2) the capture of an electron from the conduction band having chemical potential μ Fn . The ratio of the transition rates will determine the fractions of the population participating and therefore the weighting that should be given to μ Fn and μ Fp in the statistical average. and for positively charged states (q> 0) by: The value of μ F reverts back to the Fermi energy in the case of dark thermal equilibrium. The computation of the steady-state populations of defects proceeds much as in the equilibrium case by solving for charge balance: however the coupling of the defect and excess carrier populations via defect-assisted recombination must be accounted for by simultaneously solving the steady-state excess carrier populations: In Eqs 13 and 14, i indexes over all defect species, the + N D ionized donors ( − N A ionized acceptors) are determined by μ Fn (μ Fp ), n = n o + Δ n (p = p o + Δ p) are the total carrier densities including the excess electrons Δ n (holes Δ p), G is the steady-state generation rate, U BB is the band-to-band recombination rate, U Aug is the Auger recombination rate, and the U SRH terms account for Shockley-Read-Hall recombination from each of the defects' charge states herein, treated as separate defects for simplicity. The net recombination rates in Eq. 14 are expressed as:  in which β BB is the band-to-band coefficient, and C Aug,n and C Aug,p are the electron and hole Auger coefficients. All of the quantities in Eqs 13 and 14 are expressed in terms of μ Fn and μ Fp , allowing a self-consistent solution to be found for each set of conditions. The model may be extended to account for degenerate carrier densities in the bands. For the calculations, GaSb was assumed to be either undoped (N D = N A = 0) or intentionally doped with a shallow donor impurity concentration, N D , of 10 17 cm -3 and unintentionally doped with a shallow acceptor impurity concentration, N A , of 10 13 cm -3 . Defect formation energies were taken from ref. 11, in which they were calculated using density functional theory with a hybrid functional scheme. Capture cross sections for Ga Sb antisites in the 0, -1, and -2 charge states were assumed to be σ n = 3 × 10 -16 , 1 × 10 -17 cm 2 , and 5 × 10 -18 cm 2 and σ p = 3 × 10 -16 , 1 × 10 -14 , and 2 × 10 -14 cm 2 ,10 -14 cm 2 , ensuring predominant capture of holes, while the capture cross sections for the 0, + 1 and + 2 charge states for Sb Ga defects were assumed to be σ n = 3 × 10 -16 , 1 × 10 -14 , and 2 × 10 -14 cm 2 and σ p = 3 × 10 -16 , 1 × 10 -17 , and 5 × 10 -18 cm 2 , ensuring predominant capture of electrons. The capture cross sections of Ga i defects were assumed to be σ n = 10 -14 cm 2 and σ p = 10 -17 cm 2 respectively. These capture cross-section values are well within the range of those typically accepted in semiconductors, and differ by amounts appropriate for Coulomb attraction(repulsion) for opposite carrier polarities (approximately three orders of magnitude). Variation of these cross-section values did not significantly change the results of the calculation. This can be understood because emission rates depend exponentially on charge transition level but only linearly on cross sections, so only in the rare but commonly-assumed case of recombination centers at Ei with nearly equal cross sections would the model be especially sensitive. Carrier generation rates of 10 23 and 5 × 10 24 cm -3 s -1 were examined. Assuming a carrier lifetime on the order of 1 ns, these rates are expected to produce excess carrier concentrations in the range of 10 14 to 5 × 10 15 cm -3 .