Hostname: page-component-7c8c6479df-p566r Total loading time: 0 Render date: 2024-03-29T15:07:34.998Z Has data issue: false hasContentIssue false

Optical phonons and free-carrier effects in MOVPE grown AlxGa1−xN measured by Infrared Ellipsometry

Published online by Cambridge University Press:  13 June 2014

M. Schubert
Affiliation:
Center for Microelectronic and Optical Materials Research and Department of Electrical Engineering, University of Nebraska-Lincoln
A. Kasic
Affiliation:
Center for Microelectronic and Optical Materials Research and Department of Electrical Engineering, University of Nebraska-Lincoln
T.E. Tiwald
Affiliation:
Center for Microelectronic and Optical Materials Research and Department of Electrical Engineering, University of Nebraska-Lincoln
J. Off
Affiliation:
4. Physikalisches Institut, Universität Stuttgart
B. Kuhn
Affiliation:
4. Physikalisches Institut, Universität Stuttgart
Ferdinand Scholz
Affiliation:
4. Physikalisches Institut, Universität Stuttgart

Abstract

We report on the application of infrared spectroscopic ellipsometry (IR-SE) for wavenumbers from 333cm−1 to 1200cm−1 as a novel approach to non-destructive optical characterization of free-carrier and optical phonon properties of group III-nitride heterostructures. Undoped α-GaN, α-AlN, α-AlxGa1−xN (x = 0.17, 0.28, 0.5), and n-type silicon (Si) doped α-GaN layers were grown by metal-organic vapor phase epitaxy (MOVPE) on c-plane sapphire (α-Al2O3). The four-parameter semi-quantum (FPSQ) dielectric lattice-dispersion model and the Drude model for free-carrier response are employed for analysis of the IR-SE data. Model calculations for the ordinary () and extraordinary (||) dielectric functions of the heterostructure components provide sensitivity to IR-active phonon frequencies and free-carrier parameters. We observe that the α-AlxGa1−xN layers are unintentionally doped with a back ground free-carrier concentration of 1–4 1018cm−3. The ternary compounds reveal a two-mode behavior in ∈, whereas a one-mode behavior is sufficient to explain the optical response for ∈||. We further provide a precise set of model parameters for calculation of the sapphire infrared dielectric functions which are prerequisites for analysis of infrared spectra of III-nitride heterostructures grown on α-Al2O3.

Type
Research Article
Copyright
Copyright © 1999 Materials Research Society

1. Introduction

The wide band gap semiconductor systems of the group III-nitrides have gained considerable interest because of their potential application for optoelectronic devices covering the visible (VIS) spectrum with extension to near-ultraviolet (UV) wavelengths. Recent progress has enabled fabrication of GaN-based blue-light emitting diodes and laser devices, as well as high-power, high-temperature electronic devices Reference Nakamura and Fasol[1], Reference Gil[2], Reference Orton and Foxon[3], Reference Nakamura, Senoh and Mukai[4]. However, many of the physical properties of the ternary thin-film compounds Al x Ga1−x N and In x Ga1−x N, such as the lattice phonon-mode frequencies, or the anisotropy of electron and hole effective masses, are either unsecured or still unknown.

Optical phonon frequencies of Al x Ga1−x N, and coupled longitudinal optical (LO) phonon-plasmon modes in n- and p-type doped GaN have been investigated by Raman spectroscopy Reference Hayashi, Itoh, Sawaki and Akasaki[5], Reference Demangeot, Groenen, Frandon, Renucci, Briot, Clur and Aulombard[6], Reference Cros, Angerer, Stutzmann, Höpler and Metzger[7], Reference Liu, Tong, Bursill, Prawer, Nugent and Zhang[8], Reference Bergmann, Bremser, Perry, Davis, Dutta and Nemanich[9], Reference Demangeot, Frandon, Renucci, Grandjean, Beaumont, Massies and Gibart[10]. Hayashi et al. Reference Hayashi, Itoh, Sawaki and Akasaki[5] reported first Raman investigation of wurtzite AlGaN layers grown by MOVPE on c-plane sapphire. The modes observed followed the wurtzite GaN and AlN selection rules, and a one-mode behavior was stated for the mixed crystals. In most of the subsequent Raman studies, the one-mode behavior is favored over that of a two-mode behavior.

