research papers\(\def\hfill{\hskip 5em}\def\hfil{\hskip 3em}\def\eqno#1{\hfil {#1}}\)

Journal logoJOURNAL OF
APPLIED
CRYSTALLOGRAPHY
ISSN: 1600-5767

X-ray topography of subsurface crystal layers1

CROSSMARK_Color_square_no_text.svg

aInstitute of Metallurgy and Materials Science, 25 Reymonta Street, Krakow, 30-059, Poland, and bSolid State Physics Department, Yuriy Fedkovych Chernivtsi National University, Kotsyubynskiy Street 2, Chernivtsi, 58018, Ukraine
*Correspondence e-mail: ifodchuk@ukr.net

Edited by V. Chamard, Institut Fresnel, France (Received 5 December 2016; accepted 15 May 2017; online 30 May 2017)

New capabilities of the Berg–Barrett topographic method are demonstrated using a skew-asymmetric X-ray diffraction scheme for investigating structural changes near the surface of semiconductor materials. Specifying the X-ray extinction depth, the details of defects and strains are revealed with high resolution. Consequently, analysis of structural distortion of layers near the surface after various types of surface processing becomes more complete.

1. Introduction

The fundamentals of X-ray diffraction topography were formulated by Berg (1931[Berg, W. F. (1931). Naturwissenschaften, 19, 391-396.]) and supplemented in the work of Barrett (Barrett, 1945[Barrett, C. S. (1945). Trans. AIME, 161, 15-64.]; Ramachandran, 1944[Ramachandran, G. N. (1944). Proc. Indian Acad. Sci. 19, 280.]). Since then, a lot of X-ray optical schemes have been applied in various methods of X-ray topography and diffractometry (Authier, 2001[Authier, A. (2001). Dynamical Theory of X-ray Diffraction. New York: Oxford University Press.]; Bedyńska, 1973[Bedyńska, T. (1973). Phys. Status Solidi (A), 19, 365-372.]; Tanner, 1976[Tanner, B. K. (1976). X-ray Diffraction Topography. Oxford: Pergamon Press/Oxford University Press.]; Bowen & Tanner, 1998[Bowen, D. K. & Tanner, B. K. (1998). High Resolution X-ray Diffractometry and Topography. London: Taylor and Francis Ltd.]; Holy et al., 1999[Holy, V., Pietsch, U. & Baumbach, T. (1999). High-Resolution X-ray Scattering from Thin Films and Multilayers. Heilderberg: Springer.]).

The improvement of the known methods of studying the structure of materials and the creation of new methods is still relevant today. This is especially true for X-ray diffractometry and topography. In recent years, these methods have changed significantly in parallel with advances in the development of the dynamical theory of X-ray diffraction by crystals with defects and computer simulation (Härtwig, 2001[Härtwig, J. (2001). J. Phys. D Appl. Phys. 34, A70-A77.]; Holy et al., 1999[Holy, V., Pietsch, U. & Baumbach, T. (1999). High-Resolution X-ray Scattering from Thin Films and Multilayers. Heilderberg: Springer.]). High-resolution sectional and projection double-crystal topography (Authier, 2001[Authier, A. (2001). Dynamical Theory of X-ray Diffraction. New York: Oxford University Press.]; Bowen & Tanner, 1998[Bowen, D. K. & Tanner, B. K. (1998). High Resolution X-ray Diffractometry and Topography. London: Taylor and Francis Ltd.]; Indenbom & Kaganer, 1985[Indenbom, V. L. & Kaganer, V. M. (1985). Phys. Status Solidi (A), 87, 253-265.]; Novikov et al., 1995[Novikov, D., Ohler, M., Köhler, R. & Meterlik, G. (1995). J. Phys. D Appl. Phys., 28, A84-A87.]; Petroff, 1984[Petroff, J.-F. (1984). Applications of X-ray Topographic Methods to Materials Science, edited by S. Weissmann, F. Balibar & J.-F. Petroff. New York: Plenum Press.]) is sufficiently effective in the investigation of various heterostructures and multilayer systems. It presents an opportunity to study the features and mechanisms of formation of diffraction contrast of one- and two-dimensional defects on the surface and in the volume of crystals with higher resolution.

White-beam synchrotron topography (Bowen & Tanner, 1998[Bowen, D. K. & Tanner, B. K. (1998). High Resolution X-ray Diffractometry and Topography. London: Taylor and Francis Ltd.]) is highly effective for the study of defects in materials. A relatively simple implementation of the experimental scheme can easily be adapted for real-time studies of the dynamical effects on materials in in situ control mode. Such a topography scheme, because it involves the analysis of multiple reflections, allows us to obtain maps of complete long-range elastic deformations from various defects. In this case, a double-crystal scheme of topography with the use of several reflections enables the determination of the dilatations and angular misorientations (Riglet et al., 1980[Riglet, P., Sauvage, M., Petroff, J. P. & Epelboin, Y. (1980). Philos. Mag. A, 42, 339-358.]). At the same time, a direct method for plotting misorientation and dilatation maps requires an installation of the crystal analyzer after the sample (Kitano et al., 1991[Kitano, T., Ishikawa, T. & Matsui, J. (1991). Philos. Mag. A, 63, 95-109.]). The sequential movement of the analyzer makes it possible to build complete dilatation maps by the method of multiple exposures.

To characterize the three-dimensional structure of defects in a diamond single crystal, a combination of diffraction topography and microtomography is quite effective (Ludwig et al., 2001[Ludwig, W., Cloetens, P., Härtwig, J., Baruchel, J., Hamelin, B. & Bastie, P. (2001). J. Appl. Cryst. 34, 602-607.]). The use of the rocking curve method is also valid to study rotational and dilatation components of the lattice strain in silicon crystals by analyzing the dependences of the integrated intensity distribution, the half-width of the reflection curve and the peak intensity on the angular deviation from the Bragg angle (Tsoutsouva et al., 2015[Tsoutsouva, M. G., Oliveira, V. A., Baruchel, J., Camel, D., Marie, B. & Lafford, T. A. (2015). J. Appl. Cryst. 48, 645-654.]). In addition, various modifications of the topographic Berg–Barrett method are still quite effectively used today to estimate the quality of the structure of semiconductor materials (Black & Long, 2004[Black, D. R. & Long, G. G. (2004). X-ray Topography, NIST Special Publication 960-10. Washington, DC: US Government Printing Office.]; Lübbert et al., 2000[Lübbert, D., Baumbach, T., Härtwig, J., Boller, E. & Pernot, E. (2000). Nucl. Instrum. Methods Phys. Res. Sect. B, 160, 521-527.]; Lübbert & Baumbach, 2007[Lübbert, D. & Baumbach, T. (2007). J. Appl. Cryst. 40, 595-597.]; Swiatek & Fodchuk, 2016[Swiatek, Z. & Fodchuk, I. M. (2016). Arch. Metall. Mater. 61, 1931-1938.]).

In this paper new possibilities of the X-ray topography method in the skew-asymmetric Bragg diffraction scheme (which are achieved by specifying the X-ray extinction depth) are shown for studying the morphology of the surface of semiconductor materials and structural changes after various types of surface processing.

2. Grazing-incidence skew-asymmetric Bragg topography: experimental setup, geometrical conditions and principal relationships

There are a large number of different X-ray techniques that allow us to study thin subsurface layers with a thickness of about 0.1 µm (Afanas'ev & Melkonyan, 1983[Afanas'ev, A. M. & Melkonyan, M. K. (1983). Acta Cryst. A39, 207-210.]; Afanas'ev et al., 1989[Afanas'ev, A. M., Aleksandrov, P. A. & Imamov, R. M. (1989). X-ray Diffraction Diagnostics of Submicron Layers. Moscow: Nauka.]; Authier, 2001[Authier, A. (2001). Dynamical Theory of X-ray Diffraction. New York: Oxford University Press.]; Bowen & Tanner, 1998[Holy, V., Pietsch, U. & Baumbach, T. (1999). High-Resolution X-ray Scattering from Thin Films and Multilayers. Heilderberg: Springer.]; Imamov et al., 1989[Imamov, R. M., Lomov, A. A. & Novikov, D. V. (1989). Phys. Status Solidi (A), 115, K133-K134.]; Novikov et al., 1995[Novikov, D., Ohler, M., Köhler, R. & Meterlik, G. (1995). J. Phys. D Appl. Phys., 28, A84-A87.]). It is also highly efficient to use asymmetric reflections followed by rotation of the sample around the diffraction vector (Grigoriev et al., 2016[Grigoriev, D., Lazarev, S., Schroth, P., Minkevich, A. A., Köhl, M., Slobodskyy, T., Helfrich, M., Schaadt, D. M., Aschenbrenner, T., Hommel, D. & Baumbach, T. (2016). J. Appl. Cryst. 49, 961-967.]). But no less effective is the diffraction scheme for which the angle ψ between the reflecting plane and the crystal surface is slightly larger than the value of the Bragg angle θ (Fodchuk & Kshevetsky, 1992[Fodchuk, I. M. & Kshevetsky, O. S. (1992). Phys. Met. 14(5), 57-62.]; Fodchuk et al., 1995a[Fodchuk, I. M., Raransky, A. M. & Evdokimenko, A. V. (1995a). Proc. SPIE, 2647, 37-44.], 2009[Fodchuk, I., Zaplitnyy, R., Kazemirskiy, T., Litvinchuk, I. & Swiatek, Z. (2009). Phys. Status Solidi (A), 206, 1804-1808.]; Swiatek & Fodchuk, 2016[Swiatek, Z. & Fodchuk, I. M. (2016). Arch. Metall. Mater. 61, 1931-1938.]):

