An optical functional analysis of amorphous semiconductors by Sanjida Begum Minar B.Sc., Bangladesh University of Engineering and Technology, 2009 A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF MASTER OF APPLIED SCIENCE in The College of Graduate Studies (Electrical Engineering) THE UNIVERSITY OF BRITISH COLUMBIA (Okanagan) December 2012 c Sanjida Begum Minar 2012 Abstract Models for the spectral dependence of the real and imaginary components of the dielectric function, appropriate for the case of amorphous semiconductors, are considered. In the first phase of this analysis, from an empirical expression for the imaginary part of the dielectric function, this expression corresponding to that of Jellison and Modine [G. E. Jellison, Jr. and F. A. Modine, “Parameterization of the optical functions of amorphous materials in the interband region,” Applied Physics Letters, vol. 69, pp. 371-373, 1996], a closed-form expression for the real part of the dielectric function is determined using a Kramers-Kronig transformation. The resultant expression for the real component of the dielectric function corresponds with that of the model of Jellison and Modine. The subsequent comparison with experiment is found to be satisfactory. Then, in the latter stage of this analysis, through the application of a Kramers-Kronig transformation on an empirical model for the imaginary part of the dielectric function, this model stemming from ii Abstract a model for the distributions of electronic states, the spectral dependence of the real part of the dielectric function is determined. Fits with the results of experiment, taken over the near-infrared, visible, and near-ultraviolet, are also found to be satisfactory. iii Table of Contents Abstract . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ii Table of Contents . . . . . . . . . . . . . . . . . . . . . . . . . . . . iv List of Tables . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . viii List of Figures . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . List of Symbols List of Acronyms ix . . . . . . . . . . . . . . . . . . . . . . . . . . . . . xx . . . . . . . . . . . . . . . . . . . . . . . . . . . . xxiii Acknowledgements . . . . . . . . . . . . . . . . . . . . . . . . . . . xxiv Dedication . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . xxv 1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 Theoretical basis for analysis 1 . . . . . . . . . . . . . . . . . . 15 iv Table of Contents 2.1 The optical response of materials . . . . . . . . . . . . . . . . 15 2.2 The spectral dependence of the imaginary part of the dielectric function . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18 2.3 Three critical assumptions . . . . . . . . . . . . . . . . . . . . 19 2.4 A momentum matrix element formulation . . . . . . . . . . . 24 2.5 A dipole matrix element formulation . . . . . . . . . . . . . . 28 2.6 Free electron model of the DOS functions 2.7 The spectral dependence of the JDOS function 2.8 Kramers-Kronig relations 2.9 Review of models . . . . . . . . . . . . . . . . . . . . . . . . . 46 . . . . . . . . . . . 31 . . . . . . . . 34 . . . . . . . . . . . . . . . . . . . . 39 2.10 Experimental results . . . . . . . . . . . . . . . . . . . . . . . 49 3 The model of Jellison and Modine 3.1 . . . . . . . . . . . . . . . 53 Empirical models for the optical functions associated with an amorphous semiconductor . . . . . . . . . . . . . . . . . . . . 53 3.2 The Lorentz model . . . . . . . . . . . . . . . . . . . . . . . . 55 3.3 Model for the optical functions 3.4 Derivation 3.4.1 . . . . . . . . . . . . . . . . . 57 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60 Definitions . . . . . . . . . . . . . . . . . . . . . . . . 60 v Table of Contents 3.4.2 Root calculation through expansion of the denominator . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61 3.4.3 Calculation of the coefficients . . . . . . . . . . . . . . 63 3.4.4 Grouping the terms in the partial fraction expression . 71 3.4.5 Summary of the partial fraction expansion . . . . . . . 75 3.4.6 Evaluation of the constants c3 , c4 , c5 , and c6 . . . . . 76 3.4.7 Main integration . . . . . . . . . . . . . . . . . . . . . 85 3.4.8 Summation by parts . . . . . . . . . . . . . . . . . . . 90 3.4.9 Calculation of term 3: . . . . . . . . . . . . . . . . . . 95 3.4.10 Summation of term 1, term 2 and term 3 . . . . . . . . 97 3.4.11 Final Tauc-Lorentz expression for the model of Jellison and Modine . . . . . . . . . . . . . . . . . . . . . . . . 98 3.5 Comparing the results of the Jellison and Modine model with experimental data . . . . . . . . . . . . . . . . . . . . . . . . 99 4 A physically based Kramers-Kronig consistent model for the optical functions associated with amorphous semiconductors 105 4.1 Evaluating the JDOS function from the distributions of electronic states 4.2 . . . . . . . . . . . . . . . . . . . . . . . . . . . 105 Empirical model for the DOS functions . . . . . . . . . . . . 107 vi Table of Contents 4.3 JDOS evaluation and analysis . . . . . . . . . . . . . . . . . . 115 4.4 Modeling the optical functions associated with a-Si . . . . . . 119 4.5 Results 4.6 Implications of a non-unity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 132 1∞ . . . . . . . . . . . . . . . . . 137 5 Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 138 References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 140 vii List of Tables 3.1 The “best-fit” a-Si modeling parameter selections employed for the Tauc-Lorentz model of Jellison and Modine [15] for fits to experimental data corresponding to EGE a-Si, PECVD a-Si I, and PECVD a-Si II. . . . . . . . . . . . . . . . . . . . . . . 101 4.1 The “best-fit” a-Si DOS modeling parameter selections employed for the purposes of this analysis. These modeling parameters relate to Eqs. (2.17), (2.45), and (4.8) . . . . . . . . 133 viii List of Figures 1.1 Semiconductor and flat panel display shipments plotted as a function of date. This data was obtained from the 2005 Information Society Technologies proposal for advancement [10]. . . 1.2 4 A two-dimensional schematic representation of the atomic distributions associated with hypothetical crystalline, microcrystalline, and amorphous semiconductors. This figure is from Nguyen [12]. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1.3 6 A schematic representation of the light intensity as a function of the penetration depth, z, with light normally incident upon the material propagating from the left. This figure is from Thevaril [13]. . . . . . . . . . . . . . . . . . . . . . . . . . . . 1.4 8 Variations in the relative intensity, I(z)/Io , with the penetration depth, z, for different selections of the optical absorption coefficient, α. This figure is after Orapunt [14]. . . . . . . . . . 9 ix List of Figures 1.5 A representative optical transition, from an occupied electronic state in the valence band to an unoccupied electronic state in the conduction band. . . . . . . . . . . . . . . . . . . 11 2.1 The possible optical transitions within c-Si and a-Si. This figure is after Jackson et al. [16]. . . . . . . . . . . . . . . . . 22 2.2 The occupancy function, f (E), as a function of energy, E, for various selections of temperature. The Fermi energy level, EF , is the reference energy level, i.e., EF = 0 eV. . . . . . . . . . . 23 2.3 A schematic representation of a three-dimensional potential energy well. This figure is from Thevaril [13]. . . . . . . . . . . 32 2.4 The valence band DOS function, Nv (E), and the conduction band DOS function, Nc (E), plotted as a function of energy, E, determined through Eqs. (2.23) and (2.24), respectively, assuming the nominal DOS modeling parameter selections, Nvo = Nco = 2 × 1022 cm−3 eV−3/2 , Ev = 0 eV, and Ec = 2 eV. These values are representative of the case of a-Si. . . . . 36 x List of Figures 2.5 The JDOS function, J(E), plotted as a function of the photon energy, E, determined through Eq. (2.25), assuming the nominal DOS modeling parameter selections, Nvo = Nco = 2 × 1022 cm−3 eV−3/2 , Ev = 0 eV, and Ec = 2 eV. These values are representative of the case of a-Si. This plot is depicted on a linear scale. . . . . . . . . . . . . . . . . . . . . . . . . . 37 2.6 The JDOS function, J(E), plotted as a function of the photon energy, E, determined through Eq. (2.25), assuming the nominal DOS modeling parameter selections, Nvo = Nco = 2 × 1022 cm−3 eV−3/2 , Ev = 0 eV, and Ec = 2 eV. These values are representative of the case of a-Si. This plot is depicted on a logarithmic scale. . . . . . . . . . . . . . . . . . . . . . . 38 2.7 Contours of integration for Eq. (2.30). . . . . . . . . . . . . . 42 2.8 The three experimental data sets, i.e., EGE a-Si [30], PECVD a-Si I [31], and PECVD a-Si II [32], for the real part of the dielectric function, 1 (E), as a function of the photon energy, E. 51 xi List of Figures 2.9 The three experimental data sets, i.e., EGE a-Si [30], PECVD a-Si I [31], and PECVD a-Si II [32], for the imaginary part of the dielectric function, 2 (E), as a function of the photon energy, E. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52 3.1 A “best-fit” of the Jellison and Modine model with the results of experiment corresponding to the EGE a-Si expermental data set [30]. These experimental results are from Piller [30]. The fits to experiment are shown with the solid lines. The experimental data, EGE a-Si [30], is depicted with the solid points. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 102 3.2 A “best-fit” of the Jellison and Modine model with the results of experiment corresponding to the PECVD a-Si I experimental data set [31]. These experimental results are from Synowicki [31]. The fits to experiment are shown with the solid lines. The experimental data, PECVD a-Si I [31], is depicted with the solid points. . . . . . . . . . . . . . . . . . . . . . . . . . . 103 xii List of Figures 3.3 A “best-fit” of the Jellison and Modine model with the results of experiment corresponding to the PECVD a-Si II experimental data set [32]. These experimental results are from Ferluato et al. [32]. The fits to experiment are shown with the solid lines. The experimental data, PECVD a-Si II [32], is depicted with the solid points. . . . . . . . . . . . . . . . . . . . . . . . 104 4.1 A linear plot of the valence band DOS function as a function of energy, E, for various selections of γv . The DOS modeling parameters, Nvo and Ev , are nominally set to 2 ×1022 cm−3 eV−3/2 and 0.0 eV, respectively, for the purposes of this plot. This figure is after O’Leary [17]. . . . . . . . . . . . . . . . . . 111 4.2 A logarithmic plot of the valence band DOS function as a function of energy, E, for various selections of γv . The DOS modeling parameters, Nvo and Ev , are nominally set to 2 ×1022 cm−3 eV−3/2 and 0.0 eV, respectively, for the purposes of this plot. This figure is after O’Leary [17]. . . . . . . . . . . . . . . . . . 112 xiii List of Figures 4.3 A linear plot of the conduction band DOS function as a function of energy, E, for various selections of γc . The DOS modeling parameters, Nco and Ec , are nominally set to 2 ×1022 cm−3 eV−3/2 and 0.0 eV, respectively, for the purposes of this plot. This figure is after O’Leary [17]. . . . . . . . . . . . . . . . . . 113 4.4 A logarithmic plot of the conduction band DOS function as a function of energy, E, for various selections of γc . The DOS modeling parameters, Nco and Ec , are nominally set to 2 ×1022 cm−3 eV−3/2 and 0.0 eV, respectively, for the purposes of this plot. This figure is after O’Leary [17]. . . . . . . . . . . 114 4.5 A linear plot of the JDOS function as a function of the photon energy, E, for various selections of γv and γc . The DOS modeling parameters are nominally set to Nvo = Nco = 2 ×1022 cm−3 eV−3/2 , Ev = 0.0 eV, and Ec = 1.7 eV for the purposes of this plot. This figure is after O’Leary [17]. . . . . . . . . . . 116 xiv List of Figures 4.6 A logarithmic plot of the JDOS function as a function of the photon energy, E, for various selections of γv and γc . The DOS modeling parameters are nominally set to Nvo = Nco = 2 ×1022 cm−3 eV−3/2 , Ev = 0.0 eV, and Ec = 1.7 eV for the purposes of this plot. This figure is after O’Leary [17]. . . 117 4.7 A linear plot of the functional dependence of the square-root of the JDOS function on the photon energy, E, for various selections of γv and γc . The DOS modeling parameters are nominally set to Nvo = Nco = 2 ×1022 cm−3 eV−3/2 , Ev = 0.0 eV, and Ec = 1.7 eV for this plot. This figure is after O’Leary [17].118 4.8 A plot of the dependence of γv on γc . Experimental results from Sherman et al. [23], Rerbal et al. [34], Teidje et al. [35], and Winer and Ley [36] are depicted. A modeling result, obtained through a fit with experimental data, from O’Leary [17], is also depicted on this plot. This figure is after Thevaril and O’Leary [33]. . . . . . . . . . . . . . . . . . . . . . . . . . . . 121 xv List of Figures 4.9 The spectral dependence of the JDOS function. The a-Si JDOS results of Jackson et al. [16] are depicted with the solid points. The JDOS fit result, obtained by setting Nvo = Nco = 2.48 ×1022 cm−3 eV−3/2 , Ev = 0.0 eV, Ec = 1.68 eV, and γv = 40 meV, is depicted with the solid line. . . . . . . . . . . 122 4.10 R2 (E) as a function of the photon energy, E, determined using Eq. (4.7), depicted with the solid lines. For the purposes of this fit, R2o is set to 10 ˚ A2 and Ed is set to 3.4 eV. The experimental results of Jackson et al. [16] are depicted with the solid points. . . . . . . . . . . . . . . . . . . . . . . . . . . 124 4.11 The spectral dependence for the imaginary part of the dielectric function, 2 (E), determined using Eq. (2.17), where R2 (E) is determined using Eq. (4.7), on a logarithmic scale. The DOS modeling parameters are set to Nvo = Nco = 2.48 ×1022 cm−3 eV−3/2 , Ev = 0.0 eV, Ec = 1.68 eV, and γv = 40 meV for the purposes of this plot. The experimental results of Jackson et al. [16] are depicted with the solid points . . . . . . . . . . 126 xvi List of Figures 4.12 The spectral dependence of the imaginary part of the dielectric function, 2 (E), determined using Eq. (2.17), where R2 (E) is determined using Eq. (4.7), on a linear scale. The DOS modeling parameters are set to Nvo = Nco = 2.48 ×1022 cm−3 eV−3/2 , Ev = 0.0 eV, Ec = 1.68 eV, and γv = 40 meV for the purposes of this plot. The experimental results of Jackson et al. [16] are depicted with the solid points . . . . . . . . . . 127 4.13 R2 (E) as a function of the photon energy, E, determined using Eqs. (4.7) and (4.8), are depicted with the solid lines. For the A2 and Ed is set to 3.4 eV. purposes of this fit, R2o is set to 10 ˚ The asymptotic polynomial rate of attenuation, x, is set to 6.6 for the purposes of the Eq. (4.8) plot. The experimental results of Jackson et al. [16] are depicted with the solid points. 129 xvii List of Figures 4.14 The spectral dependence of the imaginary part of the dielectric function, 2 (E), determined using Eq. (2.17), where R2 (E) is determined using Eq. (4.8), on a logarithmic scale. The DOS modeling parameters are set to Nvo = Nco = 2.48 ×1022 cm−3 eV−3/2 , Ev = 0.0 eV, Ec = 1.68 eV, and γv = 40 meV for the purposes of this plot. The asymptotic polynomial rate of attenuation, x, is set to 6.6 for the purposes of this plot. The experimental results of Jackson et al. [16] are depicted with the solid points . . . . . . . . . . . . . . . . . . . . . . . 130 4.15 The spectral dependence of the imaginary part of the dielectric function, 2 (E), determined using Eq. (2.17), where R2 (E) is determined using Eq. (4.8), on a linear scale. The DOS modeling parameters are set to Nvo = Nco = 2.48 ×1022 cm−3 eV−3/2 , Ev = 0.0 eV, Ec = 1.68 eV, and γv = 40 meV for the purposes of this plot. The asymptotic polynomial rate of attenuation, x, is set to 6.6 for the purposes of this plot. The experimental results of Jackson et al. [16] are depicted with the solid points . . . . . . . . . . . . . . . . . . . . . . . 131 xviii List of Figures 4.16 A “best-fit” of the model results with the results of experiment corresponding to EGE a-Si [30]. These experimental results are from Piller [30]. The fits to experiment are shown with the solid lines. The experimental data, EGE a-Si [30], is depicted with the solid points. . . . . . . . . . . . . . . . . . . . . . . . 134 4.17 A “best-fit” of the model results with the results of experiment corresponding to PECVD a-Si I [31]. These experimental results are from Synowicki [31]. The fits to experiment are shown with the solid lines. The experimental data, PECVD a-Si I [31], is depicted with the solid points. . . . . . . . . . . 135 4.18 A “best-fit” of the model results with the results of experiment corresponding to PECVD a-Si II [32]. These experimental results are from Ferluato et al. [32]. The fits to experiment are shown with the solid lines. The experimental data, PECVD a-Si II [32], is depicted with the solid points. . . . . . . . . . . 136 xix List of Symbols 1 real part of the dielectric function 2 imaginary part of the dielectric function α optical absorption coefficient I(z) intensity of light at position z Io intensity of light at position z = 0 n(E) spectral dependence of the refractive index k(E) spectral dependence of the extinction coefficient q electron charge m electron mass V illuminated volume η polarization vector of the incident light E photon energy xx List of Symbols Ev valence band band edge EvT critical energy at which the exponential and square-root distributions interface in the valence band Ec conduction band band edge EcT critical energy at which the exponential and square-root distributions interface in the conduction band Eg energy gap γv breadth of the valence band tail γc breadth of the conduction band tail |v single-spin electronic state associated with the valence band |c single-spin electronic state associated with the conduction band P momentum operator P 2 (E) aggregate momentum matrix element squared average R2 (E) aggregate dipole matrix element squared average R the dipole operator ρA density of silicon atoms xxi List of Symbols EF Fermi energy level Nv (E) valence band density of states function Nc (E) conduction band density of states function Nvo valence band density of states prefactor Nco conduction band density of states prefactor J(E) joint density of states function E(t) time-dependent electric field J(t) time-dependent current density Σ(t) time-dependent conductivity Eo resonance energy of a bound electron at which the transition of the electron occurs between two atomic states C breadth of the region of anomalous dispersion A Lorentz oscillator strength Ed characteristic energy xxii List of Acronyms a-Si amorphous silicon c-Si crystalline silicon µc-Si microcrystalline silicon DOS density of states JDOS joint density of states EGE electron-gun evaporated PECVD plasma enhanced chemical vapor deposition xxiii Acknowledgements I am very thankful to my supervisor, Dr. Stephen Karrer O’Leary, for his guidance, full support, and constant encouragement throughout my thesis. I am grateful to him for keeping faith in me. His valuable advice always helped me to progress. I owe particular thanks Dr. Mario Beaudoin, Research Associate, Manager of the AMPELs Advanced Nanofabrication Facility, UBC, for guiding me so wonderfully during the lab training period and helping me for enlarging my vision of science. xxiv To my family xxv Chapter 1 Introduction For the past half-century, electronics has played an important role in shaping the course of human development. In advanced industrial economies, such as Canada’s, electronics has become a critical adjunct to modern life. Electronic devices are ubiquitous. They are found in houses, offices, recreational settings, vehicles, and appliances. They are used in communication systems, control systems, sensor systems, and in a plethora of other systems. With each passing year, new applications for electronic devices are conceived. Thus, the role that electronics plays in modern society continues to broaden. The transistor is the fundamental building block upon which the electronics revolution has been built. While the first microprocessors, developed in the early 1970s, were comprised of about a thousand transistors, at the time of writing, 2012, the standard personal computer, available in most offices and homes, has in excess of a billion transistors. In the near term future, this 1 Chapter 1. Introduction is only expected to exponentially increase, as it has for the past four decades; Moore [1] prophesized this exponential increase in the number of transistors on a processor in 1965, and, thus far, Moore’s prophesy, or Moore’s law, as this projection is now commonly referred to as, has held true. The electronics revolution found its genesis with the development of the first transistor in 1947 [2, 3]. Transistors are comprised of conductors, insulators, and semiconductors. The fabrication of the first transistor was achieved as a result of a detailed and quantitative understanding of the material properties of these different materials. While the properties of conductors and insulators were relatively well understood by 1947, the understanding of the material properties of semiconductors was more rudimentary; at the time, semiconductors were considered more of a laboratory curiosity. Since that time, however, fundamental progress in electronics has been achieved through continued improvements in the understanding of the material properties of the materials used within electronic devices. The next generation of electronic device will undoubtedly rely upon an even greater understanding of these material properties. Accordingly, interest in the material properties of conductors, insulators, and particularly semiconductors, remains intense. In conventional electronics, technological improvements have been achieved 2 Chapter 1. Introduction as a consequence of making the constituent electronic devices, i.e., the transistors, smaller, faster, cheaper, and more reliable. There is, however, another class of electronic device that requires size in order to properly function. Flat panel displays [4, 5] and scanners [6], acting at the interface between the human and electronic worlds, are two such examples. Digital x-ray imagers [7] and photovoltaic solar cells [8, 9] are two other examples. These electronic devices are commonly referred to as examples of large area electronic devices. Starting out as a niche field in the late 1960s, large area electronics now constitutes a sizeable and growing fraction of the overall electronics enterprise; see Figure 1.1 [10]. In conventional electronics, the focus is on producing electronic devices with sub-micron feature sizes. In contrast, for large area electronics, the focus is on depositing materials as thin-films over large substrates, i.e., of the order of a square-meter [11]. For large area electronics, the uniformity of the resultant thin-films and the expense of the deposition process are of paramount concern. Crystalline silicon (c-Si), the workhorse of conventional electronics, can not be deposited as a thin-film. Accordingly, alternative electronic materials must be employed instead. Microcrystalline silicon (µcSi) and amorphous silicon (a-Si) are two alternate silicon-based materials 3 Chapter 1. Introduction 300 Market ($ US Billions) 250 Semiconductors 200 150 100 Flat Panel Displays 50 0 1985 1990 1995 2000 2005 2010 Year Figure 1.1: Semiconductor and flat panel display shipments plotted as a function of date. This data was obtained from the 2005 Information Society Technologies proposal for advancement [10]. 4 Chapter 1. Introduction commonly used in the large area electronics field. The distinction between c-Si, µc-Si, and a-Si relates to the underlying distribution of atoms. In c-Si, the atoms are arranged in a periodic and ordered manner. In µc-Si, however, small crystallites, randomly oriented and separated by grain boundaries, are distributed throughout the volume. In a-Si, all residual crystalline order is absent, there being variations in the shortrange bonding lengths and bonding angles. In Figure 1.2 [12], a schematic two-dimensional representation of the distribution of atoms within hypothetical crystalline, microcrystalline, and amorphous semiconductors is presented. The disorder that is present, in both µc-Si and a-Si, plays an important role in shaping the resultant properties of these materials. While disorder is present within the atomic distributions of µc-Si and a-Si, there is still a considerable amount of order that is present in both of these materials. For both µc-Si and a-Si, most silicon atoms are bonded to four other silicon atoms, as in the case of c-Si. The bonding lengths and bonding angles, while exhibiting variations, are, on average, similar to those exhibited by c-Si. Accordingly, as the properties of a given material are primarily determined by its short-range order, many of the properties of these materials are similar to those of c-Si. The disorder that is present, 5 Chapter 1. Introduction Crystalline semiconductor structure Polycrystalline semiconductor structure Amorphous semiconductor structure Figure 1.2: A two-dimensional schematic representation of the atomic distributions associated with hypothetical crystalline, microcrystalline, and amorphous semiconductors. This figure is from Nguyen [12]. 6 Chapter 1. Introduction however, will lead to some important differences. These will be the focus of the subsequent analysis. When a material is exposed to light, it will respond. Consider, for example, a beam of light incident upon a semi-infinite slab of material, as seen in Figure 1.3 [13]. Ignoring reflections from the slab, if the material comprising the slab is absorbing, assuming that the material within the slab is homogeneous, it is found that the intensity of the light exponentially attenuates as it propagates through the absorbing slab. In particular, it can be shown that the intensity of the light within the slab I(z) = Io exp(−αz), (1.1) where I(z) denotes the intensity of the light at position z, Io represents the incident light intensity, and α denotes the optical absorption coefficient; it is assumed that z = 0 corresponds to the interface between the vacuum and the absorbing slab. Representative attenuations of the light intensity, for a number of selections of the optical absorption coefficient, α, are depicted in Figure 1.4 [14]. The absorption that occurs within a semiconductor arises as a consequence of optical transitions from the occupied electronic states in the valence band to the unoccupied electronic states in the conduction band. In 7 Chapter 1. Introduction Figure 1.3: A schematic representation of the light intensity as a function of the penetration depth, z, with light normally incident upon the material propagating from the left. This figure is from Thevaril [13]. 8 Chapter 1. Introduction 1 0.9 α=104 cm−1 Relative Intensity 0.8 0.7 0.6 α=105 cm−1 0.5 0.4 0.3 0.2 α=106 cm−1 0.1 0 0 0.2 0.4 0.6 0.8 1 Penetration Depth (µm) Figure 1.4: Variations in the relative intensity, I(z)/Io , with the penetration depth, z, for different selections of the optical absorption coefficient, α. This figure is after Orapunt [14]. 9 Chapter 1. Introduction Figure 1.5, a representative optical transition is depicted. The optical response of a material depends critically upon the availability of electronic states and on the probability of the optical transitions between such electronic states. The distributions of valence band and conduction band electronic states determines the range of possible optical transitions. The optical transition matrix elements, that couple the electronic states between which optical transitions occur, determine the probability of such optical transitions occurring. The optical response of a material may be characterized in terms of the spectral dependence of the optical absorption coefficient, α(E), where E denotes the photon energy of the photons within the beam of incident light. In a defect-free, disorderless semiconductor, this optical absorption coefficient terminates abruptly at the energy gap. Accordingly, measurements of the optical absorption spectrum allow one to determine the corresponding energy gap. In addition, the spectral dependence of the optical absorption coefficient beyond the energy gap provides information as to how the electronic states are distributed and on what the probability of the optical transitions between such states are. When disorder is introduced, however, the situation is more complicated, as there are now electronic states encroaching into the otherwise 10 Chapter 1. Introduction conduction band incoming photon ¯hω energy gap valence band Figure 1.5: A representative optical transition, from an occupied electronic state in the valence band to an unoccupied electronic state in the conduction band. 11 Chapter 1. Introduction empty energy gap region. Accordingly, the optical absorption spectrum no longer terminates abruptly at the energy gap. In fact, the definition of the energy gap itself becomes the subject of some controversy. While α(E) is often experimentally measured, it is the real and imaginary parts of the dielectric function, 1 (E) and 2 (E), respectively, that are the more usual focus of theoretical attention. The complex dielectric function may be expressed as the sum over its real and imaginary components, i.e., ˜(E) = where 1 (E) and 2 (E) 1 (E) + i 2 (E), (1.2) denote the real and imaginary components of this dielectric function, respectively, and i = √ −1. The optical absorption coeffi- cient, α(E), may be related to the imaginary part of the dielectric function, 2 (E), through the relationship α(E) = where E cn(E) 2 (E), (1.3) and c represent the reduced Planck’s constant and the speed of light in a vacuum, respectively, and where n(E) denotes the spectral dependence of the refractive index, where the complex refractive index, n ˜ (E) = n(E) + ik(E), (1.4) 12 Chapter 1. Introduction k(E) being the corresponding extinction coefficient. It is noted that the complex refractive index is related to the complex dielectric function through the relationship ˜(E) = n ˜ 2 (E). (1.5) From Eqs. (1.2), (1.4), and (1.5), one can thus conclude that 1 (E) = n2 (E) − k 2 (E), (1.6) = 2 n(E) k(E). (1.7) and 2 (E) These functions, α(E), 1 (E), 2 (E), n(E), and k(E), are collectively referred to as the optical functions, as they quantify the nature of the optical response of a given material. In this thesis, the optical response of amorphous semiconductors will be the focus of analysis. Following a presentation of the theoretical background required of this work, the analysis begins with a theoretical justification for the empirical model of Jellison and Modine [15], in which Kramers-Kronig consistent empirical expressions for 1 (E) and 2 (E) are formulated. Then, a model for the optical functions, i.e., 1 (E) and 2 (E), stemming from a model for the distributions of valence band and conduction 13 Chapter 1. Introduction band electronic states, is developed. Finally, these models will be employed in order to fit some relevant experimental data. Conclusions, regarding the suitability of these models, will then be drawn from these fits. The main new results presented in this thesis pertain to the theoretical justification of the Kramers-Kronig consistent empirical expressions for 1 (E) and 2 (E) corresponding to the empirical model of Jellison and Modine [15] and the new model for these optical functions stemming from a model for the distributions of valence band and conduction band electronic states. This thesis is organized in the following manner. In Chapter 2, the theoretical foundations for this analysis are provided. Then, in Chapter 3, a theoretical justification for the Kramers-Kronig consistent model of Jellison and Modine [15] is formulated. An alternate Kramers-Kronig model for the spectral dependence of 1 (E) and 2 (E), stemming from a model for the un- derlying distributions of valence band and conduction band electronic states, is then developed in Chapter 4. Finally, the conclusions of this analysis are presented in Chapter 5. 14 Chapter 2 Theoretical basis for analysis 2.1 The optical response of materials The optical response of a material is primarily shaped by the distributions of valence band and conduction band electronic states and by the magnitude of the optical transition matrix elements which couple the electronic states between which optical transitions occur. The distributions of valence band and conduction band electronic states determines what optical transitions can occur in a given material. The magnitude of the optical transition matrix elements, however, determine the probability of such optical transitions occurring. The primary goal of this analysis is the development of Kramers-Kronig consistent models for the spectral dependence of the real and imaginary components of the dielectric function, 1 (E) and 2 (E), re- spectively, for the case of amorphous semiconductors. In order to achieve this goal, the role that each of these factors plays in shaping the optical 15 2.1. The optical response of materials response of an amorphous semiconductor must be well understood. In this chapter, the theoretical foundations for the work presented in this thesis are provided. Initially, an expression for the spectral dependence of the imaginary part of the dielectric function, 2 (E), in terms of the valence band and conduction band energy levels and the magnitude of the corresponding optical transition matrix elements, is presented. Then, through the introduction of a number of suitable assumptions, a relationship between the valence band and conduction band distributions of electronic states and the spectral dependence of 2 (E) is developed. Two particular formulations for this re- lationship are developed: (1) a momentum matrix element formulation, and (2) a dipole matrix element formulation, these formulations being related to each other. For each of these formulations, an aggregate matrix element is employed. The joint density of states (JDOS) function, J(E), is introduced as a corollary to this analysis. The Kramers-Kronig relations are then presented, these relating the spectral dependence of the spectral dependence of 1 (E) and 2 (E), 2 (E) with 1 (E). 1 (E) with 2 (E) and A number of models, for both are then presented. Finally, three sets of experimental data, corresponding to different samples of a-Si, are presented, these experimental data sets being used in the subsequent analysis. 16 2.1. The optical response of materials This chapter is organized in the following manner. In Section 2.2, an expression for the spectral dependence of the imaginary part of the dielectric function, 2 (E), in terms of the valence band and conduction band energy levels and the magnitude of the corresponding optical transition matrix elements, is presented. Then, in Section 2.3, three critical assumptions are introduced, these allowing one to relate the spectral dependence of 2 (E) with that of J(E). In Section 2.4, a momentum matrix element formulation for the relationship between the imaginary part of the dielectric function, 2 (E), and the distributions of the valence band and conduction band elec- tronic states is developed, the JDOS function and the aggregate momentum matrix element being introduced as a byproduct of this particular formulation. A dipole matrix element formulation for the relationship between the imaginary part of the dielectric function, 2 (E), and the distributions of the valence band and conduction band electronic states is then developed in Section 2.5, the aggregate dipole matrix element being proposed as a byproduct of this analysis. The free electron model, for the distributions of electronic states, is then presented in Section 2.6. The spectral dependence of this JDOS function, assuming free electron distributions of valence band and conduction band electronic states, is then presented in Section 2.7. 17 2.2. The spectral dependence of the imaginary part of the dielectric function The Kramers-Kronig relationships between the real and imaginary parts of the dielectric function, 1 (E) and 2 (E), respectively, are discussed in Sec- tion 2.8, the time-symmetry requirements demanded of these relations being featured. A number of models for the spectral dependence of the optical functions, 1 (E) and 2 (E), are then presented in Section 2.9. Finally, ex- perimental a-Si data sets, corresponding to 1 (E) and 2 (E), are presented in Section 2.10, these data sets being used in the subsequent analysis. 2.2 The spectral dependence of the imaginary part of the dielectric function Optical transitions occur from occupied electronic states in the valence band (initial states) to unoccupied electronic states in the conduction band (final states). Within the framework of a one-electron picture, assuming linear optical response, Jackson et al. [16] assert that the imaginary part of the dielectric function, 2 (E) = 2πq mE 2 2 V v,c |η · P v,c |2 δ(Ec − Ev − E), (2.1) where q denotes the electron charge, m represents the electron mass, V is the illuminated volume, η is the polarization vector of the incident light, 18 2.3. Three critical assumptions and E is the photon energy; the δ(·) terms refer to the Dirac delta function. Ev and Ec denote the energy levels of representative valence band and conduction band electronic states, |v and |c , respectively. P v,c represents the momentum matrix element that couples these particular electronic states in an optical transition process; in quantum mechanics, this matrix element, P v,c = c|P |v , where P denotes the momentum operator. The |η · P v,c |2 operation refers to taking the complex modulus squared of the complex number η · P v,c . The summation in Eq. (2.1) is performed over all of the occupied valence band and unoccupied conduction band single-spin electronic states. 2.3 Three critical assumptions Eq. (2.1) expresses the spectral dependence of the imaginary part of the dielectric function, 2 (E), in terms of the valence band and conduction band energy levels and the magnitude of the corresponding momentum matrix elements. This expression applies to both the crystalline and amorphous cases. In order to relate 2 (E) with the underlying distributions of valence band and conduction band electronic states for the specific case of amorphous semiconductors, three critical assumptions are introduced into the analysis: 19 2.3. Three critical assumptions (1) it is assumed that the incident light is unpolarized, (2) it is assumed that the disorder that is present in an amorphous semiconductor completely relaxes the momentum conservation rules, and (3) it is assumed that zerotemperature statistics apply. The consequences of each approximation is outlined subsequently. The absence of polarization simplifies Eq. (2.1). With random polarization, assuming that all of the significant momentum matrix elements are of the same order of magnitude, the average value of |η · P v,c |2 reduces to 1 |Pv,c |2 , 3 where |Pv,c | denotes the amplitude of the momentum matrix element. As a result, Eq. (2.1) reduces to 2 2 (E) = (2πq) 2 2 m2 E 2 3V v,c |Pv,c |2 δ(Ec − Ev − E). (2.2) In a crystalline semiconductor, momentum must be conserved in an optical transition. In an amorphous semiconductor, however, the momentum conservation rules are relaxed. That is, optical transitions can occur from every occupied state (of a particular spin) in the valence band to every unoccupied state (of the same spin) in the conduction band; spin-flips are not permitted during an optical transition. As a consequence, in a crystalline semiconductor, the |Pv,c | values vary greatly, depending upon whether momentum is conserved or not, while in an amorphous semiconductor, |Pv,c | is 20 2.3. Three critical assumptions essentially constant for all cases. This greatly simplifies the analysis. The contrast between c-Si and a-Si is instructive. Consider the case of N silicon atoms. For the case of c-Si, from a given single-spin state in the valence band, four optical transitions are possible; there are N 2 possible values of crystal momentum for the case of c-Si. For the case of a-Si, however, 2N optical transitions are possible. These possible optical transitions are depicted in Figure 2.1. If the overall “optical strength” of c-Si and a-Si is similar, it seems likely that the |Pv,c | matrix elements associated with the small number of possible c-Si optical transitions will be of a much greater magnitude than that of their more numerous a-Si counterparts. Later on in this analysis, a normalized matrix element is introduced in order to allow for a direct comparison between the c-Si and a-Si matrix elements, where the weighting factor is directly proportional to the ratio between the number of allowed optical transitions for the c-Si and a-Si cases. The final assumption that is introduced into this analysis is that of zerotemperature statistics. In Figure 2.2, the occupancy function, f (E) = 1 , F 1 + exp( E−E ) KBT (2.3) is depicted as a function of an energy, E, for a variety of different temperature selections, EF denoting the Fermi energy level. It is noted that as the 21 Energy 2.3. Three critical assumptions EF a−Si 2ρA V 4 c−Si k Figure 2.1: The possible optical transitions within c-Si and a-Si. This figure is after Jackson et al. [16]. 22 2.3. Three critical assumptions 1 Probability of Occupation 300 K 0K 0.75 800 K 0.5 0.25 EF 0 −0.3 −0.2 −0.1 0 0.1 0.2 0.3 Electron Energy (eV) Figure 2.2: The occupancy function, f (E), as a function of energy, E, for various selections of temperature. The Fermi energy level, EF , is the reference energy level, i.e., EF = 0 eV. 23 2.4. A momentum matrix element formulation temperature goes to zero that this function, sometimes referred to as the Fermi-Dirac function, reduces to a step function, where the step occurs at the Fermi energy level, EF . For the case of semiconductors, considering that the Fermi energy level, EF , is typically in the middle of the energy gap, this implies that at zero temperature, all of the valence band electronic states are fully occupied and all of the conduction band electronic states are completely unoccupied. This assumption simplifies the subsequent analysis. 2.4 A momentum matrix element formulation It will be assumed that the overall “optical strength” of c-Si is identical to that of a-Si. As was discussed earlier, if this is the case it would seem reasonable to expect that the small number of non-zero c-Si momentum matrix elements will be larger in magnitude than their much more numerous a-Si counterparts. In order to allow for a direct contrast between these matrix elements, one could conceive of normalizing the a-Si momentum matrix elements by a factor that is proportional to the ratio between the number of allowed optical transitions, i.e., 2N . 4 Letting ρA denote the density of silicon 24 2.4. A momentum matrix element formulation atoms within a-Si, this normalization factor becomes 2ρA V 4 In order to allow for a direct relationship between . 2 (E) and the dis- tributions of valence band and conduction band electronic states, one can introduce an aggregate momentum matrix element P 2 (E) ≡ where the ρA V 2 ρA V 2 v,c |Pv,c |2 δ(Ec − Ev − E) , v,c δ(Ec − Ev − E) (2.4) prefactor allows for a direct comparison between the momen- tum matrix elements associated with c-Si and a-Si. From Eqs. (2.2) and (2.4), it can now be shown that 2 2 (E) = (2πq) The term v,c 2 m2 E 2 2 3V 2 ρA V P 2 (E) v,c δ(Ec − Ev − E). (2.5) δ(Ec − Ev − E) provides for the number of possible optical transitions between the occupied valence band and unoccupied conduction band electronic states separated by the energy E. Defining the JDOS function, J(E) ≡ 4 V2 v,c δ(Ec − Ev − E), (2.6) from Eq. (2.5), it can now be seen that 2 (E) = (2πq)2 2 2 V2 m2 E 2 3V 4 2 ρA V P 2 (E) J(E). (2.7) 25 2.4. A momentum matrix element formulation Thus, the spectral dependence of 2 (E) reduces to a product of a prefactor, the reciprocal of E 2 , the aggregate momentum matrix element, and the JDOS function, J(E). The valence band density of states (DOS) function, Nv (E), where Nv (E)∆E represents the number of valence band electronic states between [E, E +∆E], per unit volume, may be expressed as the sum over all of the valence band single-spin electronic states divided by the volume, i.e., Nv (E) = 2 V v δ(E − Ev ), (2.8) where the factor of 2 arises as a consequence of the fact that each electronic state is a single-spin state, i.e., there are two possible spins. Similarly, the conduction band DOS function, Nc (E), where Nc (E)∆E represents the number of conduction band electronic states between [E, E + ∆E], per unit volume, may be expressed as the sum over all of the conduction band singlespin electronic states divided by the volume, i.e., Nc (E) = 2 V c δ(E − Ec ), (2.9) where the factor of 2 arises as a consequence of the fact that each electronic state is a single-spin state. Assuming that all valence band states are fully occupied and that all 26 2.4. A momentum matrix element formulation conduction band states are completely unoccupied, i.e., the zero-temperature assumption, it may be shown that the JDOS function may be expressed as an integral over the product of the valence band and conduction band DOS functions, i.e., ∞ Nv (ξ)Nc (E + ξ)dξ. J(E) = (2.10) −∞ From Eqs. (2.5) and (2.6), one can thus conclude that 2 (E) = 2 (2πq)2 P 2 (E) J(E). 2 2 3ρA m E (2.11) For a wide variety of amorphous semiconductor experimental analyzes, it has been assumed that the aggregate momentum matrix element, P 2 (E), exhibits a constant spectral dependence with respect to energy, E. This is referred to as the constant momentum matrix element assumption. This assumption, which found its origins in the analysis at Tauc et al. [18], was widely accepted by researchers in the field until the early 1980s. Despite this acceptance, it lacks a theoretical justification. Moreover, Cody et al. [19] demonstrated that experimental results interpreted within the framework of this assumption are fundamentally inconsistent. Accordingly, other models for the spectral dependence of 2 (E) are sought. 27 2.5. A dipole matrix element formulation 2.5 A dipole matrix element formulation Cody et al. [19] provided experimental data which challenged the validity of the constant momentum matrix element assumption. For the purposes of their study, Cody et al. [19] prepared a large number of a-Si samples under identical conditions, but with different film thicknesses. They then measured the spectral dependence of the optical functions using transmission and reflectance measurements on these samples. Cody et al. [19] noted that while fits to this data, determined assuming a constant momentum matrix element, are satisfactory, they become greatly improved if one instead assumes a constant dipole matrix element. Moreover, such an assumption removes a film thickness dependence artifact that was noted by Cody et al. [19] within the framework of the constant momentum matrix element assumption. Accordingly, it is instructive to recast our formulation in terms of the spectral dependence of the dipole matrix element. Taking into account the experimental evidence provided by Cody et al. [19], Jackson et al. [16] recast Eq. (2.1) into a form in which the matrix elements are instead expressed in terms of the dipole matrix elements as opposed to the momentum matrix elements, i.e., Rv,c = c|R|v , where R denotes the 28 2.5. A dipole matrix element formulation dipole operator. From the commutator relations, Jackson et al. [16] find that v,c mE |η · P v,c |2 δ(Ec − Ev − E) = 2 v,c |η · Rv,c |2 δ(Ec − Ev − E). (2.12) As before, if the incident light is unpolarized, the average value of |η · Rv,c |2 reduces to 31 |Rv,c |2 , where |Rv,c | denotes the amplitude of the dipole matrix element. Accordingly, Jackson et al. [16] conclude that 2 (E) = (2πq)2 2 3V v,c |Rv,c |2 δ(Ec − Ev − E). (2.13) Introducing an aggregate dipole matrix element, defined and normalized in a manner similar to that introduced earlier for the case of the momentum matrix formalism, Jackson et al. [16] proposed that, R2 (E) ≡ ρA V 2 v,c |Rv,c |2 δ(Ec − Ev − E) , v,c δ(Ec − Ev − E) (2.14) where, as before, ρA denotes the atomic density. For this definition, it is noted that Eq. (2.13) may be represented as 2 (E) = (2πq)2 2 3V 2 ρA V R2 (E) v,c δ(Ec − Ev − E). (2.15) From Eqs. (2.6) and (2.15), one can thus conclude that 2 (E) = (2πq)2 2 R (E) J(E). 3ρA (2.16) 29 2.5. A dipole matrix element formulation That is, the spectral dependence of the imaginary part of the dielectric function, 2 (E), is proportional to the product of a prefactor, R2 (E), and J(E). Given that the density of the silicon atoms within a-Si is typically around 4.4 ×1022 cm−3 [17], it can thus be seen that, for the specific case of a-Si, Eq. (2.16) may be re-expressed as 2 (E) = 4.3 × 10−45 R2 (E) J(E), (2.17) where R2 (E) is in units of ˚ A2 and J(E) is in the units of cm−6 eV−1 [16, 17]. In light of the experimental results of Cody et al. [19], the spectral dependence of R2 (E) is often assumed to be constant. This is referred to as the constant dipole matrix element assumption. It is noted, from Eqs. (2.11) and (2.16), that the relationship between the aggregate momentum matrix element squared average and the aggregate dipole matrix element squared average is found to be P 2 (E) = m2 E 2 2 R2 (E). (2.18) Eq. (2.18) allows one to relate the one formulation for the spectral dependence of the imaginary part of the dielectric function, 2 (E), with the other. 30 2.6. Free electron model of the DOS functions 2.6 Free electron model of the DOS functions The spectral dependence of the JDOS function depends upon a con- volution over the valence band and conduction band DOS functions; recall Eq. (2.10). Accordingly, the determination of these functions is of paramount concern. While the exact form of these DOS functions depends critically upon the exact potential distribution, the free electron DOS function provides a useful limiting result that is often employed for theoretical benchmarking purposes. Consider the electronic states within a cubic box, of dimensions L×L×L, surrounded by an infinitely high potential barrier, as seen in Figure 2.3. It will be assumed that the potential within the box itself is nil. From quantum mechanics, for steady-state conditions, the wavefunctions associated with the corresponding bound electron states may be determined from Schr¨odinger’s equation, i.e., 2 2 2m Ψ(r) + V (r)Ψ(r) = EΨ(r), (2.19) where , m, V (r), and E represent the reduced Planck’s constant, the mass of the electron, the potential energy, and the electron energy, respectively. 31 2.6. Free electron model of the DOS functions Figure 2.3: A schematic representation of a three-dimensional potential energy well. This figure is from Thevaril [13]. 32 2.6. Free electron model of the DOS functions 2 denotes the mathematical operator ∂ ∂x2 + ∂ ∂y 2 + ∂ ∂z 2 in three-dimensions. It may be shown that Ψ(r) = Ψnx ,ny ,nz (x, y, z) = 2 L 3/2 sin πny πnz πnx x sin y sin z , L L L (2.20) where nx , ny , and nz are positive integers, associated with the electron motion in the different directions; these are the quantum numbers. The energy of the corresponding electronic states may be found by substituting the obtained wavefunction, i.e., Eq. (2.20), back into the Schr¨odinger’s equation, i.e., 2 Enx ,ny ,nz = 2m πnx L 2 + πny L 2 + πnz L 2 , (2.21) where Enx ,ny ,nz represents the energy associated with the nx , ny , and nz wavefunction, Ψnx ,ny ,nz (x, y, z), the labeling integers, nx , ny , and nz , exactly corresponding to the aforementioned quantum numbers. In the continuum limit, i.e., for E sufficiently large, the corresponding DOS function may be shown to be √ m3/2 √ 2 2 3 E, E ≥ 0 π . N (E) = 0, E < 0 (2.22) 33 2.7. The spectral dependence of the JDOS function Eq. (2.22) is often referred to the free electron DOS model; within the well, the potential is nil, i.e., the electrons are “free”. This form for the DOS function is often assumed in analyzes of amorphous semiconductors. 2.7 The spectral dependence of the JDOS function It is often useful to employ the free-electron model for the valence band and conduction band DOS functions in the evaluation of the JDOS function. This model suggests that the distribution of electronic states exhibits a square-root functional dependence on the energy above the conduction band minimum and below the valence band maximum. Following O’Leary [17], one can set the valence band DOS function 0, E > Ev , Nv (E) = Nvo Ev − E, E ≤ Ev (2.23) 34 2.7. The spectral dependence of the JDOS function and the conduction band DOS function Nco E − Ec , E ≥ Ec , Nc (E) = 0, E < Ec (2.24) where Nvo and Nco represent the valence band and conduction band DOS prefactors, respectively, and Ev and Ec represent the valence band and conduction band band edges; these DOS functions are also known as the Tauc DOS functions. The resultant distributions of electronic states, for this specific case, are depicted in Figure 2.4, nominal values of Nvo , Nco , Ev , and Ec being employed for the purposes of this analysis; in particular, Nvo and Nco are both set to 2 ×1022 cm−3 eV−3/2 , Ev is set to 0 eV, and Ec is set to 2 eV, these values being representative of the case of a-Si. Through the use of Eq. (2.10), it may be shown that J(E) = Nvo Nco π (E − Eg )2 , E ≥ Eg 8 , (2.25) 0, E < Eg where Eg ≡ Ec − Ev represents the energy gap of the material. The spectral dependence of the corresponding JDOS function is depicted in Figures 2.5 and 2.6. This is known as the Tauc JDOS function after the pioneering work 35 2.7. The spectral dependence of the JDOS function 15 DOS (× 1021 cm−3eV−1) Nv(E) Nc(E) 10 5 0 −0.5 0 0.5 1 1.5 Energy (eV) 2 2.5 Figure 2.4: The valence band DOS function, Nv (E), and the conduction band DOS function, Nc (E), plotted as a function of energy, E, determined through Eqs. (2.23) and (2.24), respectively, assuming the nominal DOS modeling parameter selections, Nvo = Nco = 2 × 1022 cm−3 eV−3/2 , Ev = 0 eV, and Ec = 2 eV. These values are representative of the case of a-Si. 36 2.7. The spectral dependence of the JDOS function JDOS × 1044 (cm−6eV−1) 1.6 1.2 Nvo = Nco = 2 × 1022 cm−3eV−3/2 Eg = 2 eV 0.8 0.4 0 1.2 1.6 2 2.4 Photon Energy (eV) 2.8 Figure 2.5: The JDOS function, J(E), plotted as a function of the photon energy, E, determined through Eq. (2.25), assuming the nominal DOS modeling parameter selections, Nvo = Nco = 2 × 1022 cm−3 eV−3/2 , Ev = 0 eV, and Ec = 2 eV. These values are representative of the case of a-Si. This plot is depicted on a linear scale. 37 2.7. The spectral dependence of the JDOS function 44 10 42 Nvo = Nco = 2 × 1022 cm−3eV−3/2 Eg = 2 eV JDOS (cm−6eV−1) 10 40 10 38 10 36 10 34 10 0.5 1 1.5 2 Photon Energy (eV) 2.5 3 Figure 2.6: The JDOS function, J(E), plotted as a function of the photon energy, E, determined through Eq. (2.25), assuming the nominal DOS modeling parameter selections, Nvo = Nco = 2 × 1022 cm−3 eV−3/2 , Ev = 0 eV, and Ec = 2 eV. These values are representative of the case of a-Si. This plot is depicted on a logarithmic scale. 38 2.8. Kramers-Kronig relations of Tauc et al. [18]. 2.8 Kramers-Kronig relations As a consequence of causality, there is a relationship between the spectral dependence of the real and imaginary parts of the dielectric function, and 2 (E), 1 (E) respectively. Consider a time-dependent electric field, E(t), ap- plied to a solid that produces a current density, J(t), in response. The variations with respect to position of these quantities are neglected. If only weak fields are considered, then a general linear relationship between E and J exists, where ∞ J(t) = −∞ Σ(t − t ) E(t )dt , (2.26) Σ(t) denoting the time-dependent conductivity. Through a Fourier transform on Eq. (2.26), one obtains J (ω) = σ(ω)ξ(ω), (2.27) 39 2.8. Kramers-Kronig relations where J (ω), ξ(ω), and σ(ω) are the Fourier transforms of J, E, and Σ, respectively, i.e., ∞ J (ω) = J(t)eiωt dt, (2.28a) E(t)eiωt dt, (2.28b) Σ(t)eiωt dt. (2.28c) −∞ ∞ ξ(ω) = −∞ and ∞ σ(ω) = −∞ It is also noted that the inverse Fourier transforms become J(t) = E(t) = 1 2π 1 2π ∞ −∞ ∞ J (ω)e−iωt dω, (2.29a) ξ(ω)e−iωt dω, (2.29b) σ(ω)e−iωt dω. (2.29c) −∞ and Σ(t) = 1 2π ∞ −∞ In order to further the analysis, two critical assumptions have to be made: (1) Σ(t) is zero for negative values of its argument, and (2) J(t) = 0 for t < t0 , where t0 is the earliest value of t for which E(t ) = 0. The first assumption leads to the fact that σ(ω) is an analytic function in the upper half of the complex ω plane. According to Cauchy’s theorem, the integral in Eq. (2.29c) should vanish for σ(ω) being analytic in the upper half plane. It may thus 40 2.8. Kramers-Kronig relations be shown that σ(ω0 ) = 1 2πi C1 σ(ω) dω. ω − ω0 (2.30) At first, the contour of integration of Eq. (2.30) may be chosen to be a small circle, C1 , about ω0 , as shown in Figure 2.7. It is known that the distortion of C1 into a larger contour of integration, C2 , does not change the result; see Figure 2.7, i.e., σ(ω0 ) = 1 2πi C1 σ(ω) 1 dω = ω − ω0 2πi C2 σ(ω) dω. ω − ω0 (2.31) Assuming that σ(ω) → 0 for |ω| → ∞, then only contributions from the real axis play a role in determining the integral in Eq. (2.30). Considering that ωr is the frequency along the real axis, then it may be shown that 1 σ(ω0 ) = 2πi ∞ −∞ σ(ωr ) dωr . ωr − ω0 (2.32) It is noted that ω0 may be expressed as ω0 = ω + iω . (2.33) Since lim ω →0 1 ωr − ω − iω = P 1 + iπδ(ωr − ω), ωr − ω (2.34) 41 2.8. Kramers-Kronig relations Imaginary ωi C2 ωo C1 Real ωr Figure 2.7: Contours of integration for Eq. (2.30). 42 2.8. Kramers-Kronig relations where P denotes the principal value, from Eqs. (2.32) and (2.34), it is seen that σ(ω) = ∞ 1 P πi −∞ σ(ωr ) dωr . ωr − ω (2.35) From Eq. (2.35), it is seen that σ(ω) is a complex quantity, i.e., σ(ω) = σ1 (ω) + iσ2 (ω). (2.36) Thus, it can be seen that σ1 (ω) = 1 P π ∞ −∞ σ2 (ωr ) dωr , ωr − ω (2.37a) and 1 σ2 (ω) = − P π ∞ −∞ σ1 (ωr ) dωr . ωr − ω (2.37b) That is, σ1 (ω) and σ2 (ω) are Hilbert transforms of each other. As a real E(t) produces a real J(t), then the time-dependent conductivity, σ(t), must also be real. This requires that [σ1 (ωr ) + iσ2 (ωr )]∗ = σ(−ωr ). (2.38) As a result, σ1 (ωr ) = σ1 (−ωr ), (2.39a) σ2 (ωr ) = −σ2 (−ωr ). (2.39b) and 43 2.8. Kramers-Kronig relations So, from Eq. (2.37a), it is seen that 1 σ1 (ω) = P π 1 P = π 2 = P π ∞ σ2 (ωr ) σ2 (ωr ) dωr + dωr , ωr − ω −∞ ωr − ω 0 ∞ 1 1 σ2 (ωr ) + dωr , ωr − ω ωr + ω 0 ∞ ωr σ2 (ωr ) dωr . ωr2 − ω 2 0 0 (2.40) Similarly, σ2 (ω) = − ∞ 2ω P π 0 σ1 (ωr ) dωr . ωr2 − ω 2 (2.41) Eqs. (2.40) and (2.41) constitute the Kramers-Kronig relations, applied to conductivity. It is also possible to develop Kramers-Kronig relations for the real and imaginary parts of the dielectric function. Since (ω) = 0 +i σ(ω) , ω (2.42) it may be shown that 1 (ω) =1− 2 (ω) = σ2 (ω) , 0ω (2.43a) and σ1 (ω) . 0ω (2.43b) 44 2.8. Kramers-Kronig relations Thus, 1 (ω) = 1 + ∞ 2 P π ωr 2 (ωr ) dωr , ωr2 − ω 2 (2.44a) ωr2 [1 − 1 (ωr )] dωr . ωr2 − ω 2 (2.44b) 0 and 2 (ω) = ∞ 2 P πω 0 So, the real part of the dielectric function, part of the dielectric function, 2 (E), 1 (E), is related to the imaginary through the Kramers-Kronig transfor- mation, 1 (E) where 1∞ ∞ Eg u 2 (u) du, u2 − E 2 (2.45) should ideally be set to unity. Similarly, the imaginary part of the dielectric function, 1 (E), = 2 P 1∞ + π 2 (E), is related to the real part of the dielectric function, through the Kramers-Kronig transformation 2 (E) = 2 P πE ∞ 0 u2 [1 − 1 (u)] du. u2 − E 2 (2.46) The conditions developed in Eqs. (2.39a) and (2.39b) mandate that 1 (E) = 1 (−E), 2 (E) = − 2 (−E). (2.47a) and (2.47b) 45 2.9. Review of models These are referred to as the time-symmetry requirements mandated of the Kramers-Kronig relations. 2.9 Review of models Within the framework of the momentum matrix element and the dipole matrix element formulations, the imaginary part of the dielectric function, 2 (E), may be expressed as the product of the corresponding aggregate op- tical transition matrix element, the JDOS function, and some other function of the photon energy. The optical transition matrix elements play an important role in shaping the spectral dependence of the imaginary part of the dielectric function, 2 (E). Models for the spectral dependence of the imag- inary part of the dielectric function, 2 (E), have been developed assuming a constant momentum matrix element, i.e., a momentum matrix element that is independent of the photon energy, E. Examples include the model of Forouhi and Bloomer [20], the model of McGahan et al. [21], and the model of Jellison and Modine [15]. This assumption of a constant momentum matrix element was first suggested by Tauc et al. [18]. In fact, Tauc et al. [18] suggested that the zero- 46 2.9. Review of models intercept of the plot of √ αE as a function of E, where α represents the corresponding optical absorption coefficient, allows one to define the “optical bandgap” of these materials if the aggregate momentum matrix element is constant and the valence band and conduction band DOS functions are square-root in form, i.e., of the form presented in Eqs. (2.23) and (2.24). Based upon a linear response formalism for an expression for 2 (E) 2 (E), Tauc et al. [18] derived that characterizes the absorption edge of amorphous semiconductors. They concluded that 2 (E) ∝ (E − Eg )2 , E2 (2.48) for E > Eg , where Eg denotes the corresponding “optical bandgap”. Eq. (2.48) may be used to describe the interband optical transition mechanism for the optical response of such a material, and thus, Recalling Eq. (1.7), i.e., 2 (E) 2 (E) = 0 for E < Eg . = 2n(E)k(E), where n(E) represents the spectral dependence of the refractive index and k(E) denotes the spectral dependence of the extinction coefficient, Forouhi and Bloomer [20] suggest an expression for the extinction coefficient kF B (E) = A(E − Eg )2 , E 2 − BE + C (2.49) where A, B, and C are treated as empirical fitting parameters. The term 47 2.9. Review of models (E − Eg )2 in Eq. (2.49) represents the assumed form for the JDOS function of the material itself, and arises from the assumed square-root valence band and conduction band DOS functions; recall Eq. (2.25). This formulation was then used in order to fit experimental optical functions corresponding to amorphous diamond-like carbon films [15, 21, 22, 24]. Unfortunately, as was pointed out by McGahan and Woollam [21, 22], the Forouhi and Bloomer [20] model is unable to describe the spectral dependence of these optical functions correctly. In particular, the Forouhi and Bloomer [20] formulation provides for negative band gaps, a clearly unphysical result [15, 22]. McGahan et al. [21] made some modifications to the Forouhi and Bloomer [20] model, and this was shown to successfully fit several n and k experimental data sets. Unfortunately, the Forouhi and Bloomer [20] and McGahan et al. [21] formulations have several fundamental problems, these being nicely summarized by Jellison and Modine [15]. They are as listed below: • While experiment indicates that 2 (E) → 0 as E → ∞, both the models of Forouhi and Bloomer [20] and McGahan et al. [21] instead approach non-zero values for 2 (E) as E → ∞. • Both the models of Forouhi and Bloomer [20] and McGahan et al. [21] fail to satisfy the Kramers-Kronig time-symmetry requirements, 48 2.10. Experimental results i.e., Eqs. (2.47a) and (2.47b), which mandate that k(−E) = −k(E). (2.50) Jellison and Modine [15] developed a new parameterized model for the optical functions associated with an amorphous semiconductor that does not have the problems listed above. In particular, they developed a Tauc-Lorentz expression for the imaginary part of the dielectric function, 2 (E), based on the Tauc JDOS function [18, 25], i.e., Eq. (2.25), and the Lorentz oscillator [26]. They then determined the corresponding spectral dependence of the real part of the dielectric function, 1 (E), using a Kramers-Kronig trans- formation, i.e., Eq. (2.45). As with the models of Forouhi and Bloomer [20] and McGahan et al. [21], the model of Jellison and Modine [15] is cast within the framework of the constant momentum matrix element assumption. 2.10 Experimental results In the subsequent analysis, models for the spectral dependence of the real and imaginary components of the dielectric function, 1 (E) and 2 (E), respec- tively, are fit to a-Si optical function experimental data, and from these fits, the “best” model parameters are determined, i.e., the modeling parameters 49 2.10. Experimental results that produce the “best” fits with the corresponding experimental data. The experimental a-Si optical function data considered is from Piller [30], Synowicki [31], and Ferluato et al. [32], these experimental results corresponding to samples of electron-gun evaporated (EGE) a-Si, and plasma enhanced chemical vapor deposition (PECVD) prepared a-Si; the experimental data from Piller [30] was obtained from a sample of EGE a-Si, the experimental data from Synowicki [31] was obtained from a sample of PECVD a-Si, and the experimental data from Ferluato et al. [32] was obtained from a sample of PECVD a-Si, these samples being referred to as EGE a-Si, PECVD a-Si I, and PECVD a-Si II, respectively. These three sets of experimental data are shown in Figures 2.8 and 2.9, the spectral dependence of the real part of the dielectric function being depicted in Figure 2.8 and the spectral dependence of the imaginary part of the dielectric function being depicted in Figure 2.9. The EGE a-Si [30] experimental data set was collected from an electron-gun evaporated a-Si sample, with silicon as the substrate [30]. For the case of the PECVD a-Si I [31] experimental data sets, fused silica was used as the substrate [31]. For the case of the PECVD a-Si II [32] experimental data sets, glass was used as the substrate [32]. 50 2.10. Experimental results 50 EGE a−Si [30] PECVD a−Si I [31] PECVD a−Si II [32] 40 ε1 30 20 10 0 −10 1 2 3 4 5 Photon Energy (eV) 6 7 Figure 2.8: The three experimental data sets, i.e., EGE a-Si [30], PECVD aSi I [31], and PECVD a-Si II [32], for the real part of the dielectric function, 1 (E), as a function of the photon energy, E. 51 2.10. Experimental results 50 EGE a−Si [30] PECVD a−Si I [31] PECVD a−Si II [32] 40 ε2 30 20 10 0 −10 1 2 3 4 5 Photon Energy (eV) 6 7 Figure 2.9: The three experimental data sets, i.e., EGE a-Si [30], PECVD a-Si I [31], and PECVD a-Si II [32], for the imaginary part of the dielectric function, 2 (E), as a function of the photon energy, E. 52 Chapter 3 The model of Jellison and Modine 3.1 Empirical models for the optical functions associated with an amorphous semiconductor For the purposes of materials characterization and device design, empirical models for the optical functions of an amorphous semiconductor have been sought. Forouhi and Bloomer [20] and McGahan et al. [21] have proposed empirical models for such functions. Unfortunately, these models have a number of shortcomings. While experiment indicates that as E → ∞, 2 (E) → 0, both the models of Forouhi and Bloomer [20] and McGahan et 53 3.1. Empirical models for the optical functions associated with an amorphous semiconductor al. [21] yield non-zero values for 2 (E) as E → ∞. The models of Forouhi and Bloomer [20] and McGahan et al. [21] have the additional problem of being obtained using unphysical expressions for posed models for 2 (E) 2 (E). In particular, the pro- do not terminate exact at the energy gap. Finally, both the models of Forouhi and Bloomer [20] and McGahan et al. [21] fail to satisfy the time-symmetry requirements mandated of the Kramers-Kronig relations; recall Eqs. (2.47a) and (2.47b). In 1996, Jellison and Modine [15] developed a new parameterized model for the optical functions associated with an amorphous semiconductor that remedies some of the problems found with these previous models. In this approach, a physically-based model for the imaginary part of the dielectric function, 2 (E), is proposed. This model includes a Tauc JDOS function. i.e., Eq. (2.25), and a Lorentzian momentum matrix element. Then, through the use of a Kramers-Kronig transformation on 2 (E), i.e., Eq. (2.45), the spectral dependence of the real part of the dielectric function, 1 (E), is de- termined. This chapter is organized in the following manner. In Section 3.2, the Lorentz model is described. In Section 3.3, Jellison and Modine’s model for the spectral dependence of the real and imaginary components of the 54 3.2. The Lorentz model dielectric function, 1 (E) and 2 (E), respectively, is presented [15]. Then, in Section 3.4, an explicit determination of the spectral dependence of from the use of a Kramers-Kronig transformation on 2 (E), 1 (E) i.e., Eq. (2.45), is provided. Finally, a fit with the results of experiment is presented in Section 3.5. 3.2 The Lorentz model According to the Lorentz model, an atom with electrons bounded by the nucleus may be considered in the same way as a small mass bounded to a large mass by a spring. The motion of an electron bound to the nucleus may be described by m dr d2 r + mC + mω02 r = −qEloc , 2 dt dt (3.1) where m denotes the electron mass, q represents the electron charge, and Eloc is the driving local electric field acting on the electron. It may be shown that r˜ = − 1 eEloc , 2 m (ω0 − ω 2 ) − iCω (3.2) 55 3.2. The Lorentz model and the induced dipole moment is ˜ = p 1 e2 Eloc . 2 m (ω0 − ω 2 ) − iCω (3.3) ˜ and Eloc if the displacement r is There is a linear relationship between p assumed to be sufficiently small. That is, ˜ = α p ˜ (ω)Eloc , (3.4) where α ˜ (ω) is the frequency-dependent atomic polarizability. From Eqs. (3.3) and (3.4), it may be shown that α ˜ (ω) = e2 1 . 2 m (ω0 − ω 2 ) − iCω (3.5) The complex dielectric function, ˜, may be expressed in terms of the complex polarizability, α ˜ , by ˜ = 1 + 4πN α ˜, (3.6) where N denotes the number of atoms per unit volume. From Eqs. (3.5) and (3.6), it may be shown that ˜ = 1+ 1 4πN e2 . 2 m (ω0 − ω 2 ) − iCω (3.7) Since the complex dielectric function, ˜, may be expressed as the sum over its real and imaginary components, i.e., ˜ = 1 + i 2, (3.8) 56 3.3. Model for the optical functions where 1 and 2 denote the real and imaginary components of the dielectric function, respectively, it can be seen that 1 =1+ 2 = ω02 − ω 2 4πN e2 , m (ω02 − ω 2 )2 + C 2 ω 2 (3.9a) and 4πN e2 Cω , 2 m (ω0 − ω 2 )2 + C 2 ω 2 (3.9b) Jellison and Modine [15] used this simple Lorentz model in their proposed parameterized model for the spectral dependence of the imaginary part of the dielectric function, 3.3 2, shown in the subsequent analysis. Model for the optical functions Within the framework of the momentum matrix element formulation, i.e., Eq. (2.1), Jellison and Modine [15] propose a Tauc-Lorentz expression for the imaginary part of the dielectric function, 2 (E), i.e., 2 (E) is essentially proportional to the product of a Tauc JDOS function and a Lorentz oscillator. Specifically, Jellison and Modine [15] suggest that 0, E ≤ Eg (E) = , 2T L AEo CE (E − Eg )2 , E > Eg (E 2 − Eo2 )2 + C 2 E 2 E2 (3.10) 57 3.3. Model for the optical functions where the four modeling parameters, Eo , C, Eg , and A, are in units of energy. Here, Eo denotes the resonance energy of a bound electron at which the transition of the electron occurs between two atomic states, A represents the oscillator strength, C is the breadth of the region of anomalous dispersion, and Eg is the energy gap. Jellison and Modine [15] determine the real part of the dielectric function, 1 (E), using a Kramers-Kronig transformation on 1 (E) = 2 P 1∞ + π ∞ Eg 2 (E). u 2 (u) du, u2 − E 2 That is, (3.11) where P represents the Cauchy principal part of the integral. An additional fitting parameter, 1∞ , is included in order to obtain the desired fits; in the actual Kramers-Kronig relations, 1∞ is unity. This integral may be solved, in closed-form, to yield Eg2 + Eo2 + αEg 1 A C aln ln 1 (E) = 1∞ + 2 π ζ 4 αEo Eg2 + Eo2 − αEg 2Eg + α −2Eg + α A aa tan − 4 π − tan−1 + tan−1 πζ Eo C C 2 2 2(γ − Eg ) AEo +2 Eg E 2 − γ 2 π + 2 tan−1 4 παζ αC 2 2 AEo C E + Eg |Eg − E| − ln πζ 4 E (Eg + E) 2AEo C [(Eg + E)|Eg − E|] , + E ln (3.12) g πζ 4 2 2 2 2 2 (E − E ) + E C ) o g g 58 3.3. Model for the optical functions where aln = (Eg2 − Eo2 )E 2 + Eg2 C 2 − Eo2 (Eo2 + 3Eg2 ), (3.13a) aa tan = (E 2 − Eo2 )(Eo2 + Eg2 ) + Eg2 C 2 , (3.13b) ζ 4 = [(E 2 − Eo2 )2 + C 2 E 2 ], (3.13c) χ= 4Eo2 − C 2 , (3.13d) γ= Eo2 − C2 , 2 (3.13e) and Eo , C, Eg , and A being as defined earlier. The two closed-form expressions for the optical functions, 2 (E), 1 (E) and respectively, i.e., Eqs. (3.10) and (3.12), provide a Kramers-Kronig consistent model for the optical functions of an amorphous semiconductor. This model avoids many of the problems associated with the models of Foroughi and Bloomer [20] and McGahan et al. [21]. That is, zero as E → ∞, the model for 2 (E) 2 (E) approaches is physically based, and finally, these optical function expressions obey the time-symmetry requirements mandated of the Kramers-Kronig relations, i.e., Eqs. (2.47a) and (2.47b). Jellison and Modine [15] then employed this model in order to fit a variety of experimental optical function data for the spectral dependence of 1 (E) and 2 (E), and 59 3.4. Derivation drew a number of conclusions into the character of the materials considered. 3.4 Derivation 3.4.1 Definitions From Eqs. (3.10) and (3.11), it is seen that 1 (E) = 2AEo C P 1∞ + π ∞ Eg (u − Eg )2 du.(3.14) (u2 − E 2 )[(u2 − Eo2 )2 + C 2 u2 ] Letting f (u) = (u − Eg )2 , (u2 − E 2 )[(u2 − Eo2 )2 + C 2 u2 ] (3.15) the analysis starts with a partial fraction expansion of this polynomial expression in terms of the roots of its denominator. As the denominator is a 6th -order polynomial, it has six roots, two of which may be seen immediately, from inspection, to be E and −E, respectively. Accordingly, the corresponding partial fraction expansion may be expressed as f (u) = A1 A2 A3 A4 A5 A6 + + + + + ,(3.16) u + E u − E u − r1 u − r2 u − r3 u − r4 where A1 , A2 , A3 , A4 , A5 , and A6 are the co-efficients of the partial fraction expansions, r1 , r2 , r3 , and r4 being the other four, as yet unknown, roots of the denominator of f (u). 60 3.4. Derivation 3.4.2 Root calculation through expansion of the denominator The roots of f (u) may be determined by setting (u2 − Eo2 )2 + C 2 u2 = 0. (3.17) Expanding Eq. (3.17) fully, one can see that u4 − u2 (2Eo2 − C 2 ) + Eo4 = 0. Factoring the u2 term, this equation reduces to u 2 2Eo2 − C 2 ± = (2Eo2 − C 2 )2 − 4Eo4 , 2 which leads to u2 = Eo2 − C2 C ±i 2 2 4Eo2 − C 2 . Noting that 2 u = 4Eo2 − C 2 C ±i 2 2 2 , it is seen that u = ± 4Eo2 − C 2 C ±i 2 2 . 61 3.4. Derivation Letting a= 4Eo2 − C 2 , 2 and b= C , 2 it is seen that the four roots of f (u) are r1 = a + ib, (3.18a) r2 = a − ib, (3.18b) r3 = −a + ib, (3.18c) r4 = −a − ib. (3.18d) and Through inspection, the following relationships are seen amongst these roots, i.e., r2 = r1∗ , r3 = −r1∗ , r4 = −r1 , r1 = −r3∗ , r2 = −r3 , 62 3.4. Derivation and r4 = r3∗ . 3.4.3 Calculation of the coefficients The coefficients, A1 , A2 , A3 , A4 , A5 , and A6 , may be determined using the usual approach for a partial fraction expansion. That is, for the partial fraction expansion, f (z) = i Ai , z − zi (3.19) if all the roots, zi , are distinct, then it may be seen that Ai = (z − zi )f (z)|z=zi . (3.20) Noting that all of the roots associated with the denominator of f (u) are distinct, employing the approach suggested in Eq. (3.20), the coefficients A1 , A2 , A3 , A4 , A5 , and A6 will be determined. For the case of these coefficients, a polar notation is employed, i.e., the imaginary number p + iq = p2 + q 2 exp i tan−1 q p , (3.21) may instead be expressed as p + iq = p2 + q 2 ∠ tan−1 q , p (3.22) 63 3.4. Derivation where this last expression corresponds to polar notation. That is, every complex number is expressed in this particular polar form. An alternate form for this polar notation abs (p + iq) = p2 + q 2 , (3.23) q , p (3.24) and arg (p + iq) = tan−1 may be employed. Determination of A1 A1 = f (u)(u + E)|u=−E , (−E − Eg )2 , (−E − E)[(E 2 − Eo2 )2 + C 2 E 2 ] (E + Eg )2 = − . 2E[(E 2 − Eo2 )2 + C 2 E 2 ] = (3.25) Determination of A2 A2 = f (u)(u − E)|u=E , (E − Eg )2 , (E + E)[(E 2 − Eo2 )2 + C 2 E 2 ] (E − Eg )2 = . 2E[(E 2 − Eo2 )2 + C 2 E 2 ] = (3.26) 64 3.4. Derivation Determination of A3 A3 = f (u)(u − r1 )|u=r1 , (r1 − Eg )2 , (r12 − E 2 )(r1 − r2 )(r1 − r3 )(r1 − r4 ) (r1 − Eg )2 , = (r12 − E 2 )(r1 − r1∗ )(r1 + r1∗ )(r1 + r1 ) (r1 − Eg )2 = . (r12 − E 2 )(r1 − r1∗ )(r1 + r1∗ )2r1 = (3.27) Now, by substituting r1 = a + ib, the factors in the numerator and the denominator of Eq. (3.27) may be determined and expressed in polar form. For the numerator, this yields (r1 − Eg )2 = (a + ib − Eg )2 , = (a − Eg )2 − b2 + i2b(a − Eg ), = ((a − Eg )2 − b2 )2 + 4b2 (a − Eg )2 ∠ tan−1 = (a − Eg )2 + b2 ∠ tan−1 2b(a − Eg ) . (a − Eg )2 − b2 2b(a − Eg ) , (a − Eg )2 − b2 The polar form for the other factors in the denominator of Eq. (3.27) may now be expressed as r12 − E 2 = (a2 − b2 − E 2 )2 + 4a2 b2 ∠ tan−1 2(r1 − r1∗ )(r1 + r1∗ ) = 8ab∠ π , 2 2ab , a2 − b 2 − E 2 and 65 3.4. Derivation r1 = √ b a2 + b2 ∠tan−1 . a So, the coefficient, A3 , may be expressed in polar form as [(a − Eg )2 + b2 ]∠ tan−1 2b(a−Eg ) (a−Eg )2 −b2 × (a2 − b2 − E 2 )2 + 4a2 b2 ∠ tan−1 a2 −b2ab 2 −E 2 1 √ , 2 2 a + b ∠ tan−1 ab 8ab∠ π2 [(a − Eg )2 + b2 ] √ = 8ab[ (a2 − b2 − E 2 )2 + 4a2 b2 ][ a2 + b2 ] 2b(a − Eg ) 2ab b π ×∠ tan−1 − tan−1 2 − tan−1 − . 2 2 2 2 (a − Eg ) − b a −b −E a 2 A3 = Thus, it is seen that abs(A3 ) = [(a − Eg )2 + b2 ] √ , 8ab[ (a2 − b2 − E 2 )2 + 4a2 b2 ][ a2 + b2 ] (3.28a) and that arg(A3 ) = − π 2b(a − Eg ) 2ab + tan−1 − tan−1 2 2 2 2 (a − Eg ) − b a − b2 − E 2 − tan−1 b . a (3.28b) 66 3.4. Derivation Determination of A4 A4 = f (u)(u − r2 )|u=r2 , (r2 − Eg )2 , (r22 − E 2 )(r2 − r1 )(r2 − r3 )(r2 − r4 ) (r1∗ − Eg )2 , = (r1∗2 − E 2 )(r1∗ − r1 )(r1 + r1∗ )(r1∗ + r1∗ ) (r1∗ − Eg )2 = . 2r1∗ (r1∗2 − E 2 )(r1∗ − r1 )(r1 + r1∗ ) = (3.29) Now, by substituting r1∗ = a − ib, the factors in the numerator and the denominator of Eq. (3.29) may be determined and expressed in polar form. For the numerator, this yields (r1∗ − Eg )2 = (a − ib − Eg )2 , = ((a − Eg )2 − b2 )2 + 4b2 (a − Eg )2 ∠ − tan−1 = [(a − Eg )2 + b2 ]∠ − tan−1 2b(a − Eg ) . (a − Eg )2 − b2 2b(a − Eg ) , (a − Eg )2 − b2 The polar form for the other factors in the denominator of Eq. (3.29) may be expressed as r1∗2 − E 2 = (a2 − b2 − E 2 )2 + 4a2 b2 ∠ − tan−1 2(r1∗ − r1 )(r1 + r1∗ ) = 8ab∠ −π , 2 a2 2ab , − b2 − E 2 and r1∗ = √ b a2 + b2 ∠ − tan−1 . a 67 3.4. Derivation Thus, abs(A4 ) = [(a − Eg )2 + b2 ] √ , 8ab[ (a2 − b2 − E 2 )2 + 4a2 b2 ][ a2 + b2 ] (3.30a) and arg(A4 ) = 2b(a − Eg ) 2ab π − tan−1 + tan−1 2 2 2 2 (a − Eg ) − b a − b2 − E 2 + tan−1 b . a (3.30b) Through inspection, it is seen that abs(A4 ) = abs(A3 ) and that arg(A4 ) = − arg(A3 ), i.e., A3 = A∗4 . Determination of A5 A5 = f (u)(u − r3 )|u=r3 , (r3 − Eg )2 , = (r32 − E 2 )(r3 − r1 )(r3 − r2 )(r3 − r4 ) (r3 − Eg )2 , = (r32 − E 2 )(r3 − r3∗ )(r3 + r3∗ )(r3 + r3 ) (r3 − Eg )2 = . (r32 − E 2 )(r3 − r3∗ )(r3 + r3∗ )2r3 (3.31) For this equation, for the numerator, (r3 − Eg )2 = (−a + ib − Eg )2 , = ((a + Eg )2 − b2 )2 + 4b2 (a + Eg )2 ∠ − tan−1 = [(a + Eg )2 + b2 ]∠ − tan−1 2b(a + Eg ) . (a + Eg )2 − b2 2b(a + Eg ) , (a + Eg )2 − b2 68 3.4. Derivation Similarly, for the denominator, r32 − E 2 = (a2 − b2 − E 2 )2 + 4a2 b2 ∠ − tan−1 2(r3 − r3∗ )(r3 + r3∗ ) = 8ab∠ −π , 2 a2 2ab , − b2 − E 2 and r3 = √ a2 + b2 ∠ tan−1 b . −a Thus, abs(A5 ) = [(a + Eg )2 + b2 ] √ , 8ab[ (a2 − b2 − E 2 )2 + 4a2 b2 ][ a2 + b2 ] (3.32a) and arg(A5 ) = π 2b(a + Eg ) 2ab − tan−1 + tan−1 2 2 2 2 (a + Eg ) − b a − b2 − E 2 − tan−1 b . −a (3.32b) Determination of A6 A6 = f (u)(u − r4 )|u=r4 , (r4 − Eg )2 , (r42 − E 2 )(r4 − r1 )(r4 − r2 )(r4 − r3 ) (r3∗ − Eg )2 , = (r3∗2 − E 2 )(r3∗ − r3 )(r3∗ + r3∗ )(r3∗ + r3 ) (r3∗ − Eg )2 = . (r3∗2 − E 2 )(r3∗ − r3 )(r3 + r3∗ )2r3∗ = (3.33) 69 3.4. Derivation Here, the numerator of Eq. (3.33) yields (r3∗ − Eg )2 = (−a − ib − Eg )2 , = [(a + Eg )2 + b2 ]∠ tan−1 2b(a + Eg ) . (a + Eg )2 − b2 Similarly, the terms in the denominator of Eq. (3.33) yield r3∗2 − E 2 = (a2 − b2 − E 2 )2 + 4a2 b2 ∠ tan−1 2(r3∗ − r3 )(r3 + r3∗ ) = 8ab∠ π , 2 a2 2ab , − b2 − E 2 and r3∗ = √ a2 + b2 ∠ − tan−1 b . −a Thus, abs(A6 ) = [(a + Eg )2 + b2 ] √ , 8ab[ (a2 − b2 − E 2 )2 + 4a2 b2 ][ a2 + b2 ] (3.34a) and arg(A6 ) = − π 2b(a + Eg ) 2ab + tan−1 − tan−1 2 2 2 2 (a + Eg ) − b a − b2 − E 2 + tan−1 b . −a (3.34b) Through inspection, it is also seen that abs(A5 ) = abs(A6 ) and that arg(A5 ) = − arg(A6 ), i.e., A5 = A∗6 . 70 3.4. Derivation 3.4.4 Grouping the terms in the partial fraction expression In Section 3.4.3, it was seen that A3 = A∗4 and that A5 = A∗6 . These relations may be used in order to group the terms associated with the coefficients in the partial fraction expansion. Grouping the terms associated with A3 and A4 The terms associated with A3 and A4 in the partial fraction expansion may be combined as follows: A4 abs(A3 )∠(arg(A3 )) abs(A4 )∠arg(A4 )) A3 + = + , u − r1 u − r2 u − a − ib u − a + ib abs(A3 )∠(arg(A3 )) = b (u − a)2 + b2 ∠ − tan−1 u−a abs(A3 )∠ − [arg(A3 )] + , b (u − a)2 + b2 ∠ tan−1 u−a abs(A3 ) = + (u − a)2 + b2 abs(A3 ) (u − a)2 + ∠ arg(A3 ) + tan−1 b2 b u−a ∠ − arg(A3 ) + tan−1 b u−a . 71 3.4. Derivation Letting, x = arg(A3 ) + tan−1 b u−a , and abs(A3 ) k= (u − a)2 + b2 , it may be seen that A4 A3 + = k∠x + k∠ − x, u − r1 u − r2 = k(eix + e−ix ), = 2k cos x, = 2 = 2 −2 abs(A3 ) (u − a)2 + b2 abs(A3 ) (u − + abs(A3 ) a)2 b2 cos arg(A3 ) + tan−1 cos [arg(A3 )] (u − a)2 + b2 sin [arg(A3 )] b u−a , u−a (u − a)2 + b2 b (u − a)2 + b2 , u−a (u − a)2 + b2 b −2abs(A3 ) sin [arg(A3 )] . (u − a)2 + b2 = 2abs(A3 ) cos [arg(A3 )] 72 3.4. Derivation Letting c3 = 2 abs(A3 ) cos[arg(A3 )], and c4 = −2 abs(A3 ) sin[arg(A3 )], it is also seen that, A4 A3 u−a b + = c3 + c4 . 2 2 u − r1 u − r2 (u − a) + b (u − a)2 + b2 (3.35) 73 3.4. Derivation Grouping the terms associated with A5 and A6 The terms associated with A5 and A6 in the partial fraction expansion may be combined as follows: A6 abs(A5 )∠ [arg(A5 )] abs(A5 )∠ − [arg(A5 )] A5 + = + , u − r3 u − r4 u + a − ib u + a + ib abs(A5 )∠ [arg(A5 )] = b (u + a)2 + b2 ∠ − tan−1 u+a abs(A5 )∠ − [arg(A5 )] + , b (u − a)2 + b2 ∠ tan−1 u+a abs(A5 ) = (u + a)2 + b2 abs(A5 ) a)2 = 2 (u + + abs(A5 ) b2 −2 + ∠ − arg(A5 ) + tan−1 b u+a , cos arg(A5 ) + tan−1 b u+a , (u + a)2 + b2 abs(A5 ) = 2 b u+a ∠ arg(A5 ) + tan−1 (u + a)2 + b2 abs(A5 ) cos [arg(A5 )] (u + a)2 + b2 sin [arg(A5 )] u+a (u + a)2 + b2 b (u + a)2 + b2 , u+a (u + a)2 + b2 b −2abs(A5 ) sin[arg(A5 )] . (u + a)2 + b2 = 2abs(A5 ) cos [arg(A5 )] Letting c5 = 2 abs(A5 ) cos[arg(A5 )], and 74 3.4. Derivation c6 = −2 abs(A5 ) sin[arg(A5 )]. it is seen that, A5 A6 u+a b + = c5 + c6 . 2 2 u − r3 u − r4 (u + a) + b (u + a)2 + b2 3.4.5 (3.36) Summary of the partial fraction expansion From Eqs. (3.16), (3.35), and (3.36), it is thus seen that A2 A3 A4 A5 A6 A1 + + + + + , u + E u − E u − r1 u − r2 u − r3 u − r4 A2 u−a b A1 + + c3 + c = 4 u+E u−E (u − a)2 + b2 (u − a)2 + b2 u+a b +c5 + c6 , (3.37) 2 2 (u + a) + b (u + a)2 + b2 f (u) = where A1 = −(E + Eg )2 , 2E[(E 2 − Eo2 )2 + C 2 E 2 ] (E − Eg )2 A2 = , 2E[(E 2 − Eo2 )2 + C 2 E 2 ] (3.38a) (3.38b) c3 = 2 abs(A3 ) cos[arg(A3 )], (3.38c) c4 = −2 abs(A3 ) sin[arg(A3 )], (3.38d) c5 = 2 abs(A5 ) cos[arg(A5 )], (3.38e) 75 3.4. Derivation and c6 = −2 abs(A5 ) sin[arg(A5 )]. 3.4.6 (3.38f) Evaluation of the constants c3 , c4 , c5 , and c6 Calculation of c3 From Eq. (3.38c), it is seen that c3 = 2abs(A3 ) cos[arg(A3 )], (3.39) where, from Eq. (3.28b), π 2b(a − Eg ) 2ab + tan−1 − tan−1 2 2 2 2 (a − Eg ) − b a − b2 − E 2 b . − tan−1 a arg(A3 ) = − Letting x = tan−1 y = tan−1 2b(a − Eg ) , (a − Eg )2 − b2 a2 −2ab , − b2 − E 2 and z = tan−1 −b , a 76 3.4. Derivation it is seen that arg(A3 ) = − π + x + y + z. 2 Noting that cos[arg(A3 )] = cos(− π + x + y + z) = sin(x + y + z), 2 it may be seen that cos[arg(A3 )] = cos x cos y sin z − sin x sin y sin z + sin x cos y cos z + cos x sin y cos z. Letting x = tan−1 2b(a − Eg ) , (a − Eg )2 − b2 one observes that cos x = (a − Eg )2 − b2 , (a − Eg )2 + b2 sin x = 2b(a − Eg ) . (a − Eg )2 + b2 and Similarly, letting y = tan−1 a2 −2ab , − b2 − E 2 77 3.4. Derivation one observes that cos y = a2 − b 2 − E 2 , −2ab . (a2 − b2 − E 2 )2 + (2ab)2 and sin y = (a2 − b2 − E 2 )2 + (2ab)2 Letting −b , a z = tan−1 it is seen that cos z = a (a2 + b2 , and sin z = −b (a2 + b2 . Substituting the results for abs(A3 ), as well as the cosine and sine of x, y, and z, respectively, into Eq. (3.39) yields c3 = 2 abs(A3 ) cos[arg(A3 )], 78 3.4. Derivation c3 = 2 × 8ab(a2 + b2 )[(a2 − b2 − E 2 )2 + 4a2 b2 ] 2 2 2 2 2 −b × [((a − Eg ) − b )(a − b − E ) − 2b(a − Eg )(−2ab)] + a × [2b(a − Eg )(a2 − b2 − E 2 ) + ((a − Eg )2 − b2 )(−2ab)] = 8ab(a2 + b2 )[(a2 2 × − b2 − E 2 )2 + 4a2 b2 ] 2 2 2 , 2 2b(a − Eg )[a(a − b − E ) − 2ab ] − ((a − Eg )2 − b2 )[b(a2 − b2 − E 2 ) + 2a2 b] . Calculation of c4 From Eq. (3.38d), it is seen that c4 = −2 abs(A3 ) sin[arg(A3 )], (3.40) where, from Eq. (3.28b), π 2b(a − Eg ) 2ab + tan−1 − tan−1 2 2 2 2 (a − Eg ) − b a − b2 − E 2 b − tan−1 . a arg(A3 ) = − 79 3.4. Derivation Similarly, it is seen that sin[arg(A3 )] = sin(− π + x + y + z) = − cos(x + y + z) 2 = −[cos x cos y cos z − sin x sin y cos z − sin x cos y sin z − cos x sin y sin z]. Thus, c4 = −2 abs(A3 ) sin[arg(A3 )], = 8ab(a2 + b2 )[(a2 2 × − b2 − E 2 )2 + 4a2 b2 ] 2 2 2 2 2 a × [((a − Eg ) − b )(a − b − E ) − 2b(a − Eg )(−2ab)] − b × [−2b(a − Eg )(a2 − b2 − E 2 ) − ((a − Eg )2 − b2 )(−2ab)] = 1 × 4ab(a2 + b2 )[(a2 − b2 − E 2 )2 + 4a2 b2 ] 2b(a − Eg )[b(a2 − b2 − E 2 ) + 2a2 b] + +((a − Eg )2 − b2 )[a(a2 − b2 − E 2 ) − 2ab2 ] , . 80 3.4. Derivation Calculation of c5 From Eq. (3.38e), c5 = 2 abs(A5 ) cos[arg(A5 )], (3.41) where, from Eq. (3.32b), arg(A5 ) = 2b(a + Eg ) 2ab b π − tan−1 + tan−1 2 − tan−1 . 2 2 2 2 2 (a + Eg ) − b a −b −E −a Letting x = tan−1 y = tan−1 −2b(a + Eg ) , (a + Eg )2 − b2 a2 2ab , − b2 − E 2 and z = tan−1 −b , −a it is seen that arg(A5 ) = π + x + y + z. 2 Noting that cos [arg(A5 )] = cos( π + x + y + z) = − sin(x + y + z), 2 81 3.4. Derivation it may also be seen that cos [arg(A5 )] = −[cos x cos y sin z − sin x sin y sin z + sin x cos y cos z + cos x sin y cos z]. Letting x = tan−1 −2b(a + Eg ) , (a + Eg )2 − b2 one observes that cos x = (a + Eg )2 − b2 , (a + Eg )2 + b2 sin x = −2b(a + Eg ) . (a + Eg )2 + b2 and Similarly, letting y = tan−1 a2 2ab , − b2 − E 2 one observes that cos y = a2 − b 2 − E 2 (a2 − b2 − E 2 )2 + (2ab)2 , and sin y = 2ab (a2 − b2 − E 2 )2 + (2ab)2 . 82 3.4. Derivation Letting −b , −a z = tan−1 it is seen that cos z = −a , −b . (a2 + b2 and sin z = (a2 + b2 Substituting the results of abs(A5 ), as well as the cosine and sine of x, y, and z into Eq. (3.41), yields c5 = 2 abs(A5 ) cos[arg(A5 )], = 8ab(a2 + b2 )[(a2 −2 × − b2 − E 2 )2 + 4a2 b2 ] 2 2 2 2 2 −b × [((a + Eg ) − b )(a − b − E ) − (−2b(a + Eg )) 2ab] − a × [−2b(a + Eg )(a2 − b2 − E 2 ) + ((a + Eg )2 − b2 ) 2ab] = 8ab(a2 + b2 )[(a2 2 × − b2 − E 2 )2 + 4a2 b2 ] 2 2 2 , 2 −2b(a + Eg )[a(a − b − E ) − 2ab ] + ((a + Eg )2 − b2 )[b(a2 − b2 − E 2 ) + 2a2 b] . 83 3.4. Derivation Calculation of c6 From Eq. (3.38f), c6 = −2 abs(A5 ) sin[arg(A5 )], (3.42) where, from Eq. (3.32b), arg(A5 ) = 2b(a + Eg ) 2ab π − tan−1 + tan−1 2 2 2 2 (a + Eg ) − b a − b2 − E 2 b − tan−1 . −a Similarly, sin[arg(A5 )] = sin( π + x + y + z) = cos(x + y + z), 2 = cos x cos y cos z − sin x sin y cos z − sin x cos y sin z − cos x sin y sin z. Thus, it may be seen that c6 = −2 abs(A5 ) sin[arg(A5 )], = 8ab(a2 + b2 )[(a2 −2 × − b2 − E 2 )2 + 4a2 b2 ] 2 2 2 2 2 −a × [((a + Eg ) − b )(a − b − E ) − (−2b(a + Eg )) 2ab] − b × [−2b(a + Eg )(a2 − b2 − E 2 ) − ((a + Eg )2 − b2 ) 2ab] , 84 3.4. Derivation 1 × 4ab(a2 + b2 )[(a2 − b2 − E 2 )2 + 4a2 b2 ] = 3.4.7 2b(a + Eg )[b(a2 − b2 − E 2 ) + 2a2 b] + ((a + Eg )2 − b2 )[a(a2 − b2 − E 2 ) − 2ab2 ] . Main integration So, from Eq.(3.14), it is seen that 2AEo C P π 2AEo C P 1∞ + π 1 (E) = ∞ 1∞ + = Eg ∞ (u − Eg )2 du, (u2 − E 2 )[(u2 − Eo2 )2 + C 2 u2 ] f (u) du. (3.43) Eg From the partial fraction expansion of Eq. (3.37), this yields ∞ ∞ f (u) du = Eg ∞ A1 A2 du + du Eg u − E Eg u + E ∞ u−a + c3 du (u − a)2 + b2 Eg ∞ b + c4 du (u − a)2 + b2 Eg ∞ u+a + c5 du (u + a)2 + b2 Eg ∞ b + c6 du, (u + a)2 + b2 Eg 85 3.4. Derivation In order to avoid improper integrals, the integrals will all be taken to t, and then t will go to infinity, i.e., t → ∞. Thus, it can be seen that t t A1 du + u+E f (u) du = =⇒ Eg Eg t c3 + Eg t c4 + Eg t + Eg Eg A2 du u−E u−a du (u − a)2 + b2 b du (u − a)2 + b2 c5 u+a du (u + a)2 + b2 c6 b du. (u + a)2 + b2 Eg t + t (3.44) Sequentially integrating the terms in Eq. (3.44), it is thus seen that t Eg A1 du = A1 u+E t Eg 1 du, u+E = A1 [ln(u + E)]tEg , = A1 [ln(t + E)] − A1 ln(Eg + E), (3.45) and t Eg A2 du = A2 lim →0 u−E E− Eg 1 du + lim →0 u−E t E+ 1 du , u−E t = A2 lim ln |u − E|E− Eg + lim ln |u − E|E+ , →0 →0 = A2 lim ln | − | − A2 [ln |Eg − E| + ln |t − E|] →0 −A2 lim ln | | , →0 = A2 [ln |t − E|] − A2 ln |Eg − E|. (3.46) 86 3.4. Derivation Now, for t c3 Eg u−a du, (u − a)2 + b2 letting w = (u − a)2 + b2 , =⇒ dw = 2(u − a) du, it may be seen that (t−a)2 +b2 t c3 u−a du = c3 (u − a)2 + b2 2 Eg = (Eg −a)2 +b2 c3 ln[(t − a)2 + b2 ] 2 c3 − ln[(Eg − a)2 + b2 ]. 2 1 dw, w (3.47) (3.48) Again, for t c4 Eg b du, (u − a)2 + b2 letting (u − a) = b tan θ, =⇒ du = b sec2 θ dθ, 87 3.4. Derivation it is seen that t c4 Eg b du = c4 (u − a)2 + b2 t Eg b sec2 θ dθ, b2 tan2 θ + b2 = c4 [θ]tEg , −1 = c4 tan = c4 u−a b π − tan−1 2 t , Eg Eg − a b . [as t → ∞] (3.49) Similarly, t u+a c5 du = ln[(t + a)2 + b2 ] 2 2 (u + a) + b 2 Eg c5 − ln[(Eg + a)2 + b2 ], 2 t b π Eg + a c6 du = c6 − tan−1 2 2 (u + a) + b 2 b Eg c5 (3.50) . (3.51) From Eq. (3.44), it is seen that t Eg f (u)du = A1 [ln(t + E)] − A1 ln(Eg + E) +A2 [ln |t − E|] − A2 ln |Eg − E| 88 3.4. Derivation c3 c3 ln[(t − a)2 + b2 ] − ln[(Eg − a)2 + b2 ] 2 2 π Eg − a +c4 − tan−1 2 b c5 c5 + ln[(t + a)2 + b2 ] − ln[(Eg + a)2 + b2 ] 2 2 π Eg + a +c6 − tan−1 . 2 b + (3.52) Letting φ1 = ln(t + E)] − ln(Eg + E), φ2 = ln |t − E| − ln |Eg − E|, φ3 = ln[(t − a)2 + b2 ] − ln[(Eg − a)2 + b2 ], φ4 = π − tan−1 2 Eg − a b , φ5 = ln[(t + a)2 + b2 ] − ln[(Eg + a)2 + b2 ], and φ6 = π − tan−1 2 Eg + a b , Eq. (3.52) leads to ∞ f (u)du = A1 φ1 + A2 φ2 + Eg c3 c5 φ3 + φ5 + c4 φ4 + c6 φ6 . (3.53) 2 2 89 3.4. Derivation 3.4.8 Summation by parts For simplicity, the summation, depicted in Eq. (3.53), is performed in parts such that term 1 = A1 φ1 + A2 φ2 , term 2 = c3 c5 φ3 + φ5 , 2 2 term 3 = c4 φ4 + c6 φ6 . From Eq. (3.53), it is thus seen that ∞ f (u)du = A1 φ1 + A2 φ2 + Eg term 1 c3 c5 φ3 + φ5 + c4 φ4 + c6 φ6 . 2 2 term 2 term 3 Calculation of term 1 term 1 = A1 φ1 + A2 φ2 , −(E + Eg )2 (E − Eg )2 φ + φ2 , 1 2E[(E 2 − Eo2 )2 + C 2 E 2 ] 2E[(E 2 − Eo2 )2 + C 2 E 2 ] E 2 + Eg2 2EEg = [φ2 − φ1 ] − [φ2 + φ1 ], 2 2 2 2 2 2 2E[(E − Eo ) + C E ] 2E[(E − Eo2 )2 + C 2 E 2 ] E 2 + Eg2 Eg = [φ2 − φ1 ] − [φ2 + φ1 ], 2 2 2 2 2 2 2 2E[(E − Eo ) + C E ] [(E − Eo )2 + C 2 E 2 ] E 2 + Eg2 Eg = [φ2 − φ1 ] − 4 [φ2 + φ1 ], 4 2Eζ ζ = 90 3.4. Derivation where φ2 − φ1 = ln |t − E| − ln |Eg − E| − ln(t + E)] + ln(Eg + E), |t − E| (Eg + E) + ln , t+E |Eg − E| (Eg + E) = ln , [as t → ∞] |Eg − E| = ln φ2 + φ1 = ln |t − E| − ln |Eg − E| + ln(t + E)] − ln(Eg + E), = ln[(t + E)|t − E|] − ln[(Eg + E)|Eg − E|]. 91 3.4. Derivation Calculation of term 2: term 2 = = c5 c3 φ3 + φ5 , 2 2 8ab(a2 + b2 )[(a2 1 × − b2 − E 2 )2 + 4a2 b2 ] 2b(a − Eg )[a(a2 − b2 − E 2 ) − 2ab2 ]φ3 − ((a − E )2 − b2 )[b(a2 − b2 − E 2 ) + 2a2 b]φ 3 g + 2 2 2 2 −2b(a + Eg )[a(a − b − E ) − 2ab ]φ5 + ((a + Eg )2 − b2 )[b(a2 − b2 − E 2 ) + 2a2 b]φ5 = 8ab(a2 + b2 )[(a2 , 1 × − b2 − E 2 )2 + 4a2 b2 ] 2ba[a(a2 − b2 − E 2 ) − 2ab2 ][φ3 − φ5 ] − (a2 + E 2 − b2 )[b(a2 − b2 − E 2 ) + 2a2 b][φ − φ ] 3 5 g + −2bEg [a(a2 − b2 − E 2 ) − 2ab2 ][φ3 + φ5 ] + 2aEg [b(a2 − b2 − E 2 ) + 2a2 b][φ3 + φ5 ] , 92 3.4. Derivation = 8a(Eo2 − C2 4 + C2 4 )[(Eo2 − C2 4 1 − C2 4 − E 2 )2 + 4(Eo2 − 2a2 (a2 − 3b2 − E 2 )[φ3 − φ5 ] − (a2 − b2 + Eg2 )(3a2 − b2 − E 2 )[φ3 − φ5 ] + [−2abE (a2 − b2 − E 2 ) + 2abE 2b2 ][φ + φ ] g g 3 5 + [2abEg (a2 − b2 − E 2 ) + 2abEg 2a2 ][φ3 + φ5 ] = 8a Eo2 [(E 2 1 × − Eo2 )2 + C 2 E 2 ] C2 [2(Eo2 − )(Eo2 − C 2 − E 2 )][φ3 − φ5 ] 4 − C2 C2 )4] 4 × , C2 + Eg2 )(3Eo2 − C 2 − E 2 )][φ3 − φ5 ] − 2 Eg + [φ3 + φ5 ], 2 2[(E − Eo2 )2 + C 2 E 2 ] −Eo4 + Eg2 C 2 + Eg2 E 2 − Eo2 E 2 − 3Eo2 Eg2 [φ3 − φ5 ] = 8a Eo2 [(E 2 − Eo2 )2 + C 2 E 2 ] Eg [φ3 + φ5 ], + 2 2[(E − Eo2 )2 + C 2 E 2 ] (Eg2 − Eo2 )E 2 + Eg2 C 2 − Eo2 (Eo2 + 3Eg2 ) = [φ3 − φ5 ] 8a Eo2 [(E 2 − Eo2 )2 + C 2 E 2 ] Eg + [φ3 + φ5 ], 2 2[(E − Eo2 )2 + C 2 E 2 ] aln Eg = [φ3 − φ5 ] + [φ3 + φ5 ], 2 4 4χ Eo ζ 2 ζ4 [(Eo2 93 3.4. Derivation where φ3 − φ5 = ln[(t − a)2 + b2 ] − ln[(Eg − a)2 + b2 ] − ln[(t + a)2 + b2 ] + ln[(Eg + a)2 + b2 ], (t − a)2 + b2 [(Eg + a)2 + b2 ] = ln , + ln (t + a)2 + b2 [(Eg − a)2 + b2 ] [Eg2 + a2 + b2 + 2aEg ] = ln 2 , [as t → ∞] [Eg + a2 + b2 − 2aEg ] = ln [Eg2 + Eo2 + 2aEg ] , [Eg2 + Eo2 − 2aEg ] [Eg2 + Eo2 + χEg ] = ln 2 , [recall that, χ = [Eg + Eo2 − χEg ] 4Eo2 − C 2 = 2a] φ3 + φ5 = ln[(t − a)2 + b2 ] − ln[(Eg − a)2 + b2 ] + ln[(t + a)2 + b2 ] − ln[(Eg + a)2 + b2 ], = ln[[(t − a)2 + b2 ][(t + a)2 + b2 ]] − ln[[(Eg − a)2 + b2 ][(Eg + a)2 + b2 ]], = ln[[(t − a)2 + b2 ][(t + a)2 + b2 ]] − ln[(Eg2 + Eo2 + χEg )(Eg2 + Eo2 − χEg )], = ln[[(t − a)2 + b2 ][(t + a)2 + b2 ]] − ln[(Eo2 − Eg2 )2 + Eg2 C 2 )]. 94 3.4. Derivation 3.4.9 Calculation of term 3: term 3 = c4 φ4 + c6 φ6 , = 4ab(a2 + b2 )[(a2 1 × − b2 − E 2 )2 + 4a2 b2 ] 2 2 2 2 2b(a − Eg )[b(a − b − E ) + 2a b]φ4 + ((a − E )2 − b2 )[a(a2 − b2 − E 2 ) − 2ab2 ]φ g 4 + 2b(a + Eg )[b(a2 − b2 − E 2 ) + 2a2 b]φ6 + ((a + Eg )2 − b2 )[a(a2 − b2 − E 2 ) − 2ab2 ]φ6 , 2b2 a(3a2 − b2 − E 2 ) + a(a2 + Eg2 − b2 )(a2 − 3b2 − E 2 ) = [φ4 + φ6 ] 4ab(a2 + b2 )[(a2 − b2 − E 2 )2 + 4a2 b2 ] 2b2 Eg (3a2 − b2 − E 2 ) + 2a2 Eg (a2 − 3b2 − E 2 ) [φ6 − φ4 ], + 4ab(a2 + b2 )[(a2 − b2 − E 2 )2 + 4a2 b2 ] 2b2 (3a2 − b2 − E 2 ) + (a2 + Eg2 − b2 )(a2 − 3b2 − E 2 ) [φ4 + φ6 ] = 2C Eo2 [(E 2 − Eo2 )2 + C 2 E 2 ] 2Eg [b2 (3a2 − b2 − E 2 ) + a2 (a2 − 3b2 − E 2 ) + [φ6 − φ4 ], χ C Eo2 [(E 2 − Eo2 )2 + C 2 E 2 ] = C2 (3Eo2 2 2 − C 2 − E 2 ) + (Eo2 − C2 + Eg2 )(Eo2 − C 2 − E 2 ) [φ4 + φ6 ] 2C Eo2 [(E 2 − Eo2 )2 + C 2 E 2 ] 2 2 2Eg [ C4 (3Eo2 − C 2 − E 2 ) + (Eo2 − C4 )(Eo2 − C 2 − E 2 )] + [φ6 − φ4 ], χ C Eo2 [(E 2 − Eo2 )2 + C 2 E 2 ] 95 3.4. Derivation = Eo4 − E 2 Eo2 + Eo2 Eg2 − Eg2 E 2 − Eg2 C 2 [φ4 + φ6 ] 2C Eo2 [(E 2 − Eo2 )2 + C 2 E 2 ] 2 2 2Eg [− Eo2C + Eo4 − Eo2 E 2 ] + [φ6 − φ4 ], χ C Eo2 [(E 2 − Eo2 )2 + C 2 E 2 ] −[(E 2 − Eo2 )(Eo2 + Eg2 ) + Eg2 C 2 ] = [φ4 + φ6 ] 2C Eo2 [(E 2 − Eo2 )2 + C 2 E 2 ] 2 2Eg [E 2 − Eo2 + C2 ] [φ4 − φ6 ], + χ C [(E 2 − Eo2 )2 + C 2 E 2 ] −aa tan 2Eg [E 2 − γ 2 ] = [φ + φ ] + [φ4 − φ6 ]. 4 6 2C Eo2 ζ 4 χCζ 4 where φ4 + φ6 = = = = φ4 − φ6 , = = = = Eg − a Eg + a π + − tan−1 , b 2 b −Eg + a Eg + a + tan−1 , π − tan−1 b b −2Eg + 2a 2Eg + 2a π − tan−1 + tan−1 , C C 2Eg + χ −2Eg + χ π − tan−1 + tan−1 , C C π Eg − a π Eg + a − tan−1 − + tan−1 , 2 b 2 b 2ab tan−1 2 , b − a2 + Eg2 χC tan−1 , 2 −2(a − b2 − Eg2 ) π − tan−1 2 π + tan−1 2 recall that, tan−1 2(a2 − b2 − Eg2 ) χC y π = + tan−1 x 2 , −x y 96 3.4. Derivation π + tan−1 = 2 = 3.4.10 π + tan−1 2 2 2(Eo2 − C2 − Eg2 ) χC 2(γ 2 − Eg2 ) χC , . Summation of term 1, term 2 and term 3 term 1+term 2+term 3 E 2 + Eg2 Eg = [φ2 − φ1 ] − 4 [φ2 + φ1 ] 4 2Eζ ζ aln Eg + [φ3 − φ5 ] + [φ3 + φ5 ] 2 4 4χEo ζ 2.ζ 4 −aa tan 2Eg [E 2 − γ 2 ] + [φ + φ ] + [φ4 − φ6 ], 4 6 2CEo2 ζ 4 χCζ 4 E 2 + Eg2 (Eg + E) Eg = ln − 4 ln [(t + E)|t − E|] 4 2Eζ |Eg − E| ζ Eg + 4 ln[(Eg + E)|Eg − E|] ζ [Eg2 + Eo2 + χEg ] aln + ln 8aEo2 ζ 4 [Eg2 + Eo2 − χEg ] Eg + 4 ln[[(t − a)2 + b2 ][(t + a)2 + b2 ]] 2ζ Eg − 4 ln[(Eo2 − Eg2 )2 + Eg2 C 2 )] 2ζ 2Eg + χ aa tan −2Eg + χ −1 −1 π − tan − + tan 2CEo2 ζ 4 C C 2 2 2(γ − Eg ) 2Eg [E 2 − γ 2 ] π + + tan−1 , [here t → ∞] 4 χCζ 2 χC 97 3.4. Derivation = Eg2 + Eo2 + χEg aln ln 4χEo2 ζ 4 Eg2 + Eo2 − χEg aa tan −2Eg + χ 2Eg + χ − + tan−1 π − tan−1 2 4 2CEo ζ C C 2 2 2 2 2(γ − Eg ) Eg [E − γ ] + π + 2 tan−1 4 χCζ χC 2 2 E + Eg Eg |Eg − E| [(Eg + E)|Eg − E|] + 4 ln . − ln 4 2Eζ (Eg + E) ζ (E 2 − E 2 )2 + E 2 C 2 ) o 3.4.11 g g Final Tauc-Lorentz expression for the model of Jellison and Modine Now, from Eq. (3.43), ∞ 2AEo C P f (u)du, 1 (E) = 1∞ + π Eg 2AEo C = 1∞ + [term 1 + term 2 + term 3], π Eg2 + Eo2 + χEg 2AEo C aln = 1∞ + ln π 4χEo2 ζ 4 Eg2 + Eo2 − χEg 2AEo C aa tan 2Eg + χ − π − tan−1 + tan−1 2 4 π 2CEo ζ C 2(γ 2 − Eg2 ) 2AEo C Eg [E 2 − γ 2 ] −1 π + 2 tan + π χCζ 4 χC 2 2 2AEo C E + Eg |Eg − E| − ln 4 π 2Eζ (Eg + E) 2AEo C Eg [(Eg + E)|Eg − E|] . + ln π ζ4 (E 2 − E 2 )2 + E 2 C 2 ) o g −2Eg + χ C g 98 3.5. Comparing the results of the Jellison and Modine model with experimental data =⇒ 1 (E) = Eg2 + Eo2 + χEg 1 A C aln ln 2 π ζ 4 χEo Eg2 + Eo2 − χEg A aa tan 2Eg + χ −2Eg + χ − 4 π − tan−1 + tan−1 πζ Eo C C 2 2 2(γ − Eg ) AEo +2 Eg [E 2 − γ 2 ] π + 2 tan−1 4 πχζ χC 2 2 |Eg − E| AEo C E + Eg ln − 4 πζ E (Eg + E) 2AEo C [(Eg + E)|Eg − E|] + Eg ln , 4 πζ (E 2 − E 2 )2 + E 2 C 2 ) 1∞ + o g g this being the same result as that in Eq. (3.12), all terms being defined in the same manner. 3.5 Comparing the results of the Jellison and Modine model with experimental data The parameterized model for the optical functions, 1 (E) and 2 (E), of Jellison and Modine [15] is now fit with the results of experiment. The experimental results that are considered include that of Piller [30], Synowicki [31], and Ferluato et al. [32], these experimental results corresponding to the EGE a-Si, PECVD a-Si I, and PECVD a-Si II experimental data sets mentioned 99 3.5. Comparing the results of the Jellison and Modine model with experimental data in Section 2.10, respectively; these experimental data sets are further described in Section 2.10. The approach that is adopted for each experimental fit involves the selection of the modeling parameters that “best-fits” the corresponding experimental data set. This is achieved through a selection of modeling parameters, a determination of the corresponding spectral dependencies of the optical functions, 1 (E) and 2 (E), and a comparison with the results of experiments. The “best-fit” is obtained through visual inspection. Through a sweep over a range of parameter values, the parameter selections that “best-fit” the corresponding experimental optical functions is determined. It is noted that, in all cases, the results of the model of Jellison and Modine [15] agree with that of experiment over most of the spectrum considered. The “best-fit” parameter values, corresponding to the different experimental fits, are tabulated in Table 3.1. The resultant fits, and the corresponding experimental data sets, are depicted in Figures 3.1, 3.2, and 3.3, these corresponding to the experimental results of EGE a-Si, PECVD a-Si I, and PECVD a-Si II, respectively. 100 3.5. Comparing the results of the Jellison and Modine model with experimental data Table 3.1: The “best-fit” a-Si modeling parameter selections employed for the Tauc-Lorentz model of Jellison and Modine [15] for fits to experimental data corresponding to EGE a-Si, PECVD a-Si I, and PECVD a-Si II. parameter (units) EGE a-Si [30] PECVD a-Si I [31] PECVD a-Si II [32] A (eV) 125 196 203 Eo (eV) 3.43 3.64 3.72 C (eV) 2.54 2.3 2.1 Eg (eV) 1.21 1.51 1.68 1∞ 1.0 1.0 0.65 101 3.5. Comparing the results of the Jellison and Modine model with experimental data 30 ε1 EGE a−Si [30] ε2 ε1, ε2 20 10 0 −10 1 2 3 4 Photon Energy (eV) 5 Figure 3.1: A “best-fit” of the Jellison and Modine model with the results of experiment corresponding to the EGE a-Si expermental data set [30]. These experimental results are from Piller [30]. The fits to experiment are shown with the solid lines. The experimental data, EGE a-Si [30], is depicted with the solid points. 102 3.5. Comparing the results of the Jellison and Modine model with experimental data 40 ε1, ε2 30 PECVD a−Si I [31] ε2 ε1 20 10 0 −10 1 2 3 4 Photon Energy (eV) 5 Figure 3.2: A “best-fit” of the Jellison and Modine model with the results of experiment corresponding to the PECVD a-Si I experimental data set [31]. These experimental results are from Synowicki [31]. The fits to experiment are shown with the solid lines. The experimental data, PECVD a-Si I [31], is depicted with the solid points. 103 3.5. Comparing the results of the Jellison and Modine model with experimental data 40 ε1, ε2 30 PECVD a−Si II [32] ε2 ε1 20 10 0 −10 1 2 3 4 Photon Energy (eV) 5 Figure 3.3: A “best-fit” of the Jellison and Modine model with the results of experiment corresponding to the PECVD a-Si II experimental data set [32]. These experimental results are from Ferluato et al. [32]. The fits to experiment are shown with the solid lines. The experimental data, PECVD a-Si II [32], is depicted with the solid points. 104 Chapter 4 A physically based Kramers-Kronig consistent model for the optical functions associated with amorphous semiconductors 4.1 Evaluating the JDOS function from the distributions of electronic states One weakness associated with the models of Forouhi and Bloomer [20], McGahan et al. [21], and Jellison and Modine [15], is that they do not pro105 4.1. Evaluating the JDOS function from the distributions of electronic states vide for a direct connection with the underlying distributions of electronic states. Even the model of Jellison and Modine [15], while being considerably improved over that of Forouhi and Bloomer [20] and McGahan et al. [21], fails to provide for a direct relationship between the underlying valence band and conduction band DOS functions, Nv (E) and Nc (E), respectively, with the spectral dependence of the optical functions, 1 (E) and 2 (E); Jellison and Modine [15] take a physically based model for the spectral dependence of the imaginary part of the dielectric function, 2 (E), and employ a Kramers- Kronig transformation in order to obtain the spectral dependence of the real part of the dielectric function, 1 (E). A direct relationship between the func- tional dependencies of these optical functions, 1 (E) and 2 (E), with the form of the distributions of electronic states, would allow for the gleaning of insights into the underlying electronic properties from experimental measurements of the spectral dependence of these optical functions. In this chapter, an empirical model for the valence band and conduction band DOS functions, Nv (E) and Nc (E), is proposed, this model capturing the fundamental essence of an amorphous semiconductor, i.e., exponential distributions of tail states and square-root distributions of band states. Then, the corresponding joint density of states (JDOS) function is evaluated through 106 4.2. Empirical model for the DOS functions the use of Eq. (2.10). Through the application of a model for the aggregate optical transition matrix element, the corresponding spectral dependence of the imaginary part of the dielectric function, 2 (E), is determined. A Kramers-Kronig transformation is then employed in order to determine the corresponding spectral dependence of the real part of the dielectric function, 1 (E). A comparison with the results of experiment is then shown. This chapter is organized in the following manner. In Section 4.2, an empirical model for the DOS functions is presented. Then, in Section 4.3, the corresponding JDOS function, J(E), is determined. In Section 4.4, a model for the spectral dependence of the dipole optical matrix element is presented. Finally, in Section 4.5, a fit with the results of experiment, for the spectral dependence of real and imaginary parts of the dielectric function, 1 (E) unity 4.2 and 1∞ 2 (E), respectively, is shown. Finally, the implications of a non- are described in Section 4.6. Empirical model for the DOS functions There is general consensus that the DOS functions associated with an amorphous semiconductor, i.e., Nv (E) and Nc (E), exhibit square-root func- 107 4.2. Empirical model for the DOS functions tional dependencies in the band regions and exponential functional dependencies in the tail regions. Following O’Leary [29], one can thus set EvT − Ev Ev − E Ev − EvT exp exp , E > EvT γ γ v v , Nv (E) = Nvo Ev − E, E ≤ EvT (4.1) and Nc (E) = Nco E − Ec , E ≥ EcT , EcT − Ec exp Ec − EcT γc exp E − Ec γc , E < EcT (4.2) where Nvo and Nco denote the valence band and conduction band DOS prefactors, respectively, Ev and Ec represent the valence band and conduction band band edges, γv and γc are the breadths of the valence band and conduction band tails, EvT and EcT being the critical energies at which the exponential and square-root distributions interface. It is clear, from Eqs. (4.1) and (4.2), that both Nv (E) and Nc (E) are continuous functions of energy. In order to further simplify these DOS functions, it will be further assumed that the derivatives of these DOS functions, Nv (E) and Nc (E), are also continuous functions of energy, i.e., these functions are smooth func108 4.2. Empirical model for the DOS functions tions. For this to be the case, it may be shown that EvT = Ev − γv /2 and that EcT = Ec + γc /2 [28]. With these simplifications, Eq. (4.1) reduces to γv 1 Ev − E exp − exp , E > EvT 2 2 γ v , (4.3) Nv (E) = Nvo Ev − E, E ≤ EvT and Eq. (4.2) reduces to Nc (E) = Nco γc 1 exp − 2 2 E − Ec , E ≥ EcT , exp E − Ec γc (4.4) , E < EcT It is noted that in the disorderless limit, i.e., when γv → 0 and γc → 0, these DOS functions reduce to Nv (E) → Nvo and Nc (E) → Nco 0, E > Ev , (4.5) , (4.6) Ev − E, E ≤ Ev E − Ec , E ≥ Ec 0, E < Ec respectively, these being the Tauc DOS functions, i.e., Eqs. (2.23) and (2.24), respectively. 109 4.2. Empirical model for the DOS functions The valence band DOS function, Nv (E), is depicted on a linear scale in Figure 4.1, Nvo and Ev being nominally set to 2 ×1022 cm−3 eV−3/2 and 0.0 eV, respectively, for the purposes of this analysis, i.e., the valence band band edge forms the reference energy for this plot. A variety of selections of γv are made, these spanning over the range of values representative of amorphous semiconductors. It is clearly seen that for energies below the valence band band edge, Ev , with the tail breadth set to zero, the valence DOS function terminates abruptly, i.e., there are no electronic states within the gap. For finite γv , however, a distribution of tail states encroaches into the gap region, the amount of encroachment increasing with γv . In Figure 4.2, this plot is depicted on a logarithmic scale, the same parameter selections as that employed in Figure 4.1 also being used. The conduction band DOS function, Nc (E), is depicted on linear and logarithmic scales in Figures 4.3 and 4.4, respectively, Nco and Ec being nominally set to 2 ×1022 cm−3 eV−3/2 and 0.0 eV, respectively, for the purposes of this analysis. As with the valence band case, greater γc corresponds to greater encroachment into the energy gap. 110 4.2. Empirical model for the DOS functions Nvo = 2×1022 cm−3eV−3/2 γv = 0 meV 6 γ = 25 meV v DOS (× 10 21 −3 −1 cm eV ) 8 γ = 50 meV 4 v γv = 100 meV 2 0 −0.2 −0.1 0 Energy (eV) 0.1 0.2 Figure 4.1: A linear plot of the valence band DOS function as a function of energy, E, for various selections of γv . The DOS modeling parameters, Nvo and Ev , are nominally set to 2 ×1022 cm−3 eV−3/2 and 0.0 eV, respectively, for the purposes of this plot. This figure is after O’Leary [17]. 111 4.2. Empirical model for the DOS functions 22 Nvo = 2×1022 cm−3eV−3/2 10 DOS (cm−3eV−1) 21 10 20 10 19 10 γv= 100 meV γv= 50 meV γv= 25 meV γv= 0 meV 18 10 17 10 −0.2 −0.1 0 Energy (eV) 0.1 0.2 Figure 4.2: A logarithmic plot of the valence band DOS function as a function of energy, E, for various selections of γv . The DOS modeling parameters, Nvo and Ev , are nominally set to 2 ×1022 cm−3 eV−3/2 and 0.0 eV, respectively, for the purposes of this plot. This figure is after O’Leary [17]. 112 4.2. Empirical model for the DOS functions 10 DOS (× 1021 cm−3eV−1) 8 Nco = 2×1022 cm−3eV−3/2 6 γc = 25 meV 4 2 γc = 50 meV γc = 100 meV γc = 0 meV 0 −0.5 −0.3 −0.1 0 0.1 Energy (eV) 0.3 0.5 Figure 4.3: A linear plot of the conduction band DOS function as a function of energy, E, for various selections of γc . The DOS modeling parameters, Nco and Ec , are nominally set to 2 ×1022 cm−3 eV−3/2 and 0.0 eV, respectively, for the purposes of this plot. This figure is after O’Leary [17]. 113 4.2. Empirical model for the DOS functions Nco = 2×1022 cm−3eV−3/2 20 DOS (cm−3eV−1) 10 15 10 10 10 γc= 25 meV γc= 50 meV γc= 0 meV γc= 100 meV 5 10 −0.5 −0.3 −0.1 0 0.1 Energy (eV) 0.3 0.5 Figure 4.4: A logarithmic plot of the conduction band DOS function as a function of energy, E, for various selections of γc . The DOS modeling parameters, Nco and Ec , are nominally set to 2 ×1022 cm−3 eV−3/2 and 0.0 eV, respectively, for the purposes of this plot. This figure is after O’Leary [17]. 114 4.3. JDOS evaluation and analysis 4.3 JDOS evaluation and analysis The JDOS function, J(E), corresponding to the empirical model for the DOS functions, i.e., with Nv (E) and Nc (E) as set in Eqs. (4.3) and (4.4), respectively, is now evaluated using Eq. (2.10). This JDOS function may be used in order to determine the number of possible optical transitions between the valence band and the conduction band electronic states. In Figure 4.5, the JDOS function is depicted on a linear scale for a variety of selections of γv and γc , the DOS modeling parameters being nominally set to Nvo = Nco = 2 ×1022 cm−3 eV−3/2 , Ev = 0.0 eV, and Ec = 1.7 eV; for all cases, γv is set equal to γc . This plot is shown on a logarithmic scale in Figure 4.6. Eq. (2.25) suggests that the functional dependence of J(E) on the photon energy, E, for energies well above the band gap, allows one to determine the corresponding energy gap, Eg . In Figure 4.7, a linear plot of the functional dependence of J(E) on the photon energy, E, is depicted. It is clearly seen that for γv = γc = 0, that the J(E) function termi- nates abruptly at the band gap edge. For finite γv and γc , however, a tail in the JDOS function encroaches into the gap region, this tail representing the amount of disorder that is present. 115 4.3. JDOS evaluation and analysis 1.2 −6 44 JDOS × 10 Nvo = Nco = 2 × 1022 cm−3eV−3/2 Eg = 1.7 eV −1 (cm eV ) 1 0.8 0.6 γv = γc = 100 meV 0.4 γv = γc = 50 meV γv = γc = 0 meV 0.2 0 1.8 2 2.2 Photon Energy (eV) 2.4 Figure 4.5: A linear plot of the JDOS function as a function of the photon energy, E, for various selections of γv and γc . The DOS modeling parameters are nominally set to Nvo = Nco = 2 ×1022 cm−3 eV−3/2 , Ev = 0.0 eV, and Ec = 1.7 eV for the purposes of this plot. This figure is after O’Leary [17]. 116 4.3. JDOS evaluation and analysis 44 10 Nvo = Nco = 2 × 1022 cm−3eV−3/2 E = 1.7 eV 42 g JDOS (cm−6eV−1) 10 γv = γc = 100 meV γv = γc = 50 meV 40 10 γv = γc = 0 meV 38 10 36 10 34 10 0.5 1 1.5 2 Photon Energy (eV) 2.5 Figure 4.6: A logarithmic plot of the JDOS function as a function of the photon energy, E, for various selections of γv and γc . The DOS modeling parameters are nominally set to Nvo = Nco = 2 ×1022 cm−3 eV−3/2 , Ev = 0.0 eV, and Ec = 1.7 eV for the purposes of this plot. This figure is after O’Leary [17]. 117 4.3. JDOS evaluation and analysis 1.8 J DOS × 1022(cm−3 eV −1/2) 1.6 1.4 1.2 1 0.8 Nvo = Nco = 2 × 1022 cm−3eV−3/2 Eg = 1.7 eV γv = γc = 100 meV γv = γc = 50 meV √ 0.6 0.4 γv = γc = 0 meV 0.2 0 0.5 1 1.5 2 Photon Energy (eV) 2.5 3 Figure 4.7: A linear plot of the functional dependence of the squareroot of the JDOS function on the photon energy, E, for various selections of γv and γc . The DOS modeling parameters are nominally set to Nvo = Nco = 2 ×1022 cm−3 eV−3/2 , Ev = 0.0 eV, and Ec = 1.7 eV for this plot. This figure is after O’Leary [17]. 118 4.4. Modeling the optical functions associated with a-Si 4.4 Modeling the optical functions associated with a-Si In order to determine the spectral dependence of the a-Si optical functions, 1 (E) and 2 (E), the spectral dependence of the JDOS function as- sociated with this material must be known. With a model for the JDOS function, in conjunction with a model for the spectral dependence of the aggregate dipole matrix element, R2 (E), the spectral dependence of the imaginary part of the dielectric function, 2 (E), may be determined through the use of Eq. (2.17). Then, through the use of a Kramers-Kronig transformation, i.e., Eq. (2.45), the spectral dependence of the real part of the dielectric function, 1 (E), may be determined. For the purposes of this analysis, the spectral dependence of J(E) and R2 (E) associated with a-Si will be independently modeled. Each of these models will then be contrasted with the corresponding results of experiment, these a-Si experimental results being from Jackson et al. [16]. From the spectral dependence of J(E) and R2 (E), the spectral dependence of 2 (E) will be determined, and compared with that of experiment. Following this analysis, fits with the results of experiment corresponding to other samples of a-Si will 119 4.4. Modeling the optical functions associated with a-Si be performed, it being expected that some variations in the DOS modeling parameters will occur from sample-to-sample, these variations representing the changes in the physics of each sample. Thevaril and O’Leary [33] plotted the experimentally determined a-Si valence band tail breadth parameter, γv , as a function of the conduction band tail breadth parameter, γc , for the experimental data of Sherman et al. [23], Rerbal et al. [34], Teidje et al. [35], and Winer and Ley [36]. Fit values, obtained by O’Leary [17], were also considered. This experimental data, and the fit data, is depicted in Figure 4.8. It is noted that, in all cases, the valence band tail breadth parameter, γv , exceeds that associated with the conduction band, γc . Thevaril and O’Leary [28] demonstrated that when γv exceeds γc , the optical absorption coefficient is primarily determined by γv . Thus, many of the optical properties may be determined neglecting the conduction band tail states, i.e., assuming that γc is equal to zero. For the purposes of this analysis, γc will be set to zero for all cases. In Figure 4.9, a plot of the experimental a-Si JDOS results of Jackson et al. [16] is depicted. This result is directly contrasted with that determined using the numerical evaluation procedure for the JDOS function for the DOS modeling parameters set to Nvo = Nco = 2.48 ×1022 cm−3 eV−3/2 , 120 4.4. Modeling the optical functions associated with a-Si 70 60 γv (meV) 50 γv=γc 40 O’Leary [17] 30 Sherman et al. [23] Rerbal et al. [34] 20 Teidje et al. [35] Winer and Ley [36] 10 0 0 10 20 30 40 γc (meV) 50 60 70 Figure 4.8: A plot of the dependence of γv on γc . Experimental results from Sherman et al. [23], Rerbal et al. [34], Teidje et al. [35], and Winer and Ley [36] are depicted. A modeling result, obtained through a fit with experimental data, from O’Leary [17], is also depicted on this plot. This figure is after Thevaril and O’Leary [33]. 121 4.4. Modeling the optical functions associated with a-Si 44 JDOS (cm−6eV−1) 10 Jackson et al. [16] 42 10 40 10 38 10 0 1 2 3 Photon Energy (eV) 4 Figure 4.9: The spectral dependence of the JDOS function. 5 The a-Si JDOS results of Jackson et al. [16] are depicted with the solid points. The JDOS fit result, obtained by setting Nvo = Nco = 2.48 ×1022 cm−3 eV−3/2 , Ev = 0.0 eV, Ec = 1.68 eV, and γv = 40 meV, is depicted with the solid line. 122 4.4. Modeling the optical functions associated with a-Si Ev = 0.0 eV, Ec = 1.68 eV, and γv = 40 meV. It is noted that the fit is reasonably satisfactory, except for low JDOS values, where defect absorption influences the results and high photon energies, where the non-parabolicity of the bands plays a role influencing the form of the JDOS function. This a-Si JDOS function will be used in the subsequent analysis. For the case of a-Si, Jackson et al. [16] experimentally determined the spectral dependence of the aggregate dipole matrix element squared average, R2 (E). They found that for photon energies below 3.4 eV, this matrix element remains essentially constant, while for photon energies in excess of 3.4 eV, it falls off with the algebraic dependence E −5 . Based on these observations, Thevaril and O’Leary [28] modeled R2 (E) as 5 Ed , E ≥ Ed E R2 (E) = R2o , (4.7) 1, E < Ed where Ed is the characteristic energy to which R2 (E) is constant, i.e., R2 (E) = R2o , where, for the case of the a-Si experimental data of Jackson et al. [16], R2o = 10˚ A2 . This spectral dependence for R2 (E), proposed by Thevaril and O’Leary [28], is contrasted with the experimental R2 (E) results of Jackson et al. [16] in Figure 4.10. The corresponding spectral dependence of 2 (E), 123 4.4. Modeling the optical functions associated with a-Si 3 10 Jackson et al. [16] 2 R2 (˚ A2 ) 10 Ed 1 10 10 0 10 −1 10 1 2 3 4 5 Photon Energy (eV) 6 7 Figure 4.10: R2 (E) as a function of the photon energy, E, determined using Eq. (4.7), depicted with the solid lines. For the purposes of this fit, R2o is ˚2 and Ed is set to 3.4 eV. The experimental results of Jackson et set to 10 A al. [16] are depicted with the solid points. 124 4.4. Modeling the optical functions associated with a-Si determined using Eq. (2.17), is depicted on a logarithmic scale in Figure 4.11 and on a linear scale in Figure 4.12. The DOS modeling parameters are set to Nvo = Nco = 2.48 ×1022 cm−3 eV−3/2 , Ev = 0.0 eV, Ec = 1.68 eV, and γv = 40 meV for the purposes of this plot. It is also compared with the experimental 2 (E) results of Jackson et al. [16]. While the comparison with the results of experiment is reasonably satisfactory for the case of the logarithmic plot, i.e., Figures 4.11, for the linear plot, i.e., Figure 4.12, it is seen that the spectral dependence of 2 (E) does not follow the a-Si experimental results for photon energies in excess of 3.4 eV. From Figure 4.12, it is also seen that the spectral dependence of 2 (E) exhibits a sharp peak at Ed , which does not accord with experimental observation. The formalism for the dipole matrix element squared average, R2 (E), expressed in Eq. (4.7), is responsible for this unphysical sharp peak in the spectral dependence of the imaginary part of the dielectric function, 2 (E). In order to avoid this problem created with an abrupt transition model for the spectral dependence of R2 (E), R2 (E) is instead remodeled as R2o R2 (E) = 1+ E Ed 2x 1 2 , (4.8) 125 4.4. Modeling the optical functions associated with a-Si 2 10 Ed 0 ε2 10 Jackson et al. [16] −2 10 −4 10 1 2 3 4 5 Photon Energy (eV) 6 Figure 4.11: The spectral dependence for the imaginary part of the dielectric function, 2 (E), determined using Eq. (2.17), where R2 (E) is determined using Eq. (4.7), on a logarithmic scale. The DOS modeling parameters are set to Nvo = Nco = 2.48 ×1022 cm−3 eV−3/2 , Ev = 0.0 eV, Ec = 1.68 eV, and γv = 40 meV for the purposes of this plot. The experimental results of Jackson et al. [16] are depicted with the solid points 126 4.4. Modeling the optical functions associated with a-Si 35 Ed 30 Jackson et al. [16] 25 ε2 20 15 10 5 0 1 2 3 4 5 Photon Energy (eV) 6 Figure 4.12: The spectral dependence of the imaginary part of the dielectric function, 2 (E), determined using Eq. (2.17), where R2 (E) is determined using Eq. (4.7), on a linear scale. The DOS modeling parameters are set to Nvo = Nco = 2.48 ×1022 cm−3 eV−3/2 , Ev = 0.0 eV, Ec = 1.68 eV, and γv = 40 meV for the purposes of this plot. The experimental results of Jackson et al. [16] are depicted with the solid points 127 4.4. Modeling the optical functions associated with a-Si where x represents the asymptotic polynomial rate of attenuation, i.e., for the case of the experimental observations of Jackson et al. [16], it is found that x = 5. It is noted that for photon energies E << Ed , R2 (E) = R2o . For photon energies greater than Ed , however, R2 (E) asymptotically approaches R2o ( EEd )x . The spectral dependence of this remodeled R2 (E) is contrasted with the experimental R2 (E) results of Jackson et al. [16] and the spectral dependence of R2 (E), proposed by Thevaril and O’Leary [28], i.e., Eq. (4.7), in Figure 4.13. The spectral dependence of 2 associated with the a-Si, evaluated using Eq. (2.17), where R2 (E) is determined using Eq. (4.8), is depicted on a logarithmic scale in Figure 4.14 and on a linear scale in Figure 4.15. The DOS modeling parameters are set to Nvo = Nco = 2.48 ×1022 cm−3 eV−3/2 , Ev = 0.0 eV, Ec = 1.68 eV, and γv = 40 meV for the purposes of this analysis. The asymptotic polynomial rate of attenuation, x, is set to 6.6 for the purposes of this plot. The experimental 2 (E) results of Jackson et al. [16] are depicted with the solid points. It is noted that the fit is reasonably satisfactory for both cases. The evaluation of the real part of dielectric function, 1 (E), may be directly determined from a Kramers-Kronig transformation, i.e., Eq. (2.45). Numerical integration is employed for the 128 4.4. Modeling the optical functions associated with a-Si 3 10 Jackson et al. [16] 2 R2 (˚ A2 ) 10 Ed 1 10 10 0 10 Eq. (4.8) (x = 6.6) Eq. (4.7) (x = 5) −1 10 1 2 3 4 5 Photon Energy (eV) 6 7 Figure 4.13: R2 (E) as a function of the photon energy, E, determined using Eqs. (4.7) and (4.8), are depicted with the solid lines. For the purposes of this fit, R2o is set to 10 ˚ A2 and Ed is set to 3.4 eV. The asymptotic polynomial rate of attenuation, x, is set to 6.6 for the purposes of the Eq. (4.8) plot. The experimental results of Jackson et al. [16] are depicted with the solid points. 129 4.4. Modeling the optical functions associated with a-Si 2 10 0 ε2 10 Jackson et al. [16] −2 10 −4 10 1 2 3 4 5 Photon Energy (eV) 6 Figure 4.14: The spectral dependence of the imaginary part of the dielectric function, 2 (E), determined using Eq. (2.17), where R2 (E) is determined using Eq. (4.8), on a logarithmic scale. The DOS modeling parameters are set to Nvo = Nco = 2.48 ×1022 cm−3 eV−3/2 , Ev = 0.0 eV, Ec = 1.68 eV, and γv = 40 meV for the purposes of this plot. The asymptotic polynomial rate of attenuation, x, is set to 6.6 for the purposes of this plot. The experimental results of Jackson et al. [16] are depicted with the solid points 130 4.4. Modeling the optical functions associated with a-Si 35 Ed 30 Jackson et al. [16] ε2 25 20 15 10 5 0 1 2 3 4 5 Photon Energy (eV) 6 Figure 4.15: The spectral dependence of the imaginary part of the dielectric function, 2 (E), determined using Eq. (2.17), where R2 (E) is determined using Eq. (4.8), on a linear scale. The DOS modeling parameters are set to Nvo = Nco = 2.48 ×1022 cm−3 eV−3/2 , Ev = 0.0 eV, Ec = 1.68 eV, and γv = 40 meV for the purposes of this plot. The asymptotic polynomial rate of attenuation, x, is set to 6.6 for the purposes of this plot. The experimental results of Jackson et al. [16] are depicted with the solid points 131 4.5. Results purposes of this analysis. 4.5 Results Three sets of experimental results are used in order to validate the aforementioned model for the spectral dependence of dielectric functions; these are the same experimental data sets as that considered earlier. In Figures 4.16, 4.17, and 4.18, the model results are compared with the experimental EGE a-Si data set [30], the PECVD a-Si I data set [31], and the PECVD a-Si II data set [32], respectively. The “best-fit” is obtained through visual inspection. For every case, this model perfectly fits with these three different sets of experimental data. The fitting parameters values are tabulated in Table 4.1. As this model is directly related to the DOS modeling parameters, it is possible to extract the underlying fundamental properties related to the underlying distributions of electronic states from knowledge of the spectral dependence of these optical functions. This information could help to predict the behavior of the optical response over a broad range of energies and can guide the fabrication of optical devices. 132 4.5. Results Table 4.1: The “best-fit” a-Si DOS modeling parameter selections employed for the purposes of this analysis. These modeling parameters relate to Eqs. (2.17), (2.45), and (4.8) parameter (units) EGE a-Si [30] PECVD a-Si I [31] PECVD a-Si II [32] Nvo (cm−3 eV−3/2 ) 1.88 × 1022 2.33 × 1022 2.42 × 1022 Nco (cm−3 eV−3/2 ) 1.88 × 1022 2.33 × 1022 2.42 × 1022 Ev (eV) 0.0 0.0 0.0 Ec (eV) 1.21 1.51 1.68 Eg (eV) 1.21 1.51 1.68 Ed (eV) 3.42 3.64 3.72 γv (meV) 83 40 50 1∞ 1.0 1.0 0.65 Ro2 (˚ A2 ) 10 10 10 x 5.87 6.6 6.6 133 4.5. Results 25 ε1 ε2 EGE a−Si [30] 20 ε1, ε2 15 10 5 0 −5 1 2 3 4 5 Photon Energy (eV) 6 7 Figure 4.16: A “best-fit” of the model results with the results of experiment corresponding to EGE a-Si [30]. These experimental results are from Piller [30]. The fits to experiment are shown with the solid lines. The experimental data, EGE a-Si [30], is depicted with the solid points. 134 4.5. Results 40 PECVD a−Si I [31] ε1, ε2 30 ε2 ε1 20 10 0 −10 1 2 3 4 Photon Energy (eV) 5 6 Figure 4.17: A “best-fit” of the model results with the results of experiment corresponding to PECVD a-Si I [31]. These experimental results are from Synowicki [31]. The fits to experiment are shown with the solid lines. The experimental data, PECVD a-Si I [31], is depicted with the solid points. 135 4.5. Results 40 PECVD a−Si II [32] ε1, ε2 30 ε2 ε1 20 10 0 −10 1 2 3 4 Photon Energy (eV) 5 Figure 4.18: A “best-fit” of the model results with the results of experiment corresponding to PECVD a-Si II [32]. These experimental results are from Ferluato et al. [32]. The fits to experiment are shown with the solid lines. The experimental data, PECVD a-Si II [32], is depicted with the solid points. 136 4.6. Implications of a non-unity 4.6 Implications of a non-unity Ideally, 1∞ 1∞ 1∞ should be equal to unity when all the electronic transitions are considered in the model for the spectral dependence of the imaginary part of the dielectric function, 2 (E) [32]. Ferlauto et al. [32] experienced that 1∞ departs from unity depending upon the band gap differences amongst the aSi samples. They also commented that 1∞ is less than unity for pure a-Si samples. From Tables 3.1 and 4.1, it is seen that for the EGE a-Si [30] and the PECVD a-Si I [31] samples, 1∞ is unity, whereas for the PECVD a-Si II [32] sample, it is less than unity, which is itself a pure a-Si material. 137 Chapter 5 Conclusions Models for the spectral dependence of the real and imaginary components of the dielectric function, 1 (E) and 2 (E), respectively, appropriate for the case of amorphous semiconductors, are considered. In the first phase of this analysis, from an empirical expression for imaginary part of the dielectric function, this expression corresponding with that of the model of Jellison and Modine [15], a closed-form expression for the real part of the dielectric function is determined using a Kramer-Kronig transformation. The resultant expression for the real and imaginary components of the dielectric function corresponds with the model of Jellison and Modine [15]. The comparison with experiment is found to be satisfactory. Then, in the latter stage of the analysis, through the application of a Kramers-Kronig transformation on an empirical model for the imaginary part of the dielectric function, this model stemming from a model for the distributions of electronic states, the spectral 138 Chapter 5. Conclusions dependence of the real part of the dielectric function. Fits with the results of experiment are also found to be satisfactory. There are a number of matters related to these models for the real and imaginary components of the dielectric function, 1 (E) and 2 (E), that could be pursued in the future. A critical comparison between the fits obtained using the model of Jellison and Modine [15] and that obtained using the new model would be instructive. The gleaning of insights into the underlying distributions of electronic states from experimental measurements of these optical functions, 1 (E) and 2 (E), would also be worthy of further inves- tigation. Finally, the use of these models for the purposes of device design and device optimization would be of interest. These topics will have to be addressed in the future. 139 References [1] G. E. Moore, “Cramming more components onto integrated circuits,” Electronics, pp. 114-117, 1965. [2] M. Riordan and H. Lillian, “The Moses of silicon valley,” Physics Today, vol. 50, pp. 42-47, 1997. [3] I. M. Ross, “The foundation of the silicon age,” Physics Today, vol. 50, pp. 34-39, 1997. [4] A. J. Snell, K. D. Mackenzie, W. E. Spear, P. G. Le Comber, and A. J. Hughes, “Application of amorphous silicon field effect transistors in addressable liquid crystal display panels,” Applied Physics Letters, vol. 24, pp. 357-362, 1981. [5] P. G. Le Comber, W. E. Spear, and A. Ghaith, “Amorphous-silicon field-effect device and possible application,” Electronics Letters, vol. 15, pp. 179-181, 1979. 140 References [6] M. G. Clark, “Current status and future prospects of poly-Si devices,” IEEE Proceedings-Circuits, Devices and Systems, vol. 141, pp. 3-8, 1994. [7] S. O. Kasap and J. A. Rowlands, “X-ray photoconductors and stabilized a-Se for direct conversion digital at-panel X-ray image-detectors,” Journal of Materials Science: Materials in Electronics, vol. 11, pp. 179198, 2000. [8] H. Fujioka, M. Oshima, C. Hu, M. Sumiya, N. Matsuki, K. Miyazaki, and H. Koinuma, “Characteristics of field effect a-Si:H solar cell,” Journal of Non-Crystalline Solids, vol. 227-230, pt. 2, pp. 1287-1290, 1998. [9] H. Pawlikiewicz and S. Guha, “Performance comparison of triple and tandem multi-junction a-Si:H solar cells: A numerical study,” IEEE Transactions on Electron Devices, vol. 37, pp. 1758-1762, 1990. [10] “Displays as key driver for large area electronics in intelligent environments: A vision for Europe 2007+”, A proposal for the advancement of the IST thematic priority, Information Society Technologies, May 2005. 141 References [11] R. A. Street, Editor, Technology and Applications of Amorphous Silicon. New York, New York: Springer-Verlag, 2000. [12] T. H. Nguyen, “A theory of occupation statistics for disordered semiconductors with applications,” Master’s Thesis, University of Regina, Regina, Saskatchewan, May 2000. [13] J. J. Thevaril, “The optical response of hydrogenated amorphous silicon,” Ph. D. Thesis, University of Windsor, Windsor, Ontario, Canada, 2011. [14] F. Orapunt, “The optical response of hydrogenated amorphous silicon,” Ph. D. Thesis, University of Regina, Regina, Saskatchewan, 2012. [15] G. E. Jellison, Jr. and F. A. Modine, “Parameterization of the optical functions of amorphous materials in the interband region,” Applied Physics Letters, vol. 69, pp. 371-373, 1996. [16] W. B. Jackson, S. M. Kelso, C. C. Tsai, J. W. Allen, and S.-J. Oh, “Energy dependence of the optical matrix element in hydrogenated amorphous and crystalline silicon,” Physical Review B, vol. 31, pp. 51875198, 1985. 142 References [17] S. K. O’Leary, “An empircal density of states and joint density of states analysis of hydrogenated amorphous silicon: a review,” Journal of Material Science: Materials in Electronics, vol. 15, pp. 401-410, 2004. [18] J. Tauc, R. Grigorovici, and A. Vancu, “Optical properties and electronic structure of amorphous germanium, Physica Status Solidi, vol. 15, pp. 627-637, 1966. [19] G. D. Cody, B. G. Brooks, and B. Abeles, “Optical absorption above the optical gap of amorphous silicon hydride,” Solar Energy Materials, vol. 8, pp. 231-240, 1982. [20] A. R. Forouhi and I. Bloomer, “Optical dispersion relations for amorphous semiconductors and amorphous dielectrics,” Physical Review B, vol. 34, pp. 7018-7026, 1986. [21] W. A. McGahan, T. Makovicka, J. Hale, and J. A. Woollam, “Modified Forouhi and Bloomer dispersion model for the optical constants of amorphous hydrogenated carbon thin films,” Thin Solid Films , vol. 253, pp. 57-61, 1994. [22] W. A. McGahan and J. A. Woollam, “Optical characterization and 143 References modeling of amorphous hydrogenated carbon films,” Materials Research Society Symposium Proceedings, vol. 349, pp. 453, 1994. [23] S. Sherman, S. Wagner, and R. A. Gottscho, “Correlation between the valence and conduction band tail energies in hydrogenated amorphous silicon,” Applied Physics Letters, vol. 69, pp. 3242-3244, 1996. [24] S. A. Alterovitz, N. Savvides, F. W. Smith, and J. A. Woollam, Handbook of Optical Constants of Solids II, edited by E. D. Palik (Academic Press, New York), 1991, pp. 837-840. [25] J. Tauc, Optical Properties of Solids, edited by S. Nudelman and S. S. Mitra (Plenum, New York), 1969, pp. 123. [26] F. Wooton, Optical Properties of Solids, (Academic press, New York), 1972, pp. 42-49. [27] S. K. O’Leary and S. M. Malik, “A simplified joint density of states analysis of hydrogenated amorphous silicon,” Journal of Applied Physics, vol. 92, pp. 4276-4282, 2002. [28] J. J. Thevaril and S. K. O’Leary, “A dimensionless joint density of states formalism for the quantitative characterization of the opti- 144 References cal response of hydrogenated amorphous silicon,” Journal of Applied Physics, vol. 107, pp. 083105-1-6, 2010. [29] S. K. O’Leary, “An analytical density of states and joint density of states analysis of amorphous semiconductors,” Journal of Applied Physics, vol. 96, pp. 3680-3686, 2004. [30] H. Piller, in Handbook of Optical Constants of Solids I, edited by E. D. Palik (Academic, New York, 1985), pp. 571-586. [31] R. Synowicki, J.A. Woollam Co., Inc.Ellipsometry Services Laboratory. [32] A. S. Ferlauto, G. M. Ferreira, J. M. Pearce,C. R. Wronski, R. W. Collins, X. Deng, and G. Ganguly, “Analytical model for the optical functions of amorphous semiconductors from the near-infrared to ultraviolet: Applications in thin film photovoltaics,” Journal of Applied Physics, vol. 92, pp. 2424-2436, 2002. [33] J. J. Thevaril and S. K. O’Leary, “The role that conduction band tail states play in determining the optical response of hydrogenated amorphous silicon,” Solid State Communications, vol. 151, pp. 730733, 2011. 145 References [34] K. Rerbal, J.-N. Chazalviel, F. Ozanam, and I. Solomon, “Measurement of band tail widths in hydrogenated amorphous silicon,” Journal of Non-Crystalline Solids, vol. 299-302, pp. 585-588, 2002. [35] T. Tiedje, J.M. Cebulka, D.L. Morel, B. Abeles, “Evidence for exponential band tails in amorphous silicon hydride,” Physical Review Letters, vol. 46, pp. 1425-1428, 1981. [36] K. Winer and L. Ley, “Surface states and the exponential valence-band tail in a-Si:H,” Physical Review B, vol. 36, pp. 6072-6078, 1987. 146
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- An optical functional analysis of amorphous semiconductors
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
An optical functional analysis of amorphous semiconductors Minar, Sanjida Begum 2012
pdf
Page Metadata
Item Metadata
Title | An optical functional analysis of amorphous semiconductors |
Creator |
Minar, Sanjida Begum |
Publisher | University of British Columbia |
Date | 2012 |
Date Issued | 2012-12-07 |
Description | Models for the spectral dependence of the real and imaginary components of the dielectric function, appropriate for the case of amorphous semiconductors,are considered. In the first phase of this analysis, from an empirical expression for the imaginary part of the dielectric function, this expression corresponding to that of Jellison and Modine [G. E. Jellison, Jr. and F. A. Modine, "Parameterization of the optical functions of amorphous materials in the interband region," Applied Physics Letters, vol. 69, pp. 371-373, 1996],a closed-form expression for the real part of the dielectric function is determined using a Kramers-Kronig transformation. The resultant expression for the real component of the dielectric function corresponds with that of the model of Jellison and Modine. The subsequent comparison with experiment is found to be satisfactory. Then, in the latter stage of this analysis, through the application of a Kramers-Kronig transformation on an empirical model for the imaginary part of the dielectric function, this model stemming from a model for the distributions of electronic states, the spectral dependence of the real part of the dielectric function is determined. Fits with the results of experiment, taken over the near-infrared, visible, and near-ultraviolet, are also found to be satisfactory. |
Genre |
Thesis/Dissertation |
Type |
Text |
Language | Eng |
Collection |
Electronic Theses and Dissertations (ETDs) 2008+ |
Date Available | 2012-12-07 |
Provider | Vancouver : University of British Columbia Library |
DOI | 10.14288/1.0073417 |
Degree |
Master of Applied Science - MASc |
Program |
Electrical and Computer Engineering |
Affiliation |
Applied Science, Faculty of (Okanagan) |
Degree Grantor | University of British Columbia |
Graduation Date | 2013-05 |
Campus |
UBCO |
Scholarly Level | Graduate |
URI | http://hdl.handle.net/2429/43665 |
Aggregated Source Repository | DSpace |
Download
- Media
- [if-you-see-this-DO-NOT-CLICK]
- ubc_2013_spring_minar_sanjida.pdf [ 1.14MB ]
- Metadata
- JSON: 1.0073417.json
- JSON-LD: 1.0073417+ld.json
- RDF/XML (Pretty): 1.0073417.xml
- RDF/JSON: 1.0073417+rdf.json
- Turtle: 1.0073417+rdf-turtle.txt
- N-Triples: 1.0073417+rdf-ntriples.txt
- Original Record: 1.0073417 +original-record.json
- Full Text
- 1.0073417.txt
- Citation
- 1.0073417.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Country | Views | Downloads |
---|---|---|
United States | 23 | 0 |
China | 19 | 2 |
Canada | 7 | 0 |
France | 4 | 0 |
Spain | 1 | 1 |
Japan | 1 | 0 |
City | Views | Downloads |
---|---|---|
Beijing | 19 | 0 |
Ashburn | 13 | 0 |
Unknown | 6 | 3 |
Kelowna | 6 | 0 |
Akron | 6 | 0 |
Denver | 1 | 0 |
Madrid | 1 | 0 |
Sunnyvale | 1 | 0 |
Paris | 1 | 0 |
Mississauga | 1 | 0 |
{[{ mDataHeader[type] }]} | {[{ month[type] }]} | {[{ tData[type] }]} |
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-0073417/manifest