Infrared (IR) optical investigations of III-nitride heterostructures, especially on heavily anisotropic substrates such as α-Al2O3, have not been exhaustive, and more detailed studies are needed Reference Wetzel, Haller, Amano and Akasaki[11], Reference Perlin, Litwin-Staszewska, Suchanek, Knap, Camassel, Suski, Piotrzkowski, Grzegory, Porowski, Kaminska and Chervin[12], Reference Yu, Ishikawa, Umeno, Egawa, Watanabe, Soga and Jimbo[13], Reference Wisniewski, Knapp, Malzac, Camassel, Bremser, Davis and Suski[14]. Nevertheless, very recently Wisniewski et al. presented experimental evidence of a two-mode behavior from near-normal IR-reflectometry data for hexagonal Al x Ga1−x N layers grown on 6H-SiC substrates Reference Wisniewski, Knapp, Malzac, Camassel, Bremser, Davis and Suski[14]. The two transverse modes observed were assigned to the AlN-like and GaN-like E 1(TO) phonon frequencies. The authors extrapolate the Al (Ga) impurity mode in GaN (AlN), and obtained 643cm−1 (622cm−1), respectively. Yu et al. presented near-normal IR-reflectometry data of α-Al x Ga1−x N grown on sapphire that also indicated a two-mode phonon behavior Reference Yu, Ishikawa, Umeno, Egawa, Watanabe, Soga and Jimbo[13]. A disorder-activated Raman-active E 2 mode with polarization parallel to the layer interfaces was assigned to the second TO mode Reference Yu, Ishikawa, Umeno, Egawa, Watanabe, Soga and Jimbo[13]. In view of the results presented by Wisniewski et al., and those reported in this work, the second mode observed by Yu et al.is most likely the second E 1(TO) phonon frequency in α-Al x Ga1−x N rather than a disorder-activated Raman mode. Wetzel et al. Reference Wetzel, Haller, Amano and Akasaki[11] investigated a Al0.15Ga0.85N/GaN/AlN heterostructure on c-plane sapphire by near-normal IR-reflectometry, and reported on several unidentified features within the reststrahlen band of sapphire. In particular, a transverse resonance occurred near 651cm−1, which we now would interpret as the second AlN-like E 1(TO) mode. The authors there were indecisive whether the resonance would belong to the A 1(TO) or E 1(TO) frequency, but the A 1(TO) mode is forbidden for the c-plane oriented wurtzite sample structure.

Many authors have treated the IR optical response of the group-III nitride layers isotropic, despite the wurtzite structure and the fact that the Berreman resonance Reference Berreman[15] at the A 1(LO) frequency was observed generally Reference Wetzel, Haller, Amano and Akasaki[11], Reference Yu, Ishikawa, Umeno, Egawa, Watanabe, Soga and Jimbo[13], Reference Wisniewski, Knapp, Malzac, Camassel, Bremser, Davis and Suski[14]. The observation of the Berreman effect in near-normal incidence reflectometry data is most indicative of the fact that the AlGaN layers need to be considered anisotropic. Wisniewski et al. have this structure mistakenly interpreted as the E 1(LO) mode which cannot be excited by the IR light beam incident on c-axis oriented uniaxial layers Reference Schubert, Rheinländer, Franke, Neumann, Tiwald, Woollam, Hahn and Richter[16]. Instead, the dielectric loss for polarization parallel to the film normal is probed by the p-polarized components of the incident light, which is known as the Berreman effect (The reflected p-polarized IR beam changes its amplitude and phase across the wavenumber range where the film dielectric function parallel to the growth direction is near to zero. Reference Berreman[15], Reference Schubert, Rheinländer, Franke, Neumann, Tiwald, Woollam, Hahn and Richter[16]) Wetzel et al. also used isotropic treatment of their AlGaN crystals due to their small angle-of-incidence IR-reflectometry setup, but assigned E 1(TO) and A 1(LO) modes to their spectra Reference Wetzel, Haller, Amano and Akasaki[11]. Because the E 1 and A 1 modes belong to different crystal directions, proper assignment of these modes requires anisotropic data treatment, and data analysis in Reference Wetzel, Haller, Amano and Akasaki[11] is entangled.

Spectroscopic Ellipsometry (SE) is known as excellent technique for measurement of thin-film optical properties. In contrast to reflectometry, ellipsometry is not or only little influenced by intensity fluctuations and the spectral intensity distribution of the light source. Hence, as long as some of the IR-beam reaches the detector, IR-ellipsometry is not influenced by atmospheric absorption lines. In general, ellipsometry is superior to reflectometry because two parameters, instead of one, are determined independently at each wavelength and angle of incidence. The ellipsometry parameters can be measured with high accuracy, and small systematic errors. Both real and imaginary parts of the complex dielectric response function of materials can be measured without the necessity to perform further Kramers-Kronig analysis. The two parameters also place tighter constrains to model calculations than the single rather unsecured intensity parameter. Therefore, determination of the optical properties from very thin layers can be achieved by ellipsometry Reference Jellison[17]. Recent development of polarization sensitive spectrophotometers in the middle- and far-IR now makes feasible the application of SE to photon energies that match the phonon energies of group III-V semiconducting compounds Reference Thompson, DeVries, Tiwald and Woollam[18]. Within the IR spectral range, resonant excitation of phonons by the SE probe beam strongly affects the state of polarization of the reflected beam, thereby providing high sensitivity to lattice properties of thin films. Raman scattering of light by optical phonon or plasmon excitations tends to be weaker than reflection and absorption processes, because higher-order optical interactions are involved in the scattering process whereas reflection and absorption involve the lowest order of interaction between electromagnetic waves and elementary excitations Reference Yu and Cardona[19]. Light scattering in thin layers, and especially in multilayered semiconductor structures is extremely weak, and difficult to detect. This figure changes for resonant Raman scattering, but strong light sources with wavelengths that match the band gaps of the materials are required. Sometimes, information from buried layers may still not be accessible when the embedding materials absorb the Raman beam.