[\psi - \theta \le {1^\circ}. \eqno(1)]

The position of the crystal when both the normal [\bf n] to the crystal surface and the diffraction vector [\bf h] lie in the same plane of diffraction is typical when the condition of extremely asymmetric Laue diffraction is realized. During rotation of the crystal by φ around [\bf h], the vector [\bf n] circumscribes a cone (Fig. 1[link]). In this way the Bragg condition for the incident beam can be achieved; so, [\bf n] does not lie in the plane of diffraction and various angles between the incident beam and the surface can be realized (Fig. 1[link] and 2[link]). In previous work (Swiatek & Fodchuk, 2016[Swiatek, Z. & Fodchuk, I. M. (2016). Arch. Metall. Mater. 61, 1931-1938.]) it has been noted that this technique of topography opens up new possibilities for layer-by-layer visualization of structural changes in the near-surface layers of a crystal and can also be used in a double-crystal arrangement for determining the strain values in near-surface layers (Fodchuk & Kshevetsky, 1992[Fodchuk, I. M. & Kshevetsky, O. S. (1992). Phys. Met. 14(5), 57-62.]; Fodchuk et al., 1995b[Fodchuk, I. M., Raransky, A. M. & Evdokimenko, A. V. (1995b). Inorg. Mater. 31, 1669-1675.]).

