# Effects of Electron Microscope Parameters and Sample Thickness on High Angle Annular Dark Field Imaging

Authors: Pucheng Yang; Zheng Li; Yi Yang; Rui Li; Lufei Qin; Yunhao Zou
Journal: Scanning
 (2022)

---

## Abstract

Scanning transmission electron microscopy (STEM) developed into a very important characterization tool for atomic analysis of crystalline specimens. High-angle annular dark field (HAADF) scanning transmission electron microscopy (STEM) has become one of the most powerful tools to visualize material structures at atomic resolution. However, the parameter of electron microscope and sample thickness is the important influence factors on HAADF-STEM imaging. The effect of convergence angle, spherical aberration, and defocus to HAADF imaging process has been analyzed through simulation. The applicability of two HAADF simulation software has been compared, and suggestions for their usage have been given.


---

## Body

## 1. Introduction

In recent years, dynamical diffraction is the major limitation to structure determination by electron methods, scanning transmission electron microscopy (STEM) which can effectively overcome this limitation by providing an incoherent image with electrons [1]. It has become a very popular and widespread technique and has been developed with different imaging modes (i.e., bright field (BF), annular bright field (ABF), and annular dark field (ADF)). High-resolution STEM using higher-angle scattered electrons (high-angle annular dark field, HAADF-STEM), where scattered electrons at higher angles are collected by an annular detector for STEM imaging, now attracts material scientists and semiconductor researchers [2] and becomes a powerful tool in characterizing nanoparticles [3, 4]), lithium ion batteries [5], quasicrystal [6–8], and alloy [9, 10]. HAADF has high resolution up to 40.5 pm [11, 12], image intuitive with the sensitive to chemical composition [13], and less damage to the sample, etc. It has become an important approach of material analysis at the atomic level due to the fact that its contrast depends highly on atomic number Z in a form of Znn=1.6–1.9 [14, 15].There are two steps in STEM imaging: first, the parallel electron wave is emitted by the electron gun and converged by the lens to form the convergent electron wave; second, the convergent electron wave scans the surface of sample point by point, and each sample point will get an exit wave. The effect of imaging process in these steps should be considered: the influence of lens on image including the convergence angle [16], spherical aberration and defocus [17], and detector angle [18] and the influence of sample including the sample thickness and the atomic number [15, 19]. In STEM imaging mode, there are two theories that have been proposed to reveal the image formation process, such as multislice method [20, 21] and Bloch wave method [1, 22]. The Bloch wave method has the advantage of high calculation accuracy, but it can only suit for perfect crystal [23] and has a long calculation time [24, 25]. Compared to the Bloch wave method, the multislice method has the advantages of fast calculation speed and wide applicability [26]. However, the images calculated by the multislice method are greatly influenced by the thickness of slice [27]. Thus, it is important to choose an appropriate theory according to the requirement. In recent years, there are several approximation theories applied in simulated calculation for speeding up the STEM imaging, such as the frozen phonon model [28] and the absorptive potential approximation [29]. The frozen phonon model divides the wave function into an average and a fluctuating part, which is suitable for carrying out the incoherent averages [28]. And the absorptive potential approximation was proposed based on the fast Fourier transform (FFT) multislice approach [29].In general, HAADF provides incoherent images without any phase problem and can be directly inverted to object without additional image simulations [30, 31]. However, in a series of HAADF investigations [31–34] under STEM imaging condition, crystal tilt, probe convergence angle unavoidable factor (e.g., specimen thickness [35, 36], and collection angles of the detector [36], etc.) will impact on the quality of high resolution images. For example, sample bending and deviating slightly from the zone axis will result in remarkable contrast reduction and could cause atoms to be considerably displaced in HAADF image. Intensities of atomic-resolution HAADF images of zone-axis-oriented specimens change with defocus at rates that depend on lattice spacing, thickness, and strain which also effect on the intensities [27].In this paper, the effect of electron microscope parameters and sample thickness on high angle annular dark field imaging was discussed in detail by simulation. In addition, simulation software QSTEM and Dr. Probe have been chosen for comparing their convenience in different simulation conditions.


