Undoped and Si-doped GaN films were grown by low pressure MOCVD on (0001) sapphire substrates. The angular distribution of the X-ray diffraction corresponding to the (0002), (0004), (100), (200), and (114) reflections has been measured by means of double- and triple -crystal diffractometry with Mo Ka₁ and Cu Ka₁ radiation under conditions of symmetrical and asymmetrical Bragg- and Laue-geometry. In our experiments a non-coplanar geometry was also applied. On the basis of the performed studies, five independent components of the tensor of microdistortion were evaluated and the average grain-size in two directions was determined. The type, position, and density of dislocations were established as well. The role of dislocations in strain relaxation and their influence on the optical and electrical properties are discussed.

1. Introduction

In recent years remarkable progress in the fabrication of bright UV and blue-green light emitting diodes and laser diodes has been achieved for III-N epilayers grown mainly on (0001) sapphire substrates [1]. The peculiarities specific to the growth of GaN films on lattice-mismatched sapphire substrates cause remarkable macro- and micro deformations in the epilayers for different thermal expansion coefficients of the film and the substrate. X-ray diffraction is used to study the material structural properties as well as macro- and micro-strains in GaN and related compounds. Defects in epitaxial layers are usually inspected using measurement of the angular full width at half maximum (FWHM) of the X-ray diffraction rocking curves. However, it is insufficient to establish a type of defect and to clarify the mechanism of stress relaxation and its influence on the optical and electrical properties of the epilayers. In this paper both undoped and Si-doped GaN epilayers have been studied comprehensively by means of double- and triple-crystal diffractometry.

2. Experiment

The GaN single epilayers with specular surfaces, both undoped and doped with Si, were grown in a conventional MOCVD system at a low pressure, 200 mbar. Ammonia and trimethylgallium were applied as component precursors. Purified hydrogen was used as a carrier gas. Silane was used for doping. The substrates were (0001) optical-grade polished sapphire. The growth procedure included GaN nucleation layer deposition at a low temperature of 500 ¡C followed by epilayer growth at a higher temperature of 1040 ¡C The typical layer thickness was 2-4 µm. The two-dimensional growth mode appeared to be dominant on the specular surfaces. Hexagonal growth features were lacking, whereas growth steps of ~0.3 nm height were well defined by AFM.

The angular distribution of the X-ray diffraction corresponding to the (0002), (0004), (100), (200), and (114) reflections was measured by means of triple - crystal differential diffractometry with Mo K_a1 and Cu K_a1 radiation under conditions of symmetrical and asymmetrical Bragg- and Laue-geometry. A rough drawing of measurement schemes used for the different measurement geometries is given in Figure 1. Diffraction curves were measured for q- and q-2q scanning modes. The double- crystal scheme for diffraction in non-coplanar geometry was also applied.

Schema for the diffracted beam intensity distribution in the scattering plane for symmetrical and asymmetrical Bragg- and Laue-geometry are given in Figure 2 and Figure 3. The measured values of FWHM are designated as w_h^B, w_h^L and w_a^B, w_a^L for q- and q-2q scanning mode, in Bragg- and Laue-geometry respectively. In nonplanar geometry w_nc was measured.

On the basis of these parameters, five independent components of the microdistortion tensor e_zz , e_zx , e_xx ,e_xz , and e_xy were evaluated and the average sizes of grains (average size of the coherent scattering areas ) both in the layer plane t_x, and along the c-axis t_z, were determined. This approach is discussed in detail in another paper. [2] These components of the tensor are related to the type, position, and density of dislocations. The microdistortion tensor components and the average size of the coherent scattering areas (domain size), t_x and t_z, are given in the Table 1 for undoped and Si-doped GaN epilayers. It should be noted that the Si concentration in the undoped GaN layer was taken to be 5á10¹⁶ cm^-3, corresponding to the sensitivity limit of SIMS.

The values of the sample curvature radii R and lattice parameters a and c were utilized to evaluate the biaxial stress, s_a, and the macrodeformations, e_a and e_c, using the expressions:

sa = -(1/6R){E/(1-n)}s (ts2/tl)

ea = (a - ao)/ao

ec = (c-co)/co ,

