Print version ISSN 0103-9733
Braz. J. Phys. vol.29 no.4 São Paulo Dec. 1999
A Tight-Binding Study of Acceptor Levels in Semiconductors
J. G. Menchero
Instituto de Física, Universidade Federal do Rio de Janeiro
Cx. Postal 68.528, 21945-970 Rio de Janeiro, Brazil
Received February 8, 1999
Acceptor binding energies in zinc-blende semiconductors are determined within the tight-binding formalism. The importance of fitting the valence-band masses in the (100) as well as (111) directions is discussed, and parametrizations that specifically fit the valence-band anisotropy are used to calculate Ge acceptor levels in AlxGa1-x As alloys. The sensitivity of the calculated energies to the parameters that determine bulk masses is investigated, as well as the effect of varying the on-site energy of the impurity. A comparison is made between first-neighbor and second-neighbor hopping models. For shallow levels, both approaches give the same results. For deeper levels, however, important differences arise. Experimental evidence suggests that first-neighbor models are better suited for describing intermediate to deep levels.
The effective mass theory (EMT) has long been a principal tool for investigating shallow impurity states in semiconductors [1, 2]. In EMT, it is assumed that the impurity wave function is highly delocalized in real space, which in turn implies a strong localization in k space. As a result, the electronic properties of the host material may be described by only a few parameters related to the dispersion near the k point in consideration. For instance, in the Luttinger-Kohn version of k · p theory, the electronic structure of the host is completely defined in terms of the spin-orbit energy and the three Luttinger parameters (LP's), which also determine the impurity energy level. Because each host material has a different set of LP's, and acceptor binding energies are known to vary widely from host to host, it is clear that the impurity level within EMT must depend sensitively on these bulk parameters.
Another result of EMT is that, for a given material, any singly ionized acceptor should have the same binding energy, regardless of the species. However, experimentally it is known that such energies can vary widely from impurity to impurity. For example, an In acceptor in Si has a binding energy of 157 meV, compared to 45 meV for a B acceptor in the same host . Such deviations from EMT are attributed to central-cell effects, and demonstrate that the binding energies can also be very sensitive to the details of the potential in the immediate neighborhood of the impurity.
Recently, an approach was presented  for calculating impurity states in semiconductors based on the tight-binding (TB) formalism. The localized basis set of TB provides a natural description for deep levels, while the highly delocalized shallow levels are treated by means of very large unit cells together with a scaling law that allows extrapolation to the bulk limit . In other words, the TB approach is not intrinsically restricted to either the shallow or deep limits. For the TB approach to be useful, however, a clear understanding of how the calculated energies depend on the TB parameters is required. In particular, the sensitivity of EMT on the LP's and the the role of central-cell effects must be understood in the context of the current model. Another important question particular to TB concerns differences between first-nearest neighbor (1nn) and second-nearest neighbor (2nn) descriptions of the impurity state.
The aim of this paper is to carefully investigate how the impurity binding energies depend on the TB parameters that are most physically relevant to the impurity problem. The paper is organized as follows: In Sec. II the TB approach for calculating impurity states in semiconductors is briefly reviewed. In Sec. III the dependency of the impurity binding energy on the bulk effective masses and the on-site energy of the impurity is investigated. Comparisons are also made between 1nn and 2nn TB approaches, with important differences found between the two. To illustrate these effects, Ge acceptors in GaAs, AlAs, and AlxGa1-xAs alloys are considered. The summary and conclusions are given in Sec. IV.
II Calculational Details
Our TB Hamiltonian contains terms describing the bulk material as well as the impurity, and is given by 
where the roman indices denote the site and the greek indices label the spin orbitals. The define all the on-site energies and hopping for the bulk material. In this work, both 1nn and 2nn hopping are considered. The strength of the spin-orbit interaction for atom i is defined by li, and the on-site energy for atom i due to the Coulomb potential of the impurity is given by U(ri). For our basis set the sp3s* orbitals proposed by Vogl  is used, but spin is included here leading to a total of 10 basis states per site.
The perturbation potential U(ri) is described by an isotropic q-dependent screening,
where the screening parameters A, a, b, and g are taken from Bernholc . Near the origin, the potential looks like a bare Coulomb potential Ubare = e2/ri, but far away the potential looks like a bulk-screened potential Ubulk = e2/0ri, with 0 being the static dielectric constant for the host material. Precisely at the impurity site (ri = 0), Eq. (2) is undefined and the perturbation potential is assigned a value U0, where U0 is a parameter describing central-cell effects.
In order to determine the acceptor energy, first the energy Ev at the top of the valence band is calculated for the pure system, setting U(ri) = 0 for all i. This energy is easily found with the aid of Bloch's theorem. Next, a single impurity is placed in a very large cubic supercell containing 8L3 atoms arranged in the zinc-blende structure and subject to periodic boundary conditions, with L being the length of the supercell side in units of the conventional lattice parameter. The presence of the impurity breaks translational symmetry and introduces states in the gap region with energy above Ev. The energy v(L) of the highest of these states determines the acceptor energy for supercell size L via the relation E(L) = v(L) - Ev. The main task therefore is to calculate v(L), which is accomplished here using a variational algorithm that minimizes the expectation value of áY| (H-Eref)2|Yñ, with Eref being a reference energy suitably chosen near v(L) . The computation time using this scheme scales linearly with the number of states, making solutions possible even for very large supercells. For instance, system sizes ranging up to L = 20 (64000 atoms) are handled routinely . However, using a recently proposed scaling law , the impurity binding energy can be extrapolated to the bulk limit using far smaller supercell sizes. In this scaling law, the energy E(L) as a function of supercell size L is given by
with Ea being the acceptor energy in the bulk limit (L ® ¥), and P being a constant independent of L. The acceptor energy for the infinite system can therefore be found by calculating E(L) for three relatively small supercells and then solving for Ea, P, and a in Eq. (3). When using Eq. (3) to find the binding energy in the bulk limit, caution must be exercised to verify not only that the supercell sizes are in the scaling regime, but also that the scaled energies have converged to the limiting value.
In EMT, the acceptor binding energy depends on the three Luttinger parameters g1, g2, and g3. These, in turn, can be related to the heavy-hole (HH) and light-hole (LH) masses along (100) and (111) ,
where the masses in these expressions and throughout this work are in units of the free-electron mass. The fact that the LP's depend on the valence-band masses in both the (100) and (111) directions suggests that both are very important to fit properly in the TB parametrizations. On the other hand, due to the high density of states associated with the HH band, these masses are expected to be much more important than the LH masses for the purposes of determining the acceptor binding energies. Tight-binding parametrizations specifically chosen to fit the HH masses have been determined  both for 1nn and 2nn models, and will be used in this work. In Table 1 the valence-band masses resulting from these parametrizations are shown for GaAs and AlAs, together with the experimental  and theoretical  values.
|Table 1 - Heavy- and light-hole effective masses along  and , given in units of the free-electron mass.|
The Ge acceptor in AlxGa1-xAs alloys acts as a simple substitutional impurity at the As site. The binding energy has been measured experimentally  using photoluminescence in the direct-gap range 0 < x < 0.4, and the results are given by the dotted line in Fig.1. For pure GaAs, the binding energy is ~ 40 meV, but increases rapidly with an upward curvature to a relatively deep ~ 120 meV by x = 0.40. The binding energy curves Ea(x) calculated using the 1nn and 2nn parametrizations are given in Fig.1 by the solid and broken lines, respectively. The VCA was used to obtain the Hamiltonian matrix elements for the alloy. Note, of course, that the VCA has no effect on the calculated binding energies at x = 0 or x = 1, which correspond to pure GaAs and pure AlAs respectively. The impurity perturbation potential U0 was chosen in order that the calculated acceptor energy match the experimental value for pure GaAs. This led to a value of 3.42 eV for 2nn and 2.80 eV for 1nn. These values of U0 were then used for all x, so that any variation in the acceptor energy with respect to x is attributable to the electronic response of the host material, and not to central-cell effects. The result for 2nn is in qualitative agreement with experiment, showing an increasing binding energy with increasing Al content. However, it significantly underestimates the binding energy at x = 0.4, with an energy of 71 meV compared to the experimental value of 120 meV. The result is improved considerably with 1nn, giving excellent agreement up to x = 0.2. However, the results diverge beyond that, and the binding energy of 87 meV at x = 0.4 is still substantially below the experimental value.
|Figure 1. Acceptor energies Ea for Ge impurities in AlxGa1-xAs alloys, calculated with 1nn (solid line) and 2nn parametrizations (dashed line). The experimental result is given by the dotted line. Inset: acceptor energy at x = 0.4 as a function of the mass scaling factor for 1nn and 2nn parametrizations.|
In order to investigate the possibility that small uncertainties in the effective masses for AlAs might produce better agreement with experiment, the sensitivity of the binding energy to the HH masses is studied. The HH masses for AlAs are scaled by a factor ranging from 0.60 and 1.40 in such a way as to not modify any of the band energies at G . Even though scaling the HH masses by such a large fraction may have an adverse effect on some conduction-band features, these are of little importance in determining acceptor binding energies. In the inset of Fig.1, the acceptor energy at x = 0.4 is plotted as a function of , for both 1nn and 2nn cases. A remarkably linear behavior is observed over a wide range of HH masses. Interestingly, the 1nn model is much more sensitive to variations in the HH mass, with even a modest increase greatly improving agreement with experiment. For the 2nn model, however, even increasing the HH mass by 40% leads to binding energy of only 93 meV, compared to 120 meV in experiment. It is worthwhile here to comment that in order to obtain good agreement with experiment (say 110 meV at x = 0.4) a binding energy of ~ 400 meV is required for x = 1.0 (i.e., pure AlAs). This energy serves as a useful reference for the analysis considered below.
Varying the TB parameters that determine the HH masses in effect changes the bulk properties of the host. However, it is also possible to vary the TB parameters that pertain only to the impurity, without modifying the bulk properties. These parameters would therefore describe central-cell effects, and can be incorporated within the present theory through the on-site energy U0 of the impurity itself. Up to now, it has been assumed that U0 was independent of the host. Nevertheless, it is instructive to consider what effect variations in U0 have on the impurity energies. In Fig.2(a) the acceptor energy calculated for a range of U0 is presented for GaAs and AlAs using the 2nn parametrizations. The curves are characterized by three distinct regimes: an energetically flat regime for small U0, a linearly increasing regime for large U0, and a transition region for intermediate U0. In the flat region, different impurity species (characterized by different U0) will have nearly identical binding energies, meaning that central-cell corrections will be very small. Such behavior is typical of shallow levels, which implies that EMT is expected to work well in this regime.
|Figure 2. Acceptor binding energies for Ge impurities as a function of the on-site impurity potential U0 (a) Using 2nn parametrization for AlAs and GaAs. (b) Using 1nn parametrization for AlAs and GaAs. The flat region for small U0 indicates the regime in which EMT is expected to work well.|
It is interesting to note that in going from GaAs to AlAs (i.e., increasing HH mass), the energetically flat regime gets pushed to ever smaller U0. In fact, for AlAs, the flat region is never quite reached, except perhaps for U0 < 1 eV using 2nn, and so central-cell corrections are expected to always be important for this host. It should also be kept in mind that the Coulomb potential due to the impurity as given by Eq. (2) is roughly 0.6 eV at the nearest-neighbor distance (2.45 Å). Therefore, it is unphysical to consider acceptors with U0 smaller than this value. The linear regime observed for large U0 marks the breakdown of the effective mass approach. In EMT, the impurity potential is assumed to be slowly varying on the atomic length scale. A potential of say 5.0 eV at the impurity site, falling to 0.6 eV at the nearest neighbor, clearly violates this approximation. In the linear regime, central-cell effects are expected to be very important, because even small variations in U0 (due to different impurities) will lead to large changes in binding energy.
In Fig.2(b) the impurity energy is plotted as a function of U0 for AlAs and GaAs using the 1nn parametrization. Qualitatively, the results are similar to the 2nn case shown in Fig. (a). However, the 1nn case is once again found to be far more sensitive than the 2nn case. For instance, increasing U0 to just 3.3 eV is sufficient to give a binding energy of 400 meV in AlAs with 1nn, whereas a value of 5.5 eV is required with a 2nn model. Recalling that for GaAs the value of U0 was 3.42 eV for 2nn and 2.80 eV for 1nn, and considering the chemical similarity between GaAs and AlAs, a value of 5.5 eV seems unrealistic. Therefore, the interpretation once more is that the 1nn results are more consisitent with experiment.
In order to obtain a better understanding of the connection between 1nn and 2nn models, it is necessary to compare the two cases directly, using the same value of U0. According to EMT, the acceptor binding energy should depend only on the effective masses near the k point under consideration. From this point of view, no difference should exist between 1nn and 2nn parametrizations, as long as the effective masses were the same. In order to quantitatively study the validity of this assertion, however, the 1nn and 2nn TB parameters must be adjusted so as to give the same HH masses in the (100) as well as (111) directions . Next, a value of U0 sufficiently small (0.7 eV) is chosen so that the EMT should be expected to hold, and the binding energy calculated as a function of x, where as usual x is the Al concentration of the AlxGa1-xAs alloy. The results are plotted in Fig.3, and show that the 1nn and 2nn results are in almost exact agreement over the entire range of x. This is reassuring, because for shallow levels the energy can only depend on the effective mass if the TB result is to be consisitent with EMT. The second case to be considered is for larger U0, in which the EMT can be expected to break down. A value of U0 = 3.5 eV is used and the energy Ea(x) recalculated. From Fig.3, a radical difference is now observed between 1nn and 2nn parametrizations. For GaAs, the 1nn result gives 70 meV compared with only 42 meV using 2nn. In passing, it is instructive to note that the EMT result breaks down strongly for the U0 = 3.5 eV case in GaAs but not for the U0 = 0.7 eV case in AlAs, even though the binding energies are comparable. This demonstrates that binding energy considerations alone are not enough to determine if a level can or cannot be described by EMT. In any case, the breakdown of EMT is even more dramatic for the U0 = 3.5 eV case in pure AlAs, with the 1nn approach yielding an energy of 454 meV compared with 182 meV using 2nn. The experimental data suggest, therefore, that 1nn descriptions are more appropriate than 2nn for the case of deep levels. The likely reason is that the 2nn hopping terms permit the electron to escape more easily from the impurity potential, thereby delocalizing the wave function and reducing the binding energy. To confirm this conjecture, the radial charge distribution is calculated for an AlAs host using both 1nn and 2nn with an on-site energy U0 of 3.5 eV. It is found indeed that the electron is far more delocalized using the 2nn approach. For instance, using 1nn, there is a 90% probability of localizing the electron within 7 Å of the impurity, whereas using the 2nn parametrization, the same probability is reached at roughly 14 Å.
|Figure 3. Acceptor energies calculated for shallow levels (U0 = 0.7 eV) and deep levels (U0 = 3.5 eV) using 1nn (solid line) and 2nn (dotted line) parametrizations, with the same HH effective masses. The U0 = 0.7 results are consistent with EMT, whereas major deviations are observed for U0 = 3.5.|
IV Summary and Conclusions
Binding energies have been calculated for acceptors in GaAs, AlAs, and AlxGa1-xAs hosts. It was found that the 1nn model is much more sensitive than the 2nn model to small changes in the parameters. Insight was gained into why and when central-cell corrections become important. For shallow levels in which central-cell corrections are unimportant, the 1nn and 2nn models give the same results. However, for intermediate to deep levels, the two results deviate, and experimental evidence suggests that the 1nn model is more appropriate. The reason is that the 2nn model is not as effective in binding the wave function to the impurity. The 2nn parametrizations are expected to be most useful therefore for describing shallow levels, due to their superior ability to fit the effective masses. A particularly interesting application of the 2nn approach may be for the case of donors in indirect gap materials, due to intrinsic limitations in the 1nn ability to fit conduction- band dispersion near the zone boundary.
Helpful discussions with Belita Koiller, Rodrigo Capaz, Helio Chacham, and Timothy Boykin are gratefully acknowledged. The financial support of CNPq, PRONEX, and FINEP is also greatly appreciated.
 C. Kittel and A. H. Mitchell, Phys. Rev. 96, 1488 (1954). [ Links ]
 J. M. Luttinger, and W. Kohn, Phys. Rev. 97, 869 (1955). [ Links ]
 C. Kittel, Introduction to Solid State Physics, 7th Ed., (Wiley, 1996) p. 225. [ Links ]
 J. G. Menchero, R. B. Capaz, Belita Koiller, and H. Chacham, Phys. Rev. B 59, 2722 (1999). [ Links ]
 P. Vogl, H. P. Hjalmarson, and J. D. Dow, J. Phys. Chem. Solids 44, 365 (1983). [ Links ]
 J. Bernholc and S. T. Pantelides, Phys. Rev. B 15, 4935 (1977). [ Links ]
 J. K. L. MacDonald, Phys. Rev. 46, 828 (1934); [ Links ]R. B. Capaz, G. C. de Araujo, Belita Koiller, and J. P. von der Weid, J. Appl. Phys. 74 5531 (1993); [ Links ]L. W. Wang and A. Zunger, J. Chem. Phys. 100, 2394 (1994). [ Links ]
 R. Enderlein, G. M. Sipahi, L. M. R. Scolfaro, and J. R. Leite, Phys. Stat. Sol. (b) 206, 623 (1998). [ Links ]
 J. G. Menchero and T. B. Boykin, Impurity states in semiconductors calculated via tight-binding: a parameter sensitivity study, Phys. Rev. B 59, to appear (1999). [ Links ]
 R. V. Shanabrook, O. J. Glembocki, D. A. Broido, and W. I. Wang, Phys. Rev. B 39, 3411 (1989). [ Links ]
 L. Pavesi and M. Guzzi, J. Appl. Phys. 75, 4779 (1994). [ Links ]
 G. Oelgart, B. Lippold, R. Heilmann, H. Neumann, and B. Jacobs, Phys. Stat. Sol. A 115, 257 (1989). [ Links ]