## 2. The Frozen Phonon Model

In HAADF-STEM imaging, most of the signals received by annular detector come from phonon scattering [18].While frozen phonon model provides a simple method that the crystal potential is time-dependent under the assumption of independent atomic motion, under this assumption, the intensity of image IDZ can be described as [37]
(1)IDz=∫0zW2z′,tψz′2dz′.TheW2z′,t is related to the scattering factors, and the ψz′2 is the elastic intensity. The image intensity is essentially obtained from the signal received by the detector. Under the screening of big angle annular detector in HAADF-STEM imaging mode, the final signal gxp from the probe position xp should be
(2)gxp=∫ψtk,xp2Dkd2k,where the ψtk,xp is the wave function diffracted onto the detector plane, and Dk is the detector function [38].
(3)Dk=1forkDmin≤k≤kDmax0otherwise.The inner and outer angles of detector are the product of wavelengthλ and the maximum and the minimum of the wave vector k. Therefore, kDmaxand kDmin are determined by the angle of detector. Thus, the key to analyze the intensity of HAADF image at the same detector angle is the exit wave function. While the exit wave function is the product of the incident wave function and sample transfer function which are determined by electron microscope and sample parameters, therefore, in the subsequent sections, we will analyze the effect of these parameters on HAADF image.


## 3. Exit Wave

When the parameters of the incident electron beam and the sample are known, the equation of the exit wave can be obtained.It started from the stationary Schrödinger equation:(4)−ℏ22μ∇2ψ+Uψ=Eψ,where the ψ is the wave function, which presents the electron trajectory, ћ is the reduced Planck constant, U is the potential field function, E is the energy of the electron, and μ is the electron mass under relativistic correction. Due to the equation of wave-particle duality, λ=h/p, we know that K2=2meE/h2. Therefore, the Schrödinger equation can be written as
(5)∇2ψr+4π2K2ψr+4π2Urψr=0.While under the condition of high energy approximate, the electron wave function can be written as the form of modulation wave:(6)ψr=ϕrexp2πik⋅r.Substituting it into the Schrödinger equation, then(7)Δxyϕr+∂2∂z2ϕr+4πikxy∇xyϕr+4πikzdϕrdz+4π2Urϕr=0.Under the condition of high energy approximate,∂2/∂z2φr can be ignored, and we can obtain
(8)ϕr,z=expiz4πkzΔxy+4πikxy∇xy+4π2Uϕr,0.Equation (8) shows that the relationship of the incident wave and the exit wave is determined by following factors: z is the sample thickness, kz is the reciprocal of wavelength λ, Δxy is the direction of the scattering, and U is the crystal potential. These factors, effect on imaging process, will be analyzed in subsequent sections.


## 4. The Effect of Electron Microscope Parameters on HAADF

In STEM imaging mode, the parallel electron waves from the electron gun pass through electron microscope before they are incident to the sample and become convergent electron waves. The incident waves will be affected by convergence angle, spherical aberration, and defocus. In order to illustrate the effects, BaTiO3 was taken as an example for simulation, and its structure model is shown in Figure 1. The simulation parameters are as follows: accelerating voltage is 300 kV (wavelength is 1.97 pm); the sample thickness is 82.7 Å, the tilt angle is 0, the convergence angle is 25 mrad, the spherical aberration and defocus are 0, and the detector angle is 50-250 mrad. The simulation software is QSTEM [39] which uses the frozen phonon model based on the multislice method [40].Figure 1
Projection structure model of BaTiO3: (a) at [001] direction and (b) at [011] direction.


### 4.1. Convergence Angle

