Dependence of Electron Density on Fermi Energy in N-Type Gallium Antimonide

The majority electron density as a function of the Fermi energy is calculated in zinc blende, n-type GaSb for donor densities between 1016 cm−3 and 1019 cm−3. These calculations solve the charge neutrality equation self-consistently for a four-band model (three conduction sub-bands at Γ, L, and X and one equivalent valence band at Γ) of GaSb. Our calculations assume parabolic densities of states and thus do not treat the density-of-states modifications due to high concentrations of dopants, many body effects, and non-parabolicity of the bands. Even with these assumptions, the results are important for interpreting optical measurements such as Raman measurements that are proposed as a nondestructive method for wafer acceptance tests.

The majority electron density as a function of the Fermi energy is calculated in zinc blende, n-type GaSb for donor densities between 10 16 cm -3 and 10 19 cm -3 . These calculations solve the charge neutrality equation self-consistently for a four-band model (three conduction sub-bands at Γ, L, and X and one equivalent valence band at Γ ) of GaSb. Our calculations assume parabolic densities of states and thus do not treat the density-of-states modifications due to high concentrations of dopants, many body effects, and non-parabolicity of the bands. Even with these assumptions, the results are important for interpreting optical measurements such as Raman measurements that are proposed as a nondestructive method for wafer acceptance tests.
and the charge neutrality condition to compute self-consistently the Fermi energy E F and the screening radius r s for given values of the ionized dopant concentration N I and temperature T. The static dielectric constant is ε and the permittivity of free space is ε 0 . The ionized dopant concentration is positive for ntype material (donor ions) and negative for p-type material (acceptor ions). The results reported here are for uncompensated n-type material. The results for the screening radius r s are not reported here because they are not needed to extract carrier concentrations from Raman measurements.
In this paper, we use the four-band model that has three conduction sub-bands centered at the Γ, L, and X symmetry points in the Brillouin zone and one equivalent valence band centered at the Γ symmetry point. We do not include the detailed nonparabolicity of the GaSb energy bands at Γ. Unlike GaAs, the GaSb conduction Γ, L, and X sub-band masses and energy spacings are such that for donor densities of technological interest, the conduction sub-band at L is the one that is most populated. The non-parabolicity of the conduction Γ sub-band in GaAs is discussed in Ref. [2]. If we were to use the Kane three level k · p model [2], which does not include the conduction sub-bands at L and X, we would be able to include the non-parabolicity of the conduction Γ sub-band. However, because the conduction Γ sub-band band in GaSb is not the dominant band for determining the Fermi energy, its non-parabolicity correction may not have a significant effect on the results given below and may lie within the uncertainties associated with the band masses quoted in the literature for GaSb.
The heavy hole mass m hh and light hole mass m lh for the two degenerate sub-bands at the top of the valence band are combined to give an effective mass m vΓ = (m hh 3/2 + m lh for the valence topmost sub-band. The values of these parameters are given in Table 1.
The zero of energy is at the bottom of the conduction Γ sub-band. The bottoms of the conduction L and X sub-bands are, respectively, at E cL and E cX . The top of the degenerate valence Γ sub-band is at -E G , where E G is the bandgap of GaSb. The split-off valence sub-band at Γ due to spin-orbit coupling and the non-parabolicity factor of the conduction Γ sub-band are neglected. The probabilities for typical carriers in equilibrium to occupy appreciably these states are low. This means that the Fermi energies should be sufficiently above the valence sub-band maximum at Γ. Placing exact limits on the Fermi energies for which the four-band model is valid would be tenuous, because knowledge of how the various sub-bands move relative to one another due to the dopant concentrations considered here and due to many body effects is not adequate.
The general expression [3] for the temperature dependence of conduction sub-band minima relative to the top of the valence band at Γ is in units of eV, where i = Γ, L, or X. The values for the coefficients E i0 , A i , and B i are listed in Table 2. The general expression for the parabolic densities of states for electrons and holes per band extrema and per spin direction is given by (6) where N e is the number of equivalent ellipsoids in the first Brillouin zone, the volume of the unit cell is , a L is the lattice constant, m* is one of the effective masses listed in Table 1 for the appropriate band extrema, and m 0 is the free electron mass. Because eight permutations of the wave vector in the (111) direction exist, there are eight L sub-band ellipsoids with centers located near the boundary of the first Brillouin zone. Also, because six permutations of the wave vector in the (100) direction exist, there are six X sub-band ellipsoids with centers located near the boundary of the first Brillouin zone. Since about half of each ellipsoid is in the neighboring zone, the number of equivalent subbands N cL for the L sub-band is four and the number of equivalent sub-bands N cX for the X sub-band is three.
In terms of a four-band model for room temperature n-type GaSb, the total density of states ρ c (E) for the majority carrier electrons in n-type GaSb then becomes where ρ cΓ (E), ρ cL (E), and ρ cX (E) are the sub-band densities of states for the conduction Γ, L, and X sub-bands with effective masses of m cΓ , m cL , and m cX , respectively. The density of states for the minority carrier holes is with an effective mass of m vΓ .

Results
Tables 1 and 2 contain the input parameters for the calculations of the Fermi energy as a function of the dopant donor density. We solve self-consistently, by means of an iterative procedure, Eq. (3) with Eqs. (6), (7) and (8). The independent variable is the temperature T. The Fermi energy is varied for a given temperature until Eq. (3) is satisfied. Figure 1 presents the calculated data graphically for 28 values of donor densities between 10 16 cm -3 and 10 19 cm -3 . Figure 2 gives the electron densities in the conduction sub-bands at Γ and L and the total electron density as functions of the Fermi energy. Figure 2 does not show the electron density in the conduction sub-band at X, because it is less than 10 -3 times the total electron density. Because m cΓ << m cL and E cL is much closer to E cΓ than it is to E cX , the electron density in the conduction L sub-band exceeds the electron density in the conduction Γ subband at room temperature. The solid curve in Fig. 2 is the same curve as given in Fig. 1. Figure 2 shows that the majority of electrons is in the conduction L subband and that the density of electrons in the L sub-band approaches the total density of electrons as the donor density approaches 10 19 cm -3 . Hence, even though GaSb is intrinsically a direct semiconductor, the results from Fig. 2 suggest that electrons for n-type GaSb in  the vicinity of the Fermi surface will behave as though they have many characteristics of electrons in an indirect semiconductor.
For illustrative purposes, we give here only the results for fitting the logarithm to the base 10 of the total electron density n in units of cm -3 to a polynomial in E F , namely, The analytic fits for the electron densities in the Γ, L, and X sub-bands are available by sending an email to herbert.bennett@nist.gov. During the fitting analyses, we rely substantially on graphics and keep the number of fitting parameters to a minimum, subject to the constraint that the residual standard deviation S res is acceptably small, i.e., S res ≤ 0.01. The standard deviation is a measure of the "average" error in a fitted model and thereby is a metric for assessing the quality of the fit. A smaller S res indicates a better fit. The residual standard deviation for a model where Y j are the calculated data values, the are the predicted values from the fitted model, N is the total number of data points (here N = 28), and P is the total number of parameters to be fitted in the model. We use the NIST-developed DATAPLOT [4] software for both the exploratory graphics and for the statistical analyses. Table 3 gives the four fitting parameters for the cubic l = 3 polynomial fit to the log 10 (n) as shown in Eq. (9) and the associated residual standard deviation S res = 0.0066. In general, the values of S res decrease monotonically with increasing number l of terms in the polynomial given in Eq. (9). But, care must be taken to avoid fitting noise in data sets. The general guideline for    Table 3. The four fitting parameters for a cubic polynomial fit of the theoretical calculations for the total electron density in n-type, zinc blende GaSb at 300 K as a function of the Fermi energy relative to the bottom of the conduction Γ sub-band. The ratio is the estimated value divided by the estimated standard deviation. The residual standard deviation is S res = 0.0066 many data sets is that when the absolute value of the ratio R of the estimated parameter value divided by its estimated standard deviation is less than about 2, then the rate of decrease in S res with increasing l tends to decrease. For the data given in Fig. 1, when l = 2 or P = 3, S res = 0.0644; when l = 3 or P = 4, S res = 0.0066; and when l = 4 or P = 5, S res = 0.0063. Because the change in values of S res between l = 3 and l = 4 is not significant, we use the fitting parameters for the cubic l = 3 case in this paper. Also, the ratio R for the parameter a 4 when l = 4 is -1.917, and such a value for R means that proceeding with higher l values probably is not warranted. A figure that compares the calculated total electron density as a function of the Fermi energy with the fitting results from Eq. (9) for a cubic polynomial is not given because the two curves lie on top of one another within the line widths of each curve. Fits to the calculated electron densities for each of the conduction sub-bands are available from the author upon request. Also, since the screening radii for the carriers from Eq. (2) are not needed when interpreting the proposed measurements considered here, the corresponding screening radii are not presented in this paper.

Conclusions
The results from Sec. 3 are consistent with the findings of experimental work reported in the literature such as Refs. [5] and [6]. Interpreting experiments for GaSb requires at least a three-band model and under some conditions may require a four-band model. And finally, even though GaSb is intrinsically a direct semiconductor, our results show that electrons for n-type GaSb in the vicinity of the Fermi surface will have some characteristics that are similar to those for electrons in an indirect semiconductor.