In this paper we report on recent results of IR-SE investigations on III-nitride heterostructures grown by MOVPE on c-plane sapphire. We demonstrate that IR-ellipsometry allows non-destructive characterization of free-carrier and phonon properties of multilayer samples, even if the optical response is heavily anisotropic as in the case of the wurtzite nitrides investigated here. We discuss issues of data analysis and parameter accessibility, and compare our results to recent findings by other authors. In particular, we address the two-phonon behavior problem in AlGaN, and the electron effective mass anisotropy in silicon-doped wurtzite GaN.

2. Theory

The SE parameters Ψ and Δ are defined by the complex ratio of the p- and s-polarized reflectance coefficients r p and r s, respectively Reference Jellison[17]

(1)

and depend on the angle of incidence Φ a , the thickness d of each layer, on the anisotropic dielectric functions j of the substrate, and on the anisotropic dielectric functions of all layers of the heterostructure Reference Jellison[17] Reference Schubert[20] Reference Schubert[21]. In the infrared spectral range lattice dynamics ( L), and free-carrier properties ( FC) contribute to j = j L + j FC . The four-parameter semi-quantum (FPSQ) model for multiple (l) polar phonon mode materials is used here to calculate j for electric field polarization E along j = “x”, “y”, “zReference Gervais and Piriou[22], Reference Schubert, Tiwald and Herzinger[23]

(2)

where ω TOi , ω LOi , γ TOi , and γ LOi are the transverse optical (TO) and LO-phonon frequencies and broadening parameters, respectively. Free carrier absorption affects the high-frequency conductivity, and contributes to j Reference Wolfe, Holonyak and Stillmann[24]

(3)

where n is the carrier concentration, and m j * their effective mass. (e is the electron charge. 0 is the vacuum dielectric constant.) We assume that the complex average momentum lifetime <τ m *> = <τ m (1−iωτ m )−1> equals <τ m >(1 − i ω<τ m >)−1 Reference Wolfe, Holonyak and Stillmann[24]. We substitute <τ m > with the carrier mobility μ j

(4)

For the wurtzite III-nitrides, different lattice modes [(A 1): E || c; (E 1): Ec], as well as different values m j * and μ j for E || c ( z = ||), and Ec ( x = y = ) are considered. Model calculations are employed for analysis of the IR-SE data Reference Jellison[17], and the model parameters (d and ω TOi , ω LOi , γ TOi , γ LOi , ∈j , n, m j *, μ j ; of each layer [j = “⊥”, “||”]), are varied until calculated and measured data match as closely as possible. The free carrier model requires 5 input parameters (n, m || *, μ ||, m *, μ ), but provides information, at best, to four of those. In particular, only the ratios m ⊥,|| */N, and 1/(μ ⊥,|| N) affect the dielectric response ⊥,||. If the carrier concentration is known, the remaining parameters may be obtained. Likewise, if the effective masses are known, carrier mobility and concentration may be determined. Note that the lattice-mode parameters and the high-frequency dielectric function also affect the optical free-carrier response. In particular, the high-frequency dielectric constants screen the plasma excitations and LO modes.