where s and l denote the substrate and layer, respectively. The accuracy of the values of compressive stress, dislocation density and domain size determined from the measurements did not exceed 30%. The parameters a₀ and c₀ were taken from reference [3].

The concentration of Si atoms was determined by secondary-ion mass spectrometry (SIMS). The electrical characteristics of the layers were determined using Hall- effect measurements.

3. Results

The relationship between biaxial strain and hydrostatic strain in the grown epilayers is given in Figure 4. It should be noticed that both undoped and Si-doped GaN epilayers were compressed (s_a < 0), but the stress was higher in the undoped epilayer. Moreover, the undoped GaN films tolerated practically pure biaxial strain, whereas heavily Si-doped GaN films (n> 2´10¹⁹cm^-3) tolerated hydrostatic strain mainly. In the epilayers with moderate Si-doping there was coexistence of biaxial and hydrostatic strain fields. Biaxial strain was reduced as Si-doping increased. Comparison of the measured dependence of e_a/e_c on the biaxial stress s_a with the calculated dependence suggested a type of native defect, such as interstitial Si in the GaN lattice, in these epilayers. The results of SIMS and Hall- effect measurements correlated with this suggestion (see Figure 5).

The application of differential diffractometry allowed a set of structural peculiarities to be revealed in the GaN epilayers.

Specific features of X-ray diffraction in GaN are as follows

1. w_h >> w_a, both for Bragg- and for Laue geometry;
2. w_h⁰⁰⁰²@w_a⁰⁰⁰⁴, i.e. the influence of the size effect on the value of FWHM is minimal;
3. w_h¹¹²⁴ (g₀ < g_h) < w_h⁰⁰⁰², i.e. the contribution of grain misorientation to the value of FWHM is minimal.

An analysis of the defect structure was performed using the relationship between components of the microdistortion tensor, the types of dislocations and their spacing in the GaN layer. For example, vertical screw dislocations with Burgers vector b parallel to the surface normal n (b || n) made a major contribution to the component e_zx . Vertical edge dislocations (b ^ n) and horizontal edge dislocations (b ^ n) ( so-called misfit dislocations), made a contribution to the component e_xx . Misfit dislocations also made a major contribution to the component e_xz. In addition, vertical edge dislocations (b ^ n) made a major contribution to the component e_xx .

The density of dislocations r is connected to the components of the microdistortion tensor via the expression:

r ~ w2 /b2 ~ eij2 / b2

To compare the behavior of similar types of dislocations, the dislocation density was normalized to that of the undoped layer. The dislocation density was found to increase nonlinearly in the range of Si concentration from 10¹⁸ to 8´10¹⁹ cm^-3. The dependence of the normalized density on Si concentration is presented in Figure 6 for different types of dislocations. It should be noticed that the dependence is minimized at the smallest ratio of Si concentration to electron concentration. It corresponds to a Si concentration of 8´10¹⁸ cm^-3.

The density of misfit dislocations parallel to the interface was always 9Ð25 times (sometimes up to two orders of magnitude) larger for the Si doped layers than for the undoped layers (Figure 6, curve 3). The density of micropipes also varied nonlinearly with the Si concentration. Micropipes were scarcely observed at a Si concentration of 8´10¹⁸ Ð 2´10¹⁹ cm^-3, whereas their density abruptly increased to 10⁹ cm^-2 at a Si concentration of 5Ð8 ´10¹⁹ cm^-3. For the undoped epilayers, the density was less than 10⁵ cm^-2.

The doping concentration dependence of the compression stress was also nonlinear and passed a minimum at a Si concentration of 8´10¹⁸ cm^-3 (Figure 6), again at the smallest ratio of Si concentration to electron concentration. It means that the highest portion of the electrically active Si is in the layer. Moreover, it was found that the introduction of Si reduced the size of the coherent scattering areas. The dependence of the domain size on the Si concentration is presented in Figure 7. It should be noted that the dependence is also nonlinear.

The maximum mobility in Figure 5 was obtained for the layers with the minimum compressive stress, with the minimum size of the domains along the c axis, and at a higher density of both misfit dislocations parallel to the interface and other type of dislocations, compared to the undoped layers. The density of vertical screw and edge dislocations was also minimal, while micropipes were almost absent.