When convergent electron waves are incident on the positionxp of the sample surface, the equation of incident wave function [38] is
(9)ψpx,xp=Ap∫0kmaxexp−iχk−2πik⋅x−xpd2k.Theχk is the aberration function, Ap is the factor which comments ∫ψpx,xp2d2x=1, the influence of convergence angle on the incident wave function is the upper limit of integral kmax=α/λ, and α is the convergence semiangle.As shown in Figure2(a), when the convergence angle increases, the contrast of center atom Ti will decrease. It is more obviously in thick sample (i.e., Figure 2(b)). Besides, the spots of Ba atom are larger at small convergence angle, but smaller at large convergence angle. The same situations are obtained in reference [34] which analyzed the intensity profiles of 195 nm Si0.8Ge0.2 in different convergence angles. In the HAADF-STEM model, the electron probe with small convergence angle is more sensitive to crystal potential, the atomic brightness is larger, and the spot of atom is larger in small convergence angle.Figure 2
Simulated HAADF image (thesourcesize=0.8Å, detectorangle=50mrad~250mrad, defocus=0, sphericalaberration=0) and its intensity profile along the diagonal as the arrow shown of BaTiO3 at [001] axis zones with different convergence angles: (a) 42.3 Å thickness and (b) 82.7 Å thickness.


### 4.2. Spherical Aberration and Defocus

In fact, the spherical aberration of a real STEM is not 0. The convergence angle, spherical aberration, and defocus should be consider to choose the best condition for imaging [41, 42]. In Equation (9), χk=πλk20.5Csλ2k2−Δf, it can be found that spherical aberration has a great influence on the incident convergent electron wave. In order to correct the influence of spherical aberration on the imaging process, the Scherzer focus condition [43] had been proposed. When the convergence angle is 10 mrad and spherical aberration is less than 0.1 mm, Figure 3(a) with the Scherzer focus condition has little changes as the spherical aberration increases. As shown in Figure 3(b), when defocus is 0, the simulated HAADF image becomes more and more anamorphose as the spherical aberration increases. When the convergence angle is 25 mrad, as shown in Figure 4, the anamorphose of the simulated HAADF images becomes worse as the spherical aberration is larger than 0.1 mm. It can be concluded that small convergence angle has good HAADF image at small spherical aberration conditions.Figure 3
Simulated HAADF images (samplethickness=82.7Å, sourcesize=0.8Å, detectorangle=50mrad~250mrad, convergenceangle=10mrad) with different spherical aberrations and defocus: (a) Scherzer focus condition and (b) defocus is 0.Figure 4
Simulated HAADF images (samplethickness=82.7Å, sourcesize=0.8Å, detectorangle=50mrad~250mrad, convergenceangle=25mrad) with different spherical aberrations and defocus: (a) Scherzer focus condition and (b) defocus is 0.In HAADF-STEM mode, the Scherzer focus condition (Δf=−1.15Csλ0.5) makes the transfer function have a wide flat area which leads to simulation images that have less abnormal structure information [17]. Therefore, it is obvious that the anamorphose deformation of the simulated HAADF images with the Scherzer condition is less than those who do not satisfy with it. And when the Scherzer focus condition is satisfied, the best convergence semiangle is
(10)α0=1.41λCs0.25.Combined with the limit point resolution in STEM mode,d0=0.61λ/α0, it can be concluded that
(11)d0=0.43Cs0.25λ0.75.So, appropriate defocus and convergence angle should be chosen according to its spherical aberration. At present, the resolution of STEM image with spherical aberration correction has reached 40.5 pm [11]. Therefore, for an electron microscope with a certain spherical aberration, the key to high-quality STEM image is the appropriate defocus and convergence angle.


### 4.3. Detector Angle