[Figure 1]
Figure 1
Experimental scheme of skew-asymmetric geometry of topography. S is the X-ray source, I1 and I2 are the collimation slits, F is the film, and K is the crystal.
[Figure 2]
Figure 2
(a) Diagram of the transition from the Laue condition to the Bragg condition during rotation of the crystal by φ around [\bf h]: [\varphi \ge \varphi _0'], [{\gamma _0} \ge 0]; [\varphi _0' \le \varphi \le \varphi _{0,k}'] corresponds to the simultaneous total external reflection of the incident beam [{\bf k}_0^{\rm a}] from the surface and diffraction from the (hkl) plane. (b) Skew-symmetric case of X-ray diffraction (this example shows a CdTe crystal, Cu [{K{\alpha _{1,2}}}] radiation): [\varphi = {90}]°, [{\gamma _0} = | {{\gamma _h}} |]; surface is (111), diffracting plane is (115).

The directional cosines for the angles of incidence ([{\gamma _0}]) and reflection ([{\gamma _h}]) are determined as functions of the azimuthal scan angle φ (Fodchuk et al., 1995a[Fodchuk, I. M., Raransky, A. M. & Evdokimenko, A. V. (1995a). Proc. SPIE, 2647, 37-44.]). The asymmetric diffraction in the non-skew condition, i.e. when all vectors [{\bf k}_0], [{\bf k}_h], [\bf n] and [\bf h] lie in the diffraction plane, is a reference point for the angle φ:

[{\gamma _0} = {{{{\bf k}_0} \cdot {\bf n}} \over {| {{{\bf k}_0}} |}}, \quad {\gamma _h} = {{{{\bf k}_h} \cdot {\bf n}} \over {| {{{\bf k}_h}} |}},\quad | {{{\bf k}_0}} | = | {{{\bf k}_h}} | = {1 \over \lambda },\eqno(2)]

where λ is the X-ray wavelength.

In this case the directional cosines for the angles of incidence ([{\gamma _0}]) and reflection ([{\gamma _h}]) are determined as functions of angles θ, ψ and φ:

[{\gamma _0} = \sin {\phi _0} = \sin \theta \cos \psi - \cos \theta \sin \psi \cos \varphi,\eqno(3a)]

[{\gamma _h} = \sin {\phi _h} = - (\sin \theta \cos \psi + \cos \theta \sin \psi \cos \varphi).\eqno(3b)]

When [{\phi _0} \to 0] we can determine the boundary value [{\varphi '_0}] of the angle φ from the ratio (3a)[link] (Fig. 2[link]). This angle corresponds to the transition from the Laue condition to the Bragg condition:

[\cos \varphi _0' = \tan{\theta _{\rm B}}\cot\psi. \eqno(4)]

When [\varphi = 0]° the case of asymmetric Laue diffraction is realized (Fig. 2[link]):

[\eqalignno{{\gamma _0} & = \sin {\phi _0} = \sin \theta \cos \psi - \cos \theta \sin \psi \cr & = \sin (\theta - \psi) = \textstyle\cos ({\pi \over 2} - \psi + \theta), &(5a)}]

[\eqalignno{{\gamma _h} & = \sin {\phi _h} = - \sin \theta \cos \psi - \cos \theta \sin \psi \cr & = - \sin (\theta + \psi) = \textstyle\cos ({\pi \over 2} + \psi + \theta).&(5b)}]

It is easy to determine from relation (2)[link] that the transition from the Laue condition to the Bragg condition ([{\gamma _0} \ge 0]) will occur at angles greater than [\varphi _0'] (Fig. 2[link]):

[\varphi _0' \ge \tan\theta \cot\psi. \eqno(6)]

When [\varphi = 90]°, the case of skew-symmetric Bragg diffraction occurs (Fig. 2[link]) and

[{\gamma _0} = \sin {\phi _0} = \sin \theta \cos \psi. \eqno(7)]

Moreover, at angles ranging from [\varphi _0'] to [\varphi _{0,k}'] the beam reflected by total external reflection can appear in addition to the diffracted beam (Brümmer et al., 1976[Brümmer, O., Höche, H. R. & Nieber, J. (1976). Phys. Status Solidi (A), 33, 587-593.]).

To implement the geometrical condition of skew-asymmetric diffraction we should choose the diffraction planes (hkl) for which the difference between the values of θ and ψ is negligible. The diffraction planes satisfying this condition for the standard X-ray setups and wavelengths of characteristic X-ray radiation are given in Table 1[link].

Table 1
The wavelengths and diffracted planes (hkl) for a skew-asymmetric diffraction geometry for Si, Ge, CdTe and InSb single crystals

No. Crystal Surface hkl λ (Å) θ (°) θψ (°) [{\varphi '_{0,h}}] (°) Lext (µm at φ = 90°) φ (°) (Lext = 0.1 µm)
1 Si (100) 331 1.7889, Co Kα1 45.88 −0.62 11.93 13.11 11.95
(111) 311 1.5405, Cu Kα1 28.06 −1.41 19.55 10.69 19.56
 
2 Ge (111) 220 2.2896, Cr Kα1 34.91 −0.35 8.5 2.34 8.85
(111) 113 1.6578, Ni Kα1 29.07 −0.42 10.08 4.42 10.2
 
3 GaAs (111) 311 1.6578, Ni Kα1 29.10 −0.50 10.31 2.91 10.52
 
4 InSb (111) 511 1.5405, Cu Kα1 38.20 −0.74 13.59 6.21 13.63
 
5 CdTe (111) 511 1.5405, Cu Kα1 38.14 −0.80 13.75 6.74 13.78
 
6 CdHgTe (111) 511 1.5405, Cu Kα1 38.27 −0.67 14.32 5.82 14.36

The range of azimuthal angles φ for which [{\gamma _0} \simeq {\phi _0}] ([\varphi \to {\varphi _{0,k}}]) is within several degrees. This allows us to use a simple technical adjustment to change the incident angles ([{\gamma _0} \simeq {\phi _0}]) gradually and accurately. Of greatest interest for topography are the angle values near the angle of total external reflection, [{\gamma _0} \simeq {\phi _0} \to {\phi _{0,k}}] ([\varphi \to \varphi _{0,k}']), when the incident beam is almost grazing over the surface (Figs. 1[link] and 2[link]). Then, we can achieve the following:

(a) An increase of the extent of the wavefront for the diffracted beam in comparison with the incident beam ([b = | {{\gamma _h}} |/{\gamma _0}]) gives the possibility to obtain topograms from thin near-surface layers of comparatively large area in a single crystal (Swiatek & Fodchuk, 2016[Swiatek, Z. & Fodchuk, I. M. (2016). Arch. Metall. Mater. 61, 1931-1938.]). It is important that the angle of diffraction ([{\gamma _h}]) is rather wide ([{\phi _h}] ≥ 80°). Consequently, we can obtain images virtually without distortion if the photographic material or detector is close to the sample.

(b) A decrease of [{\gamma _0}] leads to a decrease of the penetration depth of X-rays into the near-surface layers of the crystal, which is determined by the extinction length (Hartwig, 1978[Hartwig, J. (1978). Exp. Tech. Phys. 26, 535-546.]; Rustichelli, 1975[Rustichelli, F. (1975). Philos. Mag. 31, 1-12.])

[{L_{\rm ext}} = {\lambda \over {C| {{\chi _{hkl}}} |}} \left ({{\gamma _0}| {{\gamma _h}} |}\right)^{1/2} .\eqno(8)]

Here, λ is the wavelength; [{\chi _{hkl}}] is the Fourier component of the crystal polarizability, calculated according to the algorithm proposed by Stepanov (2004[Stepanov, S. (2004). Proc. SPIE, 5536, 16-26.]); and [C = \cos 2\theta] or 1 for π and σ polarization, respectively. For an unpolarized wave (the case of an X-ray tube) [C = {({1 + \cos 2\theta } )/ 2}].

3. Experimental conditions

The effectiveness of this method depends not only on the geometry and features of X-ray scattering by crystals but also on the experimental conditions: the spectral composition of the radiation, the size and shape of the source S, the distances from the source to the crystal SK and from the crystal to the photographic film (or detector) KF, and the resolution of the photographic material.

In our case, the standard X-ray topography setup was used (Fig. 1[link]), the area of the sharp focus source was S = 250–300 µm2, SK = 0.6 m, KF = 0.01 m, the resolution of the photoemulsion was 300 dashes per millimetre and the emulsion thickness was 25 µm.

To reduce the angular divergence of the beam in the plane of incidence (in the horizontal plane), two systems of collimating slits I1 and I2 (width of 2 mm and height of 50 mm) were used, both in front of the source and in front of the crystal. The total spatial resolution of the image (taking into account all instrumental errors) of the photographic material was 5–7 µm in the direction normal to the crystal surface.

In some cases work was carried out in a double-crystal arrangement in non-dispersive mode in order to obtain and analyze rocking curves. Crystals of Si with asymmetric 220 and InSb with symmetric 333 reflections were used as monochromators.

This geometry of reflection topography makes it possible to selectively visualize structural changes in the near-surface layers of crystalline compounds after different external influences with a sufficiently small step (0.05–0.1 µm) of depth penetration.

4. Results

The potential of X-ray topography to record defects in the crystal lattice is associated with different types of contrast: adsorption, orientation and extinction. These contrast types can be separated from each other in most cases (Authier, 2001[Authier, A. (2001). Dynamical Theory of X-ray Diffraction. New York: Oxford University Press.]; Petroff, 1984[Petroff, J.-F. (1984). Applications of X-ray Topographic Methods to Materials Science, edited by S. Weissmann, F. Balibar & J.-F. Petroff. New York: Plenum Press.]).

The absorption contrast is due to differences between the absorption coefficients in different areas of the sample.

The orientational contrast is determined by the difference between reflected intensities from neighboring local ranges with different orientation of the crystal lattice. To describe the orientational contrast a very good approximation applicable in a very wide range has been proposed by Bonse (1962[Bonse, U. (1962). Direct Observation of Imperfections in Crystals. New York: Willey.]):

[\Delta {I_{hkl}} = {{{\rm d}{I_{hkl}}} \over {{\rm d}\theta }}{\alpha _h}. \eqno(9)]

Here Ihkl is the intensity of the diffracted beam and [{\alpha _h}] is the lattice misorientation. An increase of the sensitivity of topography with orientational contrast to defects can be achieved by improving the collimation of the characteristic X-ray beam or its monochromatization (Tanner, 1976[Tanner, B. K. (1976). X-ray Diffraction Topography. Oxford: Pergamon Press/Oxford University Press.]).

The formation of the extinction contrast is associated with a different degree of attenuation of the incident and reflected beams due to a change in their interaction in the perfect and deformed regions of the crystal. In this case, the region with the strain gradient can be considered as a set of mosaic elements, i.e. as a mosaic crystal in which the dynamic interaction of waves is reduced, and hence the intensity of the beam reflected by this region differs from the intensity reflected by the surrounding (perfect) crystal (Kuriyama & Long, 1984[Kuriyama, M. & Long, G. G. (1984). Applications of X-ray Topographic Methods to Materials Science. New York: Plenum Press.]).

It should be noted that in topographic studies the absorption coefficient of the radiation is of great importance, since it determines the depth of X-ray penetration and the features of image formation on topograms. In our case, soft X-ray radiation is more effective (Table 1[link]).

Figs. 3–5 show examples of the application of the skew-asymmetric scheme of reflection topography for different crystals of CdTe and its compounds.

4.1. Mercury ion etching of variable-gap structure

The method of plasma-chemical etching is based on the chemical reactions between atoms or molecules of materials and ions of chemically active gases or vapor (Savitsky et al., 1998[Savitsky, V., Mansurov, L., Fodchuk, I., Izhnin, I. I., Virt, I., Lozynska, M. & Evdokimenko, A. (1998). Proc. SPIE, 3725, 299-303.]). This method has a high degree of selectivity of etching materials and is not inferior to liquid-chemical etching (Tennant et al., 1992[Tennant, W. E., Cockrum, C. A., Gilpin, J. B., Kinch, M. A., Reine, M. B. & Ruth, R. P. (1992). J. Vac. Sci. Technol. B, 10, 1359-1369.]). After plasma-chemical etching, as well as after ion etching, additional cleaning, washing and drying of the samples are not required.

Fig. 3[link] shows the results of topographic studies of the effect of mercury ion etching from a high-frequency mercury glow discharge on the morphology of the (111) Cd1−xHgxTe/CdTe surface. The following values of the accelerating voltage of Hg ions were used: I – 1200 V; II – 2400 V; III – 600 V.

[Figure 3]
Figure 3
Mercury ion etching of variable-gap Cd1−xHgxTe/CdTe structure (x = 0.19). The accelerating voltage of the Hg ions is (I) E = 1200 V, (II) E = 2400 V and (III) E = 600 V. Crystal surface is (111), reflection is 511, Cu [{K{\alpha_{ 1,2}}}]. (a) Lext = 2.61 µm (φ = 35°); (b) Lext = 2.17 µm (φ = 25°); (c) Lext = 1.11 µm (φ = 18°); (d) Lext = 1.42 µm (φ = 20°). Numbers 1 and 2 indicate cavities of hexagonal shape, and 3 indicates hill-like inclusions of another phase above the surface.

Two systems of equivalent `skew-symmetric' planes (511) and (151) were used for topographic studies, for which [\theta - \psi = - 0.8]° ([{\varphi '_0} = 13.69]°).

According to relation (8)[link], at the transition from a skew-symmetric to an extremely skew-asymmetric scheme, Lext is reduced several times from 6.74 to 1.12 µm. In this case, the Cu [{K\alpha _1}] and Cu [{K\alpha _2}] reflections are broadened, and the values of their maximum intensities are decreased. The surface morphology in irradiated regions is different. Mechanical damage (micro-scratches) is observed in region III; at the same time it is absent in regions I and II.

The difference in the reflectivity of the irradiated regions increases in the topograms with decreasing Lext (Fig. 3[link]b). Cu [{K\alpha _1}] and Cu [K{\alpha _2}] reflections suffer minor displacements and blurring on boundaries between the regions. In region II, where the ion energy is the largest, ion etching has a pronounced selectivity. This is evident in the significant rise of hill-like inclusions of another phase (indicated by number 3) above the etched matrix and the appearance of cavities of hexagonal shape (indicated by numbers 1 and 2) with a light edge. This is especially well illustrated in Fig. 3[link](c) at Lext = 1.1 µm. In the same part of the crystal, the presence of structural micro-misorientation between local regions is indicated by the changing angular distances between the Cu [K{\alpha _1}] and Cu [K{\alpha _2}] reflections and their bending.

The height s of a single hill can be estimated knowing the length of its shadow l on the surface by the ratio

[s = l \tan {\phi _0} \simeq l {\phi _0}. \eqno(10)]

In our case, the height range of the hills changes from 0.1 to 12 µm. The minimum hill size detected by this technique depends on the resolution of the photographic material or CCD camera. If, for example, the shadow length is l = 1 mm, then in the range near to angles of total eternal reflection ([{\phi _{0,k}} \simeq {15'}]) s = 0.29 µm.

Analysis of the rocking curve characteristics (half-widths, maximum values and integral intensity) allows evaluation of the average strain values in the layer, the thickness of which corresponds to the extinction length (Fodchuk et al., 2007[Fodchuk, I., Zaplitnyy, R., Kazemirskiy, T. & Swiatek, Z. (2007). Phys. Status Solidi (A), 204, 2714-2720.]). Starting with Lext ≃ 2.6 µm (Fig. 3[link]), the values of the half-widths W of the rocking curves are significantly increased in regions I and II in comparison with region III. Values of strains were averaged over the layer thickness. The ratio between the strains in the three regions is 1:2.3:1.5 at Lext ≃ 1.4 µm and 1:4:1.8 at Lext = 1.1 µm. In general, this corresponds to the values of the accelerating voltage for the Hg ions in these regions. Differences in the strain magnitudes () become notable at Lext ≤ 0.5 µm: at Lext ≃ 0.3 µm the strain values in the surface layers are ≃ 10−3 in region III and ≃ 5 × 10−4 in region I.

In region II, at layer thickness Lext ≤ 1 µm (Fig. 3[link]) the diffracted radiation is detected as background. This indicates a strong disorder of the structure in this region of the surface layer. Estimates show that the layer thickness is of the order of 0.2 µm.

The contribution of more components of the strain tensor (due to structural and mechanical defects) and strain relaxation near the surface is determinative in the diffraction pattern formation if the X-ray penetration depth is reduced (Lerche et al., 1991[Lerche, V., Dörnfelder, P. & Härtwig, J. (1991). Phys. Status Solidi (A), 128, 269-283.]). This can be used to analyze the set of X-ray topograms from one or equivalent reflections in order to improve the efficiency of X-ray topography in quantifying the inhomogeneities in the near-surface layers of crystals.

4.2. Ion irradiation of the CdTe surface

The next research objects were single crystals of CdTe(111), grown by the Bridgman method, with a dislocation density of ∼5 × 104 cm−2. Ion implantation was carried out by As ions with energy E = 100 keV and dose D = 1 × 1015 cm−2. A part of the sample surface was protected from ion influence by covering it with a mask of chemically resistant lacquer.

Under ion implantation significant disordering and damage of the surface layers are caused by the formation of excess concentrations of point defects in the two sublattices, as well as cluster formation, surface dislocations and stacking faults (Vlasov et al., 2008[Vlasov, A. P., Bonchyk, O. Yu., Kiyak, S. G., Fodchuk, I. M., Zaplitnyy, R. M., Kazemirskiy, T., Barcz, A., Zieba, P. S., Swiatek, Z. & Maziarz, W. (2008). Thin Solid Films, 516, 8106-8111.]). Implanted arsenic ions occupy anion sites of the crystal lattice. This leads to migration of ions towards the surface and, as a result, the stresses in the near-surface layers of the crystal increase.

The topogram shown in Fig. 4(a)[link] was obtained in a conventional asymmetric scheme using 440 reflections. Bending and blurring of the Cu [{K{\alpha _1}}] reflection indicate the presence of defects typical for such compounds in the crystal bulk (the part of region I delineated by a rectangle): structure fragmentation due to microtwinning or small-angle boundaries between crystalline grains, an inhomogeneous distribution of impurities, a deviation from the stoichiometry etc. For these reflections in topograms, in fact, there are no noticeable surface changes caused by implantation of As ions (region I). This is probably because the formation of the image in the topograms occurs from relatively thick layers, on which the influence of surface distortions is not significant.

[Figure 4]
Figure 4
CdTe(111). Region I contains implanted As ions (E = 100 keV, D = 1015 cm–2). Topograms were obtained with Cu [{K{\alpha _{1,2}}}] radiation. (a) Asymmetric 440 reflection, Lext = 3.68 µm; (b) skew-asymmetric 511 reflection, Lext = 2.95 µm (φ = 30°); (c) skew-asymmetric 511 reflection, Lext = 1.65 µm (φ = 20°). Numbers 1 indicates structure fragmentation, 2 cavities, 3 micro-scratches and 4 hill-like inclusions.

At the same time it is not just mechanical damage of the sample surface that is observed in topograms (Fig. 4[link]b and 4[link]c): micro-scratches as well as the boundaries between the implanted and not-implanted regions and their different reflectivity can be seen, too. Individual local regions where the X-ray beam does not diffract are detected in the absorption and partly in the extinction contrast (hill-like inclusions 3 in Fig. 4[link]). The sizes of these regions are from 60 to 300 µm and their density is ∼103 cm−2. In some cases these regions can be crystalline or amorphous inclusions of another phase or pores. At Lext ≃ 1.65 µm (Fig. 4[link]c), in the not-implanted region, blurring and fusion of Cu [{K{\alpha _1}}] and Cu [{K{\alpha _2}}] reflections are observed.

The well known software package SRIM-2008 (https://www.srim.org/#SRIM) was used to specify the model, adequately describing the defect formation due to the implantation of As ions in CdTe. During simulation it was assumed that point defects are generated only by elastic collisions of implanted ions with lattice atoms, when the energy is transferred from the implant to the CdTe nuclear subsystem. Calculations showed that energy losses of As ions on the nuclei of lattice atoms predominate during the implantation. The average free length of the ion implant is ∼0.6 µm. The maximum concentration of implanted As ions corresponds to a depth of about 0.4 µm. Moreover, the maxima of the depth dependences of the matrix ion (Cd and Te) displacements and nuclear energy losses coincide and correspond to a depth of ∼0.2 µm. This was confirmed by the data obtained from the analysis of profiles of As distribution using the SIMS method (Vlasov et al., 2008[Vlasov, A. P., Bonchyk, O. Yu., Kiyak, S. G., Fodchuk, I. M., Zaplitnyy, R. M., Kazemirskiy, T., Barcz, A., Zieba, P. S., Swiatek, Z. & Maziarz, W. (2008). Thin Solid Films, 516, 8106-8111.]).

The strain profile was obtained by the numerical solution of the Takagi equations (Takagi, 1962[Takagi, S. (1962). Acta Cryst. 15, 1311-1312.]), and the maximum strain value was found to be Δd/d ≃ 2.5 × 10−3, which corresponds to a depth of ∼0.15 µm. The range of significant surface distortions extends to a depth of 0.3–0.35 µm, and the total thickness of the distorted layer is up to 0.9 µm

4.3. Implantation of CdHgTe/CdTe by As ions to different doses

Cd1−xHgxTe/CdTe epitaxial layers are inhomogeneous and characterized by a fragmented fine-grained surface structure (Fig. 5[link]). They contain small-angle boundaries and single inclusions (usually Te). There is an insignificant bending of the planes. The average dislocation density is ∼5 × 105 cm−2 (Fodchuk et al., 2007[Fodchuk, I., Zaplitnyy, R., Kazemirskiy, T. & Swiatek, Z. (2007). Phys. Status Solidi (A), 204, 2714-2720.]).

[Figure 5]
Figure 5
Topograms of Cd1−xHgxTe/CdTe (111) epitaxial layers, obtained with the skew-asymmetric 115 reflection of Cu [{K{\alpha_{ 1,2}}}]. Region I is not implanted, region II is implanted with dose D1, and region III is implanted with dose D1 + D2. (a) [{L_{\rm ext}} \simeq 5.82 ] µm (φ = 90°); (b) [{L_{\rm ext}} \simeq 2.61] µm (φ = 30°); (c) [{L_{\rm ext}} \simeq 1.42] µm (φ = 20°).

Two stages of implantation were used to obtain a more homogeneous profile of the implanted layer, which substantially affects the electrophysical properties and structural perfection of Cd1−xHgxTe/CdTe epitaxial layers. Such implantation is carried out by changing the type of implanted ions, their doses and their energies. In our case the first dose of As ions is D1 = 1015 cm−2 (Figs. 5[link]b and 5[link]c, region II), and then part of the crystal is irradiated with an additional dose D2 = 5 × 1014 cm−2 (Figs. 5[link]b and 5[link]c, region III). The energy of the implanted ions is 100 keV.

The reflectivity of the surface layers decreases when the dose of As ion implantation increases in the different parts of the sample owing to significant structural disorder and damage of the surface layers. The thickness of such layers can be determined from a complex analysis of the rocking curve characteristics and topogram features (Fig. 5[link]): ∼0.29 µm for the first stage of implantation (D1) and 0.38 µm for the second stage (D1 + D2). The thickness of the layer that corresponds to structural disorder after both stages of ion implantation is ∼1.5 µm.

As in previous cases, the intensity of the reflected radiation originates from the relatively thin near-surface layers, so here distortions and local strains from defects significantly affect the formation of the kinematic image of defects on topograms.

5. Conclusions

The skew-asymmetric Bragg diffraction scheme gives new possibilities for studying the morphology of the surface of semiconductor materials and structural changes after various surface processing scenarios.

With a decrease of the X-ray penetration depth, the contribution to the diffraction pattern of a larger number of components of the elastic distortion tensor and the relaxation field of deformations from structural and mechanical defects increases. This increases the sensitivity and efficiency of X-ray topography and rocking curves in the quantitative characterization of inhomogeneities in the near-surface layers of crystals.

The use of the technical and logistic capabilities of modern synchrotron topography and diffractometry in a skew-asymmetric scheme would allow in situ selective investigation of structural changes with high resolution and their quantitative values in thin near-surface layers of crystalline materials from a series of topograms and rocking curves under different influences.

Footnotes

1This article will form part of a virtual special issue of the journal, presenting some highlights of the 13th Biennial Conference on High-Resolution X-ray Diffraction and Imaging (XTOP2016).

References

First citationAfanas'ev, A. M., Aleksandrov, P. A. & Imamov, R. M. (1989). X-ray Diffraction Diagnostics of Submicron Layers. Moscow: Nauka.  Google Scholar
First citationAfanas'ev, A. M. & Melkonyan, M. K. (1983). Acta Cryst. A39, 207–210.  CrossRef CAS Web of Science IUCr Journals Google Scholar
First citationAuthier, A. (2001). Dynamical Theory of X-ray Diffraction. New York: Oxford University Press.  Google Scholar
First citationBarrett, C. S. (1945). Trans. AIME, 161, 15–64.  Google Scholar
First citationBedyńska, T. (1973). Phys. Status Solidi (A), 19, 365–372.  Google Scholar
First citationBerg, W. F. (1931). Naturwissenschaften, 19, 391–396.  CrossRef CAS Google Scholar
First citationBlack, D. R. & Long, G. G. (2004). X-ray Topography, NIST Special Publication 960-10. Washington, DC: US Government Printing Office.  Google Scholar
First citationBonse, U. (1962). Direct Observation of Imperfections in Crystals. New York: Willey.  Google Scholar
First citationBowen, D. K. & Tanner, B. K. (1998). High Resolution X-ray Diffractometry and Topography. London: Taylor and Francis Ltd.  Google Scholar
First citationBrümmer, O., Höche, H. R. & Nieber, J. (1976). Phys. Status Solidi (A), 33, 587–593.  Google Scholar
First citationFodchuk, I. M. & Kshevetsky, O. S. (1992). Phys. Met. 14(5), 57–62.  Google Scholar
First citationFodchuk, I. M., Raransky, A. M. & Evdokimenko, A. V. (1995b). Inorg. Mater. 31, 1669–1675.  Google Scholar
First citationFodchuk, I. M., Raransky, A. M. & Evdokimenko, A. V. (1995a). Proc. SPIE, 2647, 37–44.  Google Scholar
First citationFodchuk, I., Zaplitnyy, R., Kazemirskiy, T., Litvinchuk, I. & Swiatek, Z. (2009). Phys. Status Solidi (A), 206, 1804–1808.  CrossRef CAS Google Scholar
First citationFodchuk, I., Zaplitnyy, R., Kazemirskiy, T. & Swiatek, Z. (2007). Phys. Status Solidi (A), 204, 2714–2720.  CrossRef CAS Google Scholar
First citationGrigoriev, D., Lazarev, S., Schroth, P., Minkevich, A. A., Köhl, M., Slobodskyy, T., Helfrich, M., Schaadt, D. M., Aschenbrenner, T., Hommel, D. & Baumbach, T. (2016). J. Appl. Cryst. 49, 961–967.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationHartwig, J. (1978). Exp. Tech. Phys. 26, 535–546.  Google Scholar
First citationHärtwig, J. (2001). J. Phys. D Appl. Phys. 34, A70–A77.  Web of Science CrossRef CAS Google Scholar
First citationHoly, V., Pietsch, U. & Baumbach, T. (1999). High-Resolution X-ray Scattering from Thin Films and Multilayers. Heilderberg: Springer.  Google Scholar
First citationImamov, R. M., Lomov, A. A. & Novikov, D. V. (1989). Phys. Status Solidi (A), 115, K133–K134.  CrossRef Google Scholar
First citationIndenbom, V. L. & Kaganer, V. M. (1985). Phys. Status Solidi (A), 87, 253–265.  CrossRef CAS Google Scholar
First citationKitano, T., Ishikawa, T. & Matsui, J. (1991). Philos. Mag. A, 63, 95–109.  CrossRef CAS Google Scholar
First citationKuriyama, M. & Long, G. G. (1984). Applications of X-ray Topographic Methods to Materials Science. New York: Plenum Press.  Google Scholar
First citationLerche, V., Dörnfelder, P. & Härtwig, J. (1991). Phys. Status Solidi (A), 128, 269–283.  CrossRef CAS Google Scholar
First citationLübbert, D. & Baumbach, T. (2007). J. Appl. Cryst. 40, 595–597.  Web of Science CrossRef IUCr Journals Google Scholar
First citationLübbert, D., Baumbach, T., Härtwig, J., Boller, E. & Pernot, E. (2000). Nucl. Instrum. Methods Phys. Res. Sect. B, 160, 521–527.  Web of Science CrossRef CAS Google Scholar
First citationLudwig, W., Cloetens, P., Härtwig, J., Baruchel, J., Hamelin, B. & Bastie, P. (2001). J. Appl. Cryst. 34, 602–607.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationNovikov, D., Ohler, M., Köhler, R. & Meterlik, G. (1995). J. Phys. D Appl. Phys., 28, A84–A87.  CrossRef CAS Google Scholar
First citationPetroff, J.-F. (1984). Applications of X-ray Topographic Methods to Materials Science, edited by S. Weissmann, F. Balibar & J.-F. Petroff. New York: Plenum Press.  Google Scholar
First citationRamachandran, G. N. (1944). Proc. Indian Acad. Sci. 19, 280.  Google Scholar
First citationRiglet, P., Sauvage, M., Petroff, J. P. & Epelboin, Y. (1980). Philos. Mag. A, 42, 339–358.  CrossRef CAS Google Scholar
First citationRustichelli, F. (1975). Philos. Mag. 31, 1–12.  CrossRef CAS Web of Science Google Scholar
First citationSavitsky, V., Mansurov, L., Fodchuk, I., Izhnin, I. I., Virt, I., Lozynska, M. & Evdokimenko, A. (1998). Proc. SPIE, 3725, 299–303.  CrossRef Google Scholar
First citationStepanov, S. (2004). Proc. SPIE, 5536, 16–26.  CrossRef Google Scholar
First citationSwiatek, Z. & Fodchuk, I. M. (2016). Arch. Metall. Mater. 61, 1931–1938.  Google Scholar
First citationTakagi, S. (1962). Acta Cryst. 15, 1311–1312.  CrossRef CAS IUCr Journals Web of Science Google Scholar
First citationTanner, B. K. (1976). X-ray Diffraction Topography. Oxford: Pergamon Press/Oxford University Press.  Google Scholar
First citationTennant, W. E., Cockrum, C. A., Gilpin, J. B., Kinch, M. A., Reine, M. B. & Ruth, R. P. (1992). J. Vac. Sci. Technol. B, 10, 1359–1369.  CrossRef CAS Google Scholar
First citationTsoutsouva, M. G., Oliveira, V. A., Baruchel, J., Camel, D., Marie, B. & Lafford, T. A. (2015). J. Appl. Cryst. 48, 645–654.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationVlasov, A. P., Bonchyk, O. Yu., Kiyak, S. G., Fodchuk, I. M., Zaplitnyy, R. M., Kazemirskiy, T., Barcz, A., Zieba, P. S., Swiatek, Z. & Maziarz, W. (2008). Thin Solid Films, 516, 8106–8111.  CrossRef CAS Google Scholar

This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.

Journal logoJOURNAL OF
APPLIED
CRYSTALLOGRAPHY
ISSN: 1600-5767
Follow J. Appl. Cryst.
Sign up for e-alerts
Follow J. Appl. Cryst. on Twitter
Follow us on facebook
Sign up for RSS feeds