Electrical Properties of Dislocations within the Nitride Based Semiconductors Gallium Nitride and Indium Nitride by Erfan Baghani B.Sc., Sharif University of Technology, 2005 M.Sc., Sharif University of Technology, 2007 A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY in The College of Graduate Studies (Applied Science) THE UNIVERSITY OF BRITISH COLUMBIA (Okanagan) October 2012 c Erfan Baghani, 2012Abstract Dislocation lines a ect the electrical and optical properties of semiconduc- tors. In this research, the e ect that the threading dislocation lines have on the free electron concentration and the electron mobility within gallium nitride and indium nitride is investigated. A formulation is developed for obtaining the screening space charge concentration and the corresponding electrostatic potential pro le surrounding the dislocation lines. The resultant electrostatic potential pro le has then been used to compute the associated electron mobility, limited by scattering from the charged dislocation lines. As part of this research, a Gibbs factor formalism is also developed that can readily obtain the occupation statistics of the defect sites associated with the threading dislocation lines. iiPreface This thesis, entitled \Electrical Properties of Dislocations within the Nitride Based Semiconductors Gallium Nitride and Indium Nitride," presents the research performed by Erfan Baghani. This research was supervised by Dr. Stephen K. O’Leary of The University of British Columbia. A portion of the contents of Chapters 2 and 3 is published (Baghani and O’Leary 2011) \Electron mobility limited by scattering from screened pos- itively charged dislocation lines within indium nitride", Appl. Phys. Lett., vol. 99, pp. 262106-1-3. These results were produced using an iterative pro- cedure that was developed by the author of this thesis. A version of the contents of Section 4.5 is published (Baghani and O’Leary 2011) \Occupation statistics of dislocations within uncompensated n-type wurtzite gallium nitride", J. Appl. Phys., vol. 109, pp. 113706-1-6. This section includes the original work of the thesis author, on the development and application of a Gibbs factor formalism to the problem of determining the occupation statistics of the dislocation lines with an assumed dangling bond dislocation core structure. A version of the results presented in Section 4.6 is published (Baghani and O’Leary 2011) \Occupation statistics of the VGa ON dislocations within n- type gallium nitride", J. Appl. Phys., vol. 110, pp. 033509-1-6. This section iiiincludes the original work of the thesis author, on the application of a Gibbs factor formalism to the problem of determining the occupation statistics of the dislocation lines with an assumed VGa ON dislocation core structure. ivContents Abstract ii Preface iii Table of Contents v List of Tables ix List of Figures x Acknowledgements xvi 1 Introduction 1 1.1 Motivation for the study of the nitride based semiconductors . 1 1.2 The crystalline structure of GaN and InN . . . . . . . . . . . . 8 1.3 Material parameters associated with wurtzite GaN and wurtzite InN . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10 v1.4 Atomic structure of the edge type dislocations within the wurtzite GaN and wurtzite InN crystal structures . . . . . . . . . . . . 11 1.5 Electrical properties of the threading dislocation lines within GaN and InN . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 2 Dislocation line charge screening 20 2.1 Modeling the charge build-up along the dislocation lines . . . 20 2.2 Formulating the screening space charge distribution surround- ing a dislocation line . . . . . . . . . . . . . . . . . . . . . . . 27 2.3 Threading dislocation line charge screening within GaN . . . . 31 2.4 Comparison with experiment . . . . . . . . . . . . . . . . . . . 34 2.5 Threading dislocation line charge screening within InN . . . . 35 3 Electron mobility limited by scattering from charged dislo- cation lines 41 3.1 The electron drift and Hall mobilities in semiconductors . . . . 41 3.2 Formulation . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46 3.3 Electron mobility limited by scattering from the charged thread- ing dislocation lines within GaN . . . . . . . . . . . . . . . . . 52 3.4 Electron mobility limited by scattering from the charged thread- ing dislocation lines within InN . . . . . . . . . . . . . . . . . 55 vi4 Occupation statistics of the dislocation defect sites 58 4.1 Di erent approaches for the determination of the dislocation defect site occupation statistics . . . . . . . . . . . . . . . . . 58 4.2 Gibbs factor formalism . . . . . . . . . . . . . . . . . . . . . . 62 4.3 Approximating the electrostatic potential induced by the dis- location core charge . . . . . . . . . . . . . . . . . . . . . . . . 69 4.4 Approximating the screening space charge distribution function 73 4.5 Occupation statistics of the dangling bond dislocation core structure within GaN . . . . . . . . . . . . . . . . . . . . . . . 78 4.6 Occupation statistics of the VGa ON dislocation core struc- ture within GaN . . . . . . . . . . . . . . . . . . . . . . . . . 86 4.7 Occupation statistics of the open core dislocation core struc- ture within InN . . . . . . . . . . . . . . . . . . . . . . . . . . 92 5 Conclusions 99 References 104 Appendices 119 Appendix A. Iterative procedure for solving the screening space charge distribution within GaN . . . . . . . . . . . . 119 viiAppendix B. Iterative procedure for solving the screening space charge distribution within InN . . . . . . . . . . . . 124 viiiList of Tables 1.1 The material parameters corresponding to wurtzite GaN and wurtzite InN. These material parameters are drawn from the references indicated within the square brackets. The bulk donor energy levels are provided with respect to the valence band maximum. . . . . . . . . . . . . . . . . . . . . . . . . . . 10 ixList of Figures 1.1 A transmission electron microscope image of the threading dislocation lines present within epitaxially grown InN. This image is from Cimalla et al. [37]. . . . . . . . . . . . . . . . . 7 1.2 (a) The zinc-blende and (b) wurtzite crystal structures. . . . . 9 1.3 A section of a cubic crystal lattice, including a pure edge (A) and a pure screw dislocation (B). This gure is from Darling [52]. 12 1.4 Schematic atomic con gurations of the: (a) 4-ring, (b) 5/7- ring, (c) 8-ring, (d) gallium vacancy, (e) nitrogen vacancy, and (f) VGa ON dislocation core structures within wurtzite GaN. Parts (a), (b), and (c) are after Takie and Nakayama [59]. Parts (d) and (e) are after Lee et al. [57]. Part (f) is after You et al. [60]. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13 x1.5 The atomic structure of a dislocation line within the diamond crystal structure. This gure is after Read [73]. The direction of the dislocation line is indicated with the dashed line in this gure. The symbol ? is frequently used in order to denote an edge type dislocation in the literature. . . . . . . . . . . . . . 18 2.1 Modeling the space charge build-up around a threading dis- location line. Rdc denotes the dislocation core radius. Rsc represents the screening space charge cylinder radius. Rsc will be introduced subsequently. . . . . . . . . . . . . . . . . . . . 22 2.2 Schematic diagram for the energy band bending around a threading dislocation line within GaN. . . . . . . . . . . . . . 24 2.3 Schematic diagram for the energy band bending around a threading dislocation line within InN. . . . . . . . . . . . . . . 25 2.4 The screening space charge concentration of a threading dis- location line within GaN. The bulk donor concentration, Nd, is set to 1018 cm 3 and the temperature set to 300 K for the purposes of this analysis. . . . . . . . . . . . . . . . . . . . . . 32 xi2.5 The magnitude of the electrostatic potential surrounding a threading dislocation line within GaN. The bulk donor con- centration, Nd, is set to 1018 cm 3 and the temperature set to 300 K for the purposes of this analysis. . . . . . . . . . . . . . 33 2.6 Comparison of the results from the present analysis with the experimental data of Cherns et al. [69]. For this analysis, f is set to 2 and the bulk doping concnetration, Nd, is set to 6 1017 cm 3, this being the value measured by Cherns et al. [69]. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36 2.7 Screening space charge concentration of a threading disloca- tion line within InN. The bulk donor concentration, Nd, is set to 1018 cm 3 and the temperature set to 300 K for the purposes of this analysis. . . . . . . . . . . . . . . . . . . . . . . . . . . 38 2.8 The magnitude of the electrostatic potential surrounding a threading dislocation line within InN. The bulk donor concen- tration, Nd, is set to 1018 cm 3 and the temperature set to 300 K for the purposes of this analysis. . . . . . . . . . . . . . 39 3.1 (a) Scattering of a plane wave with a wave vector ~k?, from a cylindrically symmetric potential, and (b) de ning the wave vector di erence, ~ ?. . . . . . . . . . . . . . . . . . . . . . . . 47 xii3.2 (a) The electron drift mobility and (b) the electron Hall mo- bility limited by scattering from charged threading dislocation lines within GaN. . . . . . . . . . . . . . . . . . . . . . . . . . 53 3.3 (a) The electron drift mobility and (b) the electron Hall mo- bility limited by scattering from charged threading dislocation lines within InN. . . . . . . . . . . . . . . . . . . . . . . . . . 56 4.1 The schematic diagram of energy levels in the band gap for four di erent dislocation core structures within GaN. This g- ure is after Lee et al. [57]. . . . . . . . . . . . . . . . . . . . . 65 4.2 The four possible occupation states of a dislocation defect site. 66 4.3 Comparison of the exact and the suggested approximation to the value of the summation, PN n=1 1 n , as a function of N . . . . 72 4.4 The approximation of the screening space charge distribution function surrounding a threading dislocation line within GaN by a step function. . . . . . . . . . . . . . . . . . . . . . . . . 74 4.5 The approximation of the screening space charge surrounding a dislocation line within InN by an exponential function. . . . 77 4.6 The fraction of ionized and neutral dislocation defect sites for a dangling bond core structure of the dislocation lines within n-type GaN as a function of the bulk donor concentration. . . 83 xiii4.7 The average negative charge per dislocation dangling bond as a function of the free electron concentration within n-type GaN. The calculations were all performed at 300 K for a dislocation line density of 108 cm 2. . . . . . . . . . . . . . . . . . . . . . 85 4.8 Bonding and anti-bonding energy formation from the two ni- trogen dangling bonds within the VGa ON dislocation core structure. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87 4.9 The fraction of ionized and neutral dislocation defect sites for the VGa ON core structure of the dislocation lines within n-type GaN as a function of the bulk donor concentration. . . 89 4.10 The average negative charge per VGa ON dislocation defect site as a function of the free electron concentration within n- type GaN. Results corresponding to the energy minimization approach employed by You et al. [60] for the same VGa ON dislocation core structure are also depicted. The results de- picted correspond to a temperature of 300 K and a dislocation line density of 108 cm 2. . . . . . . . . . . . . . . . . . . . . . 91 4.11 Schematic energy band diagram of the neutral state of the open core dislocation core structure within InN. . . . . . . . . 93 xiv4.12 The fraction of ionized and neutral dislocation defect sites for the open core structure of the dislocation lines within n-type InN as a function of the bulk donor concentration. . . . . . . . 96 4.13 The average positive charge per dislocation defect site for the open core structure of the dislocation lines within n-type InN as a function of the free electron concentration. The calcula- tions are performed at 300 K for a dislocation line density of 108 cm 2. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 98 A.1 An iterative construction of the screening space charge sur- rounding a threading dislocation line within n-type GaN. . . . 121 B.1 An iterative construction of the screening space charge sur- rounding a threading dislocation line within n-type InN. . . . 126 xvAcknowledgements I would like to express my sincere thanks to my supervisor, Dr. Stephen O’Leary, who gave me the chance for being involved in this interesting and challenging research project, and spent hundreds of hours of his time dis- cussing various aspects of my work. His warm attitude has always inspired and encouraged me through these years. I would also like to thank my com- mittee members, Dr. Jonathan Holzman and Dr. Kenneth Chau, for the fruitful discussions and helpful hints, and for providing me with support along the way. The nancial support for my study was provided by a research grant from the Natural Sciences and Engineering Research Council of Canada (NSERC) and from the College of Graduate Studies of the Okanagan campus of The University of British Columbia. xviChapter 1 Introduction 1.1 Motivation for the study of the nitride based semiconductors The III-V nitride semiconductors, gallium nitride (GaN), indium nitride (InN), and aluminum nitride (AlN), possess a number of interesting material properties that make them suitable for a wide variety of important electronic and optoelectronic device applications [1, 2, 3, 4]. In terms of electronics, the wide energy gap of GaN, 3.39 eV at room temperature, makes it possible to use this material in high-frequency and/or high-power electron device ap- plications [5, 6]. In terms of optoelectronics, the direct nature of the energy 1gap associated with these materials, and the fact that ternary and quater- nary alloys of these materials can be epitaxially grown, allows for their use in optoelectronic device applications, ranging from the near-infrared to the ultraviolet [7]. While initial attempts to fabricate electron devices with these materials were hindered by growth di culties, many of these problems have now been solved. As a consequence, functional electron devices have been fabricated using GaN, InN, and AlN, and alloys of these materials [7, 8, 9]. This has generated considerable interest in the III-V nitride semiconductors. Pioneering investigations into GaN were performed in the early 1930s by Johnson et al. [10]. The GaN samples considered in these early investigations were prepared through the interaction of ammonia gas with metallic gallium at elevated temperatures. Later, reactive sputtering was used in order to grow thin lms of this material. Workers, such as Hovel and Cuomo [11], Malyutin et al. [12], Vesely et al. [13], Hariu et al. [14], Lakshmi et al. [15], Matsushita et al. [16], and Kubota et al. [17], employed this approach in order to understand the physical properties of this material. Unfortunately, the materials produced were polycrystalline in nature and not of the high quality crystalline form required for the electron device applications presently envisioned for this material. A successful means of epitaxially growing high quality crystals of GaN was 2found through the use of chemical vapor deposition (CVD). This approach to growing crystalline GaN was rst employed by Maruska and Tietjen in 1969 [18]. In their approach, HCl vapor was allowed to react with metallic gallium, the resultant chemical species reacting with ammonia gas at the substrate. A number of variations on this approach have been developed through the years, including those by Pankove et al. [19], Crouch et al. [20] and Adonin et al. [21]. CVD remains a popular means whereby crystalline GaN may be fabricated. Molecular beam epitaxy (MBE) may also be employed for the deposition of crystalline GaN. Initially, it was thought that MBE would allow for growth at lower temperatures than those required of CVD. In principle, this should allow for greater control over the growth process, thereby permitting the increased incorporation of nitrogen and reducing the unwanted high back- ground electron concentration typically found in GaN. In addition, MBE was viewed as a promising technique for the growth of nitride-based het- erostructures. Unfortunately, the promise that MBE o ered failed to meet some of these lofty expectations, probably as a consequence of low reaction rates. Nevertheless, MBE remains a popular means of depositing customized devices fabricated with GaN. Juza and Hahn [22] were the rst researchers to synthesize InN. They 3grew their samples using the same deposition procedure as that employed by Johnson et al. [10], i.e., through the interaction of ammonia gas with metallic indium at elevated temperatures. It was not until the 1960s, however, that research into InN was further pursued. Following the lead of GaN researchers, reactive sputtering was used in order to grow thin lms of InN. Researchers, such as Pastrnak and Souckova [23], Hovel and Cuomo [11], and Tansley and Foley [24] used this technique in order to explore the properties of InN. Later, CVD and MBE were used in order to fabricate high quality samples of crystalline InN [23, 25, 26]. A major breakthrough occurred in 2002, stemming from the much im- proved quality of InN lms grown through MBE. The bandgap of the re- sulting epitaxial lms was found to be 0.64 eV [27]; previous measurements had suggested a bandgap value of about 1.89 eV at room temperature for this material [28], although it should be noted that these previous measure- ments were performed on powdered samples of InN [11, 28, 29, 30]. These new results suggest that InN, and that alloys of InN with GaN, can be used to fabricate a family of optical devices that span over a broad swath of the electromagnetic spectrum, ranging from the near-infrared to the ultraviolet regions [7]. This nding ignited tremendous interest in the material proper- ties of the narrow-bandgap group-III nitride semiconductors. 4Unfortunately, a major obstacle for the further development of the ni- tride based semiconductor device technologies remains the absence of a lat- tice matched substrate for the epitaxial growth of these materials. Presently, substrates, such as sapphire (A12O3) and silicon carbide (SiC), are commonly used for the growth of high quality GaN lms; it should be noted, however, that there has been great recent interest in the growth of GaN-based ma- terials on silicon (Si) substrates, o ering not only low production costs but also attractive opportunities to incorporate GaN-based electron devices into Si-based technologies [31, 32]. One of the major challenges in depositing GaN-based materials on any of these substrates is the large lattice mismatch between the substrate and GaN-based materials, which inhibits high quality crystal growth [32]. Therefore, a bu er layer, such as AlN, is often employed in order to ameliorate this problem. Intensive studies have been carried out to improve matters. Recently, an enhancement of InGaN lm quality has been reported by Kurouchi et al. [33] with the insertion of a thin InN bu er layer. In another study, the lm quality of InGaN was found to be improved when zinc oxide (ZnO) is used as the substrate instead of sapphire [34]. Mismatches between the growing epilayer and the underlying substrate lead to large concentrations of threading dislocation lines, i.e., between 107 and 1011 cm 2 [35, 36], these threading dislocation lines extending from the 5substrate and/or bu er interface all the way into the deposited samples. In Figure 1.1, a typical transmission electron microscope image of the thread- ing dislocation lines within InN is shown, this gure being from Cimalla et al. [37]. Dislocations are known to a ect the electrical and the optical properties of semiconductors [38]. As a result, the study of the threading dislocation lines within group-III nitride semiconductors has been the fo- cus of some attention in recent years. Electrically, it is known that within GaN, electrons trapped by the threading dislocation lines diminish the over- all electron mobility [39, 40]. Optically, the trapping level associated with the dislocation defects can act as non-radiative recombination centers, and thus, are responsible for shaping features observed in the optical absorption and photoluminescence spectra [41]. In this thesis, the occupancy of the defect sites associated with the thread- ing dislocation lines within GaN and InN will be examined. The way this occupancy shapes the material properties of these semiconductors will also be probed. The focus will be primarily upon those properties that have a direct impact upon electron device performance. The case of AlN is not con- sidered as its energy gap is so wide that it e ectively acts as an insulator and is thus electrically inert. 6 Figure 1.1: A transmission electron microscope image of the threading dislo- cation lines present within epitaxially grown InN. This image is from Cimalla et al. [37]. 71.2 The crystalline structure of GaN and InN Depending on the epitaxial growth conditions, GaN, InN, and related alloys, have been successfully grown in both the zinc-blende and the wurtzite crystal structures. In Figures 1.2(a) and (b), the atomic con gurations associated with the zinc-blende and wurtzite crystal structures have been schematically depicted, respectively. The zinc-blende crystal structure can be thought of as being similar to the crystal structure associated with a dia- mond, with every other atom in the diamond crystal structure being replaced by an atom of the second element. The wurtzite crystal structure, however, is structurally closer to the graphite crystal structure. In order to distinguish the associated honey-comb meshes within the wurtzite crystal structure de- picted in Figure 1.2(b), this crystal structure has to viewed from the top. The majority of work on GaN-based materials has been focused on the wurtzite structure, as most III-V nitride semiconductors have been grown epitaxially on sapphire substrates, resulting in wurtzite based nitride lms [42, 43]. As a result, for the purposes of the analysis presented in this thesis, the focus will be on the electrical properties of the dislocation lines present within the wurtzite crystal structures associated with both GaN and InN. 8 (a) ( b ) Atom B Atom A Atom A Atom B Figure 1.2: (a) The zinc-blende and (b) wurtzite crystal structures. 91.3 Material parameters associated with wurtzite GaN and wurtzite InN In Table 1.1, the material parameters, corresponding to wurtzite GaN and wurtzite InN, are presented. The symbol, me, denotes the free electron mass. Table 1.1: The material parameters corresponding to wurtzite GaN and wurtzite InN. These material parameters are drawn from the references indi- cated within the square brackets. The bulk donor energy levels are provided with respect to the valence band maximum. Parameter Symbol GaN InN Electron e ective mass m e 0:2 me [44] 0:04 me [45] Energy gap "g (eV) 3:4 [44] 0:7 [46] Bulk donor energy level "d (eV) 3:39 [47] 0:69 [48] Relative dielectric constant r 8:9 [44] 15:3 [49] c-lattice constant c ( A) 5:1 [50] 5:71 [51] a-lattice constant a ( A) 3:16 [50] 3:53 [51] 101.4 Atomic structure of the edge type dis- locations within the wurtzite GaN and wurtzite InN crystal structures Dislocations are one dimensional defects embedded in the bulk of a crys- tal structure. Structurally, dislocations can be pure edge, pure screw, or of a mixed type. A schematic atomic con guration of such defects is depicted in Figure 1.3 [52]. Experimental results, presented in the literature, have shown that the majority of the dislocation lines within strained epitaxially grown samples of the nitride based semiconductors are pure edge in charac- ter [53, 54, 55, 56]. As a result, in this thesis, only edge type dislocations for the core structure of the dislocation lines will be considered. Several di erent atomic con gurations, for the core structure of edge type disloca- tions within the wurtzite crystal structure of GaN, have been suggested in the literature [57, 58]. These include the 4-ring, the 5/7-ring, the 8-ring, the gallium vacancy, and the nitrogen vacancy dislocation core structures. The 5/7-ring and the 8-ring core structures have alternatively been referred to as the open core and full core dislocations, respectively. The defect site associated with each of these dislocation core structures is depicted in Figure 1.4. Parts (a), (b), and (c) of this gure are after Takie 11 Figure 1.3: A section of a cubic crystal lattice, including a pure edge (A) and a pure screw dislocation (B). This gure is from Darling [52]. 12 (a) ( b ) ( c ) ( d ) ( e ) ( f ) Nitrogen Gallium Nitrogen Gallium Nitrogen Gallium Figure 1.4: Schematic atomic con gurations of the: (a) 4-ring, (b) 5/7- ring, (c) 8-ring, (d) gallium vacancy, (e) nitrogen vacancy, and (f) VGa ON dislocation core structures within wurtzite GaN. Parts (a), (b), and (c) are after Takie and Nakayama [59]. Parts (d) and (e) are after Lee et al. [57]. Part (f) is after You et al. [60]. 13and Nakayama [59], while parts (d) and (e) are after Lee et al. [57]. Part (f) is after You et al. [60]. The wurtzite lattice c-direction, along with the [0001] crystallographic growth direction, runs out of the page for all parts of Figure 1.4. Similar dislocation core structures have also been suggested for the threading dislocation lines within wurtzite InN [59, 61]. Of particular importance in this thesis, is another dislocation core struc- ture within the wurtzite crystal structure of GaN, depicted in Figure 1.4(f), and henceforth referred to as the VGa ON dislocation core structure. This dislocation core structure can be obtained from the gallium vacancy dislo- cation core structure through the replacement of the nitrogen atom at the center of the gallium vacancy dislocation core structure with an oxygen atom. In Section 4.6, the energy con guration, as well as the occupation statistics associated with the VGa ON dislocation core structure, will be discussed in greater detail. Regardless of the particular atomic con guration of the dislocation core structure, in this thesis, the smallest repetitive structural unit of the dislo- cation, along the direction of the dislocation line, will henceforth be referred to as a dislocation defect site. 141.5 Electrical properties of the threading dis- location lines within GaN and InN Dislocation lines within semiconductors were identi ed as possessing electrical properties associated with them as early as the 1950s [62, 63]. It was observed that after the introduction of dislocations through plastic de- formation into samples of n-doped germanium and silicon, the free electron concentration, as well as the overall electron mobility, decreased measurably within these semiconductor materials. A thorough review article on the elec- trical, as well as the optical properties of dislocations within germanium and silicon, exists in the literature [64]. Within GaN and InN, it has also been observed that the electron mobility is adversely a ected by the presence of the threading dislocation lines [39, 40, 65, 66, 67]. Studies suggest that the threading dislocation lines a ect the free electron concentration of the associated semiconductor material di er- ently within n-type GaN and n-type InN; within n-type GaN, experimental results, obtained from the electron holographic mapping of individual thread- ing dislocation lines, have shown the dislocation defect sites to be negatively charged [68, 69, 70, 71]. In other words, within GaN, the dislocation defect sites act to reduce the free electron concentration of the material by trap- 15ping some fraction of the free electrons from the conduction band. Within n-type InN, however, several investigations suggest that the threading dislo- cation lines contribute to the problem of the unintentional n-type nature of this semiconductor material [59, 61, 66, 72]. In other words, within InN, the dislocation defect sites donate electrons to the conduction band, thereby in- creasing the free electron concentration within this semiconductor material. The underlying reason for this di erence in the contribution of the threading dislocation lines to the free electron concentration within GaN and InN, is the di erence in the positioning of the dislocation defect energy level, "t, with respect to the conduction band minimum for each case. Within GaN, the dislocation defect energy level lies deep within the energy gap. In the case of InN, however, the dislocation defect energy level is suggested to lie above the conduction band minimum. Several approaches have been pursued in the literature in order to deter- mine the occupancy of the dislocation defect sites within n-type GaN. The earliest modeling of the charge build up along a dislocation line is due to Read [73]. The model of Read remains the most widely used model for the treatment of the occupancy of the dislocation lines within semiconductors to date. The model of Read [73] pictures a dislocation line within the diamond crystal structure of germanium as being comprised of a row of identical and 16equally spaced dangling bonds; see Figure 1.5. A revised electrostatic theory of charged dislocation lines within semi- conductors has been formulated in this thesis. In addition to the screening space charge distribution pro le, this procedure also gives a computational solution to the electrostatic potential surrounding the charged dislocation lines. The electrostatic potential thus obtained enables a determination of the electron mobility limited by scattering from the charged dislocation lines with greater certainty than that which has previously been achieved in the literature. In addition, the development of theoretical tools, that may be more easily applied to such problems, has been a focus of study. The nature of the screening mechanisms within these materials is also examined. Ulti- mately, fresh insights into the nature of these charged defects, and how they may in uence the performance of GaN and InN related electron devices, have been acquired. This thesis is organized in the following manner. In Chapter 2, a for- mulation for obtaining the most general screening space charge distribution surrounding a dislocation line within a semiconductor will be developed. Then this formulation will be applied to the speci c case of dislocation lines within GaN and InN. In particular, it will be shown that distinctively dif- ferent screening space charge distribution pro les exist around the threading 17 Figure 1.5: The atomic structure of a dislocation line within the diamond crystal structure. This gure is after Read [73]. The direction of the dislo- cation line is indicated with the dashed line in this gure. The symbol ? is frequently used in order to denote an edge type dislocation in the literature. 18dislocation lines within GaN and InN. In Chapter 3, the formulation for electron scattering from charged dislocation lines will be presented. Both the electron drift and Hall mobilities, limited by scattering from the thread- ing dislocation lines within GaN and InN, will be computed. In Chapter 4, the dislocation occupation statistics will be discussed. In particular, it is suggested that a Gibbs factor formalism can be used to obtain the occupa- tion statistics of the defect sites along threading dislocation lines within a semiconductor. Finally, conclusions are drawn in Chapter 5 and recommen- dations for future research are made. 19Chapter 2 Dislocation line charge screening 2.1 Modeling the charge build-up along the dislocation lines As was mentioned in the previous chapter, the defect sites along the threading dislocation lines within GaN and InN can be electrically active. Taking charge neutrality considerations into account, a negative or positive charge residing on the defect sites along threading dislocation lines should be compensated by a spatial charge distribution of the opposite polarity 20and equal in magnitude surrounding each threading dislocation line. This surrounding space charge distribution will henceforth be referred to as the screening space charge distribution. In this thesis, for computational simplicity, the screening space charge distribution will be assumed to be invariant along the direction of the dislo- cation lines. As a result, in cylindrical coordinates, a solely radial functional form will be assumed for the screening space charge distribution. This in- variance assumption is valid if (1) the dislocation lines are in nitely long, and (2) the discrete charge residing at the defect sites along these dislocation lines can be modeled with a uniform linear charge along the lines. In order to avoid the electrostatic potential surrounding the dislocation lines from becoming singular at r = 0, the core charge distribution will be considered to be spread uniformly in a cylinder of radius, Rdc, which will henceforth be referred to as the dislocation core radius. For the purposes of this analysis, the basal lattice constant, a, in the wurtzite crystal structures of GaN and InN will be taken as the dislocation core radius, Rdc, for each case. As a result, in this chapter, the charge build-up at each threading dislocation line will be modeled with two coaxial, in nitely long, space charge cylinders, as shown in Figure 2.1. The threading dislocation lines extend in the [0001] crystallographic direction of the wurtzite 21 s cR d cR r Figure 2.1: Modeling the space charge build-up around a threading disloca- tion line. Rdc denotes the dislocation core radius. Rsc represents the screening space charge cylinder radius. Rsc will be introduced subsequently. 22crystal structure for both GaN and InN. The question as to whether the defect sites along the dislocation lines within a given semiconductor material become predominantly positively or negatively charged depends upon (1) the positioning of the dislocation de- fect energy level, "t, within the energy band structure of the semiconductor material, and (2) the bulk doping conditions of the semiconductor material. For the purposes of this analysis, the primary focus is on n-type bulk dop- ing conditions of the semiconductor material; it is experimentally observed that most samples of GaN and InN are n-type doped in nature [66, 72]. It has been suggested that within GaN, the dislocation defect sites introduce a deep acceptor level within the energy gap of GaN. As a result, the dislocation defect sites within GaN have a natural tendency to accept electrons, thereby becoming negatively charged. This situation is schematically depicted in Fig- ure 2.2. Within InN, however, the dislocation defect energy level, "t, has been shown in the literature to lie above the conduction band minimum [59, 61]. As a result, the dislocation defect sites within InN have a natural tendency to donate their electrons to the conduction band and thereby become positively charged. This situation is schematically depicted in Figure 2.3. In this chapter, a quantitative analysis of the screening space charge dis- tribution surrounding a threading dislocation line is presented. The analysis 23 Radi al distance from the dislocation core line Energy 0q t cd Figure 2.2: Schematic diagram for the energy band bending around a thread- ing dislocation line within GaN. 24 Radi al distance from the dislocation core line En ergy 0q t cd Figure 2.3: Schematic diagram for the energy band bending around a thread- ing dislocation line within InN. 25is cast within the framework of the aforementioned radial spatial charge dis- tribution assumption. Following the presentation of a general formulation for this problem, the dislocation line charge screening, speci c to GaN and InN, are then computed. The analyzes presented in the subsequent chapters will be built upon the framework of the formulation developed in this chapter. This chapter is organized in the following manner. In Section 2.2, a general formulation for the screening of the dislocation lines within a semi- conductor material will be developed. Then, in Section 2.3, this formulation will be applied to the speci c case of n-type GaN in order to obtain the screening space charge distribution, as well as the associated electrostatic potential function, surrounding the threading dislocation lines within this semiconductor. In Section 2.4, the electrostatic potential obtained from this formulation will be compared with existing electron holographic experimen- tal data in the literature. In Section 2.5, the formulation will be applied to the speci c case of n-type InN. 262.2 Formulating the screening space charge distribution surrounding a dislocation line In modeling the space charge distribution associated with a dislocation line, it is useful to treat the charge accumulated on the defect sites along the core of the dislocation line with a uniform space charge of constant concen- tration, Ndc, con ned to an inner cylinder of radius Rdc, i.e., ndc(r) = ( Ndc; r Rdc 0; r > Rdc : (2.1) This space charge cylinder will henceforth be referred to as the dislocation core charge cylinder. Assuming that the dislocation defect sites are separated by the c-lattice constant associated with the relevant wurtzite crystal struc- ture, c, and denoting the charge donated per dislocation defect site along a dislocation line by f , the constant dislocation core space charge concen- tration, Ndc, in units of the fundamental electron charge, q, can thus be expressed as Ndc = f R2dcc : (2.2) Surrounding each dislocation core charge cylinder is a screening space charge cylinder of the opposite charge. For computational convenience, the radius of this space charge cylinder, Rsc, will be assumed to be nite. Un- der predominantly n-type bulk doping conditions, the contributions from the 27bulk acceptor and bulk hole concentrations to this screening space charge concentration can be neglected. Thus, the screening space charge concentra- tion, nsc(r), in units of the fundamental electron charge, q, may be written as nsc(r) = ( n(r) N+d (r); r Rsc 0; r > Rsc: : (2.3) In this equation, N+d (r) and n(r) represent the bulk ionized donor and free electron concentrations, respectively, at a radial distance, r, from the dislo- cation line core center. Poisson’s equation may be expressed as r2’ = r 0 ; (2.4) where r denotes the relative dielectric constant of the material, 0 represents the dielectric constant associated with free space, and is the space charge concentration. Within the framework of cylindrical coordinates, r2’ = 1 r @ @r r @’ @ r + 1 r2 @2’ @ 2 + @2’ @z2 ; (2.5) and thus, Eq. (2.4) can be written as 1 r @ @r r @’ @ r + 1 r2 @2’ @ 2 + @2’ @z2 = r 0 : (2.6) A cylindrical spatial charge distribution associated with a dislocation line implies a cylindrically symmetric electrostatic potential. That is, within the 28present framework of the dislocation line charge modeling, ’(r) does not depend on the and z coordinates. As a result, Eq. (2.6) further reduces to 1 r @ @r r @’(r) @ r = r 0 : (2.7) Noting that the space charge concentration associated with a given threading dislocation line may be expressed as the di erence between the core space charge concentration and the screening space charge concentration, i.e., = q(ndc nsc); (2.8) Eq. (2.7) can be integrated to obtain the resultant electrostatic potential ’(r) = q r 0 Z r 0 1 r0 Z r0 0 r00 [ndc(r 00) nsc(r 00)] dr00 dr0 + cint; (2.9) where cint is a constant of integration that will ultimately be determined through the requirement that ’(r Rsc) = 0. The e ect of the electrostatic potential, ’(r), in the neighbourhood of the dislocation core is to shift all of the electronic energy levels by the amount q’(r) relative to the Fermi energy level, "F. For the bulk ionized donor and free electron concentrations employed in Eq. (2.3), this can be achieved by expressing N+d (r) = Nd 1 + 2 exp [("F ("d q’(r))) =kT ] ; (2.10) 29and n(r) = Nc 2 p Z 1 0 p x 1 + exp [x ("F ("c q’(r)))=kT ] dx; (2.11) where "d and "c denote the bulk donor and conduction band minimum en- ergy levels, respectively. Nd represents the bulk donor concentration and Nc denotes the e ective density of states of the conduction band, i.e., Nc 2 m ckT 2 ~2 3=2 ; (2.12) m c , k, ~, and T being the electron e ective mass, Boltzmann’s constant, the reduced Planck’s constant, and the temperature, respectively. The position of the Fermi energy level, "F, can be obtained from the local charge neutrality condition at distances far from the dislocation lines, where the potential, ’(r), reduces to zero. That is Nc 2 p Z 1 0 p x 1 + exp [x ("F "c)=kT ] dx = Nd 1 + 2 exp [("F "d) =kT ] : (2.13) It is seen that Eq. (2.13) is simply equivalent to saying that nsc(r) = 0, which should be the case at distances far from the dislocation lines, i.e., for r Rsc. In the present analysis, the space charge concentration, nsc(r), and the electrostatic potential, ’(r), will be taken as the unknown variables of the for- mulation. In order to simultaneously solve for nsc(r) and ’(r) from Eqs. (2.3) 30and (2.9), an iterative approach will be proposed that involves nding suc- cessively more accurate solutions to the screening space charge distribution function, nsc(r). The details of this iterative approach, for the cases of GaN and InN, are provided in Appendices A and B, respectively. 2.3 Threading dislocation line charge screen- ing within GaN The screening space charge concentration pro le, nsc(r), and the cor- responding electrostatic potential function, ’(r), within GaN are plotted as a function of the radius, r, for ve di erent values of the negative charge donated per dislocation defect site, f , in Figures 2.4 and 2.5, respectively. The bulk donor concentration has been xed to 1018 cm 3 for all of the curves depicted in these gures. The GaN material parameters used in obtaining these results are speci ed in Table 1.1. The radius of the dislocation core charge cylinder has been set equal to the basal lattice constant of the wurtzite crystal structure of GaN, a, i.e., 3:16 A. The corresponding calculations are all performed at 300 K. The results presented in Figures 2.4 and 2.5 warrant the following obser- vations: 310 10 20 30 40 50 60 0 1 2 3 4 5 6 7 8 9 10 11 x 1017 r (nm) n sc (r) (c m− 3 ) f = 0.2 f = 0.5 f =1 f =1.5 f =2 Nd = 10 18 cm−3 T = 300 K GaN Figure 2.4: The screening space charge concentration of a threading disloca- tion line within GaN. The bulk donor concentration, Nd, is set to 1018 cm 3 and the temperature set to 300 K for the purposes of this analysis. 320 10 20 30 400 1 2 3 4 5 6 r (nm) Electrostatic Potential (V ) f = 0.2 f = 0.5 f =1 f =1.5 f =2 Nd = 10 18 cm−3 T = 300 K GaN Figure 2.5: The magnitude of the electrostatic potential surrounding a threading dislocation line within GaN. The bulk donor concentration, Nd, is set to 1018 cm 3 and the temperature set to 300 K for the purposes of this analysis. 33 it is observed that with increased values of f , more and more of the bulk donor atoms at larger radial distances from the dislocation core loose their donated electrons to the defect sites along the threading dislocation lines due to the requirement of charge neutrality around each threading dis- location line, the positive charge contained within the screening space charge distribution, nsc(r), is directly proportional to the value of f , i.e., f / Z Rsc 0 nsc(r)2 r dr; (2.14) for each curve depicted in Figure 2.4 close to the dislocation line, the screening space charge concentration reaches its maximum possible value of 1018 cm 3; however, as we start to move away from the dislocation line, it is observed that screening from free electrons gradually reduces the screening space charge con- centration to zero at large radial distances from the dislocation line. 2.4 Comparison with experiment O -axis electron holography is an interferometric technique that can measure the phase shift in an electron wave passing through an object. Since 34the phase shift in the electron wave function is sensitive to the local electro- static potential within the sample, electron holography can be used to provide quantitative information about potential variations with nanometer-scale res- olution. The electrostatic potential associated with a charged dislocation line within wurtzite GaN has been experimentally determined by Cherns et al. [69]; through the use of electron holography, Cherns et al. [69] experi- mentally determined the electrostatic potential associated with an individual threading dislocation line within n-type wurtzite GaN. Their experimental results are presented in Figure 2.6. For the bulk doping concentration found by Cherns et al. [69], i.e., Nd = 6 1017 cm 3, a satisfactory t to this ex- perimental data may be achieved through the present analysis by setting the negative charge donated per dislocation defect site to two, i.e., f = 2. This best t from the reults of the present analysis is also depicted in Figure 2.6. 2.5 Threading dislocation line charge screen- ing within InN Using the formulation developed in Section 2.2, the screening space charge concentration pro le, nsc(r), and the corresponding electrostatic po- tential function, ’(r), within InN are plotted as a function of the radius, r, 35−20 −15 −10 −5 0 5 10 15 20 −5 −4 −3 −2 −1 0 1 2 3 4 n−type GaN r (nm) Electrostatic potential (V ) n−type GaN T = 300 K Cherns et al. [69] Present analysis (Nd=6×10 17 cm−3) Figure 2.6: Comparison of the results from the present analysis with the experimental data of Cherns et al. [69]. For this analysis, f is set to 2 and the bulk doping concnetration, Nd, is set to 6 1017 cm 3, this being the value measured by Cherns et al. [69]. 36for di erent values of the negative charge donated per dislocation defect site, f , in Figures 2.7 and 2.8, respectively. Similar to the computations for the case of GaN, the bulk donor concentration has been xed to 1018 cm 3 for all of the curves depicted in these gures. The InN material parameters used for these computations are speci ed in Table 1.1. The radius of the dislo- cation core charge cylinder has been set equal to the basal lattice constant of the wurtzite crystal structure of InN, a, i.e., 3:53 A. The corresponding calculations are all performed at 300 K. The results presented in Figures 2.7 and 2.8 warrant the following obser- vations: it is observed that with increasing values of f , the screening space charge concentration increases rapidly at distances close to the dislo- cation line due to the requirement of charge neutrality around each of the thread- ing dislocation lines, the negative charge contained within the screening space charge distribution, nsc(r), is directly proportional to the value of f , i.e., f / Z Rsc 0 nsc(r)2 r dr; (2.15) for each curve depicted in Figure 2.7 370 2 4 6 8 10 12 14 16 18 20 0 2 4 6 8 10 12 x 1019 r (nm) n sc (c m− 3 ) f = 0.2 f = 0.5 f =1 f =1.5 f =2 Nd = 10 18 cm−3 T = 300 K InN Figure 2.7: Screening space charge concentration of a threading dislocation line within InN. The bulk donor concentration, Nd, is set to 1018 cm 3 and the temperature set to 300 K for the purposes of this analysis. 380 2 4 6 8 10 12 14 16 18 20 0 0.5 1 1.5 2 2.5 r (nm) Electrostatic Potential (V ) f = 0.2 f = 0.5 f =1 f =1.5 f =2 Nd = 10 18 cm−3 T = 300 K InN Figure 2.8: The magnitude of the electrostatic potential surrounding a threading dislocation line within InN. The bulk donor concentration, Nd, is set to 1018 cm 3 and the temperature set to 300 K for the purposes of this analysis. 39 unlike what was seen in the screening of the dislocation core charge within GaN, within InN the screening close to the dislocation line is provided predominantly by free carriers. This free carrier screening constitutes a much more e ective screening of the dislocation core charge when compared with the screening of the dislo- cation core charge o ered by the ionized donors present within GaN. As a result, it is seen that for the same value of f , the values for the computed electrostatic potential function, ’(r), are much lower within InN than for the case of GaN; compare Figures 2.5 and 2.8. 40Chapter 3 Electron mobility limited by scattering from charged dislocation lines 3.1 The electron drift and Hall mobilities in semiconductors The electron mobility characterizes how quickly an electron can move through a material under the action of an applied electric eld. When an electric eld, E , is applied across a piece of a semiconductor material, the 41electrons respond by moving with an average velocity, vdrift, referred to as the electron drift velocity. Within a given semiconductor material, the electron mobility depends on the magnitude of the electric eld, particularly at high elds, when velocity saturation occurs. In the present thesis, however, the focus is on low- eld electron mobilities. At low applied electric eld strengths, the electron drift velocity is linearly proportional to the applied electric eld strength. That is vdrift = driftE ; (3.1) where drift denotes the electron drift mobility and the negative sign re ects the fact that electrons are pulled in the opposite direction of the applied elec- tric eld. Electron mobilities are almost always speci ed in units of cm2=Vs. The electron mobility arises as a consequence of imperfections within a semiconductor material and as a result of the interaction of transiting elec- trons with the atoms within the material, which naturally engage in vibra- tions as a result of thermal processes. Each material system possesses crystal defects, vacancies, dislocations, impurities, and other features which will act to scatter transiting electrons. At elevated temperatures, due to their thermal energy, the atoms within the crystal structure of a semiconductor material will vibrate about their equilibrium lattice sites. The transiting electrons will collide with such vibrating atoms, and, as a result, will be scattered. 42As with the molecules of a gas within a container, in the absence of an applied electric eld, the conducting electrons within a semiconductor material move in random directions, with a wide range of di erent energies. When an electric eld is applied, a net average velocity is obtained, this velocity being in the opposite direction of the applied electric eld. The average time it takes for an electron within a semiconductor material to move between successive scattering events is denoted by h i. The bracket notation in this expression, h i, is used to denote an averaging over all possible energies of the incident electron. In the literature, the physical quantity, h i, is often referred to as the mean free time, or the relaxation time, of the electrons. Under the action of an applied electric eld, E , electrons will accelerate at a rate equal to qE=m e between the successive scattering events. Noting that the force exerted on an electron by the periodic net positive charge of the lattice atoms is already implicitly contained within the electron e ective mass term, m e, the electron drift velocity can thus be expressed as vdrift = qE m e h i : (3.2) That is, on average, an electron’s velocity in the opposite direction of the applied electric eld is proportional to the product of the applied electric eld strength and the average time between scattering events, and inversely proportional to the electron e ective mass. A comparison of Eq. (3.2) with 43Eq. (3.1), gives the following expression for the electron drift mobility drift = q h i m e : (3.3) Empirically, electron mobility is often measured through the use of a Hall e ect experiment. In the Hall e ect experiment, in addition to the applied electric eld, E , the transport of the conduction electrons is a ected by the presence of a weak magnetic eld, B, applied in a direction perpendicular to that of the applied electric eld. As a result, the experimentally measured electron Hall mobility di ers slightly from the electron drift mobility de ned earlier. Speci cally, it can be shown that the electron Hall mobility is related to the electron drift mobility through the following relationship Hall = h 2i h i2 drift; (3.4) where, once again, the bracket notation in the expressions for h 2i and h i, represents an averaging of the relaxation time squared and the relaxation time over all of the possible energies of the incident wave vectors [74]. For multiple scattering mechanisms, the e ective mean free time is derived from the individual mean free times of the representative scattering events, i.e., 1 h i = 1 h disi + 1 h impi + 1 h piezoi + 1 h latticei + :::; (3.5) 44where h disi, h impi, h piezoi, and h latticei denote the dislocation scattering limited mean free time, the impurity scattering limited mean free time, the piezoelectric scattering limited mean free time, and the lattice vibration scattering limited mean free time, respectively. Multiplying both sides of Eq. (3.5) by m e=q, one obtains the celebrated Matthiessen’s rule, relating the overall electron mobility, , to the partial electron mobilities, i.e., 1 = 1 dis + 1 imp + 1 piezo + 1 lattice + :::; (3.6) where dis, imp, piezo, and lattice denote the dislocation limited, impurity limited, piezoelectric limited, and lattice vibration limited electron mobilities, respectively. It should be mentioned that Matthiessen’s rule is an approx- imation and is not universally valid. This rule is not valid if the factors a ecting the mobility depend on each other, because individual scattering probabilities cannot be summed unless they are independent of each other. Among the di erent electron scattering mechanisms present within GaN and InN, in this thesis, only the electron drift and Hall mobility limited by scattering from the dislocation lines will be examined. Speci cally, in the present work, the dependence of the electron drift and Hall mobilities on the threading dislocation line density, and the amount of charge accumulated along the threading dislocation lines within GaN and InN, will be computed separately. 45This chapter is organized in the following manner. In Section 3.2, the formulation corresponding to the scattering of electrons from the dislocation lines within semiconductors will be presented. Results for the electron drift and Hall mobilities, limited by scattering from the threading dislocation lines within GaN and InN, are then presented in Sections 3.3 and 3.4, respectively. 3.2 Formulation The formulation presented in this section follows closely the formula- tion developed in a paper by You et al. [60]. Similar to the other scattering mechanisms present within semiconductors, scattering from the dislocation lines should be treated quantum mechanically. In Figure 3.1(a), a repre- sentative scattering event from a dislocation line is depicted. In this gure, it is assumed that an electron, of a plane wave functional dependence, i.e., (~k; ~r) = 1p V exp(i~k ~r), is incident upon the dislocation line and is scattered in a direction indicated by the wave vector, ~k0. The incident and the scattered wave vectors, i.e., ~k and ~k0, can be decomposed into components parallel and perpendicular to the dislocation line direction. The component parallel to the dislocation line direction, denoted by ~kk and ~k0k, will not change during the scattering process, i.e., ~kk = ~k0k. For the perpendicular components of 46 ?k ? ?k ? k ? ?k ?? (a) ( b ) Figure 3.1: (a) Scattering of a plane wave with a wave vector ~k?, from a cylindrically symmetric potential, and (b) de ning the wave vector di erence, ~ ?. 47the incident and scattered wave vectors, k? = k0?, i.e., elastic scattering will be assumed. In the Dirac bracket notation, the matrix element of the scattering po- tential in the incident and scattered wave functions may be expressed as M( ?) = h (~k; ~r) j’(~r)j (~k 0; ~r)i; (3.7) where ? is the magnitude of the wave vector di erence, ~ ?; see Figure 3.1(b). The Dirac bracket notation can be expanded in terms of the associated inte- gral over the volume of the semiconductor material, V , to give M( ?) = 1 V Z V exp(i~ ~r)’(~r) d3~r: (3.8) Since the scattering potential, ’(~r), has a solely radial dependence in cylindri- cal coordinates, this integration in cylindrical coordinates may be simpli ed, leading to M( ?) = 1 V Z H 0 Z 2 0 Z R 0 exp(i ?r cos( ))’(r)r dz d dr: (3.9) Performing the integrations over the z and coordinates of the cylindrical coordinate system, the matrix element may thus be expressed as M( ?) = 2 S Z R 0 J0( ?r)’(r)r dr; (3.10) where S is the surface area of the semiconductor material and J0( ?r) denotes a Bessel function of the rst kind of order 0, evaluated at ?r. Ideally, the 48integration in Eq. (3.10) should be computed for the radius, R, approaching in nity. In practice, however, since it is assumed that ’(r) = 0 for r > Rsc, one need only take the radius, R, to be equal to Rsc. In quantum physics, Fermi’s golden rule is the way to calculate the transi- tion rate (probability of transition per unit time) from one energy eigenstate of a quantum system into a continuum of energy eigenstates, due to a per- turbation. According to Fermi’s golden rule, the probability of the incident electron being scattered in the direction denoted by the wave vector, ~k0?, is proportional to the square modulus of the scattering matrix element, i.e., jM( ?)j 2. The scattering rate is thus given by W (~k; ~k0) = 2 ~ jM( ?)j 2 (Ek Ek0); (3.11) where the in the (Ek Ek0) term on the right hand side of Eq. (3.11), denotes the Dirac delta function, not to be confused with the ? de ned earlier as the di erence between the incident and scattered wave vectors in the perpendicular direction. The average time that it takes for an electron to be successively scattered between the dislocation lines depends on the energy of the electron. For an electron characterized by a plane wave with a wave vector, ~k, the kinetic 49energy of the electron can be written as E(k) = p2 2m e : (3.12) Substituting for the electron crystal momentum, p = ~k, into Eq. (3.12), the energy of the electron becomes related to its wave vector through the following relation E(k) = ~2k2 2m e : (3.13) It can be shown that the relaxation time for the elastic scattering of an electron with incident wave vector, ~k, is related to the scattering rate W (~k; ~k0), through the following relationship [74] 1 (k) = Z W (~k; ~k0) " 1 ~v~k0 ~E ~v~k ~E # d3k0 (2 )3 ; (3.14) where ~v~k and ~v~k0 denote the incident and the scattered electron velocities, respectively. As a result, in the present case of scattering from a cylindrical potential function, and for an electric eld, ~E , applied in a direction per- pendicular to the direction in which the threading dislocation lines extend, the relaxation time from one scatterer in a volume V , will be given by the following relationship 1 dis(k) = m e 2 ~3k3 Z 2k 0 jM( ?)j 2 ? 2 p 1 ( ?=2k)2 d ?: (3.15) In Eq. (3.16), the integration runs over all possible values of the wave vector di erence, ?, i.e., it runs from the lower limit, 0, corresponding to negligible 50scattering, to the higher limit, 2k, corresponding to the complete backward scattering of the incident plane wave. If the density of the threading dislocation lines is dis, then there are disS scatterers in the volume. It will be assumed that these dislocation lines are randomly placed and are well separated from each other so that the scattering is incoherent. In this case, according to quantum mechanics, one can simply sum the scattering rate from each scatterer. The total relaxation time is thus 1 dis(k) = dism e 2 ~3k3 Z 2k 0 ~M( ?) 2 ? 2 p 1 ( ?=2k)2 d ?; (3.16) where ~M( ?) SM( ?). The energy selective mean free time, dis(E), has to be integrated over all the energies of the incident wave vector in order to obtain the energy averaged mean free time, i.e., h disi = 4 3 (1=2)(kBT )5=2 Z 1 0 dis(E)E 3=2 exp( E=kBT ) dE: (3.17) In Eq. (3.17), the term E3=2 exp( E=kBT ) in the integrand corresponds to the distribution of energies of the electrons within a semiconductor, which essentially acts as a weighting factor for the integration. Similarly, the mean free time squared energy averaged, can be obtained from the following relationship [60] 2dis = 4 3 (1=2)(kBT )5=2 Z 1 0 2dis(E)E 3=2 exp( E=kBT ) dE: (3.18) 51Finally, substitution of the expressions for h disi and h 2disi, from Eqs. (3.17) and (3.18), into Eqs. (3.3) and (3.4), gives the associated dislocation scatter- ing limited electron drift and Hall mobilities, respectively. 3.3 Electron mobility limited by scattering from the charged threading dislocation lines within GaN Substituting the electrostatic potential pro les obtained earlier in Fig- ure 2.5 for ’(r), in the formulation developed in Section 3.2, the electron drift and Hall mobilities within GaN can be computed numerically. The re- sultant dislocation limited electron drift and Hall mobilities are plotted as a function of the value of the negative charge accumulated along the disloca- tion lines, f , in Figures 3.2(a) and (b), respectively. The lled blue circles in these gures correspond to the result of these computations. The same GaN material parameters from Table 1.1, used earlier in obtaining the elec- trostatic potentials of Figure 2.5, have also been utilized for the purposes of the electron mobility computations in this section. Furthermore, since the electrostatic potentials in Figure 2.5 were obtained at room temperature 520 0.5 1 1.5 2 101 102 103 104 f Dislocation limited electron drift mobility (c m2 /V sec ) (a) Nd = 1018 (cm−3) σdis = 10 10 (cm−2) T = 300 K Numerical Computation µdriftdis ∝ 1/f 2.51 µdriftdis ∝ 1/f 2 0 0.5 1 1.5 2 101 102 103 104 f Dislocation limited electron Hall mobility (c m2 /V sec ) (b) Nd = 1018 (cm−3) σdis = 10 10 (cm−2) T = 300 K Numerical Computation µHalldis ∝ 1/f 2.51 µHalldis ∝ 1/f 2 Figure 3.2: (a) The electron drift mobility and (b) the electron Hall mobility limited by scattering from charged threading dislocation lines within GaN. 53and for a bulk donor concentration of Nd = 1018 cm 3, the same values of T = 300 K and Nd = 1018 cm 3 have also been assumed in the computation of the electron mobilities in this section. The electron mobilities presented in Figures 3.2(a) and (b) have been plotted for a dislocation line density of dis = 1010 cm 2. Since, within the formulation developed in Section 3.2, the electron mobilities are inversely proportional to the dislocation line den- sity, dis (recall Eq. (3.16)), deducing the values for the associated electron drift and Hall mobilities at other dislocation line densities, from the results depicted in Figures 3.2(a) and (b) would be straightforward. If the electrostatic potentials depicted in Figure 2.5 were linearly propor- tional to the value of f , Eqs. (3.10) and (3.16) suggest that the computed electron drift and Hall mobilities would have a 1=f 2 functional dependence. This 1=f 2 functional dependence of the electron mobility is depicted with the dashed lines in Figures 3.2(a) and (b). However, from Figures 3.2(a) and (b), it is observed that a best t to the computed electron mobility results is obtained for a 1=f 2:51 functional dependence. Since the negative charge in the dislocation core cylinder scales linearly with the value of f , this overall dependence of the electron mobility on the value of f suggests that the bulk ionized donor screening occurring within GaN is less e ective than a linear screening with increasing values of f . 543.4 Electron mobility limited by scattering from the charged threading dislocation lines within InN In Figures 3.3(a) and (b), the dislocation limited electron drift and Hall mobilites are plotted as a function of the positive charge accumulated along the dislocation lines within InN, f . The electrostatic potentials obtained earlier in Figure 2.8, along with the InN material parameters given in Ta- ble 1.1, were used for the computation of these dislocation limited electron mobility results within InN. Since the electrostatic potentials in Figure 2.8 were obtained at room temperature and for a bulk donor concentration of Nd = 1018 cm 3, the same values of T = 300 K and Nd = 1018 cm 3 have also been assumed in the computation of the electron mobilities in this sec- tion. Furthermore, the electron mobilities depicted in Figures 3.3(a) and (b) have been plotted for a dislocation line density of dis = 1010 cm 2. Within the framework of the formulation of Section 2.2, utilized in obtaining the electron mobility results in this section, determining the dislocation limited electron mobilities for other threading dislocation line densities would be straightforward. For the case of InN, from Figures 3.3(a) and (b) it is observed that a best 550 0.5 1 1.5 2 102 103 104 105 f Dislocation limited electron drift mobility (c m2 /V sec ) (a) Nd = 1018 (cm−3) σdis = 10 10 (cm−2) T = 300 K Numerical Computation µdriftdis ∝ 1/f 1.7 µdriftdis ∝ 1/f 2 0 0.5 1 1.5 2 102 103 104 105 f Dislocation limited electron Hall mobility (c m2 /V sec ) (b) Nd = 1018 (cm−3) σdis = 10 10 (cm−2) T = 300 K Numerical Computation µHalldis ∝ 1/f 1.7 µHalldis ∝ 1/f 2 Figure 3.3: (a) The electron drift mobility and (b) the electron Hall mobility limited by scattering from charged threading dislocation lines within InN. 56 t to the computed electron mobility results is obtained for a 1=f 1:7 functional dependence. Since the positive charge in the dislocation core cylinder scales linearly with the value of f , this overall dependence of the electron mobilities on the value of f suggests that the predominantly free electron screening of the dislocation lines within InN is more e ective than a linear screening with increasing values of f . Furthermore, in terms of the absolute electron mobility values, under the same n-type bulk doping conditions, it is observed that the electron drift and Hall mobilities within InN are much higher than those computed for the case of GaN. This higher electron mobility within InN can be attributed to the much stronger screening of the dislocation core charge within InN as well as the lower value for the electron e ective mass of InN. 57Chapter 4 Occupation statistics of the dislocation defect sites 4.1 Di erent approaches for the determina- tion of the dislocation defect site occupa- tion statistics In the previous two chapters, computations of the screening space charge distribution functions, the scattering potentials, and the electron mobilities, were reported for ve di erent representative values of f , for the cases of both GaN and InN. In general, however, the value of f can be kept as a 58variable of the formulation and determined as a function of the bulk donor concentration. In fact, determining the value of f as a function of the free carrier concentration and/or bulk doping concentration within n-type GaN has been the subject of several studies in the literature [39, 60, 75, 76, 77]. In this chapter, theoretical means of determining the value of f from knowledge of the free carrier concentration and/or the bulk doping concen- tration within GaN and InN is developed. Initially, a survey of theoretical formulations that are available in the literature is provided. Later, a simple Gibbs factor formalism will be introduced, which allows for a determination of the occupation statistics of multiple ionization states associated with a dislocation defect site. Thus far, the treatment of multiple ionization states associated with a dislocation defect site has only been achieved in the litera- ture through the use of either detailed numerical density functional calcula- tions [76] or through tedious energy minimization computations [60, 77]. The formulation provided herein, represents a break from this tradition, allowing for the presentation of intuitively appealing results. Understanding the occupation statistics of the defect sites associated with the dislocation lines within semiconductor materials has been the focus of much attention over the years. The model of Read [73], developed in 1954, remains the dominant approach used for the characterization of the occupa- 59tion of the defect sites associated with the threading dislocation lines within semiconductor materials. The model of Read [73] pictures a dislocation line within the crystal structure (as was initially applied to germanium) as being comprised of a column of identical and equally spaced dangling bonds, as was shown in Figure 1.5. Each dangling bond can potentially trap an extra electron from a nearby donor atom. The model of Read [73] assumes that the screening of the charge trapped by the dangling bonds along a dislocation line is provided for by the complete ionization of the bulk donor atoms in the immediate vicinity of the dislocation line, up to a certain radius, referred to as the Read radius, R. The Read radius can be readily obtained in terms of the charge trapped along the dislocation lines per dislocation defect site, f , and the bulk donor concentration, Nd, i.e., R = r f cNd : (4.1) As some of the dangling bonds trap an extra electron, an electrostatic po- tential begins to build up, preventing the remaining dangling bonds from becoming negatively charged. Read [73] determined the fraction, f , of dan- gling bonds that would actually trap electrons by minimizing the electrostatic energy within the space charge region formed about each dislocation line. In 1998, Weimann et al. [39] employed the model of Read [73] as the basis for obtaining the occupation statistics of the defect sites associated with the 60threading dislocation lines within GaN. For the purposes of their computa- tions, Weimann et al. [39] assumed a simple dangling bond model for the dislocation core structure. Gurusinghe and Andersson [75] also followed the general approach presented by Read [73]. However, instead of minimizing the electrostatic potential energy within the space charge cylinder, they min- imized the corresponding Helmholtz free energy. You et al. [60] also followed the energy minimization approach of Read [73]. However, instead of pictur- ing each dislocation line as a simple column of dangling bonds, they assumed a realistic VGa ON core structure for the threading dislocation lines. Le- ung et al. [76] obtained the occupation statistics of the threading dislocation lines within GaN for a number of di erent dislocation core structures, for both p-type and n-type doping conditions. Leung et al. [76] used a simu- lated annealing approach, which is essentially a numerical density functional calculation. Results obtained through the application of the proposed Gibbs factor formulation presented in this thesis to the problem of determining the occupation statistics of threading dislocation lines within GaN will be compared with these other approaches in this chapter. This chapter is organized in the following manner. In Section 4.2, the general Gibbs factor formalism, and its formulation within semiconductor materials, will be presented. In Section 4.3, the electrostatic potential at the 61site of a dislocation defect, resulting from the charge on the other dislocation defect sites along the dislocation line, will be obtained. In Section 4.4, simple analytical expressions will be proposed to approximate the actual screening space charge distributions of the dislocation lines within GaN and InN, sep- arately. In Section 4.5, the proposed Gibbs factor formalism will be applied to the speci c case of dislocation lines within GaN with an assumed dangling bond dislocation core structure. In Section 4.6, the Gibbs factor formalism will be applied to the problem of determining the occupation statistics of the VGa ON dislocation defect sites within GaN. In Section 4.7, a 5/7-ring core structure will be assumed for the determination of the occupation statistics of the dislocation defect sites within wurtzite InN. 4.2 Gibbs factor formalism The probability of a quantum mechanical system, in thermal and dif- fusive equilibrium with a reservoir, to exist in each of its di erent ionization states is proportional to the corresponding Gibbs factor of the system in that particular state. The absolute probability that the system will be found in a state with N1 electrons, with an associated total energy of "1, is given by the Gibbs factor associated with that state, divided by the grand partition 62function, i.e., P (N1; "1) = exp [(N1 "1) =kBT ] Z ; (4.2) where denotes the chemical potential of the system/reservoir. The grand partition function, Z, is simply the summation of the Gibbs factors associated with the di erent possible states of the quantum mechanical system, i.e., Z("F; T ) = 1X Ne=0 X s(Ne) exp Ne "s(Ne) =kBT : (4.3) In Eq. (4.3), the rst summation runs over the number of electrons within the quantum mechanical system. Even though in Eq. (4.3) this summation runs from zero to in nity, in practice, only a few ionization states of the system are most probable and thus, to a good approximation, this summation need only be run over two or three electrons. The second summation, i.e., the summation over s(Ne), is a summation that runs over all the di erent possible quantum mechanical states of the system when it consists of Ne electrons. In the context of semiconductor physics, the chemical potential associ- ated with the electrons, , is referred to as the Fermi energy level, "F. Fur- thermore, for the case of the dislocation lines within semiconductors, each dislocation defect site serves the purpose of the quantum mechanical sys- tem, and the bulk of the semiconductor material serves the purpose of the reservoir with which the dislocation defect site should remain in thermal and 63di usive equilibrium. In general, it is possible that for some dislocation core structures, the associated defect sites introduce a multitude of energy levels into the energy band structure of the semiconductor material. For example, in Figure 4.1, the multitude of energy levels that the full core, the open core, the gallium vacancy, and the nitrogen vacancy dislocation core structures in- troduce into the energy band of GaN are depicted. However, under any given bulk doping condition, most of those energy levels remain either completely occupied or completely unoccupied, and thus, need not be included in the evaluation of the grand partition function. As a result, in this thesis, it will be assumed that each dislocation defect site introduces a single energy level, "t, into the energy gap of the semiconductor. This energy level can exist in any of the four di erent occupation states depicted in Figure 4.2. According to the Gibbs factor formalism, the probability of the dislocation defect site to exist in each of these states is given by fs(0) = 1 Z ; (4.4) fs(1") = fs(1#) = fs(1) 2 = exp "s(1) "F =kBT Z ; (4.5) fs(2) = exp "s(2) 2"F =kBT Z ; (4.6) with Z = 1 + 2 exp "s(1) "F =kBT + exp "s(2) 2"F =kBT ; (4.7) 64 Figure 4.1: The schematic diagram of energy levels in the band gap for four di erent dislocation core structures within GaN. This gure is after Lee et al. [57]. 65 𝜀𝑡 s ( 0 ) s ( 1 ) s ( 1 ) s ( 2 ) Figure 4.2: The four possible occupation states of a dislocation defect site. 66"s(1) and "s(2) representing the total energy of the dislocation defect site in the states labelled by s(1) and s(2) in Figure 4.2, respectively. In order to relate the energies "s(1) and "s(2) to the dislocation defect energy level, "t, it is recalled, from Section 2.2, that the e ect of the elec- trostatic potential, ’(r), in the neighbourhood of the dislocation core is to shift all of the electronic energy levels by the amount q’(r) relative to the Fermi energy level, "F. For the single electron occupancy of a dislocation defect site, i.e., in either of the s(1 ") or s(1 #) states, the energy of the single electron residing in the dislocation defect site will thus be shifted by the value of the electrostatic potential at the site of the dislocation defect, ’0, multiplied by the charge of the electron, i.e., "s(1) = "t q’0: (4.8) If the electron-electron repulsion of the two electrons in the state s(2) is neglected, the energy of the double occupancy of the dislocation defect site will simply be twice that of the single occupancy of the dislocation defect site, i.e., "s(2) = 2 "s(1): (4.9) In general, however, this electron-electron repulsion causes the dislocation defect energy level, "t, to split into two energies, " (1) t and " (2) t , much like 67the rst and second ionization energies of a helium atom. In the case of a dislocation defect site in the presence of inner electron-electron interactions, "(1)t will be the energy of the single electron in either of the states s(1 ") or s(1 #). Due to the electrostatic repulsion from the rst existing electron in the state s(2), the energy of the second electron in the state s(2), "(2)t , will be higher than that of the single electron in the state s(1), i.e., "(2)t > " (1) t . In the presence of electron-electron repulsions within each dislocation defect site, the total energies of the defect site in its di erent occupation states should thus be written as "s(1) = " (1) t q’0; (4.10) and "s(2) = "s(1) + "(2)t q’0 : (4.11) The value of the electrostatic potential at the site of a dislocation defect, ’o, is due to the presence of (1) the charge accumulated on the other defect sites situated along the dislocation line, and (2) the screening space charge distribution. The contribution from these two charge build-ups will hence- forth be denoted by ’odc and ’osc , respectively. As a result, the electrostatic potential at the site of a dislocation defect can be expressed as ’o = ’odc + ’osc : (4.12) 68Taking the discreteness of the charge residing on the dislocation defect sites into account, in the following section, an approximate expression for ’odc will be derived. Then, in Section 4.4, simple analytical expressions will be obtained for the electrostatic potential, ’osc , by virtue of approximating the actual screening space charge distributions, obtained earlier in Figures 2.4 and 2.7, by some alternative simple analytical functions. 4.3 Approximating the electrostatic potential induced by the dislocation core charge In this section, the aim is to obtain an approximate analytical expression for the electrostatic potential, ’odc . Speci cally, the discreteness of the charge on the defect sites along the dislocation line will be taken into account and the dislocation line will be assumed to have a nite height, h. Assuming the distance between the successive dislocation defect sites to be equal to the c- lattice constant associated with the wurtzite crystal structure of GaN/InN, a total of h=c defect sites exist on the dislocation line. Consequently, the number of defect sites residing in the four dislocation defect occupation states s(0), s(1 #), s(1 "), and s(2), can be obtained from the following relations Ns(0) = hfs(0) c ; (4.13) 69Ns(1#) = Ns(1") = Ns(1) 2 = hfs(1#) c = hfs(1") c = hfs(1) 2c ; (4.14) and Ns(2) = hfs(2) c ; (4.15) respectively. Furthermore, it will be assumed that, in units of the funda- mental electron charge, the dislocation defect site in each of the four states s(0), s(1 #), s(1 "), and s(2), has a net charge of is(0), is(1#), is(1"), and is(2) associated with it. The average charge per dislocation defect site, in units of the fundamental electron charge, can thus be obtained from the following relation f = is(0) Ns(0) + 2 is(1#) Ns(1#) + is(2) Ns(2) : (4.16) Recall that in Chapters 2 and 3, ve di erent representative values for the average charge per dislocation defect site, f , were considered, i.e., 0.2, 0.5, 1, 1.5, and 2. In the context of the Gibbs factor formalism developed in this chapter, however, from Eq. (4.16) and Eqs. (4.4) to (4.6), it is now seen how the average charge per dislocation defect site, f , can be related to the value of the electrostatic potential at the site of a dislocation defect. For a dislocation defect site situated in the middle of a dislocation line of height, h, one half of the other dislocation defect sites lie above this defect site and one half lie below this dislocation defect site. As a result, taking the 70discreteness of the charge on the other dislocation defect sites into account, the electrostatic potential due to the other defect sites along the dislocation line, at the site of a defect in the middle of the dislocation line, can be expressed as ’odc = 2 "0fs(0)is(0) q Ns(0)=2X n=1 1 n + 2 "0fs(1)is(1) q Ns(1)=2X n=1 1 n + 2 "0fs(2)is(2) q Ns(2)=2X n=1 1 n ; (4.17) where "0 q2 4 r oc : (4.18) In order to express the result of the summation, PN n=1 1 n , in a closed analyt- ical form, Read [73] employed the approximation, PN n=1 1 n ln (N) + 0:577, which is valid for large values of N . In the present thesis, however, in order to make a corresponding approximation, which will be valid for small values of N as well, the following series expansion approximation will be made NX n=1 1 n ln (N + 1) + 0:577; (4.19) where the value, 0.577, is an approximation to Euler’s constant. In Fig- ure 4.3, the value of the two expressions on the left and right hand side of Eq. (4.19) have been plotted as a function of N , with the red and blue curves, respectively. Substituting the approximate expression for PN n=1 1 n , i.e., Eq. (4.19), into 710 20 40 60 80 100 0 1 2 3 4 5 6 N Summation / Approximate analytical expressio n Σ n=1 N 1/n ln(N+1)+0.557 Figure 4.3: Comparison of the exact and the suggested approximation to the value of the summation, PN n=1 1 n , as a function of N . 72Eq. (4.17), the electrostatic potential, ’odc , can thus be nally written as ’odc = 2 "0 q fs(0)is(0) ln Ns(0)=2 + 1 + 0:577 + fs(1)is(1) ln Ns(1)=2 + 1 + 0:577 + fs(2)is(2) ln Ns(2)=2 + 1 + 0:577 g: (4.20) 4.4 Approximating the screening space charge distribution function In general, one would like to take the value of the charge per disloca- tion defect site, f , as a variable, in addition to the two existing variables, ns(r) and ’(r), and solve simultaneously for the three variables. However, simple analytical expressions can only be obtained for the value of the elec- trostatic potential, if the screening space charge distribution function, nsc(r), is approximated with a simple analytical functional form. In the case of the predominantly ionized bulk donor screening, the screen- ing space charge distribution function may be approximated by a step func- tion. In Figure 4.4, both the exact and the suggested approximation to the screening space charge distribution functions have been depicted for the rep- resentative case of GaN. The red curve in Figure 4.4 is identical to the exact 730 20 40 600 2 4 6 8 10 x 1017 r (nm) n sc (r) (c m− 3 ) Nd = 10 18 cm−3 T = 300 K f = 1 R = 24.98 nm Figure 4.4: The approximation of the screening space charge distribution function surrounding a threading dislocation line within GaN by a step func- tion. 74screening space charge distribution result for the value of f = 1, obtained earlier in Figure 2.4. In the suggested approximate screening space charge distribution function depicted by the blue curve in Figure 4.4, the radius at which the screening space charge distribution drops o to zero, has been set to the Read radius, R, i.e., it is assumed that nsc(r) ( Nd; r R 0; r > R : (4.21) When determining the radial functional form of the electrostatic potential resulting from such a step function space charge distribution, Read [73] as- sumed the corresponding space charge cylinder to be of a nite height, h, and determined the corresponding electrostatic potential function, at mid- point along the height of the cylinder to be ’sc(r) = f "0 q 1 + 2 ln h R r2 R2 ; (4.22) where f refers to the average negative charge per dislocation defect site. Thus, at the site of a dislocation defect, i.e., at a radial distance r = 0 from the dislocation line ’osc = f "0 q 1 + 2 ln h R : (4.23) For the case of InN, the predominantly free electron screening will be approximated by an exponential function, i.e., nsc(r) n0 exp( r= ): (4.24) 75The parameter n0 in Eq. (4.24) can be obtained through the requirement that the total negative charge subtended by the exponential function equals the positive charge of f per dislocation defect site along the dislocation line, i.e., n0 = 2f "0 q 2 : (4.25) In Figure 4.5, both the exact and the suggested approximation to the screen- ing space charge distribution functions have been depicted for the representa- tive case of InN. The red curve in Figure 4.5 is identical to the exact screening space charge distribution result for the value of f = 1, obtained earlier in Figure 2.7. It is observed that a best exponential functional t to the exact screening space charge distribution in Figure 4.5 is obtained for being set to 2:6 nm. In the case of an exponential screening space charge distribution function, the resultant radial functional form for the electrostatic potential, midpoint along a cylinder of height, h, can be obtained as ’sc(r) = 2f "0 q ( exp ( r= ) + Z r= 0 1 exp( ) d + ln h 0:577 ) ; (4.26) where "0 is the same as that de ned earlier, i.e., in Eq. (4.18). The value of the electrostatic potential resulting from the exponential charge distribution at the site of a dislocation defect, i.e., at radial distance r = 0 from the 760 5 10 15 20 25 30 0 1 2 3 4 5 6 x 10 19 r (nm) n sc (r) (c m− 3 ) Nd = 10 18 cm−3 T = 300 K f = 1 λ = 2.6 nm Figure 4.5: The approximation of the screening space charge surrounding a dislocation line within InN by an exponential function. 77dislocation line, can be readily obtained from Eq. (4.26), i.e., ’osc = 2f "0 q 1 0:577 + ln h : (4.27) In Sections 4.5 and 4.6, the expression obtained for ’osc in Eq. (4.23) will be used in order to obtain the occupation statistics of the threading dislo- cation lines within n-type GaN for two di erent dislocation core structures. In Section 4.7, the expression obtained for ’osc in Eq. (4.27) will be used in order to obtain the occupation statistics of the threading dislocation lines within n-type InN. 4.5 Occupation statistics of the dangling bond dislocation core structure within GaN In this section, the Gibbs factor formalism formulated in Section 4.2 will be applied to the problem of determining the occupation statistics of the dislocation lines within GaN. For the sake of comparison with results already published in the literature, the dislocation core structure, assumed for the computations in this section, will be the same picture of dislocations used by Weimann et al. [39], i.e., a column of identical and equally spaced dangling bonds, with the distance between the successive dangling bonds 78being taken as the c-lattice constant, c, associated with wurtzite GaN. In this and the following section, in order to take the e ect of the charge trapped by the defect sites along the threading dislocation lines on the free electron concentration of the entire semiconductor material into account, the following overall charge neatrality condition will be used Nd = n+ f dis c : (4.28) In writing Eq. (4.28), the bulk donor atoms have been assumed to be com- pletely ionized, i.e., N+d = Nd. The term, n, on the right hand side of Eq. (4.28) is then the fraction of the electrons donated by the bulk donor atoms which exist as free electrons, whereas the term, f dis=c, will be the fraction of the electrons donated by the bulk donor atoms which have become trapped by the dislocation defect sites. Furthermore, the full degenerate re- lation, n = Nc 2 p Z 1 0 p x 1 + exp [x ("F "c)=kT ] dx; (4.29) for the free electron concentration within a semiconductor material, can be used along with Eq. (4.28) in order to obtain the position of the Fermi energy level, "F, for known values of Nd and dis. For the dangling bond dislocation core structure, the neutral state of a dangling bond corresponds to either of the states, s(1 #) or s(1 "), i.e., is(1#) = 79is(1") = 0. As the dangling bond loses an electron, it becomes positively charged and thus is(0) = +1. On the other hand, as each dangling bond traps an extra electron, it becomes negatively charged and the occupation state, s(2), corresponds to a state of charge is(2) = 1. For these ionization states, the value of the electrostatic potential, ’odc , can be obtained from Eq. (4.20) as follows ’odc = 2 "0 q fs(0) ln Ns(0)=2 + 1 + 0:577 fs(2) ln Ns(2)=2 + 1 + 0:577 ; (4.30) and, from Eq. (4.16), the average charge per dangling bond is given by f = fs(2) fs(0): (4.31) The value of the electrostatic potential at the site of a dangling bond, due to the screening space charge of ionized bulk donor atoms, was obtained earlier in Section 4.4 as ’osc = f "0 q 1 + 2 ln Ns(2)c=fs(2) R ; (4.32) where in this relation, the equivalent expression, Ns(2)c=fs(2), has been re- placed for the height h in Eq. (4.23). Finally, substituting for ’odc and ’osc , from Eqs. (4.30) and (4.32) into Eq. (4.12), the value of the electrostatic potential at the site of a dangling bond within wurtzite GaN, can be written 80as ’o = 2 "0 q fs(0) ln Ns(0)=2 + 1 + 0:577 fs(2) ln Ns(2)=2 + 1 + 0:577 + fs(2) fs(0) 2 1 + 2 ln Ns(2)c=fs(2) R : (4.33) The quantities, "0, R, fs(0), fs(2), and Ns(0) in Eq. (4.33) were all given in terms of the material parameters and/or other quantities, earlier in Eqs. (4.18), (4.1), (4.4), (4.6), and (4.13), respectively. In particular, from Eqs. (4.4) and (4.6), it is seen that fs(0) and fs(2) are expressed in terms of the energies, "s(1), "s(2), and "F. As was mentioned at the outset of this section, for given values of dis and Nd, "F can be obtained from Eq. (4.28). Furthermore, for the purposes of the computations in this section, the energies, "s(1) and "s(2), will be related to ’o and "t, through the use of Eqs. (4.8) and (4.9), i.e., the inner-dangling bond electron-electron interactions will be ignored. The posi- tion of the dislocation defect energy level, "t, will be set to 1:2 eV above the valence band maximum, in agreement with the value adopted by Weimann et al. [39] and Gurusinghe and Andersson [75]. For concreteness, the value of Ns(2) will be set to 105, i.e., long dislocation lines will be assumed; as long as the value of Ns(2) is su ciently large, its value will not a ect the electro- static potential, ’o, obtained. Thus, for any speci c value of the bulk donor 81density, Nd, Eq. (4.33) can be solved for the single variable, ’o, once all the necessary substitutions have been made from the other equations. After the electrostatic potential, ’o, has been solved for from Eq. (4.33), its value can be plugged into Eqs. (4.4) through (4.7) in order to obtain the occupation probabilities, fs(0), fs(1), and fs(2), associated with a dislocation dangling bond. In Figure 4.6, the fraction of the dislocation dangling bonds in their di erent possible fs(0), fs(1), and fs(2) occupation states have been plotted as a function of the bulk donor concentration; in order to emphasize the charge residing on the dangling bonds, in Figure 4.6 the occupation fractions, fs(0), fs(1), and fs(2), have alternatively been labelled as f+1, f0, and f 1, respectively, these new subscripts representing the corresponding charge states. The material parameters used for these computations are as speci ed in Table 1.1. These calculations were performed at 300 K for a dislocation line density of 108 cm 2. It is seen from this gure that, under bulk n-type doping conditions, most of the dislocation defect sites are negatively charged. However, at lower bulk donor concentrations, a non-negligible fraction of the dislocation defect sites also become positively charged. Experimentally, it has been observed that the threading dislocation lines within GaN become positively charged under bulk p-type doping condi- tions [69]. The proposed Gibbs factor formalism allows for the dislocation 821015 1016 1017 1018 1019 1020 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 f −1 f0 f +1 Doping concentration (cm−3) Fraction of ionized and neutral dislocation defect site s Figure 4.6: The fraction of ionized and neutral dislocation defect sites for a dangling bond core structure of the dislocation lines within n-type GaN as a function of the bulk donor concentration. 83defect sites to be positively charged as well. As a result, the proposed Gibbs factor formalism can potentially be applied to p-type doping conditions as well. In Figure 4.7, the average charge per dangling bond, f , is plotted as a function of the free electron concentration, n. Also plotted in this gure are the results from Weimann et al. [39] and from Gurusinghe and Andersson [75] for the same dislocation line density of dis = 108 cm 2. From Figure 4.7, it is observed that the results obtained from the application of the Gibbs factor formalism in this work [78], clearly favour the results of Gurusinghe and Andersson [75] over the results reported by Weimann et al. [39]. The considerable di erence between the reproduced result of Weimann et al. [39] and the others depicted in Figure 4.7 are most likely a consequence of di erent parameter value selections; through a detailed critical comparative analysis, it is found that the results of Weimann et al. [39] may be exactly reproduced by selecting a dislocation energy level, "t, of 2.02 eV above the valence band edge, all other parameters being as set in Table I. 841015 1016 1017 1018 1019 1020 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Free electron concentration (cm−3) Average negative charge per dangling bon d Gurusinghe and Andersson [75] Present Gibbs factor formalism Weimann et al. [39] Figure 4.7: The average negative charge per dislocation dangling bond as a function of the free electron concentration within n-type GaN. The calcula- tions were all performed at 300 K for a dislocation line density of 108 cm 2. 854.6 Occupation statistics of the VGa ON dis- location core structure within GaN In this section, another example of the application of the proposed Gibbs factor formalism to the problem of obtaining the occupation statistics of the threading dislocation lines within GaN is provided. In this section, unlike the discussion in Section 4.5, a more realistic VGa ON core structure will be assumed for the threading dislocation lines. The particular atomic con guration of this dislocation core structure was depicted earlier in Fig- ure 1.4(f). It is noted, from Figure 1.4(f), that each of the two small black nitrogen atoms lying beneath the larger white gallium atoms at the center of the VGa ON dislocation defect site, is missing a bonding to a neighbouring atom. This missing bond is due to the gallium vacancy at the defect site, hence the notation, VGa, in the VGa ON dislocation core structure. These two nitrogen atoms in Figure 1.4(f) are close enough in order to form bonding and anti-bonding states [79]. The energy band con guration of this bonding and anti-bonding energy formation at the VGa ON dislocation defect site is schematically depicted in Figure 4.8. It is clear from Figure 4.8 that the N-N anti-bonding energy level is the en- ergy level which is electrically active in the VGa ON dislocation core struc- 86 N - N anti - bo nding bo nd N - N bonding bo nd N dangling bond N dangling bond Figure 4.8: Bonding and anti-bonding energy formation from the two nitro- gen dangling bonds within the VGa ON dislocation core structure. 87ture within GaN. In its neutral state, depicted in Figure 4.8, this energy level holds no electrons, i.e., is(0) = 0. When this energy level is occupied by one electron, it becomes singly negatively charged, i.e., is(1#) = is(1") = 1, and when it is occupied by two electrons, it becomes doubly negatively charged, i.e., is(2) = 2. For these ionization states, it can be readily shown that the value of the electrostatic potential at the site of a VGa ON dislocation defect will be given by ’o = 2 "0 q fs(1) ln Ns(1)=2 + 1 + 0:577 2 fs(2) ln Ns(2)=2 + 1 + 0:577 + 2 fs(2) + fs(1) 2 1 + 2 ln Ns(2)c=fs(2) R : (4.34) Similar to Eq. (4.33), Eq. (4.34) can be solved for the single variable, ’o, once all the necessary substitutions have been made from the other equa- tions into it. In particular, the inner dislocation defect site electron-electron interactions will be ignored and the value of the dislocation defect energy level, "t, will be set to 1:2 eV above the valence band maximum. After the electrostatic potential, ’o, has been solved for from Eq. (4.34), its value can be plugged into Eqs. (4.4) through (4.7) in order to obtain the occuption probabilities, fs(0), fs(1), and fs(2) associated with a VGa ON dis- location defect site within GaN. In Figure 4.9, the fraction of the VGa ON 881015 1016 1017 1018 1019 1020 0 0.2 0.4 0.6 0.8 1 Doping concentration (cm−3) Fraction of ionized and neutral defect site s f0 f −1 f −2 Figure 4.9: The fraction of ionized and neutral dislocation defect sites for the VGa ON core structure of the dislocation lines within n-type GaN as a function of the bulk donor concentration. 89dislocation defect sites in their di erent possible fs(0), fs(1), and fs(2) occupa- tion states have been plotted as a function of the bulk donor concentration; in order to emphasize the charge residing on the dislocation defect sites, in Figure 4.9 the occupation fractions, fs(0), fs(1), and fs(2), have alternatively been labelled as f0, f 1, and f 2, respectively, these new subscripts repre- senting the corresponding charge states. The material parameters used for these computations are as speci ed in Table 1.1. These calculations were performed at 300 K for a dislocation line density of 108 cm 2. It is found that in the case of the VGa ON dislocation core structure, the Gibbs factor formalism gives a notable computational advantage over the extension of the model of Read [73], as has been pursued by You et al. [60, 77]. Furthermore, the application of the Gibbs factor formalism gives a more physically satis- fying picture for the occupation of the dislocation lines for this dislocation core structure [79]. In Figure 4.10, the average charge per dislocation defect site, f , is plotted as a function of the free electron concentration, n. For comparison, results corresponding to You et al. [60] are also plotted in this gure. It is seen that the overall trends in the two curves depicted in Figure 4.10 are similar. More speci cally, at low free electron concentrations, when most of the N-N anti- bonding energy levels are singly occupied, the di erences between the results 901015 1016 1017 1018 1019 1020 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 Free electron concentration (cm−3) Average negative charge per dislocation defec t Present Gibbs factor formalism You et al. [60] Figure 4.10: The average negative charge per VGa ON dislocation defect site as a function of the free electron concentration within n-type GaN. Re- sults corresponding to the energy minimization approach employed by You et al. [60] for the same VGa ON dislocation core structure are also depicted. The results depicted correspond to a temperature of 300 K and a dislocation line density of 108 cm 2. 91obtained in this thesis and those of You et al. [60] are negligible. However, at higher free electron concentrations, when a larger fraction of the N-N anti- bonding energy levels are doubly occupied, di erences between the results are more pronounced. As a result, these di erences may be attributed to the fact that in the computations of this section, electron-electron interactions within individual N-N anti-bonding energy levels have been neglected when determining the energy associated with the fs(2) occupation state of the N-N anti-bonding energy level. 4.7 Occupation statistics of the open core dis- location core structure within InN In this section, an open core structure will be assumed for the threading dislocation lines within InN. The atomic con guration of this dislocation core structure was depicted earlier in Figure 1.4(b). The energy band con gura- tion corresponding to this dislocation core structure within InN is depicted in Figure 4.11. It is clear, from Figure 4.11, that the N-N anti-bonding energy level is the energy level which is electrically active in the open core dislocation core structure within InN. In its neutral state, depicted in Figure 4.11, this energy 92 In dangling bond In dangling bond N dangling bond N dangling bond N - N bonding bo nd In - In anti - bo nding bo nd In - In bo nding bo nd N - N anti - bo nding bo nd 0.25 eV Figure 4.11: Schematic energy band diagram of the neutral state of the open core dislocation core structure within InN. 93level holds two electrons, i.e., is(2) = 0. When this energy level loses one of its electrons to the conduction band, it becomes singly positively charged, i.e., is(1#) = is(1") = 1, and when this energy level holds no electrons it is doubly positively charged, is(0) = 2. For these ionization states, the value of the electrostatic potential at the site of the open core dislocation defect within InN, ’o, can be obtained as follows ’o = 2 "0 q fs(1) ln Ns(1)=2 + 1 + 0:577 + 2 fs(2) ln Ns(2)=2 + 1 + 0:577 2 fs(2) + fs(1) 2 1 0:577 + ln Ns(2)c=fs(2) : (4.35) Note that due to the free electron screening of the positive charge on the dislocation defect sites within InN, Eq. (4.27) has been used to replace ’osc when obtaining the above relation. Similar to Eq. (4.33), Eq. (4.35) can be solved for the single variable, ’o, once all the necessary substitutions have been made from the other equations into it. For the case of InN, when computing the position of the Fermi energy level, "F, the slightly di erent overall charge neutrality condition, Nd + f dis c = n; (4.36) should be used, due to the fact that within this semiconductor material the dislocation defect sites themselves donate electrons instead of trapping 94the electrons donated by the bulk donor atoms; compare Eq. (4.36) with Eq. (4.28). The inner dislocation defect site electron-electron interactions will be ignored. Density functional calculations of Takie and Nakayama [59] have shown that in the open core dislocation core structure within InN, the N-N anti-bonding energy level lies 0:25 eV above the conduction band minimum. As a result, for the purposes of the computations in this section, the dislocation defect energy level, "t, will be set equal to 0:25 eV above the conduction band minimum. After the electrostatic potential, ’o, has been solved for from Eq. (4.35), its value can be plugged into Eqs. (4.4) through (4.7) in order to obtain the occupation probabilities, fs(0), fs(1), and fs(2), associated with an open core dislocation defect site within wurtzite InN. In Figure 4.12, the fraction of the open core dislocation defect sites in their di erent possible fs(0), fs(1), and fs(2) occupation states have been plotted as a function of the bulk donor con- centration; in order to emphasize the charge residing on the dislocation defect sites, in Figure 4.12 the occupation fractions, fs(0), fs(1), and fs(2), have alter- natively been labelled as f+2, f+1, and f0, respectively, these new subscripts representing the corresponding charge states. The material parameters used for these computations are as speci ed in Table 1.1. These calculations were performed at 300 K for a dislocation line density of 108 cm 2. It is seen 951015 1016 1017 1018 1019 0 0.1 0.2 0.3 0.4 0.5 Doping concentration (cm−3) Fraction of ionized and neutral defect site s f +2 f0 f +1 Figure 4.12: The fraction of ionized and neutral dislocation defect sites for the open core structure of the dislocation lines within n-type InN as a function of the bulk donor concentration. 96from this gure that, under bulk n-type doping conditions, the dislocation defect sites are predominantly positively charged. In contrast, for the case of n-type GaN, the dislocation defect sites were primarily negatively charged; recall Figures 4.6 and 4.9. In Figure 4.13, the average positive charge per dislocation defect site for the open core dislocation core structure within wurtzite InN has been plotted as a function of the free electron concentration. In comparison to the disloca- tion defect occupation statistics results for GaN, obtained in Figures 4.7 and 4.10, it is observed that the average charge per dislocation defect site within InN decreases with increasing free electron concentration. This distinction in the general behavior of the curve depicted in Figure 4.13, when compared to the curves depicted in Figures 4.7 and 4.10, can be attributed to the fact that within n-type GaN, with increasing free electron/bulk donor concentration, the position of the Fermi energy level tends to move away from the position of the dislocation defect energy level. On the other hand, within n-type InN, as the bulk donor concentration/free electron concentration increases, the Fermi energy level moves into the conduction band, and thus, gets closer to the position of the dislocation defect energy level within this semiconductor material. 971016 1017 1018 1019 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 Free electron concentration (cm−3) Average positive charge per dislocation defect sit e Figure 4.13: The average positive charge per dislocation defect site for the open core structure of the dislocation lines within n-type InN as a function of the free electron concentration. The calculations are performed at 300 K for a dislocation line density of 108 cm 2. 98Chapter 5 Conclusions In conclusion, a generalized electrostatic theory of charged dislocation lines within semiconductor materials was developed in this thesis. Speci - cally, the generalizations involved, allow for a determination of the occupation statistics of multiple ionization states of the defect sites along a dislocation line. In addition, the formulation developed in this thesis allows for the determination of the most general screening space charge pro le associated with the dislocation lines [80]. The materials considered in this analysis were the nitride based semiconductors, GaN and InN, for which strained epitaxial crystal growth results in the presence of a high density of threading disloca- tion lines within their respective crystal structures. It was found that the distinctively di erent positionings of the disloca- 99tion defect energy levels relative to the conduction band minimum in the two semiconductor materials, GaN and InN, lead to distinctively di erent pro les for the space charge distribution, screening the threading disloca- tion line charge in each material; while within n-type GaN, the screening is dominated by ionized bulk donor atoms, within n-type InN, this screening is provided for by the free electrons donated by the defect sites situated along the threading dislocation line itself. The recognition of this distinction in the screening space charge distribution of the threading dislocation lines within InN when compared to the case of GaN was missing in the literature and its computation has thus been one of the original contributions of this research. The dislocation limited electron drift and Hall mobilities within GaN and InN were also computed as part of this research. Under the same n-type bulk doping conditions, it was observed that the electron drift and Hall mobilities within InN are much higher than those computed for the case of GaN. This higher electron mobility within InN can be attributed to the much stronger screening of the dislocation core charge within InN as well as the lower value for the electron e ective mass of InN. Current experimental data existing in the literature on the threading dislocation line density dependence of the electron Hall mobility within InN, di er by an order of magnitude [66, 67] and furthermore, within each experimental data set, are much too scattered 100to enable a conclusive comparison with the theoretically computed values in this thesis for the electron Hall mobility within InN. Nevertheless, in the future, as the quality of the epitaxially grown layers of InN improve, and more accurate experimental determinations of the electron Hall mobility within InN are performed, it is expected that the theoretical results obtained in this thesis will provide a conclusive answer as to how much charge is accumulated along the threading dislocation lines. This is important, because of the fact that the threading dislocation lines within InN have been blamed for the high unintentional n-type conductivity of InN. This high unintentional n- type conductivity is in turn hampering the successful p-type doping of this important semiconductor material and thus, the fabrication of functional electronic/optoelectronic devices from InN. The application of the Gibbs factor formalism to the problem of deter- mining the occupation statistics of the defect sites along a dislocation line has been another original contribution of this work. It was shown that in special cases the formulation essentially reduces to the model of Read [73] for the treatment of the dislocation lines. In continuation of the work presented in this thesis, several avenues for future investigation are conceivable. These topics for further investigation include the following: 101(1) In the literature, only limiting case analytical expressions have been used in order to obtain the electron drift mobility limited by scattering from the threading dislocation lines within GaN [39, 40, 75]. The scattering po- tential that has been computed in this thesis, as well as experimental data from the literature [69, 71], suggest that the actual scattering potential devi- ates substantially from the simple limiting case analytical expressions used in the literature for the purposes of computing the dislocation limited electron mobility within GaN. (2) In Section 3.2, when formulating the dislocation limited electron drift and Hall mobilities, it was assumed that the dislocation lines are well sep- arated from each other. Under high threading dislocation line density con- ditions, i.e., when the screening space charge cylinders of nearby dislocation lines start to overlap, this assumption of the formulation developed in Chap- ter 3 will no longer be valid. Furthermore, under high dislocation line den- sity conditions, the occupation statistics of the threading dislocation lines obtained in Chapter 4 relies on the implicit assumption that the dislocation lines are uniformly distributed within the semiconductor sample. Thus, in the high dislocation line density limit, a more rigorous treatment of the e ect of non-uniformity in the distribution of the threading dislocation lines, on both the free carrier concentration and the electron mobility, remains to be 102addressed. 103References [1] S. Nakamura, \InGaN/AlGaN blue-light-emitting diodes," J. Vac. Sci. Technol. A, vol. 13, no. 3, pp. 705-710, May/June 1995. [2] S. Nakamura, \The roles of structural imperfections in InGaN-based blue light-emitting diodes and laser diodes," Science, vol. 281, pp. 956- 961, August 1998. [3] S. Strite and H. Morko c, \GaN, AlN, and InN: A review," J. Vac. Sci. Technol. B, vol. 10, no. 4, pp. 1237-1266, July/August 1992. [4] H. Morko c, S. Strite, G. B. Gao, M. E. Lin, B. Sverdlov, and M. Burns, \Large-band-gap SiC, III-V nitride, and II-VI ZnSe-based semiconduc- tor device technologies," J. Appl. Phys., vol. 76, no. 3, pp. 1363-1398, August 1994. 104[5] S. N. Mohammad and H. Morko c, \Progress and prospects of group- III nitride semiconductors," Prog. Quantum Electron., vol. 20, no. 5/6, pp. 361-525, 1996. [6] S. J. Pearton, J. C. Zolper, R. J. Shul, and F. Ren, \GaN: Processing, defects, and devices," J. Appl. Phys., vol. 86, no. 1, pp. 1-78, July 1999. [7] J. Q. Wu, \When group-III nitrides go infrared: New properties and perspectives," J. Appl. Phys., vol. 106, pp. 011101-1-28, July 2009. [8] F. A. Ponce, and D. P. Bour, \Nitride-based semiconductors for blue and green light-emitting devices," Nature, vol. 386, pp. 351-359, March 1997. [9] U. K. Mishra, L. Shen, T. E. Kazior, and Y.-F. Wu, \GaN-based rf power devices and ampli ers," Proceedings of the IEEE, vol. 96, no. 2, pp. 287-305, February 2008. [10] W. C. Johnson, J. B. Parsons, and M. C. Crew, \Nitrogen compounds of gallium. III," J. Phys. Chem., vol. 36, no. 10, pp. 2651-2654, January 1932. 105[11] H. J. Hovel and J. J. Cuomo, \Electrical and optical properties of rf- sputtered GaN and InN," Appl. Phys. Lett., vol. 20, no. 2, pp. 71-73, January 1972. [12] K. V. Malyutin and A. A. Veshkin, \Preparation of GaN lms by ca- thodic sputtering," Neorg. Mater., vol. 8, no. 5, pp. 821-823, May 1972. [13] J. C. Vesely, M. Shatzkes, and P. J. Burkhard, \Space-charge-limited current ow in gallium nitride thin lms," Phys. Rev. B, vol. 10, no. 2, pp. 582-590, July 1974. [14] T. Hariu, T. Usuba, H. Adachi, and Y. Shibata, \Reactive sputtering of gallium nitride thin lms for GaAs MIS structures," Appl. Phys. Lett., vol. 32, no. 4, pp. 252-253, February 1978. [15] E. Lakshmi, B. Mathur, A. B. Bhattacharya, and V. P. Bhargava, \The growth of highly resistive gallium nitride lms," Thin Solid Films, vol. 74, no. 1, pp. 77-82, December 1980. [16] K. Matsushita, Y. Matsuno, T. Hariu, and Y. Shibata, \A comparative study of the deposition conditions in the plasma-assisted deposition of gallium nitride thin lms," Thin Solid Films, vol. 80, no. 1-3, pp. 243- 247, June 1981. 106[17] K. Kubota, Y. Kobayashi, and K. Fujimoto, \Preparation and properties of III-V nitride thin lms," J. Appl. Phys., vol. 66, no. 7, pp. 2984-2988, October 1989. [18] H. P. Maruska and J. J. Tietjen, \The preparation and properties of vapor-deposited single-crystal-line GaN," Appl. Phys. Lett., vol. 15, no. 10, pp. 327-329, November 1969. [19] J. I. Pankove, J. E. Berkeyheiser, H. P. Maruska, and J. Wittke, \Lu- minescent properties of GaN," Solid State Commun., vol. 8, no. 13, pp. 1051-1053, July 1970. [20] R. K. Crouch, W. J. Debnam, and A. L. Fripp, \Properties of GaN grown on sapphire substrates," J. Mater. Sci., vol. 13, no. 11, pp. 2358- 2364, 1978. [21] A. S. Adonin, V. A. Evmenenko, L. N. Mikhailov, and N. G. Ryabtsev, \Study of certain kinetic features of the growth and doping of heteroepi- taxial layers of gallium nitride, their cathodoluminescence and electric properties," Neorg. Mater., vol. 17, no. 9, pp. 1608-1613, September 1981. 107[22] R. Juza and H. Hahn, \Crystal structures of Cu3N, GaN and InN - Metallic amides and metallic nitrides V Announcement," Zeitschr. Anorgan. Allgem. Chem., vol. 239, no. 3, pp. 282-287, October 1938. [23] J. Pastrnak and L. Souckova, \Production of thin layers of aluminum, gallium and indium nitride under a gas discharge," Phys. Status Solidi, vol. 3, no. 2, pp. K71-K74, 1963. [24] T. L. Tansley and C. P. Foley, \Electron mobility in indium nitride," Electron. Lett., vol. 20, no. 25-26, pp. 1066-1068, December 1984. [25] L. A. Marasina, I. G. Pichugin, and M. Tlaczala, \Preparation of InN epitaxial layers in InCl3 NH3 system," Krist. Technol., vol. 12, no. 6, pp. 541-455, January 1977. [26] H. Lu, W. J. Scha , and L. F. Eastman, \Surface charge accumulation of InN lms grown by molecular-beam epitaxy," Appl. Phys. Lett., vol. 82, no. 11, pp. 1736-1768, March 2003. [27] V. Y. Davydov, A. A. Klochikhin, R. P. Seisyan, V. V. Emtsev, S. V. Ivanov, F. Bechstedt, J. Furthmuller, H. Harima, A. V. Mudryi, J. Aderhold, O. Semchinova, and J. Graul, \Absorption and emission of hexagonal InN. Evidence of narrow fundamental band gap," Phys. Status Solidi (b), vol. 229, no. 3, pp. r1-r3, February 2002. 108[28] T. L. Tansley and C. P. Foley, \Optical band gap of indium nitride," J. Appl. Phys., vol. 59, no. 9, pp. 3241-3244, May 1986. [29] V. A. Tyagai, A. M. Evstigneev, A. N. Krasiko, A. F. Andreeva, and V. Y. Malakhov, \Optical properties of indium nitride lms," Sov. Phys. Semicond., vol. 11, no. 11, pp. 1257-1259, November 1977. [30] K. L. Westra, R. P. W. Lawson, and M. J. Brett, \The e ects of oxygen contamination on the properties of reactively sputtered indium nitride lms," J. Vac. Sci. Technol. A, vol. 6, no. 3, pp. 1730-1732, May/June 1988. [31] Y. E. Romanyuk, D. Kreier, Y. Cui, K. M. Yu, J. W. Ager III, and S. R. Leone, \Molecular beam epitaxy of InGaN thin lms on Si(111): E ect of substrate nitridation," Thin Solid Films, vol. 517, no. 24, pp. 6512- 6515, October 2009. [32] K. P. Beh, F. K. Yam, C. W. Chin, S. S. Tneh, and Z. Hassan, \The growth of III-V nitrides heterostructure on Si substrate by plasma- assisted molecular beam epitaxy," J. Alloys Compd., vol. 506, no. 1, pp. 343-346, September 2010. [33] M. Kurouchi, T. Yamaguchi, H. Naoi, A. Suzuki, T. Araki, and Y. Nanishi, \Growth of In-rich InGaN on InN template by radio-frequency 109plasma assisted molecular beam epitaxy," J. Cryst. Growth, vol. 275, no. 1-2, pp. e1053-e1058, February 2005. [34] N. Li, S. -J. Wang, E. -H. Park, Z. C. Feng, H. -L. Tsai, J. -R. Yang, and I. Ferguson, \Suppression of phase separation in InGaN layers grown on lattice-matched ZnO substrates," J. Cryst. Growth, vol. 311, no. 22, pp. 4628-4631, November 2009. [35] D. Kapolnek, X. H. Wu, B. Heying, S. Keller, B. P. Keller, U. K. Mishra, S. P. Den Baars, and J. S. Speck, \Structural evolution in epitaxial metalorganic chemical vapor deposition grown GaN lms on sapphire," Appl. Phys. Lett., vol. 67, no. 11, pp. 1541-1543, September 1995. [36] F. A. Ponce, B. S. Krusor, J. S. Major, Jr., W. E. Plano, and D. F. Welch, \Microstructure of GaN epitaxy on SiC using AlN bu er layers," J. Appl. Phys., vol. 67, no. 3, pp. 410-412, July 1995. [37] V. Cimalla, V. Lebedev, F. M. Morales, M. Niebelsch utz, G. Ecke, R. Goldhahn, and O. Ambacher, \Origin of n-type conductivity in nomi- nally undoped InN," Mat.-wiss. u. Werksto tech, vol. 37, No. 11, pp. 924-928, November 2006. [38] R. Jaszek, \Carrier scattering by dislocations in semiconductors," J. Mater. Sci.: Mater. Electron., vol. 12, no. 1, pp. 1-9, 2001. 110[39] N. G. Weimann, L. F. Eastman, D. Doppalapudi, H. M. Ng, and T. D. Moustakas, \Scattering of electrons at threading dislocations in GaN," J. Appl. Phys., vol. 83, no. 7, pp. 3656-3659, April 1998. [40] D. C. Look, and J. R. Sizelove, \Dislocation scattering in GaN," Phys. Rev. Lett., vol. 82, no. 6, pp. 1237-1240, February 1999. [41] F. A. Ponce, D. P. Bour, W. Gotz, and P. J. Wright, \Spatial distribution of the luminescence in GaN thin lms," Appl. Phys. Lett., vol. 68, no. 1, pp. 57-59, January 1996. [42] T. Lei, M. Fanciulli, R. J. Molnar, T. D. Moustakas, R. J. Graham, and J. Scanlon, \Epitaxial growth of zinc blende and wurtzitic gallium nitride thin lms on (001) silicon," Appl. Phys. Lett., vol. 59, no. 8, pp. 944-946, August 1991. [43] S.-H. Park and S.-L. Chuang, \Comparison of zinc-blende and wurtzite GaN semiconductors with spontaneous polarization and piezoelectric eld e ects," J. Appl. Phys., vol. 87, no. 1, pp. 353-364, January 2000. [44] M. E. Levinshtein, S. L. Rumyantsev, and M. Shur, Properties of Ad- vanced Semiconductor Materials: GaN, AlN, InN, BN, SiC, SiGe. New York: John Wiley and Sons Inc., pp. 1-30, 2001. 111[45] T. Unishima, M. Highashiwaki, and T. Matsui, \Optical properties of Si-doped InN grown on sapphire (0001)," Phys. Rev. B, vol. 68, no. 23, pp. 235204-1-7, December 2003. [46] J. Wu, W. Walukiewicz, K. M. Yu, J. W. Ager, E. E. Haller, H. Lu, W. J. Scha , Y. Saito, and Y. Nanishi, \Unusual properties of the fundamental band gap of InN," Appl. Phys. Lett., vol. 80, no. 21, pp. 3967-3969, May 2002. [47] H. Wang, and A.-B. Chen, \Calculation of shallow donor levels in GaN," J. Appl. Phys., vol. 87, no. 11, pp. 7859-7863, June 2000. [48] A. Janotti and C. G. Van de Walle, \Sources of unintentional conduc- tivity in InN," Appl. Phys. Lett., vol. 92, pp. 032104-1-3, January 2008. [49] M. E. Levinshtein, S. L. Rumyantsev, and M. Shur, Properties of Ad- vanced Semiconductor Materials: GaN, AlN, InN, BN, SiC, SiGe. New York: John Wiley Sons Inc., pp. 49-66, 2001. [50] O. Lagerstedt and A. Monemar, \Variation of lattice parameters in GaN with stoichiometry and doping," Phys. Rev. B, vol. 19, no. 6, pp. 3064- 3070, March 1979. 112[51] I. G. Pichugin and M. Tlachala, \X-Ray analysis of indium nitride," Inorg. Mater., vol. 14, no. 1, pp. 135-136, 1978. [52] David Darling, \The encyclopaedia of science, Solid State Physics" http://www.daviddarling.info/encyclopedia/D/dislocation.html [53] X. H. Wu, L. M. Brown, K. Kapolnek, S. Keller, B. Keller, S. P. Den- Baars, and J. S. Speck, \Defect structure of metal-organic chemical vapor deposition-grown epitaxial (0001) GaN=Al2O3," J. Appl. Phys., vol. 80, no. 6, pp. 3228-3237, September 1996. [54] T. Hino, S. Tomiya, T. Miyajima, K. Yanashima, S. Hashimoto, and M. Ikeda, \Characterization of threading dislocations in GaN epitaxial layers," Appl. Phys. Lett., vol. 76, no. 23, pp. 3421-3423, June 2000. [55] J. Arvanitidis, D. Christo los, G. A. Kourouklis, A. Delimitis, M. Kat- sikini, P. Komninou, S. Ves, E. Dimakis, and A. Georgakilas, \Depth pro le of the biaxial strain in a 10 m thick InN (0001) lm," J. Appl. Phys., vol. 100, pp. 113516-1-5, December 2006. [56] C. Rauch, F. Tuomisto, A. V. Clemente, B. Lacroix, P. Ruterana, S. Kraeusel, B. Hourahine, and W. J. Scha , \Defect evolution and inter- play in n-type InN," Appl. Phys. Lett., vol. 100, pp. 091907-1-3, March 2012. 113[57] S. M. Lee, M. A. Belkhir, X. Y. Zhu, Y. H. Lee, Y. G. Hwang and T. Frauenheim, \Electronic structures of GaN edge dislocations," Phys. Rev. B, vol. 61, no. 23, pp. 16033-16039, June 2000. [58] L. Lymperakis, J. Neugebauer, M. Albrecht, T. Remmele, and H. P. Strunk, \Strain induced deep electronic states around threading dis- locations in GaN," Phys. Rev. Lett., vol. 93, no. 19, pp. 196401-1-4, November 2004. [59] Y. Takei and T. Nakayama, \Electron-carrier generation by edge dislo- cations in InN lms: First-principles study," J. Cryst. Growth, vol. 311, no. 10, pp. 2767-2771, May 2009. [60] J. H. You, J. Q. Lu, and H. T. Johnson, \Electron scattering due to threading edge dislocations in n-type wurtzite GaN," J. Appl. Phys., vol. 99, pp. 033706-1-10, February 2006. [61] E. Kalesaki, J. Kioseoglou, L. Lymperakis, P. Komninou, and T. Karakostas, \E ect of edge threading dislocations on the electronic structure of InN," Appl. Phys. Lett., vol. 98, pp. 072103-1-3, February 2011. [62] C. J. Gallagher, \Plastic deformation of germanium and silicon," Phys. Rev., vol. 88, no. 4, pp. 721-722, November 1952. 114[63] G. L. Pearson, W. T. Read, and F. J. Morin, \Dislocations in plastically deformed germanium," Phys. Rev., vol. 93, no. 4, pp. 666-667, February 1954. [64] W. Schr oter and H. Cerva, \Interaction of point defects with dislocations in silicon and germanium: Electrical and optical e ects," Solid State Phenom., vol. 85-86, pp. 67-143, 2002. [65] C. Shi, P. M. Asbeck, and E. T. Yu, \Piezoelectric polarization associ- ated with dislocations in wurtzite GaN," Appl. Phys. Lett., vol 74, no. 4, pp. 573-575, January 1999. [66] H. Wang, D. S. Jiang, L. L. Wang, X. Sun, W. B. Liu, D. G. Zhao, J. J. Zhu, Z. S. Liu, Y. T. Wang, S. M. Zhang, and H. Yang, \Investigation on the structural origin of n-type conductivity in InN lms," J. Phys. D: Appl. Phys., vol. 41, no. 13, pp. 135403-1-5, June 2008. [67] C. S. Gallinat, G. koblm uller, and J. S. Speck, \The role of threading dis- locations and unintentionally incorporated impurities on the bulk elec- tron conductivity of In-face InN," Appl. Phys. Lett., vol. 95, pp. 022103- 1-3, July 2009. 115[68] J. Cai and A. Ponce, \Determination by electron holography of the elec- tronic charge distribution at threading dislocations in epitaxial GaN," Phys. Stat. Sol. (a), vol. 192, no. 2, pp. 407-411, July 2002. [69] D. Cherns, C. G. Jiao, H. Mokhtari, J. Cai, and F. A. Ponce, \Electron holography studies of the charge on dislocations in GaN," Phys. Stat. Sol. (b), vol. 234, no. 3, pp. 924-930, December 2002. [70] B. S. Simpkins, E. T. Yu, P. Waltereit, and J. S. Speck, \Correlated scanning Kelvin probe and conductive atomic force microscopy studies of dislocations in gallium nitride," J. Appl. Phys., vol. 94, no. 3, pp. 1448- 1453, August 2003. [71] E. M uller, D. Gerthsen, P. Br uckner, F. Scholz, T. Gruber, and A. Waag, \Probing the electrostatic potential of charged dislocations in n- GaN and n-ZnO epilayers by transmission electron holography," Phys. Rev. B, vol. 73, no. 24, pp. 245316-1-9, June 2006. [72] L. F. J. Piper, T. D. Veal, C. F. McConville, H. Lu, and W. J. Scha , \Origin of the n-type conductivity of InN: The role of positively charged dislocations," Appl. Phys. Lett., vol. 88, pp. 252109-1-3, June 2006. [73] W. T. Read, Jr., \Theory of dislocations in germanium," Philos. Mag., vol. 45, no. 367, pp. 775-796, August 1954. 116[74] J. Singh, Electronic and Optoelectronic Properties of Semiconductor Structures., Cambridge: Cambridge University Press, Chapter 4, 2003. [75] M. N. Gurusinghe and T. G. Andersson, \Mobility in epitaxial GaN: Limitations of free-electron concentration due to dislocations and com- pensation," Phys. Rev. B, vol. 67, no. 23, pp. 235208-1-7, June 2003. [76] K. Leung, A. F. Wright, and E. B. Stechel, \Charge accumulation at a threading edge dislocation in gallium nitride," Appl. Phys. Lett., vol. 74, no. 17, pp. 2495-2497, April 1999. [77] J. H. You, J. Q. Lu, and H. T. Johnson, \Atomistically informed electro- static model of an edge dislocation in a complex crystalline material," Mathematics and Mechanics of Solids, vol. 13, no. 3-4, pp. 267-291, May 2008. [78] E. Baghani and S. K. O’Leary, \Occupation statistics of dislocations within uncompensated n-type wurtzite gallium nitride," J. Appl. Phys., vol. 109, pp. 113706-1-6, June 2011. [79] E. Baghani and S. K. O’Leary, \Occupation statistics of the VGa ON dislocations within n-type gallium nitride," J. Appl. Phys., vol. 110, pp. 033509-1-6, August 2011. 117[80] E. Baghani and S. K. O’Leary, \Electron mobility limited by scatter- ing from screened positively charged dislocation lines within indium ni- tride," Appl. Phys. Lett., vol. 99, pp. 262106-1-3, December 2011. 118Appendices Appendix A. Iterative procedure for solving the screening space charge distribution within GaN The iterative procedure described in this and the following appendix, is based upon nding successively more accurate solutions to the screening space charge distribution function, nsc(r). For concreteness, in the example solved in this appendix, the value of f has been xed to 1, i.e., it is assumed that on average, one electron is being trapped per dislocation defect site. The bulk donor concentration, Nd, has been xed to the representative value of 1018 cm 3. The same procedure described in this appendix can be read- ily applied in order to obtain an accurate approximation to the screening space charge distribution function, nsc(r), for other values of f and Nd as 119well. Throughout this appendix, the dislocation core charge cylinder will be assumed to be the same as that de ned earlier, i.e., in Eq. (2.1). In the case of the screening of the dislocation core charge within GaN, the model of Read [73] provides us with a rst approximation to the screening space charge distribution function, nsc(r), i.e., nsc(r) ( Nd; r R 0; r > R ; (A.1) where the Read radius, R, is de ned in Section 4.3. This approximate screen- ing space charge distribution function is depicted with the blue curve in Figure A.1(a). Substitution of the approximate screening space charge dis- tribution function from Eq. (A.1) into Eq. (2.9) and subsequently, the resul- tant electrostatic potential into Eqs. (2.10) and (2.11), yields a new screening space charge distribution function depicted by the red curve in Figure A.1(a). If the screening space charge distribution suggested in Eq. (A.1) were the ex- act simultaneous solution of Eqs. (2.3) and (2.9), the blue and red curves in Figure A.1(a) would coincide. From Figure A.1(a) it is noted that when compared with the blue curve, the red curve in this gure is de cient in the amount of positive charge con- tained within it. As a result, the space charge distribution function depicted with the red curve in Figure A.1(a) does not meet the charge neutrality con- dition, i.e., the total amount of positive charge contained within it does not 1200 20 40 60 0 2 4 6 8 10 x 1017 (a) r (nm) n sc (r) (c m− 3 ) 0 20 40 60 0 2 4 6 8 10 x 1017 r (nm) n sc (r) (c m− 3 ) (b) 0 20 40 60 0 2 4 6 8 10 x 1017 r (nm) n sc (r) (c m− 3 ) (c) 0 20 40 60 0 2 4 6 8 10 x 1017 r (nm) n sc (r) (c m− 3 ) (d) Figure A.1: An iterative construction of the screening space charge surround- ing a threading dislocation line within n-type GaN. 121exactly equal the xed negative charge within the dislocation core charge cylinder. For the next suggested approximation to the screening space charge distri- bution function, nsc(r), this de ciency in the net positive charge needs to be remedied. In order to compensate for this de ciency, in the next suggested approximation for the screening space charge distribution, part of the red curve from Figure A.1(a) is replaced by an exponential tail at radial distances above a truncation radius, rtrunc, indicated by the black line on the horizontal axis in Figure A.1(b). In other words, for 0 r < rtrunc, the new suggested approximation to the screening space charge distribution function, nsc(r), is identical to the red curve in Figure A.1(a). For rtrunc r < Rsc, this new suggested approximation is simply an exponential function, A exp( r= ). The two variables, A and , can be uniquely determined by the requirement that: (1) nsc(r) should remain continuous at rtrunc, and (2) the total negative charge subtended by the newly suggested approximation to nsc(r) should be equal to the xed negative charge in the dislocation core charge cylinder. The resultant screening space charge distribution is depicted by the blue curve in Figure A.1(b). Subsequent substitution of the screening space charge distri- bution depicted by the blue curve in Figure A.1(b) into Eqs. (2.9), (2.10), and (2.11), yields the red curve depicted in Figure A.1(b). 122In Figures A.1(c) and A.1(d), the same process of replacing an exponential tail to the red curve from the previous iteration has been repeated. The radius at which the exponential tail is replaced in each iteration has been chosen so that improved convergence between the blue and red curves in successive iterations is maintained. It is observed that su cient convergence is reached in Figure A.1(d), i.e., a self-consistent solution to Eqs. (2.3) and (2.9) has been found. 123Appendix B. Iterative procedure for solving the screening space charge distribution within InN For concreteness, in the example solved in this appendix, the value of f has been xed to 1, i.e., it is assumed that on average, each dislocation defect site donates one electron. The bulk donor concentration, Nd, has been xed to the representative value of 1018 cm 3. The same procedure described in this appendix can be readily applied in order to obtain an accurate approxi- mation to the screening space charge distribution function, nsc(r), for other values of f and Nd as well. Throughout this appendix, the dislocation core charge cylinder will be assumed to be the same as that de ned earlier, i.e., in Eq. (2.1). In the case of InN, the screening of the dislocation core charge is pro- vided for by the free electrons donated by the dislocation defect sites. As their name suggests, free electrons are free to aggregate around the positive charge on the dislocation core charge cylinder. Consequently, the free elec- tron concentration is expected to be higher near the dislocation core and to drop o at radial distances further away from the dislocation line. The process may thus be initiated by approximating such screening space charge 124distribution by a simple exponential function, i.e., nsc(r) n0 exp( r= ): (B.1) Initially, the value of is set equal to the degenerate Debye length within the bulk of InN, which for the free electron concentration of 1018 cm 3, equals 10:57 nm. The other unknown variable in Eq. (B.1), namely, n0, is uniquely determined from the requirement that the total negative charge subtended by the screening space charge distribution, nsc(r), equals the positive charge within the dislocation core charge cylinder. This suggested screening space charge distribution is depicted with the blue curve in Figure B.1(a). Substituting the screening space charge distribution function, nsc(r), from Eq. (B.1) into Eq. (2.9), the electrostatic potential resulting from this space charge distribution, screening the dislocation core charge, can be obtained. Substituting the resultant electrostatic potential function into Eqs. (2.10) and (2.11), results in a new screening space charge distribution function, depicted by the red curve in Figure B.1(a). It is seen from Figure B.1(a) that the two screening space charge functions depicted deviate substantially. This deviation is an indication that the initial screening space charge distribution function, suggested by Eq. (B.1), is not a good approximation to the actual self-consistent solution of Eqs. (2.3) and (2.9). The steepness of the resultant red curve in Figure B.1(a) suggests that improved convergence might be 1250 10 20 30 0 1 2 3 4 5 6 x 1019 (a) r (nm) n sc (r) (c m− 3 ) 0 10 20 30 40 0 1 2 3 4 5 6 x 1019 r (nm) n sc (r) (c m− 3 ) (b) 0 10 20 30 40 0 1 2 3 4 5 6 x 1019 r (nm) n sc (r) (c m− 3 ) (c) 0 10 20 30 40 0 1 2 3 4 5 6 x 1019 r (nm) n sc (r) (c m− 3 ) (d) Figure B.1: An iterative construction of the screening space charge surround- ing a threading dislocation line within n-type InN. 126achieved by choosing a steeper exponential function than the initial guess for nsc(r). In Figure B.1(b), the same procedure as that used in obtaining Fig- ure B.1(a) has been followed, however, this time, the value of has been selected to be 4 nm, instead of the Debye length. Note that once the value of is changed, the value for n0 should also change accordingly, in order to maintain charge neutrality, i.e., equality of the negative charge subtended by the newly suggested screening space charge distribution function with the xed positive charge in the dislocation core charge cylinder. The next suggested approximation to the actual screening space charge distribution function is depicted with the blue curve in Figure B.1(c). For 0 r < rtrunc, this suggested approximation is identical to the red curve in Figure B.1(b). For rtrunc r < Rsc, this suggested approximation is simply an exponential function, A exp( r= 0). The truncation radius, rtrunc, is indicated with the black line on the horizontal axis in Figure B.1(c). The two variables, A and 0, can be uniquely determined by the requirements that: (1) nsc(r) should remain continuous at rtrunc, and (2) the total negative charge subtended by the newly suggested approximation to nsc(r), should be equal to the xed positive charge in the dislocation core charge cylinder. In Figure B.1(d), the same procedure of replacing an exponential tail 127to the red curve, obtained from the previous iteration, has been followed. The radius at which the exponential function is replaced has once again been indicated by the black line on the horizontal axis. It is observed that su cient convergence has been reached in Figure B.1(d), i.e., a self-consistent solution to Eqs. (2.3) and (2.9) has been obtained. 128
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Electrical properties of dislocations within the nitride...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Electrical properties of dislocations within the nitride based semiconductors gallium nitride and indium.. Baghani, Erfan 2012-11-14
pdf
Page Metadata
Item Metadata
Title | Electrical properties of dislocations within the nitride based semiconductors gallium nitride and indium nitride |
Creator |
Baghani, Erfan |
Publisher | University of British Columbia |
Date Issued | 2012 |
Description | Dislocation lines affect the electrical and optical properties of semiconductors. In this research, the effect that the threading dislocation lines have on the free electron concentration and the electron mobility within gallium nitride and indium nitride is investigated. A formulation is developed for obtaining the screening space charge concentration and the corresponding electrostatic potential profile surrounding the dislocation lines. The resultant electrostatic potential profile has then been used to compute the associated electron mobility, limited by scattering from the charged dislocation lines. As part of this research, a Gibbs factor formalism is also developed that can readily obtain the occupation statistics of the defect sites associated with the threading dislocation lines. |
Genre |
Thesis/Dissertation |
Type |
Text |
Language | eng |
Date Available | 2012-11-14 |
Provider | Vancouver : University of British Columbia Library |
DOI | 10.14288/1.0073375 |
URI | http://hdl.handle.net/2429/43581 |
Degree |
Doctor of Philosophy - PhD |
Program |
Electrical and Computer Engineering |
Affiliation |
Applied Science, Faculty of Engineering, School of (Okanagan) |
Degree Grantor | University of British Columbia |
Graduation Date | 2013-05 |
Campus |
UBCO |
Scholarly Level | Graduate |
Aggregated Source Repository | DSpace |
Download
- Media
- 24-ubc_2013_spring_baghani_erfan.pdf [ 1.75MB ]
- Metadata
- JSON: 24-1.0073375.json
- JSON-LD: 24-1.0073375-ld.json
- RDF/XML (Pretty): 24-1.0073375-rdf.xml
- RDF/JSON: 24-1.0073375-rdf.json
- Turtle: 24-1.0073375-turtle.txt
- N-Triples: 24-1.0073375-rdf-ntriples.txt
- Original Record: 24-1.0073375-source.json
- Full Text
- 24-1.0073375-fulltext.txt
- Citation
- 24-1.0073375.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.24.1-0073375/manifest