In HAADF-STEM mode, the annular detector mainly receives high angle scattered electrons. Its scattering intensity can be expressed as the internal of inner angleθ1 to outer angle θ2 [18]:
(12)σθ1,θ2=mm0Z2λ44π3a021θ12+θ02−1θ22+θ02.In equation (12), m is the mass of high-velocity electrons, m0 is the static mass of electron, Z is the atomic number, λ is the wavelength, α0 is the Bohr radius, and θ0 is the Born characteristic scattering angle.When the thickness of a sample ist, the intensity of exit electron wave is I, the number of atoms per unit volume is N, and the scattering intensity can be expressed as equation (13):
(13)Is=σθ1,θ2⋅NtI.Equation (12) can be transformed into equation (14):
(14)σθ1,θ2=mm0Z2λ44π3a02θ22−θ12θ12+θ02θ22+θ02.It obvious in Equation (14) that the larger inner angle θ1, the weaker the image intensity. The partial derivative of outer angle θ2 can be obtained:
(15)∂σθ1,θ2∂θ2=mm0Z2λ44π3a022θ2θ22+θ022.The overall intensity of the image will be stronger with the increase of the outer angle. This also reflects that a larger detector angle range can obtain a stronger intensity image. However, the small inner angle may cause other signals to be detected (e.g., the diffraction contrast caused by Bragg reflection). Larger inner angle will also affect the relationship between the contrast of the atomic column and the atomic numberZ in the final image which will be analyzed in the following sections. The intensity of the image atomic column is not simply linear with the sample thickness and the range of the detector angle [1]. Therefore, it is necessary to consider both the detector angle and image interpretation.


### 4.4. The Effect of Sample Thickness on HAADF

In imaging process, the incident wave is converged by the electron microscope and the surface of the sample is scanned point by point. There is no doubt that, sample parameters are also important influence on HAADF image. This section will discuss the influence of sample thickness on imaging process which is based on the simulation calculation of the multislice method.It can be found in Figure5 that the contrast of center atoms (Ti) is stronger in thick sample. This phenomenon has been explained in [15]. Equation (12) shows that the intensity of atomic column is proportional to the square of the atomic number Z, while it should not consider the scattering intensity of one atomic array only in HAADF-STEM. For thin sample, if convergence semiangle α0 and the inner angle θ1 of the detector have the relationship:
(16)θ1≥3α0,then the coherent effect between different atomic arrays of thin samples can be ignored. The approximation of the image intensity to the n-th power of the atomic number Z is very accurate. The range of n is 1.6-1.9 which is related to the inner and outer angle of detector.Figure 5
Simulated HAADF-STEM image (sourcesize=0.8Å, detectorangle=50mrad~250mrad, convergenceangle=25mrad, defocus=0, sphericalaberration=0) of BaTiO3 sample at different thicknesses and its intensity profile along the diagonal as the arrow shown.


## 4.1. Convergence Angle

When convergent electron waves are incident on the positionxp of the sample surface, the equation of incident wave function [38] is
(9)ψpx,xp=Ap∫0kmaxexp−iχk−2πik⋅x−xpd2k.Theχk is the aberration function, Ap is the factor which comments ∫ψpx,xp2d2x=1, the influence of convergence angle on the incident wave function is the upper limit of integral kmax=α/λ, and α is the convergence semiangle.As shown in Figure2(a), when the convergence angle increases, the contrast of center atom Ti will decrease. It is more obviously in thick sample (i.e., Figure 2(b)). Besides, the spots of Ba atom are larger at small convergence angle, but smaller at large convergence angle. The same situations are obtained in reference [34] which analyzed the intensity profiles of 195 nm Si0.8Ge0.2 in different convergence angles. In the HAADF-STEM model, the electron probe with small convergence angle is more sensitive to crystal potential, the atomic brightness is larger, and the spot of atom is larger in small convergence angle.Figure 2
Simulated HAADF image (thesourcesize=0.8Å, detectorangle=50mrad~250mrad, defocus=0, sphericalaberration=0) and its intensity profile along the diagonal as the arrow shown of BaTiO3 at [001] axis zones with different convergence angles: (a) 42.3 Å thickness and (b) 82.7 Å thickness.


## 4.2. Spherical Aberration and Defocus