A decrease in both the size of the domains in the columnar structure along the c axis and of misfit dislocations parallel to the interface seems to be favorable for forming channels for carrier transport which promoted an increase in mobility in the Si-doped epilayers. The results suggest that the controlling structural factors influencing carrier transport are the size of the domains along the c axis and the density of misfit dislocations parallel to the interface.

4. Conclusion

The investigations showed that Si doping significantly changed relaxation of stress in the structure of GaN epilayerÐAl₂O₃. Si doping reduced biaxial stress in the layer, and varied the size of the domains in the columnar structure along the c axis, the density of micropipes , the density of vertical screw and edge dislocations, and the density of misfit dislocations parallel to the interface.

Both undoped and Si-doped GaN epilayers were compressed, but the stress was higher in the undoped epilayer. Moreover, the undoped GaN films tolerated practically pure biaxial strain, whereas heavily Si-doped GaN films (n> 2´10¹⁹cm^-3) tolerated mainly hydrostatic strain. In the epilayers with moderate Si-doping there was coexistence of biaxial and hydrostatic strain fields. Biaxial stress reduced as Si-doping increased. At the same time an increase in electron mobility from 20 to 200-260 cm² V^-1 s^-1 was observed

A comparison of the measured dependence of e_a/e_c on the biaxial stress s_a with the calculated dependence suggested a type of the native defect, such as interstitial Si, in these epilayers.

Si doping resulted in increasing both diagonal components e_zz , e_xx and off-diagonal components e_zx ,e_xz , and e_xy of the microdistortion tensor. At the same time, the average size of the coherent scattering areas (domain size) t_x and t_z decreased.

A decrease in the size of the domains in the columnar structure along the c axis and of misfit dislocations parallel to the interface seems to be favorable for forming channels for carrier transport.

References

[1]Proceedings of The Second International Conference on Nitride Semiconductors, October 27-31, 1997, Tokushima, Japan

[2] R. N. Kyutt, T. S. Argunova, Nuov. Cim. D 19, 267-275 (1997).

[3] M. Leszczynski, T. Suski, P. Perlin, M. Teisseyre, J. Grzegory, M. Bockowski, Yun, S. Porowski, J. Major, J. Phys. D 28, A149-A153 (1995).

[4] R. N. Kyutt, L. M. Sorokin, T. S. Argunova, S. S. Ruvimov, Phys. Solid State 36, 1473 (1994).

[5] C. Kisielowski, J. Krüger, S. Ruvimov, T. Suski, J. W. Ager, E. Jones, Z. Lilienthal-Weber, M. Rubin, M. D. Bremser, R. F. Davis, Phys. Rev. B 54, 17745 (1996).

Table 1

Figure 1. X-ray measurement of epitaxial GaN layers on double-crystal (DCD)- and triple-crystal (TCD)- diffractometer

Figure 2. Distribution of intensity around (1124) reflection in asymmetric Bragg geometry for g0<gh and g0>gh, k0 and kh Ð vectors of incident and diffracted X-ray waves, wh - double-crystal FWHM for q - scanning mode.

Figure 3. Distribution of intensity around (1124) reflection in asymmetric Bragg geometry for g0<gh and g0>gh, k0 and kh Ð vectors of incident and diffracted X-ray waves, wh - double-crystal FWHM for q - scanning mode.

Figure 4. Dependence of ec/ea on biaxial stress sa. + Ð experiment, solid line Ð calculated. Calculations were performed based on equation from ref. [5]. For pure biaxial strain ec/ea=-2n/(1-n). For combination of biaxial and hydrostatic strain ea=(1-bC)(1+(1-n)sa/E)-1, ec=(1-bC)(1-2nsa/E)-1, where n - Poisson ratio, E Ð Young modulus, C- impurity concentration, b Ð contraction coefficient.

Figure 5. The density of vertical edge (1), screw (2) dislocations and misfit dislocations (3) parallel to the interface in relation to Si doping

Figure 6. The density of vertical edge (1), screw (2) dislocations and misfit dislocations (3) parallel to the interface in relation to Si doping

Figure 7. Variation of the domain size of columnar structure with Si doping

© 1998-2000 The Materials Research Society

M R S Internet Journal of Nitride Semiconductor Research