Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Application of digital phase shifting moiré method in interface and dislocation location recognition and real strain characterization from HRTEM images

Open Access Open Access

Abstract

In high-resolution transmission electron microscopy (HRTEM) images of heterostructures, it is always difficult to accurately determine the interface position and identify dislocations in a large field of view at tens to hundreds of nanometers due to the small lattice differences. However, in the heterostructure, the determination of the interface position is the key to obtain the true mismatch stress/strain field of the interface. Due to the magnifying effect of the digital moiré method on small differences, digital moiré technology was applied to determine Ge/Si heterostructure interfaces and large-area identification interface dislocations in HRTEM lattice diagrams in this study. By optimizing the frequency and angle of the reference lattice, the interface and dislocation position are clearly and intuitively displayed. How to accurately determine the position of the heterostructure interface and the dislocation of the large-area recognition interface from HRTEM images are studied through simulation experiments. The results show that when the frequency of the reference lattice and the specimen lattice are close, and the angle between them is within 10°, the position of the heterostructure interface can be accurately and intuitively determined by the naked eye according to the distortion characteristics of the moiré fringe. When the frequency of the reference lattice is 0.7 to 0.9 times of the specimen lattice, and the rotation angle is within 8°, the visually clear crossover phenomenon of the moiré fringes is used for large-area identification of interface dislocations. Using the phase measurement interface position sensitivity can reach the Å level. Using the phase-shifting digital moiré method the strain field on the dislocation core at the Ge/Si heterostructure interface and the interface stress distribution were quantitatively analyzed. Compared with the Peierls-Nabarro dislocation model and the Foreman dislocation model, Foreman's variable factor α = 4 is more suitable for describing the strain field of misfit dislocations on the Ge/Si heterostructure interface.

© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Semiconductor heterostructure material is a new generation of semiconductor materials based on the molecular beam epitaxy and metal organic chemical vapor deposition [1,2]. Owing to its excellent photoelectric properties, it has been widely used in electronic, optoelectronic and energy conversion devices [35]. However, in the production of the heterostructure, due to the difference in the lattice constants, lattice mismatch is inevitable. The resulting strain energy accumulation would generate misfit dislocations or other defects at the film-substrate interface during the growth process [68]. The performance, yield strength, reliability, and degradation behavior of semiconductor devices would suffer from the presence of defects, especially dislocations, which can adversely affect its electronic performance [9]. At the same time, with semiconductor devices sizes decreasing into the nanometer scale, strain engineering in the Si channel becomes a critical part of improving performance for large scale integrated circuits. The presence of stress and strain in the interface film may also greatly improve the performance of the micro-nano electronic components. Accurate measurement of these complex interface strain and stress fields will provide a scientific basis for the optimal design of microelectronics and optoelectronic devices. The strain at more general heterophase boundaries containing interfacial dislocations has been treated within elasticity theory as a bimaterial of two joined half-spaces, each having the elastic properties of the bulk [10]. Therefore, it is especially important to determine the interface of the two-phase material during the calculation process, but the interface between two adjacent atomic phases is difficult to be accurately observed in high-resolution images. How to accurately measure the misfit dislocation strain field and interface strain/stress distribution at the interface of heterostructures has become a key issue. Moreover, there are many dislocation models for theoretical analysis of dislocations, such as the elastic theory model, the discrete model, the Peierls-Nabarro (P-N) dislocation model [11,12], the improved P-N dislocation model by Foreman [13], and so on. However, it is still not clear that which dislocation model is best suited for describing the misfit dislocation at the Ge/Si heterostructure interface.

The sampling moiré [14] and phase-shift digital moiré method [15] are widely used in deformation testing due to their fast and high precision. For example, S Ri [16] et al. used the sampling moiré method to accurately measure the small deformation distribution of the three-point bend of the steel beam. Some researchers have applied the sampling moiré method to the defect detection of the GaN atomic structure [17]. Masako Kodera [18] et al. developed a two-dimensional fast-Fourier-transform sampling moiré technique to visually and quantitatively determine the locations of minute defects in a transmission electron microscopy image. Some researchers used a phase-shifted moiré method to measure thermal deformation in a four-plane package electronic package under an atomic force microscope [19]. Some researchers applied the digital phase shift method to microscopic scale 3D deformation measurement [20] Optical moiré technique is widely used in related fields since 1948, which is still a powerful and promising technique for researchers. Some researchers have developed many new methods based on the principle of geometric moiré technology [21,22]. Some researchers use HREM crystal material images and unidirectional grating superposition to form a nano-moiré method with nanometer sensitivity and spatial resolution [23]. It can be used to measure important displacements, strains, dislocations, superposition faults and atomic bond ruptures in nano-mechanics. Some researchers used a TEM lattice image to generate a reference lattice and a sample lattice [24]. A digital rotating moiré strain measurement method based on HRTEM lattice image is proposed. The sampling moiré method, the digital phase shift method, and the digital moiré method described above are mostly used for deformation measurements that do not include an interface. Although, GPA is widely used in various deformation tests due to its high precision and strong noise immunity, such as strain field analysis close to dislocations [2527]. However, the GPA method can only measure the strain of a single type of atomic lattice structure. For multilayer structures, the measurement results are almost all relative strain fields [2832]. And during the GPA calculation process, the selection process involves the selection of the size and shape of the diffraction point and filtering methods, which may make the calculating result different by different operator [33]. However, in materials of semiconductor heterostructures, the real strain/stress at the interface is more important for mechanical characterization, especially for structural properties on both sides of the interface. How to accurately determine the interface of the heterostructure and correctly select the reference area of each part is the key problem of the actual interface strain/stress analysis.