In fact, the spherical aberration of a real STEM is not 0. The convergence angle, spherical aberration, and defocus should be consider to choose the best condition for imaging [41, 42]. In Equation (9), χk=πλk20.5Csλ2k2−Δf, it can be found that spherical aberration has a great influence on the incident convergent electron wave. In order to correct the influence of spherical aberration on the imaging process, the Scherzer focus condition [43] had been proposed. When the convergence angle is 10 mrad and spherical aberration is less than 0.1 mm, Figure 3(a) with the Scherzer focus condition has little changes as the spherical aberration increases. As shown in Figure 3(b), when defocus is 0, the simulated HAADF image becomes more and more anamorphose as the spherical aberration increases. When the convergence angle is 25 mrad, as shown in Figure 4, the anamorphose of the simulated HAADF images becomes worse as the spherical aberration is larger than 0.1 mm. It can be concluded that small convergence angle has good HAADF image at small spherical aberration conditions.Figure 3
Simulated HAADF images (samplethickness=82.7Å, sourcesize=0.8Å, detectorangle=50mrad~250mrad, convergenceangle=10mrad) with different spherical aberrations and defocus: (a) Scherzer focus condition and (b) defocus is 0.Figure 4
Simulated HAADF images (samplethickness=82.7Å, sourcesize=0.8Å, detectorangle=50mrad~250mrad, convergenceangle=25mrad) with different spherical aberrations and defocus: (a) Scherzer focus condition and (b) defocus is 0.In HAADF-STEM mode, the Scherzer focus condition (Δf=−1.15Csλ0.5) makes the transfer function have a wide flat area which leads to simulation images that have less abnormal structure information [17]. Therefore, it is obvious that the anamorphose deformation of the simulated HAADF images with the Scherzer condition is less than those who do not satisfy with it. And when the Scherzer focus condition is satisfied, the best convergence semiangle is
(10)α0=1.41λCs0.25.Combined with the limit point resolution in STEM mode,d0=0.61λ/α0, it can be concluded that
(11)d0=0.43Cs0.25λ0.75.So, appropriate defocus and convergence angle should be chosen according to its spherical aberration. At present, the resolution of STEM image with spherical aberration correction has reached 40.5 pm [11]. Therefore, for an electron microscope with a certain spherical aberration, the key to high-quality STEM image is the appropriate defocus and convergence angle.


## 4.3. Detector Angle

In HAADF-STEM mode, the annular detector mainly receives high angle scattered electrons. Its scattering intensity can be expressed as the internal of inner angleθ1 to outer angle θ2 [18]:
(12)σθ1,θ2=mm0Z2λ44π3a021θ12+θ02−1θ22+θ02.In equation (12), m is the mass of high-velocity electrons, m0 is the static mass of electron, Z is the atomic number, λ is the wavelength, α0 is the Bohr radius, and θ0 is the Born characteristic scattering angle.When the thickness of a sample ist, the intensity of exit electron wave is I, the number of atoms per unit volume is N, and the scattering intensity can be expressed as equation (13):
(13)Is=σθ1,θ2⋅NtI.Equation (12) can be transformed into equation (14):
(14)σθ1,θ2=mm0Z2λ44π3a02θ22−θ12θ12+θ02θ22+θ02.It obvious in Equation (14) that the larger inner angle θ1, the weaker the image intensity. The partial derivative of outer angle θ2 can be obtained:
(15)∂σθ1,θ2∂θ2=mm0Z2λ44π3a022θ2θ22+θ022.The overall intensity of the image will be stronger with the increase of the outer angle. This also reflects that a larger detector angle range can obtain a stronger intensity image. However, the small inner angle may cause other signals to be detected (e.g., the diffraction contrast caused by Bragg reflection). Larger inner angle will also affect the relationship between the contrast of the atomic column and the atomic numberZ in the final image which will be analyzed in the following sections. The intensity of the image atomic column is not simply linear with the sample thickness and the range of the detector angle [1]. Therefore, it is necessary to consider both the detector angle and image interpretation.


## 4.4. The Effect of Sample Thickness on HAADF