Note that although the number of unknown parameters is apparently large, the number of independent data points acquired during the ellipsometry experiment exceeds by far the number of fit parameters. During the course of the analysis of the sample presented here, we successively build a data reference basis for the different Al x Ga1−x N compounds investigated. We started to measure and to thereafter model the optical response of the sapphire substrate, and continued to determine the AlN and GaN dielectric functions. In doing so, the number of unknown parameters for analysis of the more complex heterostructures was reduced to those of the ternary layers, or in the case of the doped GaN, to the free-carrier parameters only. As a result, we obtained very stable fits. There is, however, a limitation of parameter accuracy due to limited sensitivity caused by the geometry of the experiment. Strong absorption of the incident polarized light at the A 1(LO) and E 1(TO) frequencies leads to resonant structures within the Ψ and Δ spectra. The latter modes can be determined precisely. Due to the symmetry of the samples (the wurtzite group-III nitride films have their c-axes parallel to the sample normal), the A 1(TO) and E 1(LO) modes are forbidden, and cannot be excited by the IR probe beam. Although claimed by other authors Reference Wetzel, Haller, Amano and Akasaki[11], Reference Wisniewski, Knapp, Malzac, Camassel, Bremser, Davis and Suski[14], resonance features due to these modes cannot be observed in IR reflectometry or ellipsometry experiments. Nevertheless, both ordinary and extraordinary dielectric functions influence the complex reflection coefficients (and therefore the ellipsometric parameters, see Eqs. C5 and C6 within the Appendix C in Reference Schubert[20]), and the A 1(TO) and E 1(LO) frequencies can still be obtained from the line-shape analysis of multiple angle-of-incidence data. However, the error bars for these frequencies will be much larger than for the A 1(LO) and E 1(TO) frequencies.

3. Experimental

Undoped AlN, GaN, and Al x Ga1−x N epilayers (x = 0.17, 0.28, 0.5), and Si-doped n-type GaN epilayers with wurtzite structure ([ 0 0 0 1 ] || to the substrate normal) were grown by MOVPE on c-plane sapphire Reference Scholz, Härle, Bolay, Steuber, Kaufmann, Reyher, Dörnen, Gfrörer, Im and Hangleiter[25], Reference Scholz, Härle, Steuber, Bolay, Dörnen, Kaufmann, Syganow and Hangleiter[26], Reference Im, Kollmer, Off, Sohmer, Scholz and Hangleiter[27]. TEGa, TMAl, and ammonia were used as precursors. Prior to the epilayers an AlN-buffer layer was deposited. The free-electron concentration n, and the lateral mobility μ of the n-type GaN layers were determined from electrical Hall measurements. All samples were measured at room temperature for wavenumbers from 333cm−1 to 1200cm−1 with 2cm−1 resolution, and at multiple angles of incidence (55° – 75°) with a rotating-polarizer, rotating-compensator, Fourier-transform-based variable angle-of-incidence spectroscopic ellipsometer Reference Thompson, DeVries, Tiwald and Woollam[18].

4. Results and Discussion

Figure 1 and Figure 2 show experimental data (symbols), and best-fit calculations (solid lines) of Ψ and Δ for (a) c - plane sapphire (α-Al2O3), (b) undoped α-GaN(793nm)/α-AlN(80nm)/α-Al2O3, (c) α-Al0.17Ga0.83N(556nm)/α-AlN(49nm)/α-Al2O3, (d) α-Al0.28Ga0.72N(535nm)/α-AlN(126nm)/α-Al2O3, (e) α-Al0.50Ga0.50N(316nm)/α-AlN(61nm)/α-Al2O3, and (f) α-AlN(243nm) /α-Al2O3. (The thickness information corresponds to our best-fit values. Data from 72° angle of incidence are shown only, and presentation of further angle-of-incidence data is suppressed to add clarity to the figures.) Excellent agreement is obtained between the best-fit calculation and the experimental data. Because sapphire was the substrate for the III-nitride films, we first determined || and for α-Al2O3 precisely. (The existing knowledge on the IR dielectric response had to be revised, and details will be given in Ref. Reference Schubert, Tiwald and Herzinger[23].) Following the common notation Reference Gervais and Piriou[22], Reference Barker[28], our best-fit values {ω TO, ω LO, γ TO, γ LO} for the (E u ) and (A 2u ) phonon modes of sapphire are given in Table 1. Equation (2) then reproduces || and . The reststrahlen bands of sapphire, indicated by the brackets below spectra (a) in Figure 1 and Figure 2 dominate spectra (b) - (f). However, the lattice modes of the III-nitride films affect the IR-SE data, especially at the (E 1) ω TO (solid arrows), and (A 1) ω LO (dotted arrows) frequencies.

Table 1. Room temperature transverse and longitudinal optical phonon frequencies and broadening parameters for α-Al2O3 in units of cm−3 (∈∞⊥ = 3.077, ∈∞|| = 3.070) [23]

Figure 1. Experimental data (symbols), and best-fit calculation (solid lines) of Ψ from (a) α-Al2O3 (c-plane), (b) undoped α-GaN, (c) α-Al0.17Ga0.83N, (d) α-Al0.28Ga0.72N, (e) α-Al0.50Ga0.50N, and (f) α-AlN. Vertical arrows indicate frequencies at which resonant excitation of E 1(TO) (solid arrows), and A 1(LO) phonons (dotted arrows) occur. The vertical lines denote the modes of the AlN-nucleation layer (f), as these are present in spectra (b) - (e) as well. The angle of incidence is Φ a = 72°. Spectra are to scale, but shifted for convenience by 60° each. The brackets below spectra (a) indicate the sapphire A 2u phonon frequencies (TO: solid, LO: dotted). See also Figure 2. The sapphire phonon modes are given in the text.