In this study, the application of digital moiré has an amplification feature for small deformations. By precisely adjusting the frequency and angle between the reference grid and the sample lattice, the position of the interface in the high-resolution image of the heterostructure is clearly and intuitively displayed, and the position of the defect in the dislocation can be easily recognized on a large field of view. The simulation test results show that when the frequency of the reference lattice and the specimen lattice are close, and the angle between them is within 10°, the position of the heterostructure interface can be accurately and intuitively determined according to the distortion characteristics of the moiré fringes. When the frequency of the reference lattice is 0.7 to 0.9 times of the specimen lattice, and the rotation angle is within 8°, the visually clear crossover phenomenon of the moiré fringes is used for large-area identification of interface dislocations. And combined with the four-step phase-shifting method, the phase-shifting digital moiré method is used to accurately measure the strain field around a mismatching dislocation at the interface of the Ge/Si heterostructure and the interface stress distribution were quantitatively analyzed using elastic theory. The theoretical strain fields calculated by the P-N and Foreman dislocation models were analyzed and compared. The accurate measurement of the internal strain/stress distribution of the semiconductor heterostructure and the evaluation of the overall heterogeneous lattice quality provide an important basis for judging the optoelectronic performance of the whole device and the optimization design of the heterostructure.

2. Method principle

2.1 Digital moiré formation method

Moiré method is a relatively mature method of photomechanics, which uses two groups of the grating with periodic changes to generate moiré fringes. The moiré fringes will move macroscopically with a slight change in the geometry of the grating, where the movement of the moiré fringes is directly used for the displacement and strain fields analysis of the tested object. In the moiré method, the deformed grating is the carrier of the deformation information of the object. Therefore, the appropriate grating can be utilized to quickly and easily extract the deformation information of the object. In this study, the reference lattice is edited by MATLAB software to ensure there is not any unknown deformation in the reference lattice. That is, a high-resolution image of the crystal material (specimen lattice) is superimposed with an unidirectional digital lattice (reference lattice) to form a moiré fringe, which shows the difference between the specimen lattice and the reference lattice, and this difference contains information such as displacement, strain, and atomic defects. Therefore, the digital moiré method effectively reduces the possible experimental error. Simultaneously, the reference lattice can be accurately moved in the computer, and the reference lattice frequency and the angle between the reference lattice and the specimen lattice can be easily adjusted. Moiré fringes with different frequencies and angles can be obtained. Interface position can be quickly determined in heterostructure. And defect locations such as dislocations can be conveniently identified over a large field of view.

In this study, an HRTEM image of Ge/Si heterostructure (11-1) crystal face were used as the specimen lattice. The intensity distribution of the periodic specimen lattice can be expressed as

$${I_s}(x,y) = {N_{s0}} + {N_{s1}}\cos \left( {\frac{{2\pi }}{{{p_s}}}y + {\varphi_s}} \right)$$
where Ns0 is the mean background intensity, Ns1 is the modulation amplitude, and ps is the pitch of the specimen lattice, φs is the phase component that carrying the deformations.

The reference lattice was generated by MATLAB programming. The intensity distribution of the periodic specimen lattice and reference lattice can be expressed as

$${I_r}(x,y) = {N_{r0}} + {N_{r1}}\cos \left( {\frac{{2\pi }}{{{p_r}}}y + {\varphi_r}} \right)$$
where Nr0 and Nr1 are the constants related to the mean background intensity and modulation amplitude, respectively, and pr is the pitch of the reference or specimen lattice. φr is the initial phase.

According to the moiré formation principle, when the reference lattice is superimposed with the specimen lattice, a clear moiré fringe is produced. The intensity distribution function can be calculated by