In imaging process, the incident wave is converged by the electron microscope and the surface of the sample is scanned point by point. There is no doubt that, sample parameters are also important influence on HAADF image. This section will discuss the influence of sample thickness on imaging process which is based on the simulation calculation of the multislice method.It can be found in Figure5 that the contrast of center atoms (Ti) is stronger in thick sample. This phenomenon has been explained in [15]. Equation (12) shows that the intensity of atomic column is proportional to the square of the atomic number Z, while it should not consider the scattering intensity of one atomic array only in HAADF-STEM. For thin sample, if convergence semiangle α0 and the inner angle θ1 of the detector have the relationship:
(16)θ1≥3α0,then the coherent effect between different atomic arrays of thin samples can be ignored. The approximation of the image intensity to the n-th power of the atomic number Z is very accurate. The range of n is 1.6-1.9 which is related to the inner and outer angle of detector.Figure 5
Simulated HAADF-STEM image (sourcesize=0.8Å, detectorangle=50mrad~250mrad, convergenceangle=25mrad, defocus=0, sphericalaberration=0) of BaTiO3 sample at different thicknesses and its intensity profile along the diagonal as the arrow shown.


## 5. Comparison HAADF Simulation Software

There is several software for HAADF simulation, and Dr. Probe and QSTEM are two free software which has wide used for researchers. In Dr. Probe software, the frozen-lattice approach [44] was used to simulate thermal-diffuse scattering (TDS). In this section, a horizontal comparison has been made between the two software to provide some reference for researchers. The simulation parameters of the two software are shown in Table 1.Table 1
The simulation parameters of QSTEM and Dr. Probe.
ParametersQSTEM and Dr. ProbeAccelerating voltage300 kVSources size0.8 ÅConvergence angle25 mradCc1 mmCs0 mmHADDF detector angle50-250 mrad


### 5.1. Unit Cell with Small Size

PbTiO3 was chosen as the simulated material whose lattice constants are a=3.90Å, b=3.90Å, and C=4.15Å. As shown in Figure 6, the white highlights are Pb (Z=82) atoms, and the gray and white ones are Ti atoms. The atomic number of O is too small which does not show in the image. Figure 6(c) is the structure projection of PbTiO3 at the [001] zone axis. The simulation results of the two software for small crystal cell have good quality, and both correctly reflect the crystal structure and atomic phase arrangement, while the contrast of Figure 6(b) is better. The contrast of Ti atoms in QSTEM is very fuzzy and difficult to distinguish. Moreover, the difference of the contrast and image point size of Pb atoms and Ti atoms in Figure 6(a) are too large. The image quality of Figure 6(a) is unnatural and unreal compared with Figure 6(b).Figure 6
Simulated HAADF-STEM images (samplethickness=83.0Å, sourcesize=0.8Å, detectorangle=50mrad~250mrad, convergenceangle=25mrad, defocus=0, sphericalaberration=0) along PbTiO3 [001]: (a) QSTEM, (b) Dr. Probe, and (c) structure model projection of PbTiO3 at the [001] zone axis.Figure7 shows the simulated HAADF image of PbTiO3 along [001] at thicknesses with 40 unit cells. As shown in Figure 7(a), there are black dots in the center of the white bright spots representing Pb, but the contrast of Ti atom is improved obviously as compared with Figure 6(a). However, Figure 7(b) has problems in showing the contrast of Ti atoms, and it was difficult to observe the distribution of Ti atoms.Figure 7
Simulated HAADF-STEM images (samplethickness=166.0Å, sourcesize=0.8Å, detectorangle=50mrad~250mrad, convergenceangle=25mrad, defocus=0, sphericalaberration=0) along PbTiO3 [001]: (a) QSTEM and (b) Dr. Probe.


### 5.2. Unit Cells with Moderate Size

MgAlO4 (a=b=c=8.0858Å) was chosen as the simulated material with moderate unit cell size whose lattice parameters are between PbTiO3 and Mg44Rh7. The thickness is 10 unit cells. The simulation results are shown in Figure 8. Compared with Figure 8(c), the white spots in Figure 8(a) and 8(b) represent Mg (Z=12) atom. In Figure 8(a), the brightest atom is Mg atom, and the surrounding gray atom is Al atom. Due to the Mg atoms which are too close to the Al atoms in the projection of [001], the simulated image cannot distinguish them. Obviously, when unit cell has moderate size, QSTEM does better.Figure 8
Simulated HAADF-STEM images (samplethickness=80.9Å, sourcesize=0.8Å, detectorangle=50mrad~250mrad, convergenceangle=25mrad, defocus=0, sphericalaberration=0) along MgAlO4 [001]: (a) QSTEM, (b) Dr. Probe, and (c) structure model projection of MgAlO4 at the [001] zone axis.