Figure 2. Same as Figure 1 for cosΔ. The brackets below spectra (a) indicate the E u modes of sapphire. See text for further explanation.

The Ψ and Δ spectra (f) from the 243nm-thick single layer of AlN are sufficient to calculate || and (best-fit values in Table 2 and Table 3). The AlN buffer layer in the GaN (b), and the α-Al x Ga1−x N samples (c) - (e) can be clearly identified by the resonant absorption at the AlN phonon mode frequencies. The vertical lines mark (E 1) ω TO, and (A 1) ω LO for α-AlN. Spectra (c) allow derivation of the α-GaN dielectric functions. Note the sharp resonance near (A 1) ω LO due to absorption of p-polarized light near the LO-phonon frequency for E||c, which is well-known as the Berreman-effect Reference Berreman[15]. (Dotted arrows at spectra (b) in Figure 1 and Figure 2;. See Table 2 and Table 3 for the best-fit parameters of α-GaN.) Spectra (c) - (e) of the α-Al x Ga1−x N compounds are influenced by free-carrier effects. Table 4 contains the free-carrier parameters obtained from IRSE data analysis or Hall measurements. The dotted arrows in Figure 1 and Figure 2 denote the high-frequency coupled LO-phonon-plasmon modes at ~ 811cm−1 (x = 0.17), ~ 870cm−1 (x = 0.28), and ~ 830cm−1 (x = 0.50) observed here. The free-carrier effective mass was assumed to be isotropic, and set at a fixed value of 0.22m 0. No significant anisotropy was obtained for the mobility, and hence we assumed μ ~ μ || = μ.

Table 2. IR (A 1) lattice-mode best-fit parameters for ∈|| of the Al x Ga1−x N samples investigated here. The (A 1) ω TO (ω LO) frequencies have large, <±25% {small, <±1%} uncertainty due to the c-plane orientation of the III-N films (a) Isotropically averaged between “||” and “⊥”). See also Table 3.

Table 3. IR (E 1) lattice-mode best-fit parameters for ∈ of the Al x Ga1−x N samples investigated here. The (E 1) ω TO (ω LO) frequencies have small, <±1% {large, <±25%} uncertainty due to the c-plane orientation of the III-N films (a) Isotropically averaged between “||” and “⊥”). See also Table 2.

In Figure 1 and Figure 2, the GaN-like (E 1 GaN) ω TO frequency shifts gradually toward that of AlN with increasing x. A second (E 1 AlN) TO mode occurs below the AlN (E 1) ω TO frequency (solid arrows). The best-fit lattice mode parameters for α-Al x Ga1−x N are given in Table 2 and Table 3. As mentioned above, the (A 1) TO and (E 1) LO modes cannot be directly observed in this type of IR-SE experiment. These parameters follow from the lineshape analysis of the IR-SE spectra, and contain much larger error bars. Similar to Wisniewski et al. Reference Wisniewski, Knapp, Malzac, Camassel, Bremser, Davis and Suski[14], and Yu et al. Reference Yu, Ishikawa, Umeno, Egawa, Watanabe, Soga and Jimbo[13], we observe two transverse resonances of E 1 symmetry within the ternary α-Al x Ga1−x N compounds. Two polar modes are included in Equation 2 for calculation of the AlGaN dielectric response perpendicular to the c-axis, i.e., for . In agreement with the assignment of Wisniewski et al., we identify one mode as the AlN-like E 1 transverse-optical sublattice vibration in wurtzite AlGaN (“E 1 AlN” in Table 3). The GaN-like sublattice vibration frequencies (“E 1 GaN” in Table 3) are in good agreement with the E 1-modes observed by Hayashi et al. using Raman scattering Reference Hayashi, Itoh, Sawaki and Akasaki[5]. Assuming a linear dependence of the GaN-like and AlN-like E 1 modes on the alloy composition x, we estimate the impurity modes for Al (Ga) in GaN (AlN) as I Al = 630cm−1 (I Ga = 605cm−1). These findings agree well with the results obtained by Wisniewski et al. for AlGaN films grown on 6H-SiC (I Al = 643cm−1 , I Ga = 622cm−1) Reference Wisniewski, Knapp, Malzac, Camassel, Bremser, Davis and Suski[14]. Both the different lattice mismatch and the isotropic treatment of the reflectometry data may cause the deviations between the impurity modes obtained here and those reported by Wisniewski et al. Besides the two E 1(TO) frequencies in Table 3, two E 1(LO) frequencies are given which follow from the best-fit analysis of . The LO frequencies are used as parameters to model the amount of polarity of the respective phonon branch. The low-frequency LO mode is silent according to the best-fit line shape of ( remains less than zero between the GaN-like E 1(TO) and the high-frequency E 1(LO) mode). For polarization parallel to the layer c-axis, i.e., for ||, we did not observe such a two-mode behavior. There, a single polar lattice mode was sufficient to model the ellipsometric data. To the best of our knowledge, the A 1(TO) phonon-mode behavior of wurtzite AlGaN has not been discussed within reports of IR investigations. Wisniewski et al. Reference Wisniewski, Knapp, Malzac, Camassel, Bremser, Davis and Suski[14], Yu et al. Reference Yu, Ishikawa, Umeno, Egawa, Watanabe, Soga and Jimbo[13], and Wetzel et al. Reference Wetzel, Haller, Amano and Akasaki[11] did not attempt to extract anisotropy information from the near-normal IR-reflectometry data, and treated their samples isotropic. Grille, Schnittler and Bechstedt have studied the lattice dynamics of random wurtzite-AlGaN alloys using a modified random-element isodisplacement and a rigid-iron model. There the authors obtain a two-mode behavior for the E 1(TO) mode, and a one-mode behavior for the A 1(TO) mode. As well, both E 1(LO) and A 1(LO) modes obey a one-mode behavior Reference Grille, Schnittler and Bechstedt[29]. Furthermore, a second E 1(LO) mode is always existent within this model, but silent throughout the compositional range (zero spectral weight) Reference Grille[30]. These recent theoretical results agree excellent with our observations reported here.