$$I = {I_r} - {I_s} = {N_{r0}} + {N_{r1}}\sin \left( {\frac{{2\pi }}{{{p_r}}}y + {\varphi_r}} \right) - {N_{s0}} - {N_{s1}}\sin \left( {\frac{{2\pi }}{{{p_s}}}y + {\varphi_s}} \right)$$
In Eq. (3), the information of the reference lattice and specimen lattice is clear and evident, but the Moiré information is implicit. Fortunately, in the traditional Moiré method, there is a classical relationship among the reference spacing (pr), specimen spacing (ps) and Moiré spacing (δ):
$$\frac{1}{\delta } = \frac{1}{{{p_r}}} - \frac{1}{{{p_s}}}$$
Substituting Eq. (4) into Eq. (3), we obtain
$$I \approx {N_0} + {N_1}\cos \left( {\left. {\frac{{2\pi }}{\delta }y + {\varphi_0}} \right)} \right.$$
Due to the intensity distribution of digital moiré fringes is still a trigonometric function distribution, the phase shift algorithm is equally applicable. The phase shift method obtain several fringe patterns with different initial phases by phase-shifting the phase field of a fringe pattern, and the main phase values in the range of [-π, π] can be obtained by combining multiple fringe patterns. This measurement method has the advantages of high resolution, high accuracy, simple algorithm, and fast processing speed. The most extensive four-step phase shift measurement technique with a step length of π/2 is used here. The principle is to calculate the main phase value of the surface topography information of the measured object through four fringe patterns with π/2 phase shift between each other. The phase difference between the four fringe images satisfying the four-step phase shift is 0, π/2, π, and 3π/2, respectively. According to Eq. (5), we can get
$${I_n} \approx {N_0} + {N_1}cos \left( {\left. {\frac{{2\pi }}{\delta }y + {\varphi_0} + \frac{{i\pi }}{2}} \right),({n = 1,2,3,4;i = 0,1,2,3} )} \right.$$
After calculation, the main phase values of the fringe map can be obtained
$$\varphi (x,y) = \arctan \frac{{{I_4}(x,y) - {I_2}(x,y)}}{{{I_1}(x,y) - {I_3}(x,y)}}$$
The strain components can be obtain from the displacement gradients
$$\left\{ \begin{array}{l} {u_x} = {u_r}\cos {\alpha_1} + {u_n}\cos {\alpha_2}\\ {u_y} = {u_r}\sin {\alpha_1} + {u_n}\sin{\alpha_2} \end{array} \right.$$
$$\left\{ \begin{array}{l} {\varepsilon_x} = \partial {u_x}/\partial x\\ {\varepsilon_y} = \partial {u_y}/\partial y \end{array} \right.$$
where ux and uy represent displacements in the x and y directions, respectively; ur and un represent two different crystal orientations [111] and [111] of the lattice structure, respectively. And αi (i = 1, 2) is the angle between the normal direction of the grating and the X-axis of the rectangular coordinate.

The stress is obtained by the following equation:

$${\sigma _x} = \frac{E}{{1 - {v^2}}}({\varepsilon _x} + v{\varepsilon _y})$$
where E is the Young's modulus; v is the Poisson's ratio, ${\varepsilon _x}$ and ${\varepsilon _y}$ are strains in the x and y directions, respectively.

2.2 Method for judging the position of the interface position

The two semiconductors that normally form a heterojunction have the same crystal structure and close atomic spacing. At the same time, any one of the crystal directions in the HRTEM image can be regarded as a periodic orthogonal lattice. The direction and spacing of the moiré fringe are related to the frequency of the reference lattice and the included angle (θ) between reference lattice and specimen lattice. Due to the different lattice spacing at the ends of the heterostructure interface, there are two different fringes at the ends of the interface. When the reference lattice is rotated in the plane, the moiré spacing (δ), the angle (β) between the main direction of the moiré and the horizontal direction can be expressed as

$$\left\{ \begin{array}{l} {\delta_1} = \frac{{{p_r}{p_{s1}}}}{{\sqrt {p_{s1}^2 + p_r^2 - 2{p_r}{p_{s1}}\cos {\theta_1}} }}\\ \tan {\beta_1} = \frac{{{p_r}\sin {\theta_1}}}{{{p_{s1}} - {p_r}\cos {\theta_1}}} \end{array} \right.$$
$$\left\{ \begin{array}{l} {\delta_2} = \frac{{{p_r}{p_{s2}}}}{{\sqrt {p_{s2}^2 + p_r^2 - 2{p_r}{p_{s2}}\cos {\theta_2}} }}\\ \tan {\beta_2} = \frac{{{p_r}\sin {\theta_2}}}{{{p_{s2}} - {p_r}\cos {\theta_2}}} \end{array} \right.$$
where θ1 and θ2 are the angles between the reference lattice and the specimen lattice at both ends of the interface, pr is the space of the reference lattice, ps1 and ps2 are the spacing of the specimen lattice at both ends of the interface, respectively.

The angle ($\Delta \beta $) between the moiré fringes at both ends of the interface can be obtained from Eqs. (11) and (12).

$$\Delta \beta = {\beta _2} - {\beta _1} = {\mathop{\textrm {arc}}\nolimits} \tan \frac{{{p_r}\sin {\theta _2}}}{{{p_{s2}} - {p_r}\cos {\theta _2}}} - {\mathop{\textrm {arc}}\nolimits} \tan \frac{{{p_r}\sin {\theta _1}}}{{{p_{s1}} - {p_r}\cos {\theta _1}}}$$
During the simulation, the angle between the two sets of reference lattice and the specimen lattice is the same. So,
$${\theta _2} = {\theta _1} = \theta $$
Substituting Eq. (14) into Eq. (13), we obtain
$$\Delta \beta = {\beta _2} - {\beta _1} = {\mathop{\textrm {arc}}\nolimits} \tan \frac{{{p_{s1}}\sin \theta + {p_{s2}}\sin \theta }}{{\frac{{{p_{s1}}{p_{s2}}}}{{{p_r}}} + {p_r} - {p_{s1}}\cos \theta - {p_{s2}}\cos \theta }}$$
From the Eq. (15), we know that when θ is constant, pr is close to ps1 or ps2, $\Delta \beta $ is the largest, indicating that the fringe is bent at the interface. When prps1 is very close, substituting it into Eq. (15), we obtain
$$\Delta \beta = {\beta _2} - {\beta _1} = {\mathop{\textrm {arc}}\nolimits} \tan \frac{{\sin \theta }}{{1 - \cos \theta }}$$
It is analyzed from the Eq. (16) that when pr is constant, the closer θ is to 0°, the more $\Delta \beta $ is the most largest, indicating that the fringe bending phenomenon at the interface is more obvious.

The rotation of the moiré fringes was generated by simulation experiments, which verified the derivation of the above theoretical equation. Figure 1(a) is a simulated orthogonal lattice (as a specimen lattice) image with different frequency (the frequencies are 0.136 and 0.130 respectively) using MATLAB software. In order to observe the variation of the moiré fringe more intuitively, series moiré patterns are produced, as shown in Figs. 1(b)–1(d). The red dashed line is the contact interface of two orthogonal lattice. The frequency of the reference lattice changes from 0.115 to 0.135, with an interval of 0.01. The angle between the reference lattice and the sample lattice is constant at 5° (θ12). As the frequency of the reference lattice approaches the specimen lattice (f1 = 0.136, f2 = 0.130), the moiré fringe is bent at the interface (see the red dotted line). Therefore, the simulation test results are consistent with the inference of Eq. (15). As shown in Figs. 1(f)–2(h), the angle between the reference lattice and the specimen lattice varies from 1° to 15°. The frequency of the reference lattice is constant at 0.13. As the angle decreases, the moiré fringe width becomes wider and wider, and the bending angle at the interface becomes larger and larger. The simulation results verify the results of Eq. (16). The simulation test illustrates that when the frequency of the reference lattice and the specimen lattice are close, the bending of the moiré fringe at the interface can be observed with the naked eye at an angle of 10°, which is used to judge the position of the interface.

 figure: Fig. 1.

Fig. 1. Rotating moiré, (a) simulated lattice structure with different frequency; (b) θ = 5°, fr = 0.115; (c) θ = 5°, fr = 0.125; (d) θ = 5°, fr = 0.135. (e) θ = 1°, fr = 0.130; (f) θ = 5°, fr = 0.130; (g) θ = 10°, fr = 0.130 ; (h) θ = 15°, fr = 0.130.

Download Full Size | PDF

Figure 2(a) is a four-step phase shift fringe obtained by Gaussian filtering, and the phase shift amounts are 0,π/2, π, and 3π/2, respectively. Figure 2(b) shows the phase field after wrapped using a four-step phase shifting. Figure 2(c) shows the gradient line obtained by the phase field along the stripe direction. If the fringes are not changed, the gradient lines obtained by the phase field are continuous line segments. When the moiré fringe changes in bending, the continuity of the gradient line at the bend point changes. Therefore, the interface position can be accurately determined according to the continuity of its phase field gradient line. This method improves the accuracy and sensitivity of measurement and realizes automatic measurement. The accuracy of the phase is one percent of the displacement, so the theoretical accuracy of the phase shift method used to determine the position of the interface can reach the Å level.

 figure: Fig. 2.

Fig. 2. (a) Four-step phase-shifted moiré fringe patterns extracted from Fig. 1(d); (b) wrapped phase field; (c) gradient line obtained from the phase field.

Download Full Size | PDF

2.3 Method for identifying dislocations over a large area

High-resolution transmission electron microscopy (HRTEM) is commonly used to observe dislocation information at the heterostructure interfaces. However, it is difficult to quickly identify dislocations in large areas in HRTEM images, since the lattice anomalies caused by dislocations are only within a few atomic ranges in thousands of atoms. In this paper study, the digital moiré method is applied to quickly identify dislocations in a large area in the HRTEM image of heterostructure. We used software simulation to generate three common dislocations in HRTEM images, as indicated by the arrows in the white boxes in Fig. 3. It is found through experiments that when there are dislocation defects in the crystal, the misalignment phenomenon of the moiré fringes caused by the sliding of atoms in different rows (as shown in the white box of Fig. 3(b)). At the same time, the intersection of the moiré fringes is always at the dislocation position and does not change with the angle, as shown in the white box of Fig. 3(c). So, there is an intersection after the superposition and the intersection location does not change with the rotation of the reference lattice, we can determine the existence of the dislocation defect. It is found through experiments that when the frequency of the reference lattice is 0.8 ∼ 0.9 times of the specimen lattice, and the angle is within 10°, the formed moiré fringes are suitable for large-area recognition of dislocations in HRTEM images. That is, a high-resolution image of the crystal material (specimen lattice) is superimposed with an unidirectional digital lattice (reference lattice) to form a moiré fringe, which shows the difference between the specimen lattice and the reference lattice, and this difference contains information such as displacement, strain, and atomic defects.

 figure: Fig. 3.

Fig. 3. (a) Three white boxes A, B, and C are the dislocation locations generated by the software, (b) and (c) different angles of rotational moiré fringes formed by the reference lattice and the specimen lattice, the white square is the location of the moiré fringe intersection.

Download Full Size | PDF

3. Experiment and discussion

An ultra-high vacuum chemical vapor deposition system equipped with pyrolytic BN effusion cells was used to grow Ge films on Si (001) substrate. The TEM specimens were prepared for cross-sectional imaging along the [110] direction using a standard technique, which involved mechanical grinding followed by ion milling. An HRTEM experiment was performed on a JEM-2010 transmission electron microscope at 200 kV.

3.1 Application of digital moiré method in the determining Ge/Si heterostructure interface position

Figure 4(a) is an HRTEM image of a Si/Ge heterostructure in which the yellow region is the interface between Si and Ge, and the two sides are single crystal Ge and single crystal Si regions, respectively. In the HRTEM image, the position of the interface can only be roughly judged. To accurately determine the position of the interface, it is necessary to continue to develop new methods. In the paper, since the lattice constants of Si (a = 0.5428 nm) and Ge (a = 0.5658 nm) are different in the Ge/Si heterostructure, the (11-1) interplanar crystal spacing is also different. Through simulation experiments, we found that the moiré fringes formed by superimposing with the reference lattice when there are two frequencies in the specimen lattice have bending phenomenon, and the line connecting the bending points is the interface between the two frequencies in the specimen lattice. In this paper, moiré fringe images of various angles were generated, in which the best result was obtained when parallel moiré is formed between the reference lattice (fr = 0.130) and the specimen lattice (the (11-1) crystal face of the Ge/Si heterostructure). There is an obvious bending point in the moiré fringes. Therefore, the bending points of moiré fringe are sequentially connected (see the red dashed line in Fig. 4(c)), and the interface position of Ge/Si heterostructure can be determined. The red dotted line in Fig. 4(d) is the interface position determined in the HRTEM image using the above method. It is stated here that the interface determined by this method is not an interface in the strict physical sense, but is the middle position of the transition area between the two substances.

 figure: Fig. 4.

Fig. 4. (a) HRTEM image of Si/Ge heterostructure, the yellow area is the interface between Si and Ge; (b) reference lattice parallel to the Si/Ge heterostructure (11-1) crystal face; (c) interface on moiré fringe; (d) real interface on Si/Ge heterojunction HRTEM image.

Download Full Size | PDF

3.2 Digital moiré method to quickly identify misfit dislocations on the heterostructure interface

In diamond crystal structures, such as Si, Ge and Si1-xGex, dislocations mainly propagate in the {111} crystal face, so misfit dislocations are mainly located in the {111} crystal face near the interface [34]. Figure 5(a) shows the HRTEM image of Ge/Si heterostructure. Figure 5(b) shows a rotation moiré image formed by superimposing the (11-1) crystal face of the HRTEM image with the reference lattice (p = 9 pixels) after a 9° (the angle between the reference lattice and the specimen lattice) counterclockwise rotation. The intersection of fringes can be clearly seen in the four white rectangular boxes region of A, B, C, and D. This phenomenon is very similar to the presence of misfit dislocations in the Ge/Si heterostructure, so it is determined that there may be misfit dislocations. Four white squares A, B, C, D with the same size and position are also marked in the HRTEM image (see Fig. 5(a)). Four misfit dislocations on the (11-1) crystal face of the Ge/Si heterostructure was observed in the squares, as indicated by the white arrow in the Fig. 5(a). There is only one dislocation in each square of this direction, which is corresponding to the change of moiré fringes in Fig. 5(b). Therefore, it can be determined that the intersection of moiré fringes in Fig. 5(b) is caused by misfit dislocations on the interface, and the position information of misfit dislocations in the Ge/Si heterostructure can be quickly and accurately captured by the digital moiré method.

 figure: Fig. 5.

Fig. 5. (a) HRTEM images of Ge/Si heterostructure; (b) rotation moiré formed by the reference lattice and the specimen lattice.

Download Full Size | PDF

3.3 Application of sampling moiré method in determination of burgers vector

Figure 6(a) is the interface region intercepted in the HRTEM image of Si/Ge heterostructure. To facilitate the calculation, the interface is rotated to the horizontal direction where the red dashed line is the interface position of the Ge / Si heterostructure. The upper region is Ge, and the lower is Si. Figure 6(b) shows the Fourier transform of Fig. 6(a), in which the analyzed spots have been encircled. Usually, the type of these misfit dislocations can be directly determined from the HRTEM images by drawing a Burgers circuit around the dislocation. It is observed from Fig. 6(a) that the lattice atoms in the HRTEM image are not obvious, which make it difficult to determine the type of dislocation for drawing the Burgers circuit. For this reason, the sampling moiré method, with the advantages of large field of view and a significant amplification of the distortion, was used to identify the Burger circuit, where the pixels were acquired by sampling in the horizontal and vertical directions with a constant pixel interval in the HRTEM image respectively. The sampling interval approximates the interplanar crystal spacing of Ge/Si heterostructure in the horizontal (13 pixels) and vertical (11 pixels) directions, respectively. Figure 6(c) shows the sampling moiré image of Fig. 6(a), in which the misfit dislocation core in the Ge/Si interface can be seen clearly. The Burgers vector is b = 1/2 [110].

 figure: Fig. 6.

Fig. 6. (a) HRTEM image of Ge/Si heterostructure interface; (b) Fourier transform of (a); (d) burgers vector with sampling moiré.

Download Full Size | PDF

3.4 Strain field of dislocation core at Ge/Si interface

Then, the digital moiré method is used to characterize the strain field around a misfit dislocation. Figure 7(a) calculated region is shown in the white rectangle in Fig. 6(a). It is very convenient to realize the phase shift by moving the reference lattice phase. Taking the parallel moiré pattern as an example, the reference lattice is superposed with the specimen lattice (the (11-1) crystal face of the Ge/Si heterostructure). When the four reference lattices are superposed with the same specimen lattice, the four-step phase-shifted moiré pattern can be obtained, and the phase shift amounts are 0, π/2, π, and 3π/2, respectively, as shown in Fig. 7(b). Figure 7(c) is a moiré pattern extracted by filtering from Fig. 7(b).

 figure: Fig. 7.

Fig. 7. (a) HRTEM image of Ge/Si heterostructure interface; (b) digital moiré fringe patterns; (c) fringe patterns extracted from (b).

Download Full Size | PDF

Figure 8 shows the strain maps around the misfit dislocations at the Ge/Si heterostructure interface. The strain maps of the misfit dislocations are all eight-shaped, with the strain value ranging from -8% to + 8%. The experimental ɛxx strain field calculated using the four-step phase shift method is shown in Fig. 8(a). Figure 8(b) shows the strain component ɛxx calculated using the P-N dislocation model [25]. Figures 8(c)–8(k) demonstrate the strain components ɛxx given by the Foreman model [10], with the Poisson's ratio νSi = 0.278, Burgers vector bSi = 0.384 nm, Young's modulus is 130 GPa; Poisson's ratio νGe = 0.32, bGe = 0.4 nm, Young's modulus is 79.9 GPa, and the corresponding alterable factors10 ranging from α = 2 to α = 10. The eight-shaped strain maps become wider and the dislocation width increase with the increase of the factor α in the Foreman model. There exists a convergence region around the dislocation core. The maximum strain value appears in the core area of the dislocation and the strain value become smaller away from the dislocation core. The strain map of the Foreman model (α = 4) is qualitatively similar to the experimentally determined strain in Fig. 8. Thus, the Foreman model with α = 4 best describes the strain fields of the misfit dislocation core in the Ge/Si heterostructure interface.

 figure: Fig. 8.

Fig. 8. Experimental and theoretical strain fields: (a) experimental strain field; (b) strain field of P-N dislocation model; (c) - (k) strain field of Foreman dislocation model, with α ranging from 2 to 10, respectively.

Download Full Size | PDF

3.5 The strain and stress distribution near the interface of dislocation

At the center of dislocations, atoms are arranged in a disorderly manner and beyond the range of elastic deformation. Hooke's law cannot be applied in the dislocation core, but away from the center area the Hooke’s law can be still used for elastic deformation analysis. In this paper, Ge/Si heterostructure are separated at the interface and treated as two materials with two connected half-spaces, and each half-space is considered to be within the elastic deformation range so that the interface stress can be calculated through elastic deformation.

Figure 9(a) is an interface area extracted away from the dislocation core. The red dashed line is the center of the Ge/Si heterostructure interface. The corresponding parallel moiré fringes are obtained by superimposing specimen lattice (the (11-1) crystal face of the Ge/Si heterostructure) and the reference lattice. Taking the phase far away from the interface as the reference zone, the corresponding displacement ur and un can be obtained by the phase shifting method. Then, the displacement in the horizontal and vertical directions can be obtained using ux= urx + unx and uy= ury + uny, respectively. The strain fields of the Ge/Si heterostructure in the horizontal is shown in Fig. 9(b). There exists a convergence region of the strain field around the interface. The maximum strain value appears in the interface area and the strain value become smaller away from the dislocation core. The upper Ge region is compressively strained and lower Si region is tensile strained. This is due to the fact that Ge (a = 0.5658 nm) has a larger lattice constant than Si (a = 0.5428 nm). When Ge and Si are combined at the interface, due to the difference in lattice constant, the upper Ge is compressed by the Si with a smaller lattice constant. Conversely, the lower Si is stretched. Then, stress is calculated suing the elastic theory.

 figure: Fig. 9.

Fig. 9. Interfacial strain and stress of Ge/Si heterostructure: (a) intercepted Ge/Si heterostructure interface; (b) strain field in the horizontal direction; (c) stress profile corresponding to the dash-dotted line in (b).

Download Full Size | PDF

Figure 9(c) shows the stress profile of the Ge/Si heterostructure corresponding to the dash-dotted line in Fig. 9(b). The stress approaches 0 at a distance away from the interface in the Ge region, and the compressive stress become larger when getting closer to the interface, with a maximum value of 3 GPa. When passing through the interface to reach the Si region, the compressive stress becomes tensile stress and reaches a maximum of 6.6 GPa near the interface.

4. Conclusions

In this study, digital moiré techniques were used to accurately determine the interface of heterostructures and large-area recognition dislocation defects in HRTEM lattice images. The experimental results show that the digital moiré method magnifies the small deformation in the HRTEM image, and can visually observe the existence of interface misfit dislocations and mismatched strains in the moiré pattern. When the frequency of the reference lattice and the specimen lattice are close, and the angle between them is within 10°, the position of the heterostructure interface can be accurately and intuitively determined according to the distortion characteristics of the moiré fringes. When the frequency of the reference lattice is 0.7 to 0.9 times of the specimen lattice, and the rotation angle is within 8°, the visually clear crossover phenomenon of the moiré fringes is used for large-area identification of interface dislocations. The digital phase shifting method is used to characterize the strain field distribution near the dislocation. Compared with the P-N dislocation model and the Foreman dislocation model, Foreman's variable factor α = 4 is more suitable for describing the strain field of misfit dislocations on the Ge/Si heterostructure interface. Simultaneously, quantitative analysis is made on the strain and stress distribution of the interface away from the dislocation, which provides a scientific basis for the optimal design of microelectronics and optoelectronic devices.

Funding

National Natural Science Foundation of China (11232008, 11372037, 11572041).

Acknowledgments

The authors are grateful to all of those who were involved in this work

Disclosures

The authors declare no conflicts of interest.

References

1. S. V. Elshocht, M. Caymax, T. Conard, S. D. Gendt, I. Hoflijk, M. Houssa, F. Leys, R. Bonzom, B. D. Jaeger, J. V. Steenbergen, W. Vandervorst, M. Heyns, and M. Meuris, “Study of CVD high-k gate oxides on high-mobility Ge and Ge/Si substrates,” Thin Solid Films 508(1-2), 1–5 (2006). [CrossRef]  

2. L. Vescan and S. Wickenhauser, “Relaxation mechanism of low temperature SiGe/Si (0 0 1) buffer layers,” Solid-State Electron. 48(8), 1279–1284 (2004). [CrossRef]  

3. K. Brunner, “Si/Ge nanostructures,” Rep. Prog. Phys. 65(1), 27–72 (2002). [CrossRef]  

4. J. Xiang, W. Lu, Y. J. Hu, Y. Wu, H. Yan, and C. M. Lieber, “Ge/Si nanowire heterostructures as Highperformance field-effect transistors,” Nature 441(7092), 489–493 (2006). [CrossRef]  

5. Y. Hu, H. O. Churchill, D. J. Reilly, J. Xiang, C. M. Lieber, and C. M. Marcus, “A Ge/Si heterostructure nanowire-based double quantum dot with integrated charge sensor,” Nat. Nanotechnol. 2(10), 622–625 (2007). [CrossRef]  

6. P. M. J. Maree, J. C. Barbour, d. V. G. F. Van, and K. L. Kavanagh, “Generation of misfit dislocations in semiconductors,” J. Appl. Phys. 62(11), 4413–4420 (1987). [CrossRef]  

7. N. Hirashita, N. Sugiyama, E. Toyoda, and S. Takagi, “Relaxation processes in strained Si layers on silicon-germanium- on-insulator substrates,” Appl. Phys. Lett. 86(22), 221923 (2005). [CrossRef]  

8. R. Loo, R. Delhougne, M. Caymax, and M. Ries, “Formation of misfit dislocations at the thin strained Si∕strain-relaxed buffer interface,” Appl. Phys. Lett. 87(18), 182108 (2005). [CrossRef]  

9. S. Mahajan, “Defects in Semiconductors and their Effects on Devices,” Acta Mater. 48(1), 137–149 (2000). [CrossRef]  

10. R. Bonnet and M. Loubradou, “Atomic positions around misfit dislocations on a planar heterointerface,” Phys. Rev. B 49(20), 14397–14402 (1994). [CrossRef]  

11. R. E. Peierls, “The size of a dislocation,” Proc. Phys. Soc. 52(1), 34–37 (1940). [CrossRef]  

12. F. R. N. Nabarro, “Dislocations in a simple cubic lattice,” Proc. Phys. Soc. 59(2), 256–272 (1947). [CrossRef]  

13. A. J. Foreman, M. A. Jaswon, and J. K. Wood, “Factors controlling dislocation widths,” Proc. Phys. Soc. A 64(2), 156–163 (1951). [CrossRef]  

14. S. Ri and T. Muramatsu, “Theoretical error analysis of the sampling moiré method and phase compensation methodology for single-shot phase analysis,” Appl. Opt. 51(16), 3214–3223 (2012). [CrossRef]  

15. A. Asundi and K. H. Yung, “Phase-shifting and logical moiré,” J. Opt. Soc. Am. A 8(10), 1591–1600 (1991). [CrossRef]  

16. S. Ri, M. Fujigaki, and Y. Morimoto, “Sampling Moiré Method for Accurate Small Deformation Distribution Measurement,” Exp. Mech. 50(4), 501–508 (2010). [CrossRef]  

17. Q. H. Wang, S. Ri, H. Tsuda, M. Kodera, K. Suguro, and N. Miyashita, “Visualization and automatic detection of defect distribution in GaN atomic structure from sampling Moiré phase,” Nanotechnology 28(45), 455704 (2017). [CrossRef]  

18. M. Kodera, Q. Wang, S. Ri, H. Tsuda, A. Yoshioka, and T. Sugiyama, “Characterization technique for detection of atom-size crystalline defects and strains using two-dimensional fast-fourier-transform sampling moiré method,” Jpn. J. Appl. Phys. 57(4S), 04FC04 (2018). [CrossRef]  

19. H. Xie, A. Asundi, G. B. Chai, Y. Lu, Z. W. Zhong, and B. K. A. Ngoi, “High resolution AFM scanning Moiré method and its application to the micro-deformation in the BGA electronic package,” Microelectron. Reliab. 42(8), 1219–1227 (2002). [CrossRef]  

20. D. Wu, H. Xie, C. Li, and R. Wang, “Application of the digital phase-shifting method in 3D deformation measurement at micro-scale by SEM,” Meas. Sci. Technol. 25(12), 125002 (2014). [CrossRef]  

21. G. A. Geach and R. Phillip, “Moiré Patterns in Transmission Electron Micrographs of Sub-Boundaries of Aluminium,” Nature 179(4573), 1293 (1957). [CrossRef]  

22. G. A. Bassett, J. W. Menter, and D. W. Pashley, “Moiré Patterns on Electron Micrographs, and their Application to the Study of Dislocations in Metals,” Proc. R. Soc. Lond. A 246(1246), 345–368 (1958). [CrossRef]  

23. F. Dai and Y. Xing, “Nano-moire method,” Acta Mech. Sin. 15(3), 283–288 (1999). [CrossRef]  

24. H. Xing, Z. Gao, H. Wang, Z. Lei, L. Ma, and W. Qiu, “Digital rotation moirémethod for strain measurement based on high-resolution transmission electron microscope lattice image,” Opt. Laser. Eng. 122(2), 347–353 (2019). [CrossRef]  

25. J. Li, C. Zhao, Y. Xing, S. Su, and B. Cheng, “Full-Field Strain Mapping at a Ge/Si Heterostructure Interface,” Materials 6(6), 2130–2142 (2013). [CrossRef]  

26. M. Mao, A. Nie, J. Liu, H. Wang, S. X. Mao, Q. Wang, K. Li, and X.-X. Zhang, “Atomic resolution observation of conversion-type anode RuO2 during the first electrochemical lithiation,” Nanotechnology 26(12), 125404 (2015). [CrossRef]  

27. A. Nie, L.-Y. Gan, Y. Cheng, H. Asayesh-Ardakani, Q. Li, C. Dong, R. Tao, F. Mashayek, H.-T. Wang, and U. Schwingenschlögl, “Atomic-scale observation of lithiation reaction front in nanoscale SnO2 materials,” ACS Nano 7(7), 6203–6211 (2013). [CrossRef]  

28. Q. Liu, C. W. Zhao, Y. M. Xing, S. J. Su, and B. W. Cheng, “Quantitative strain analysis of misfit dislocations in a Ge/Si heterostructure interface by geometric phase analysis,” Opt. Laser. Eng. 50(5), 796–799 (2012). [CrossRef]  

29. N. Cherkashin, T. Denneulin, and M. J. Hÿtch, “Electron microscopy by specimen design: application to strain measurements,” Sci. Rep. 7(1), 12394 (2017). [CrossRef]  

30. F. Hüe, M. Hÿtch, H. Bender, F. Houdellier, and A. Claverie, “Direct mapping of strain in a strained silicon transistor by high-resolution electron microscopy,” Phys. Rev. Lett. 100(15), 156602 (2008). [CrossRef]  

31. N. Cherkashin, M. J. Hÿtch, E. Snoeck, A. Claverie, J. M. Hartmann, and Y. Bogumilowicz, “Quantitative strain and stress measurements in Ge/Si dual channels grown on a Si0.5Ge0.5 virtual substrate,” Mater. Sci. Eng., B 124-125, 118–122 (2005). [CrossRef]  

32. H. Zhang, Z. Liu, H. Wen, H. Xie, and C. Liu, “Subset geometric phase analysis method for deformation evaluation of HRTEM images,” Ultramicroscopy 171, 34–42 (2016). [CrossRef]  

33. J. J. Peters, R. Beanland, M. Alexe, J. W. Cockburn, D. G. Revin, S. Y. Zhang, and A. M. Sanchez, “Artefacts in geometric phase analysis of compound materials,” Ultramicroscopy 157, 91–97 (2015). [CrossRef]  

34. J. P. Hirth, J. Lothe, and T. Mura, Theory of Dislocations, (2 edition. Wiley, 1982).

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (9)

Fig. 1.
Fig. 1. Rotating moiré, (a) simulated lattice structure with different frequency; (b) θ = 5°, fr = 0.115; (c) θ = 5°, fr = 0.125; (d) θ = 5°, fr = 0.135. (e) θ = 1°, fr = 0.130; (f) θ = 5°, fr = 0.130; (g) θ = 10°, fr = 0.130 ; (h) θ = 15°, fr = 0.130.
Fig. 2.
Fig. 2. (a) Four-step phase-shifted moiré fringe patterns extracted from Fig. 1(d); (b) wrapped phase field; (c) gradient line obtained from the phase field.
Fig. 3.
Fig. 3. (a) Three white boxes A, B, and C are the dislocation locations generated by the software, (b) and (c) different angles of rotational moiré fringes formed by the reference lattice and the specimen lattice, the white square is the location of the moiré fringe intersection.
Fig. 4.
Fig. 4. (a) HRTEM image of Si/Ge heterostructure, the yellow area is the interface between Si and Ge; (b) reference lattice parallel to the Si/Ge heterostructure (11-1) crystal face; (c) interface on moiré fringe; (d) real interface on Si/Ge heterojunction HRTEM image.
Fig. 5.
Fig. 5. (a) HRTEM images of Ge/Si heterostructure; (b) rotation moiré formed by the reference lattice and the specimen lattice.
Fig. 6.
Fig. 6. (a) HRTEM image of Ge/Si heterostructure interface; (b) Fourier transform of (a); (d) burgers vector with sampling moiré.
Fig. 7.
Fig. 7. (a) HRTEM image of Ge/Si heterostructure interface; (b) digital moiré fringe patterns; (c) fringe patterns extracted from (b).
Fig. 8.
Fig. 8. Experimental and theoretical strain fields: (a) experimental strain field; (b) strain field of P-N dislocation model; (c) - (k) strain field of Foreman dislocation model, with α ranging from 2 to 10, respectively.
Fig. 9.
Fig. 9. Interfacial strain and stress of Ge/Si heterostructure: (a) intercepted Ge/Si heterostructure interface; (b) strain field in the horizontal direction; (c) stress profile corresponding to the dash-dotted line in (b).

Equations (16)

Equations on this page are rendered with MathJax. Learn more.

I s ( x , y ) = N s 0 + N s 1 cos ( 2 π p s y + φ s )
I r ( x , y ) = N r 0 + N r 1 cos ( 2 π p r y + φ r )
I = I r I s = N r 0 + N r 1 sin ( 2 π p r y + φ r ) N s 0 N s 1 sin ( 2 π p s y + φ s )
1 δ = 1 p r 1 p s
I N 0 + N 1 cos ( 2 π δ y + φ 0 )
I n N 0 + N 1 c o s ( 2 π δ y + φ 0 + i π 2 ) , ( n = 1 , 2 , 3 , 4 ; i = 0 , 1 , 2 , 3 )
φ ( x , y ) = arctan I 4 ( x , y ) I 2 ( x , y ) I 1 ( x , y ) I 3 ( x , y )
{ u x = u r cos α 1 + u n cos α 2 u y = u r sin α 1 + u n sin α 2
{ ε x = u x / x ε y = u y / y
σ x = E 1 v 2 ( ε x + v ε y )
{ δ 1 = p r p s 1 p s 1 2 + p r 2 2 p r p s 1 cos θ 1 tan β 1 = p r sin θ 1 p s 1 p r cos θ 1
{ δ 2 = p r p s 2 p s 2 2 + p r 2 2 p r p s 2 cos θ 2 tan β 2 = p r sin θ 2 p s 2 p r cos θ 2
Δ β = β 2 β 1 = arc tan p r sin θ 2 p s 2 p r cos θ 2 arc tan p r sin θ 1 p s 1 p r cos θ 1
θ 2 = θ 1 = θ
Δ β = β 2 β 1 = arc tan p s 1 sin θ + p s 2 sin θ p s 1 p s 2 p r + p r p s 1 cos θ p s 2 cos θ
Δ β = β 2 β 1 = arc tan sin θ 1 cos θ
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.