### 5.3. Unit Cell with Big Size

Mg44Rh7 (a=b=c=20.148Å) was chosen as the simulated material. The simulation results are shown in Figure 9. As shown in Figures 9(a) and 9(b), the white bright spots represent Rh (Z=45) atoms, and the gray bright spots around the white bright spots represent Mg (Z=12) atoms. Figure 9(c) is the projection structure model of Mg44Rh7 at the [001] zone axis. By comparing Figure 9(c), the simulated results of the two software are consistent with the projection crystal structure. However, there are differences in Mg atoms. Figure 9(b) can well display the contrast of Mg atoms, but Figure 8(a) is not obvious in showing the contrast of Mg atoms, and the contrast of Mg atoms among the four Rh atoms is insufficient.Figure 9
Simulated HAADF-STEM images (samplethickness=201.5Å, s, detectorangle=50mrad~250mrad, convergenceangle=25mrad, defocus=0, sphericalaberration=0) along Mg44Rh7 [001]: (a) QSTEM, (b) Dr. Probe, and (c) structure model projection of Mg44Rh7 at the [001] zone axis.


### 5.4. Calculate Speed

In order to control variables, Dr. Probe and QSTEM were kept consistent during calculation, and single-core calculation was performed. QSTEM and Dr. Probe are used on the same PC. According to the simulation time records, the calculate time of the two software is shown in Table2. QSTEM generally takes much longer time than Dr. Probe. In the simulation of small cell PbTiO3, the time of QSTEM was about 397% longer than that of Dr. Probe; in the simulation of medium cell MgAlO4, the time of QSTEM was about 164% longer than that of Dr. Probe; in the simulation of large cell Mg44Rh7, the time of QSTEM was about 440% longer than that of Dr. Probe. Therefore, Dr. Probe has higher computational efficiency than that of QSTEM.Table 2
Simulate time of QSTEM and Dr. probe.
CrystalQSTEMDr. ProbePbTiO34641 s934 sMgAlO434443 s13045 sMg44Rh7204573 s37861 s


## 5.1. Unit Cell with Small Size

PbTiO3 was chosen as the simulated material whose lattice constants are a=3.90Å, b=3.90Å, and C=4.15Å. As shown in Figure 6, the white highlights are Pb (Z=82) atoms, and the gray and white ones are Ti atoms. The atomic number of O is too small which does not show in the image. Figure 6(c) is the structure projection of PbTiO3 at the [001] zone axis. The simulation results of the two software for small crystal cell have good quality, and both correctly reflect the crystal structure and atomic phase arrangement, while the contrast of Figure 6(b) is better. The contrast of Ti atoms in QSTEM is very fuzzy and difficult to distinguish. Moreover, the difference of the contrast and image point size of Pb atoms and Ti atoms in Figure 6(a) are too large. The image quality of Figure 6(a) is unnatural and unreal compared with Figure 6(b).Figure 6
Simulated HAADF-STEM images (samplethickness=83.0Å, sourcesize=0.8Å, detectorangle=50mrad~250mrad, convergenceangle=25mrad, defocus=0, sphericalaberration=0) along PbTiO3 [001]: (a) QSTEM, (b) Dr. Probe, and (c) structure model projection of PbTiO3 at the [001] zone axis.Figure7 shows the simulated HAADF image of PbTiO3 along [001] at thicknesses with 40 unit cells. As shown in Figure 7(a), there are black dots in the center of the white bright spots representing Pb, but the contrast of Ti atom is improved obviously as compared with Figure 6(a). However, Figure 7(b) has problems in showing the contrast of Ti atoms, and it was difficult to observe the distribution of Ti atoms.Figure 7
Simulated HAADF-STEM images (samplethickness=166.0Å, sourcesize=0.8Å, detectorangle=50mrad~250mrad, convergenceangle=25mrad, defocus=0, sphericalaberration=0) along PbTiO3 [001]: (a) QSTEM and (b) Dr. Probe.