Figure 3 shows experimental data (symbols), and best fit calculations (solid lines) of Ψ from Si-doped n-type α-GaN/α-AlN/α-Al2O3 samples [(a) undoped, same as (c) in Figure 1, (b) n = 7 1018cm−3, d GaN = 950nm, d AlN = 80nm, (c) n = 1 1019cm−3, d GaN = 920nm, d AlN = 80nm]. The concentration of free carriers increases from top to bottom in Figure 3, where the lattice resonance bands of the substrate and the epilayers are screened gradually, as can be seen from spectra (a) - (c). The strong resonance at the GaN (A 1) ω LO frequency in (a) disappears in spectra (b), (c) due to coupling of LO phonon-plasmon modes. However, a small dip still occurs in spectra (b), (c) which indicates a surface-depletion layer (d dep), where the carrier density is too small for screening of the LO phonon resonance [(b) d dep ~ 20nm, (c) d dep ~ 10nm]. The lattice modes of the GaN and AlN layer, as well as those of the sapphire substrate, can still be identified. The vertical lines indicate the (E 1) ω TO, and (A 1) ω LO frequencies of GaN (dash-dotted line), and AlN (dotted line). The lattice parameters for spectra (b), (c) are found to be consistent with those for the undoped epilayer. The remaining free-carrier best-fit parameters are given in Table 4. The lateral mobilities μ are found to be in fair agreement with electrical mobility data extracted from Hall measurements. The vertical mobility is higher than the lateral mobility, which can be explained by the columnar growth of the GaN films. The free-carrier concentration was assumed from Hall measurements. We could not observe substantial anisotropy of the free-carrier effective mass. However, there is a small indication that the effective mass m* ~ 1.01m* ||.. 1.05m* ||. No significant change of the electron effective mass with increase in carrier concentration is obtained. These results from GaN films agree with those reported by Perlin et al. obtained from IR-reflectometry investigations on n-type bulk GaN crystals Reference Perlin, Litwin-Staszewska, Suchanek, Knap, Camassel, Suski, Piotrzkowski, Grzegory, Porowski, Kaminska and Chervin[12].

Figure 3. Experimental data (symbols), and best-fit calculation (solid lines) of Ψ from differently Si-doped n-type α-GaN films. Vertical lines indicate the E 1(TO), and A 1(LO) phonon frequencies of the GaN epilayer (dash-dotted), and the AlN buffer layer (dotted).

Table 4. Free-carrier parameters obtained from the IRSE data analysis and Hall measurements of the AlxGa1-xN films investigated here. (a) assumed values from Hall measurements; b) isotropically averaged between “⊥”, and “||”, but varied during IRSE best-fit analysis; c) assumed values, fixed during best-fit analysis; d) sample not measured)

As can be seen in Figure 3, deviations between our calculated and measured data increase with increasing concentration of silicon as dopant, and hence with the increase of the number of free carriers. We have identified additional IR active resonance structures within Si-doped MBE and MOVPE wurtzite GaN samples grown on c-plane sapphire similar to those investigated here. We will report in a follow up manuscript on further results about the IR optical properties of Mg-doped p and Si-doped n-type α-GaN films Reference Kasic, Schubert, Tiwald, Einfeldt, Böttcher, Hommel, Off, Kuhn and Scholz[31].

