Methods for enhancing P-type doping in III-V semiconductor films
09721810 · 2017-08-01
Assignee
Inventors
- Feng Liu (Salt Lake City, UT, US)
- Gerald Stringfellow (Salt Lake City, UT, US)
- Junyi Zhu (Lakewood, CO, US)
Cpc classification
Y10S117/913
GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
International classification
H01L21/20
ELECTRICITY
H01L21/02
ELECTRICITY
Abstract
Methods of doping a semiconductor film are provided. The methods comprise epitaxially growing the III-V semiconductor film in the presence of a dopant, a surfactant capable of acting as an electron reservoir, and hydrogen, under conditions that promote the formation of a III-V semiconductor film doped with the p-type dopant. In some embodiments of the methods, the epitaxial growth of the doped III-V semiconductor film is initiated at a first hydrogen partial pressure which is increased to a second hydrogen partial pressure during the epitaxial growth process.
Claims
1. A method of doping a III-V semiconductor film, the method comprising: presenting a p-type dopant which comprises cadmium to an epitaxial growth process; presenting a surfactant capable of acting as an electron reservoir to the epitaxial growth process; presenting hydrogen to the epitaxial growth process; and growing a III-V semiconductor film under conditions including the p-type dopant, the surfactant, and the hydrogen throughout the epitaxial growth process that promote the formation of a p-type doped III-V semiconductor film.
2. The method recited in claim 1, wherein the presenting a surfactant capable of acting as an electron reservoir to the epitaxial growth process comprises presenting Sb.
3. The method recited in claim 1, wherein the presenting a surfactant capable of acting as an electron reservoir to the epitaxial growth process comprises presenting Bi.
4. The method recited in claim 1, wherein the presenting hydrogen to the epitaxial growth process comprises: initiating the epitaxial growth process at a first hydrogen partial pressure; and increasing the hydrogen partial pressure from the first hydrogen partial pressure to a second hydrogen partial pressure during the epitaxial growth process.
5. The method recited in claim 1, further including annealing the p-type doped III-V semiconductor film to remove co-doped hydrogen from the III-V semiconductor film.
6. The method recited in claim 1, wherein growing the III-V semiconductor film includes using an OMVPE growth process.
7. The method recited in claim 1, wherein growing the III-V semiconductor film includes using a MBE process.
8. The method recited in claim 1, wherein growing the III-V semiconductor film includes using a CBE process.
9. The method of claim 1, wherein the III-V semiconductor film is a GaP semiconductor film.
10. A method of doping a semiconductor film, comprising: growing a semiconductor film in the presence of a p-type dopant which comprises cadmium; a surfactant; and a hydrogen source; and creating conditions during the growing of the semiconductor film that promote the formation of a semiconductor film doped with the p-type dopant, such that the p-type dopant, the surfactant and the hydrogen source are present throughout growth of the semiconductor film doped with the p-type dopant.
11. The method recited in claim 10, wherein the surfactant is capable of acting as an electron reservoir.
12. The method recited in claim 11, wherein the surfactant is Bi.
13. The method recited in claim 11, wherein the surfactant is Sb.
14. The method recited in claim 13, wherein growing a semiconductor film in the presence of a hydrogen source further includes: initiating the growth of the semiconductor film at a first hydrogen partial pressure; and increasing the hydrogen partial pressure from the first hydrogen partial pressure to a second hydrogen partial pressure during the growth of the semiconductor film.
15. The method of claim 10, wherein the semiconductor film is a GaP semiconductor film.
16. A method of doping a III-V semiconductor film, the method comprising: growing epitaxially the III-V semiconductor film entirely in the presence of a p-type dopant which is cadmium, a surfactant capable of acting as an electron reservoir, and hydrogen, under conditions that promote the formation of a p-type doped III-V semiconductor film.
17. The method recited in claim 16, wherein the growing epitaxially of the III-V semiconductor film further includes: initiating the epitaxial growth at a first hydrogen partial pressure and maintaining the first hydrogen partial pressure for a first period of time during the growing; and increasing the hydrogen partial pressure to a second hydrogen partial pressure and maintaining the second hydrogen partial pressure for a second period of time during the growing.
18. The method of claim 16, further including annealing the p-type doped III-V semiconductor film at a temperature and for a time sufficient to remove hydrogen from the p-type doped III-V semiconductor film.
19. The method of 16, wherein the III-V semiconductor film is a GaP semiconductor film.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
(1) In order to describe the manner in which the above-recited and other advantages and features of the invention can be obtained, a more particular description of the invention briefly described above will be rendered by reference to specific example embodiments thereof which are illustrated in the appended drawings. Understanding that these drawings depict only typical implementations of the invention and are not therefore to be considered to be limiting of its scope, the invention will be described and explained with additional specificity and detail through the use of the accompanying drawings.
(2)
(3)
(4)
(5)
(6)
(7)
DETAILED DESCRIPTION OF THE INVENTION
(8) As used herein, and unless otherwise specified, “a” or “an” means “one or more.” As will be understood by one skilled in the art, for any and all purposes, particularly in terms of providing a written description, all ranges disclosed herein also encompass any and all possible subranges and combinations of subranges thereof. Any listed range can be easily recognized as sufficiently describing and enabling the same range being broken down into at least equal halves, thirds, quarters, fifths, tenths, etc. As a non-limiting example, each range discussed herein can be readily broken down into a lower third, middle third and upper third, etc. As will also be understood by one skilled in the art, all language such as “up to,” “at least,” “greater than,” “less than,” and the like includes the number recited and refers to ranges which can be subsequently broken down into subranges as discussed above. Finally, as will be understood by one skilled in the art, a range includes each individual member.
(9) One example embodiment of the present invention provides one or more methods for p-type doping III-V semiconductors, such as GaP, using a surfactant capable of acting as an electron reservoir. Examples of such surfactants include such as Sb and Bi. In addition, example embodiments of the present invention provide methods for p-type doping III-V semiconductors using a surfactant in the presence of H. For example, the use of the above example surfactants in the presence of H can=create what can be called a dual-surfactant effect of Sb and H for p-type doping enhancement for dopants such as Zn, Cd, Mg and Be. The dual-surfactant effect, which relies on two surface elements, greatly broadens the scope and application of the conventional surfactant effect which relies on only one surface element. The dual-surfactant effect relies on two surface H atoms and an H in the bulk serves as a co-dopant. Specifically, in order to accommodate the p-type dopant incorporation, the role of the surfactant (e.g., Sb—a “metallic” element) is to provide an electron reservoir to redistribute electrons, the role of the two surface H atoms is to satisfy the electron counting rule (ECR) at the surface before doping, and the role of the co-doped H (a “single electron”) is to add one electron to satisfy the ECR in the bulk after the doping. Later annealing allows the co-doped H to diffuse away.
(10) For the purposes of this disclosure, a surfactant in epitaxial growth is defined as an active element that floats on top of the growing film surface. It is usually a single foreign element. Here, the term “film” doping energy is used to refer to the energy change for a dopant atom (e.g., Zn) replacing a cation atom (e.g., Ga) in the sub-surface position during the epitaxial growth process, as distinguished from the “bulk” doping energy for a dopant atom replacing a cation atom in the bulk. The methods described herein can also be incorporated in MBE and CBE processes. In addition, the term “under conditions that promote the formation of III-V semiconductor film doped with p-type dopant,” includes a variety of such conditions that are known to someone skilled in the art.
EXAMPLES
(11) First-principles calculations of Zn, Mg, Be and Cd incorporation in (001) GaP films under the influence of surface Sb and H were carried out. It was found that Sb alone has little effect on the film doping energy of all the dopants in GaP film, and it was only when H was also present that the film doping energy was substantially lowered by Sb. Also, surface H did not function as effectively alone without Sb. Without intending to be bound to any particular theory of the invention, the inventors believe that it is the combined effect of Sb and H (a dual-surfactant effect) that makes the p-type doping processes thermodynamically favorable. The role of Sb in producing the dual-surfactant effect is to serve as an electron reservoir to accommodate the redistribution of electrons, in a similar spirit to the generalized ECR in the semiconductor surface with metal elements [14]. The role of H is to supply the one electron missing from the p-type dopant, so that the system can satisfy the ECR [15]. Experimentally, it is difficult to achieve high doping levels of Mg and Zn in III-V systems [8,16]. The present invention represents an important breakthrough that provides a common strategy for improved p-type doping during epitaxial growth of III-V compounds.
(12) In order to arrive at the present invention, calculations were performed using the Vienna ab initio simulation package [17] within the local density approximation. GaP (001) films were modeled by the same technique as applied in the calculation on the dual surfactant effect of Sb and H on Zn in GaP described in Zhu et al., Phys. Rev. Lett. 101 (2008) 196103, the entire disclosure of which is incorporated herein by reference.
(13) In the calculation, the (2×2) reconstructed surface was chosen, as shown in (E
.sub.undoped) is the total energy of the doped (undoped) system, i.e., supercell; and μ.sub.Ga(μ.sub.Dopant) is the chemical potential of Ga (dopant). In general, μ.sub.Ga may vary from μ.sub.Ga[bulk]+ΔH.sub.f[GaP] (the P-rich condition) to μ.sub.Ga[bulk] (the Ga-rich condition) [21], where ΔH.sub.f[GaP] is the GaP enthalpy of formation; and μ.sub.Dopant equals μ.sub.Dopant[bulk]. In calculating the change of film doping energy due to surfactant Sb, μ.sub.Ga and μ.sub.Zn do not appear. In the case of H, the film doping energy depends on the chemical potential of H
(μ
.sub.H) if one additional H is added to the system upon doping. The chemical potential of H is a variable strongly depending upon the growth condition. Here, the typical value of μ.sub.H=−0.67 eV, one-half of the energy of an H.sub.2 molecule at T=900K, p=1 atm was chosen [22]. This is different from the value used in the analysis of Zhu et al., Phys. Rev. Lett. 101 (2008) 196103, which would correspond to a H partial pressure of 10.sup.−10 atm, unrealistically low for OMVPE growth. Also, in OMVPE growth, there is a large catalytic effect of the phosphide surface on the decomposition of precursors [23]. So the actual chemical potential of the atomic H may be difficult to estimate due to non-equilibrium conditions at the surface. The chemical potential given here in the gas phase may be considered as the lower limit, because the actual H chemical potential could be higher due to the surface adsorption and decomposition of the precursor [23]. This would lead to a larger doping energy drop due to the extra H than our estimation.
(14) First, the doping energy differences of bulk vs. film without surfactant were studied. The p-type doping energy in bulk GaP was calculated as a reference. Then, the p-type film doping energy in the surface positions of a GaP (0 0 1) film was calculated by replacing the Ga in the first cation layer (
(15) Next the surfactant effect of Sb was investigated. When Sb is introduced during OMVPE growth, it stays on top of the surface replacing the P surface dimers, due to its large atomic size and lower dangling bond energy (see
(16) Importantly, it was found that for all p-type dopants studied, the film doping energy difference, ΔE.sub.1 (a vs. b in
(17) TABLE-US-00001 TABLE 1 Doping energies of bulk and film without surfactant for different dopants. Bulk Film E.sub.Film − E.sub.Bulk (eV) Be 0.23 eV + μ.sub.Ga − μ.sub.Be 0.02 eV + μ.sub.Ga − μ.sub.Be −0.21 Zn 2.69 eV + μ.sub.Ga − μ.sub.Zn 2.46 eV + μ.sub.Ga − μ.sub.Zn −0.23 Mg 2.61 eV + μ.sub.Ga − μ.sub.Mg 2.14 eV + μ.sub.Ga − μ.sub.Mg −0.47 Cd 3.69 eV + μ.sub.Ga − μ.sub.Cd 2.74 eV + μ.sub.Ga − μ.sub.Cd −0.95
(18) TABLE-US-00002 TABLE 2 Film doping energy difference for p-type dopants with different surface configurations in GaP as shown in FIG. 2. Zn Mg Cd Be ΔE.sub.1(eV) 0.07 −0.15 0.19 0.36 ΔE.sub.2(eV) 0.22 0.07 0.55 0.41 ΔE.sub.3(eV) −0.78 −0.47 −0.49 −0.92 ΔE.sub.4(eV) −0.93 −0.69 −0.85 −0.97
(19) Next, the role of H in the doping process was studied. OMVPE growth produces atomic H on the surface [23]. A significant concentration of H (presumably in the form of Zn—H and C—H complexes) is observed in the GaP epitaxial films [8]. Surface H allows the surface to satisfy the ECR [15]. The clean GaP (0 0 1)−(2×2) model surfaces discussed above do not satisfy the ECR. This provides a thermodynamic driving force for H to incorporate into the film.
(20) In order for the GaP (0 0 1)−(2×2) surface to satisfy the ECR, one H can be added to each surface dimer on alternating sides, causing dimer buckling. This has been shown by calculations to have the lowest energy [13]. When Sb is incorporated as a surfactant, the H bonds to the Sb dimers in the same way (see
(21) Next, the combined effects of Sb and H were investigated. Assuming again a surface covered with 2H/cell before and after the p-type dopant incorporation (see
(22) The above results indicate that the surfactant effect of Sb is enhanced when surface H is also introduced. The underlying physical reason giving rise to the effect of Sb is probably due to the lower electronegativity of Sb in comparison with P, in a similar spirit to the generalized ECR in the semiconductor surface with metal elements proposed recently [14]. Antimony is more metallic than P, so that Sb can serve more effectively as an electron reservoir to accommodate the distribution of electrons when a p-type dopant is present. Since the p-type dopant will violate the ECR by having one less electron, it is easier for Sb than for P to accommodate the missing electron in order to “partially” satisfy the ECR. Also, it can be seen that the dual surfactant effect in decreasing the film doping energy is smallest for Mg, compared with the other three elements. This is probably due to the electronegativity difference among these dopants. Mg, Be, Cd and Zn have electronegativities of 1.31, 1.57, 1.69, 1.65, respectively [25]. So the electronegativity of Mg is noticeably the lowest among all these four elements. Thus, relatively more electrons of Mg can contribute to the electron redistribution, and the dopant Mg itself can be considered as an electron reservoir which might negate part of the electron reservoir effect that Sb contributes to the system. Consequently, Sb and H have the smallest dual surfactant effect for Mg incorporation.
(23) Next, the surface configuration with three surface hydrogen atoms per four Sb atoms was investigated. The ECR is violated after the dopant is incorporated if the surface H remains constant. In order to satisfy the ECR, one additional H has to be added. In a previous study of Zn, to quantify the role of the extra H, the Zn film doping energy was calculated by assuming a surface covered with 2H before, but with 3H after doping (see
(24) So far, the dual-surfactant effect of Sb and H in enhancing the p-type doping in the first cation layer, i.e., the surface position, has been shown. It is also important to investigate the configuration where the dopant replaces a Ga atom at the second cation layer, i.e., the subsurface or “bulk” position. In the previous study of Zn, two possible configurations were considered: one with Zn in the second cation layer between two surface Sb dimers (
(25) The 0.08 eV difference between the two reflects the dependence of Zn film doping energy on the “atomic-level” stress at these sites [26] and Zn is slightly favored at the tensile sites between surface dimers relative to the compressive sites directly below surface dimers. Here, this treatment has been expanded to include Mg/Cd/Be and a similar result was found. The respective film doping energies are: ΔE.sub.2ndMg,between.sup.Sb-2H.fwdarw.2H=0.3 eV+μ.sub.Ga−μ.sub.Mg and ΔE.sub.2ndMg,below.sup.Sb-2H.fwdarw.2H=0.58 eV+μ.sub.Ga−μ.sub.Mg; ΔE.sub.2ndCd,between.sup.Sb-2H.fwdarw.2H=1.1 eV+μ.sub.Ga−μ.sub.Cd and ΔE.sub.2ndCd,below.sup.Sb-2H.fwdarw.2H=1.46 eV÷μ.sub.Ga−μ.sub.Cd; ΔE.sub.2ndBe,between.sup.Sb-2H.fwdarw.2H=−2.25 eV+μ.sub.Ga−μ.sub.Be and ΔE.sub.2ndBe,below.sup.Sb-2H.fwdarw.2H=−2.31 eV+μ.sub.Ga−μ.sub.Be. The 0.28 eV energy difference between the two sites for Mg and 0.36 eV energy difference between the two sites for Cd, which are larger than that of Zn are probably due to the larger covalent radii of Mg and Cd [27]. So the tensile site is even more favored by a larger Mg and Cd. Be is slightly favored at the compressive site, which is due to the smaller covalent radius of Be.
(26) Lastly, co-doping of the p-type dopant and H into the GaP bulk was investigated. Experimental observation indicated that complexes of Zn, P and H form during doping of Zn in GaP [28]. This suggests some H goes into the bulk with the acceptor. In the previous Zn study, the energies associated with one H incorporation into the subsurface with Zn was calculated, and the Zn—P—H complex structures were determined.
(27) When a surface H goes into the subsurface, it changes the dopant bonding configuration, forming a dopant-P—H complex. Without H, the dopant bonds with four neighboring P atoms in an sp.sup.3 hybridization (a tetrahedral structure,
(28) The co-doping of H with dopant into the subsurface is also found to be energetically favorable. The film doping energy at the 2nd-cation-layer positions forming the Zn—P—H complex was found previously to be ΔE.sub.2ndZn.sup.Sb-2H.fwdarw.1bulkH=−0.19 eV+μ.sub.Ga−μ.sub.Zn [13], which is about 0.2 eV lower than the corresponding cases with all 3H atoms remaining on the surface. This suggests that there exists a thermodynamic driving force for one H atom to go into the subsurface (or “bulk”) with the dopant, i.e., co-doping of the Zn with H. The 2.sup.nd cation layer Mg/Cd/Be doping calculation was also performed. The film doping energy was ΔE.sub.2nd,Mg.sup.Sb-2H.fwdarw.1bulkH=−0.29 eV+μ.sub.Ga−μ.sub.Mg, and ΔE.sub.2nd,Cd.sup.Sb-2H.fwdarw.1bulkH=0.80 eV+μ.sub.Ga−μ.sub.Cd and ΔE.sub.2ndMg,Be.sup.Sb-2H.fwdarw.1bulkH=−2.59 eV+μ.sub.Ga−μ.sub.Be. All the energies were lower than the corresponding case with all 3H atoms on the surface. So, similar to H co-doping with Zn into the bulk, it is confirmed that there is a thermodynamic driving force for one H atom to go into the bulk with Mg/Cd/Be. On the other hand, it is known that H would compensate the p-type dopant in bulk GaP, mitigating the doping effect [29,30], which is also confirmed by electronic structure calculations. To activate the acceptor, annealing can be done after the doping process to remove the H [29,30]. In other words, the co-doped H may be removed after it serves its purpose to assist the dopant incorporation.
(29) Now, based on all the above calculations, a plausible complete p-type doping process during OMVPE growth is postulated, driven by the dual surfactant effect of Sb and H as shown in
(30) In a subsequent example step, one additional H is added to the surface to further assist the doping process by adding one more electron to satisfy the ECR (
(31) The corresponding film doping energy changes are shown in
(32) Consequently, there is an overall downhill energy landscape for the whole doping process thermodynamically. Notice that the second-step film doping energy with one H added to the system depends on the H chemical potential, while the third-step H film co-doping energy with one H moving from surface to bulk does not. This implies that one may increase the partial pressure of H and hence H chemical potential during growth to further enhance p-type doping in the second step.
(33) This theoretical picture of the dual-surfactant effect and its underlying physical mechanism based on the ECR is qualitatively consistent with the experimental observation of both the enhanced Zn doping with Sb introduction during OMVPE growth and the co-incorporation of H with Zn [7,8]. The overall trend of the calculated energy changes is expected to be correct.
(34) Accordingly,
(35) For example,
(36) In addition, the method 600 comprises an act 604 of presenting a surfactant to the epitaxial growth process. For example, a surfactant that can act as an electron reservoir can be presented to the epitaxial growth process. For example, as described above, surfactants that can act as the electron reservoir include Sb and Bi.
(37) In addition, the method 600 comprises an act 606 of presenting hydrogen to the epitaxial growth process. For example, hydrogen can be exposed to the epitaxial growth process at a variety of partial pressures.
(38) Furthermore, the method 600 comprises an act 608 of epitaxially growing a III-V semiconductor film. For example, the III-V semiconductor film can be grown under conditions that promote the formation of a p-type doped semiconductor film.
(39) Accordingly, the diagrams and figures provided in
(40) The present invention may be embodied in other specific forms without departing from its spirit or essential characteristics. The described embodiments are to be considered in all respects only as illustrative and not restrictive. The scope of the invention is, therefore, indicated by the appended claims rather than by the foregoing description. All changes which come within the meaning and range of equivalency of the claims are to be embraced within their scope.
REFERENCES
(41) The above Background and Detailed Description cite to various references. A complete list of these references is included below, and each is incorporated herein by reference in its entirety. All patents, applications, references, and publications cited herein are incorporated by reference in their entirety to the same extent as if they were individually incorporated by reference. [1] M. Copel, M. C. Reuter, E. Kaxiras, R. M. Tromp, Phys. Rev. Lett. 63 (1989) 632. [2] G. Rosenfeld, et al., Phys. Rev. Lett. 71 (1993) 895. [3] J. Meyer, J. Vrijmoeth, H. van der Vegt, R. Behm, Phys. Rev. B 51 (1995) 14790. [4] M. R. Pillai, S.-S. Kim, S. T. Ho, S. A. Barnett, J. Vac. Sci. Technol. B 18 (2000) 1232. [5] E. Rudkevich, F. Liu, D. E. Savage, T. F. Kuech, L. McCaughan, M. G. Lagally, Phys. Rev. Lett. 81 (1998) 3467. [6] C. M. Fetzer, R. T. Lee, J. K. Schurtleff, G. B. Stringfellow, S. M. Lee, T. Y. Seong, Appl. Phys. Lett. 76 (2000) 1440. [7] D. C. Chapman, A. D. Howard, G. B. Stringfellow, J. Crystal Growth 287 (2006) 647. [8] A. D. Howard, D. C. Chapman, G. B. Stringfellow, J. Appl. Phys. 100 (2006) 44904. [9] R. R. Wixom, L. W. Rieth, G. B. Stringfellow, J. Crystal Growth 265 (2004) 367-374. [10] W. Shockley, J. L. Moll, Phys. Rev. 119 (5) (1960) 1480. [11] S. J. Pearton, H. Cho, F. Ren, J. I. Chyi, J. Han, R. G. Wilson, in: R. Feenstra, T. Myers, M. S. Shur, H. Amano (Eds.), GaN and Related Alloys, vol. 595, Materials Research Society, 2000, pp. W10.6.1-W10.6.10. [12] J. K. Shurtleff, S. W. Jun, G. B. Stringfellow, Appl. Phys. Lett. 78 (2001) 3038. [13] J. Y. Zhu, F. Liu, G. B. Stringfellow, Phys. Rev. Lett. 101 (2008) 196103. [14] L. Zhang, E. G. Wang, Q. K. Xue, S. B. Zhang, Z. Zhang, Phys. Rev. Lett. 97 (2006) 126103. [15] P. K. Larsen, D. J. Chadi, Phys. Rev. B 37 (1988) 8282. [16] S. D. Burnham, G. Namkoong, W. Henderson, W. A. Doolittle, J. Crystal Growth 279 (2005) 26-30. [17] G. Kresse, J. Hafner, Phys. Rev. B 49 (1994) 14251; G. Kresse, J. Furthmuller, Comput. Mater. Sci. 6 (1996) 15. [20] R. Wixom, N. A. Modine, G. B. Stringfellow, Phys. Rev. B 67 (2003) 115309. [21] C. G. VandeWalle, J. Neugebauer, Phys. Rev. Lett. 88 (2002) 066103. [22] O. Knacke, O. Kubaschewski, K. Hesselmann (Eds.), Thermochemical Properties of Inorganic Substances, second Ed., Springer-Verlag, 1991, pp. 728, 743, 803. [23] G. B. Stringfellow, Organometallic Vapor-Phase Epitaxy, Theory and Practice, second ed., Academic Press, 1999, pp. 141, 243. [24] M. Sugiyama, S. Maeyama, F. Maeda, M. Oshima, Phys. Rev. B 52 (1995) 2678. [25] A. L. Allred, J. Inorg, Nucl. Chem. 17 (1961) 215. [26] F. Liu, M. G. Lagally, Phys. Rev. Lett. 76 (1996) 3156. [27] B. Cordero, V. Gmez, A. E. Platero-Prats, M. Revs, J. Echeverra, E. Cremades, F. Barragn, S. Alvarez, Dalton Trans. 2008 (2008) 2832-2838. [28] M. D. McCluskey, E. E. Haller, J. W. Walker, N. M. Johnson, Appl. Phys. Lett. 65 (1994) 2191. [29] E. V. K. Rao, B. Theys, Y. Gottesman, H. Bissesur, in: 11.sup.th International Conference on Indium Phosphide and Related Materials, IEEE, 0-7803-5562-8/99, 1999. [30] C. H. Chen, S. A. Stockman, M. J. Peanasky, C. P. Kuo, in: G. B. Stringfellow, M. G. Craford (Eds.), Semiconductors and Semimetals, vol. 48, Academic Press, New York, 1997, pp. 122-127. [31] M. D. Sturge, Statistical and Thermal Physics: Fundamentals and Applications, illustrated ed., A K Peters, Ltd., 2003, p. 136. [32] M. Ilegems, H. C. Montgomery, J. Phys. Chem. Solids 34 (1973) 885.