## 5.2. Unit Cells with Moderate Size

MgAlO4 (a=b=c=8.0858Å) was chosen as the simulated material with moderate unit cell size whose lattice parameters are between PbTiO3 and Mg44Rh7. The thickness is 10 unit cells. The simulation results are shown in Figure 8. Compared with Figure 8(c), the white spots in Figure 8(a) and 8(b) represent Mg (Z=12) atom. In Figure 8(a), the brightest atom is Mg atom, and the surrounding gray atom is Al atom. Due to the Mg atoms which are too close to the Al atoms in the projection of [001], the simulated image cannot distinguish them. Obviously, when unit cell has moderate size, QSTEM does better.Figure 8
Simulated HAADF-STEM images (samplethickness=80.9Å, sourcesize=0.8Å, detectorangle=50mrad~250mrad, convergenceangle=25mrad, defocus=0, sphericalaberration=0) along MgAlO4 [001]: (a) QSTEM, (b) Dr. Probe, and (c) structure model projection of MgAlO4 at the [001] zone axis.


## 5.3. Unit Cell with Big Size

Mg44Rh7 (a=b=c=20.148Å) was chosen as the simulated material. The simulation results are shown in Figure 9. As shown in Figures 9(a) and 9(b), the white bright spots represent Rh (Z=45) atoms, and the gray bright spots around the white bright spots represent Mg (Z=12) atoms. Figure 9(c) is the projection structure model of Mg44Rh7 at the [001] zone axis. By comparing Figure 9(c), the simulated results of the two software are consistent with the projection crystal structure. However, there are differences in Mg atoms. Figure 9(b) can well display the contrast of Mg atoms, but Figure 8(a) is not obvious in showing the contrast of Mg atoms, and the contrast of Mg atoms among the four Rh atoms is insufficient.Figure 9
Simulated HAADF-STEM images (samplethickness=201.5Å, s, detectorangle=50mrad~250mrad, convergenceangle=25mrad, defocus=0, sphericalaberration=0) along Mg44Rh7 [001]: (a) QSTEM, (b) Dr. Probe, and (c) structure model projection of Mg44Rh7 at the [001] zone axis.


## 5.4. Calculate Speed

In order to control variables, Dr. Probe and QSTEM were kept consistent during calculation, and single-core calculation was performed. QSTEM and Dr. Probe are used on the same PC. According to the simulation time records, the calculate time of the two software is shown in Table2. QSTEM generally takes much longer time than Dr. Probe. In the simulation of small cell PbTiO3, the time of QSTEM was about 397% longer than that of Dr. Probe; in the simulation of medium cell MgAlO4, the time of QSTEM was about 164% longer than that of Dr. Probe; in the simulation of large cell Mg44Rh7, the time of QSTEM was about 440% longer than that of Dr. Probe. Therefore, Dr. Probe has higher computational efficiency than that of QSTEM.Table 2
Simulate time of QSTEM and Dr. probe.
CrystalQSTEMDr. ProbePbTiO34641 s934 sMgAlO434443 s13045 sMg44Rh7204573 s37861 s


## 6. Conclusion

With the development of scanning transmission electron microscopy, image interpretation in HAADF-STEM mode has become particularly important. In the imaging process, the influence of the electron microscope parameters and sample parameters must be considered. In this paper, the effect of electron microscope parameters and sample thickness on high angle annular dark field imaging were discussed in detail by simulation and experiment. In addition, simulation software QSTEM and Dr. Probe have been chosen for comparing their convenience in different simulation conditions. The conclusion is as follows:(1)
Appropriate convergence angle is one of the key parameters for getting a good HAADF image. Small convergence angle has better image quality(2)
Appropriate defocus and convergence angle should be chosen according to its spherical aberration(3)
Under the condition of thin sample, sample thickness has little effect on HAADF image(4)
QSTEM and Dr. Probe both are excellent like simulation software. Dr. Probe has higher computational efficiency than that of QSTEM



---

*Source: 8503314-2022-03-20.xml*