5. Conclusion

To conclude, we have used IR-SE for non-destructive optical characterization of lattice and free-carrier properties of group-III nitride heterostructures. The IR-SE data provide sensitivity to thickness, optical lattice modes, and free-carrier properties of each layer within the sample. Model calculations allow quantification of anisotropic phonon frequencies, carrier effective masses and concentration, as well as optical mobility values. So far, simple approximations are used for parameterization of the materials dielectric response. Further studies will consider frequency dependence of carrier life-time, as well as anharmonicity in the plasmon-LO-phonon coupling mechanism.

Part of this work is supported by the DFG under SPP1032 (“Gruppe III-nitride und ihre Heterostrukturen”) grant Nr. Rh 28-3/1 (Germany). Part of this work was also supported by the CMOMR at UNL (U.S.A.). The authors thank Bernd Rheinländer at Univerität Leipzig for continuing interest. M.S. and A.K. thank Daniel W. Thompson at UNL for technical support. We also acknowledge Craig M. Herzinger at J.A.Woollam Co., Lincoln, for many stimulating discussions. We are indebted to Dr. Hartmut Grille, Technische Universität Ilmenau, and Professor Friedhelm Bechstedt, Friedrich-Schiller Universität Jena, for providing us with their results prior to publication.

References

Nakamura, Shuji, Fasol, Gerhard, The Blue Laser Diode - GaN based Light Emitters and Lasers (Springer-Verlag, Heidelberg, 1997).Google Scholar
Gil, B., “Group III Nitride Semiconductor Compounds: Physics and Applications” (Clarendon Press, Oxford, 1998)Google Scholar
Orton, J. W, Foxon, C. T., Rep. Prog. Phys. 61, 1 (1998).CrossRefGoogle Scholar
Nakamura, Shuji, Senoh, Masayuki, Mukai, Takashi, Appl. Phys. Lett. 62, 2390-2392 (1993).CrossRefGoogle Scholar
Hayashi, K., Itoh, K., Sawaki, N., Akasaki, I., Sol. St. Comm. 77, 115 (1991).CrossRefGoogle Scholar
Demangeot, F., Groenen, J., Frandon, J., Renucci, M. A., Briot, O., Clur, S., Aulombard, R. L., Appl. Phys. Lett. 72, 2674 (1998).CrossRefGoogle Scholar
Cros, A., Angerer, H., Stutzmann, M., Höpler, R., Metzger, T., Sol. St. Comm. 104, 35 (1997).CrossRefGoogle Scholar
Liu, M. S., Tong, Y. Z., Bursill, L. A., Prawer, S., Nugent, K. W., Zhang, G. Y., Sol. St. Comm. 108, 765 (1998).CrossRefGoogle Scholar
Bergmann, L., Bremser, M. D., Perry, W. G., Davis, R. F., Dutta, M., Nemanich, R. J., Appl. Phys. Lett. 71, 2157 (1997).CrossRefGoogle Scholar
Demangeot, F, Frandon, J, Renucci, MA, Grandjean, N, Beaumont, B, Massies, J, Gibart, P, Sol. St. Comm. 106, 491-494 (1998).CrossRefGoogle Scholar
Wetzel, C., Haller, E. E., Amano, H., Akasaki, I., Appl. Phys. Lett. 68, 2547-2549 (1996).CrossRefGoogle Scholar
Perlin, P., Litwin-Staszewska, E., Suchanek, B., Knap, W., Camassel, J., Suski, T., Piotrzkowski, R., Grzegory, I., Porowski, S., Kaminska, E., Chervin, J. C., Appl. Phys. Lett. 68, 1114-1116 (1996).CrossRefGoogle Scholar
Yu, G., Ishikawa, H., Umeno, M., Egawa, T., Watanabe, J., Soga, T., Jimbo, T., Appl. Phys. Lett. 73, 1472 (1998).CrossRefGoogle Scholar
Wisniewski, P., Knapp, W., Malzac, J. P., Camassel, J., Bremser, M. D., Davis, R. F., Suski, T., Appl. Phys. Lett. 73, 1760 (1998).CrossRefGoogle Scholar
Berreman, D. W., Phys. Rev. 130, 2193 (1963).CrossRefGoogle Scholar
Schubert, M., Rheinländer, B., Franke, E., Neumann, H., Tiwald, T. E., Woollam, J. A., Hahn, J., Richter, F., Phys. Rev. B 56, 13306 (1997).CrossRefGoogle Scholar
Jellison, G. E., Thin Sol. Films 313, 33 (1998).CrossRefGoogle Scholar
Thompson, D. W., DeVries, M. J., Tiwald, T. E., Woollam, J. A., Thin Sol. Films 313, 341 (1998).CrossRefGoogle Scholar
Yu, P. Y., Cardona, M., “Fundamentals of Semiconductors”, Springer-Verlag, Berlin-Heidelberg-New York, 1996 Google Scholar
Schubert, M., Phys. Rev. B 53, 4265 (1996).CrossRefGoogle Scholar
Schubert, M., Thin Sol. Films 313, 323 (1998).CrossRefGoogle Scholar
Gervais, F., Piriou, B., J. Phys. (Paris) C7, 2374 (1974).Google Scholar
Schubert, M., Tiwald, T. E., Herzinger, C. M., unpublished (2000).Google Scholar
Wolfe, C. M., Holonyak, N., Stillmann, G. E., “Physical Properties of Semiconductors” (Prentice Hall, New Jersey, 1989)Google Scholar
Scholz, F., Härle, V., Bolay, H., Steuber, F., Kaufmann, B., Reyher, G., Dörnen, A., Gfrörer, O., Im, S. -J., Hangleiter, A., Sol. St. Electr. 41, 141 (1997).CrossRefGoogle Scholar
Scholz, F., Härle, V., Steuber, F., Bolay, H., Dörnen, A., Kaufmann, B., Syganow, V., Hangleiter, A., J. Cryst. Growth 170, 321-324 (1997).CrossRefGoogle Scholar
Im, J. S., Kollmer, H., Off, J., Sohmer, A., Scholz, F., Hangleiter, A., Phys. Rev. B 57, R9435 (1998).Google Scholar
Barker, A. S., Phys. Rev. 132, 1474 (1963).CrossRefGoogle Scholar
Grille, H., Schnittler, Ch., Bechstedt, F., unpublished (1999).Google Scholar
Grille, H., private communicationsGoogle Scholar
Kasic, A., Schubert, M., Tiwald, T. E. Einfeldt, S., Böttcher, T., Hommel, D., Off, J., Kuhn, B., Scholz, F., unpublished (1999).Google Scholar
Figure 0

Table 1. Room temperature transverse and longitudinal optical phonon frequencies and broadening parameters for α-Al2O3 in units of cm−3 (∈∞⊥ = 3.077, ∈∞|| = 3.070) [23]

Figure 1

Figure 1. Experimental data (symbols), and best-fit calculation (solid lines) of Ψ from (a) α-Al2O3 (c-plane), (b) undoped α-GaN, (c) α-Al0.17Ga0.83N, (d) α-Al0.28Ga0.72N, (e) α-Al0.50Ga0.50N, and (f) α-AlN. Vertical arrows indicate frequencies at which resonant excitation of E1(TO) (solid arrows), and A1(LO) phonons (dotted arrows) occur. The vertical lines denote the modes of the AlN-nucleation layer (f), as these are present in spectra (b) - (e) as well. The angle of incidence is Φa= 72°. Spectra are to scale, but shifted for convenience by 60° each. The brackets below spectra (a) indicate the sapphire A2u phonon frequencies (TO: solid, LO: dotted). See also Figure 2. The sapphire phonon modes are given in the text.

Figure 2

Figure 2. Same as Figure 1 for cosΔ. The brackets below spectra (a) indicate the Eu modes of sapphire. See text for further explanation.

Figure 3

Table 2. IR (A1) lattice-mode best-fit parameters for ∈|| of the AlxGa1−xN samples investigated here. The (A1) ωTO (ωLO) frequencies have large, <±25% {small, <±1%} uncertainty due to the c-plane orientation of the III-N films (a) Isotropically averaged between “||” and “⊥”). See also Table 3.

Figure 4

Table 3. IR (E1) lattice-mode best-fit parameters for ∈ of the AlxGa1−xN samples investigated here. The (E1) ωTO (ωLO) frequencies have small, <±1% {large, <±25%} uncertainty due to the c-plane orientation of the III-N films (a) Isotropically averaged between “||” and “⊥”). See also Table 2.

Figure 5

Figure 3. Experimental data (symbols), and best-fit calculation (solid lines) of Ψ from differently Si-doped n-type α-GaN films. Vertical lines indicate the E1(TO), and A1(LO) phonon frequencies of the GaN epilayer (dash-dotted), and the AlN buffer layer (dotted).

Figure 6

Table 4. Free-carrier parameters obtained from the IRSE data analysis and Hall measurements of the AlxGa1-xN films investigated here. (a) assumed values from Hall measurements; b) isotropically averaged between “⊥”, and “||”, but varied during IRSE best-fit analysis; c) assumed values, fixed during best-fit analysis; d) sample not measured)