Development of Advanced Electron Tomography in Materials Science Based On TEM and STEM
Development of Advanced Electron Tomography in Materials Science Based On TEM and STEM
         Mao-hua LI1, Yan-qing YANG1, Bin HUANG1, Xian LUO1, Wei ZHANG1, Ming HAN1, Ji-gang RU2
  1. State Key Laboratory of Solidification Processing, Northwestern Polytechnical University, Xi’an 710072, China;
                        2. Beijing Institute of Aeronautical Materials, Beijing 100095, China
                                       Received 24 November 2013; accepted 12 May 2014
Abstract: The recent developments of electron tomography (ET) based on transmission electron microscopy (TEM) and scanning
transmission electron microscopy (STEM) in the field of materials science were introduced. The various types of ET based on TEM
as well as STEM were described in detail, which included bright-field (BF)-TEM tomography, dark-field (DF)-TEM tomography,
weak-beam dark-field (WBDF)-TEM tomography, annular dark-field (ADF)-TEM tomography, energy-filtered transmission electron
microscopy (EFTEM) tomography, high-angle annular dark-field (HAADF)-STEM tomography, ADF-STEM tomography,
incoherent bright field (IBF)-STEM tomography, electron energy loss spectroscopy (EELS)-STEM tomography and X-ray energy
dispersive spectrometry (XEDS)-STEM tomography, and so on. The optimized tilt series such as dual-axis tilt tomography, on-axis
tilt tomography, conical tilt tomography and equally-sloped tomography (EST) were reported. The advanced reconstruction
algorithms, such as discrete algebraic reconstruction technique (DART), compressed sensing (CS) algorithm and EST were
overviewed. At last, the development tendency of ET in materials science was presented.
Key words: electron tomography; materials science; transmission electron microscopy; scanning transmission electron microscopy
 Foundation item: Projects (51071125, 51201135) supported by the National Natural Science Foundation of China; Project (B08040) supported by the
                  Program of Introducing Talents of Discipline to Universities, China
 Corresponding author: Yan-qing YANG; Tel/Fax: +86-29-88460499; E-mail: yqyang@nwpu.edu.cn
 DOI: 10.1016/S1003-6326(14)63441-5
3032                        Mao-hua LI, et al/Trans. Nonferrous Met. Soc. China 24(2014) 3031í3050
physical properties of the object.                                  2.1.2 DF-TEM tomography
     The aim of this review is mainly firstly to                         In general, DF-TEM images do not fulfill the
summarize the novel tomographic imaging modes, then                 projection requirement because the image intensity
to introduce the optimized tilt series and the advanced             varies rapidly in a complicated manner with sample
reconstruction algorithms, and finally to present the               orientation. It is very difficult to acquire a tilt series of
further developments of ET.                                         DF-TEM images of a crystalline specimen for
                                                                    tomography. Fortunately, KIMURA et al [9] have
2 Novel tomographic imaging modes                                   successfully reconstructed the D1a-ordered Ni4Mo
                                                                    precipitates in Ni–Mo alloy by DF-TEM tomography,
2.1 Based on TEM                                                    and introduced how to obtain a tilt series of DF-TEM
2.1.1 BF-TEM tomography                                             images of the Ni4Mo precipitates. Firstly, a systematic
     BF-TEM tomography is suitable for amorphous                    row containing the D1a superlattice reflection at (4/5, 2/5,
materials, where mass thickness contrast is dominant. A             0) was parallel to the tilt axis of the holder by placing
few successful examples were reported, such as studying             carefully the specimen on the specimen holder, which
the location and distribution of metal (oxide) particles in         made the D1a superlattice reflection exist in the tilt series
zeolites and catalyst materials in 3D, and the location             from í60° to +60° (see Fig. 2(a)). Secondly, by adopting
and the connectivity of pores in mesoporous materials               the (4/5, 2/5, 0) superlattice reflection, a tilt series of
[4í8]. Figure 1 shows the 3D reconstruction of an                   DF-TEM images of the Ni4Mo precipitates were
Au/SBAí15 model catalyst particle, revealing the size               recorded (see Fig. 2(b)). Finally, the 3D shape and
and location of Au particles inside the support material            position of the Ni4Mo precipitates were reconstructed
clearly [8]. However, for crystalline materials, BF-TEM             (see Fig. 2(c)).
images are mainly dominated by diffraction contrast that            2.1.3 WBDF-TEM tomography
is highly sensitive to the direction of the incident beam.               Although coherent diffraction contrast seems to
Therefore, it is very difficult to obtain BF-TEM images             violate the projection requirement, WBDF-TEM
with clear contrast through the entire tilt range. In               tomography has recently been proposed to investigate the
general, BF-TEM tomography is unsuitable for                        3D distribution of dislocation network in GaN film
crystalline materials.                                              [10,11] and dislocation-precipitate interactions in an
Fig. 1 Au catalyst nanoparticles inside SBA-15 mesoporous support (a), surface rendering of slice (b) from (a), surface rendering of
Au nanoparticles only (c), BF images at í55° (d), 0° (e) and +55° (f) [8]
                          Mao-hua LI, et al/Trans. Nonferrous Met. Soc. China 24(2014) 3031í3050                         3033
AlíMgíSc alloy [12]. Figure 3(a) shows the complex              reconstruction. To meet the projection requirement, the
3D dislocation network in a GaN epilayer, revealing             WBDF-TEM condition must maintain constant through
threading dislocations at small-angle grain boundaries of       the entire tilt range, so that the dislocations are visible in
individual domains (D), a dislocation bundle associated         the entire tilt series. Unfortunately, data acquisition in
with a crack tip (B), in-plane dislocations caused by           WBDF-TEM tomography is very difficult because the
threading dislocations turning over (T), a jog of an            diffraction condition must be lined up exactly along the
in-plane dislocation by a threading dislocation (J), and        entire tilt range. Besides, during tilting, general
mixed dislocation (M) [10]. Figures 3(b) and (c) show           acquisition software cannot automatically complete due
3D relationships between the dislocations and the Al3Sc         to the unstable image contrast. Consequently, the tilt
particle, among the twist boundary, the dislocations, and       series of WBDF-TEM images have to be obtained
the Al3Sc particle, respectively. The WBDF-TEM                  manually.
technique is suitable for observation of dislocations           2.1.4 ADF-TEM tomography
because a WBDF-TEM image of dislocation is narrow in                  A schematic diagram of ADF-TEM is shown in
width and close to the dislocation core. Therefore,             Fig. 4(a). An annular objective aperture is inserted in the
WBDF-TEM tomography is probable for 3D dislocation              back focal plane of the objective lens, which shuts a
3034                        Mao-hua LI, et al/Trans. Nonferrous Met. Soc. China 24(2014) 3031í3050
Fig. 3 3D reconstructions based on WBDF-TEM imaging: (a) Complex 3D dislocation network in GaN epilayer; (b), (c) 3D
relationships between dislocations and Al3Sc particles, among twist boundary, dislocations, and Al3Sc particle, respectively [12]
Fig. 4 Schematic diagram of ADF-TEM (left) and secondary electron image of annular aperture (right) (a), ADF-TEM image of two
CdTe tetrapods (b) and isosurface visualization of 3D CdTe tetrapods by ADF-TEM tomography (c)
Fig. 5 EFTEM tomographic reconstructions: (a) Magnetite crystal chain with Fe (red) and O (green) (One section in the upper right
corner is perpendicular to z axis. Other sections in A–E are perpendicular to the chain axis from five of the crystals [16]); (b) C
nanohorn with Cu core from elemental maps [17]
Fig. 6 3D reconstructions using STEM-HAADF tomography: (a) Ge precipitates with various morphologies such as plates (blue),
tetrahedral (green), octahedral (orange), rod-shaped (yellow) and irregular shapes (white) [20]; (b) Distribution of S precipitates
along helical dislocation [21]; (c) Morphology and distribution of S precipitates from direction perpendicular to central part of grain
boundary [21,22]; (d) Configuration of S precipitates (yellow) around dispersoid (blue), with both needle-like and granular
morphologies [21]
in AlíCuíMg alloy (see Figs. 6(b)í(d)). Besides, this                 Exposure (scanning) time per image generally requires
technology has been used to reconstruct GuinieríPreston               10í30 s. Long time scanning not only produces more
zones in the AlíAg system [23], Au nanoparticles                      scanning noise but also makes it unsuitable for
supported on TiO2 catalyst [24], Sn-rich quantum dots                 beam-sensitive samples. HAADF-STEM tomography is
embedded in a Si matrix [25], CeO2 nanoparticles [26],                unfit for some nanostructures such as lattice defects,
TiO2 nanotubes [27], GaPíGaAs core-shell nanowires                    grain boundaries, orientation variants of ordered
[28], and so on.                                                      domains, because these nanostructures are only related to
     Although HAADF-STEM tomography is a very                         diffraction contrast images. HAADF-STEM tomography
powerful tool to reconstruct actual 3D structure of                   is unsuitable for multiphase materials which consist of
nanoscale materials with complex morphology and                       the elements with neighboring atomic number because of
chemistry, it also has some unavoidable disadvantages.                low-contrast.
                            Mao-hua LI, et al/Trans. Nonferrous Met. Soc. China 24(2014) 3031í3050                        3037
2.2.2 ADF-STEM tomography                                          The ADF-STEM tomography allows large deviations
      Both the inner collection angle of an annular                from the exact diffraction vector to image dislocations
detector and the convergence angle for ADF-STEM                    along a tilt series, and the dislocation contrast is much
imaging are smaller than that for HAADF-STEM                       more consistent by ADF-STEM imaging than by WBDF
imaging. TANAKA et al [29,30] have successfully                    imaging in the entire tilt range. As a result, automation of
reconstructed dislocations in a Si crystal by ADF-STEM             data acquisition in ADF-STEM tomography is
tomography. Figure 7(a) shows a schematic diagram                  practicable for reconstruction.
showing the diffraction vector of g=220 maintained in              2.2.3 IBF-STEM tomography
the annular detector during the entire tilt series. Figure               For HAADF-STEM technique, the acceptable
7(b) shows 3D dislocations viewed from different                   sample thickness is 50í100 nm. When samples are
orientations. The dislocations are colored according to            thicker, the image contrast reversal would arise [31]. The
their slip systems.                                                reason for this phenomenon is that the thick and/or high
                                                                   atomic number materials lead to increased backscatter
                                                                   and multiple scattering events to high angles. The angle
                                                                   of the electron beam detector in IBF-STEM is generally
                                                                   set to be 0í100 mrad to suppress diffraction contrast, and
                                                                   avoid contrast reversal at all thick materials. ERCIUS et
                                                                   al [31] have investigated thick copper microelectronic
                                                                   structures with a stress void by IBF-STEM tomography
                                                                   and HAADF-STEM tomography, respectively (see
                                                                   Fig. 8). The results indicate that the 3D reconstruction
                                                                   using IBF-STEM imaging shows the precise location of
                                                                   the stress void without artifacts.
                                                                   2.2.4 EELS-STEM tomography
                                                                         EELS-STEM spectrum imaging is successfully
                                                                   utilized to analyze the elemental, physical and chemical
                                                                   state information in nanoscale materials [44í46]. In
                                                                   Refs. [32í34], researchers have used EELS-STEM
                                                                   tomography to probe a W-to-Si contact from
                                                                   semiconductor device, a ZnO thin film, mesoporous a
                                                                   Co3O4 particles filled with FexCo3íxO4 and multiwalled
                                                                   carbon nanotubes (MWCNT)–nylon nanocomposite in
                                                                   3D. Figure 9(a) shows a schematic diagram revealing the
                                                                   collection of an energy-loss series at each tilt angle and
                                                                   each energy window with a specific voxel P in the 3D
                                                                   volume. Figure 9(b) shows surface renders of the
                                                                   MWCNT (purple) within the nylon (gray) viewed from
                                                                   different angles using plasmon-loss electrons. The hole
                                                                   in the nylon and the voids in the MWCNT are obviously
                                                                   visible in 3D [34].
                                                                   2.2.5 XEDS-STEM tomography
                                                                         Recently, XEDS spectrum imaging technique has
                                                                   been utilized to resolve 3D chemical studies in nanoscale
Fig. 7 Schematic diagram of diffraction vector of g=220            materials [35í37]. Compared with HAADF-STEM
maintained in annular detector [29] (a), 3D reconstruction of      tomography, XEDS-STEM tomography has obvious
dislocations in Si crystal by tilt series of ADF-STEM images (b)   superiority in enhancing chemical contrast especially
[30]                                                               between neighboring atomic numbers in the multiphase
                                                                   materials. Compared with EFTEM and EELS spectrum
     For dislocation, ADF-STEM tomography technique                imaging, XEDS imaging can obtain the maps of a large
is more feasible than WBDF-TEM tomography. The                     number of elements simultaneously. GENC et al [35]
reasons include the following aspects: 1) The                      have demonstrated 3D distribution of elements in a
ADF-STEM technique suppresses diffraction contrast                 Li1.2Ni0.2Mn0.6O2 (LNMO) nanoparticle based on
because of multiple beams and the convergent beam,                 XEDS-STEM imaging. Figure 10(a) shows 3D
while the WBDF imaging utilizes only one reflection; 2)            elemental distribution of Mn (red) indicating relatively
3038                         Mao-hua LI, et al/Trans. Nonferrous Met. Soc. China 24(2014) 3031í3050
Fig. 11 Schematic diagrams of geometrical configuration between pillar-shaped specimen and XEDS detector [36] (a), and
symmetrically arranged four XEDS detectors around specimen (b) [35]
3040                         Mao-hua LI, et al/Trans. Nonferrous Met. Soc. China 24(2014) 3031í3050
Fig. 12 Schematic diagram of 4D ET with time resolution (a), a series of 2D images at various projection angles and time steps (b),
snapshots of 4D visualization of nanotubes for two angles at relatively early time (c) and 4D tomography of bracelet at longer time (d)
(White arrows indicate the direction of motion) [39]
Fig. 13 Schematic diagram of missing wedge generated due to restricted tilt range in Fourier space (a) (ș is the sampling angle and Į
is the maximum tilt angle) and missing wedge (gray) reduced to missing pyramid (green) by dual-axis tilt tomography (b) [47,48]
Fig. 14 3D reconstructions of CdTe tetrapods: (a) Reconstructions from same volume using one single-axis tilt series; (b) Another
perpendicular single-axis tilt series; (c) Dual-axis tilt series; (d) Magnification of tetrapod boxed shown in (c) [47]
3042                         Mao-hua LI, et al/Trans. Nonferrous Met. Soc. China 24(2014) 3031í3050
although the missing legs in Fig. 14(a) are present, some             (Figs. 15(b), (d)). Long white arrows indicate fanning
legs from another direction are missing (indicated by the             effects due to the missing wedge in Figs. 15(b) and (d). A
arrows). Figure 14(c) illustrates that no legs are missing.           CNT (indicated by short white arrow) is clearly observed
Figure 14(d) shows the magnification of the tetrapod                  in Fig. 15(a), but is faint in Fig. 15(b). Another short
boxed shown in Fig. 14(c).                                            white arrow (Fig. 15(d)) indicates a fanning effect, which
                                                                      might be misinterpreted as a CNT. The results indicate
3.2 On-axis tilt tomography                                           that on-axis tilt series can improve the fidelity of
     On-axis tilt is a full 360° single-axis tilt. A                  reconstructions. In addition, pillar-shaped specimens can
micropillar sample mounted on a dedicated tomography                  be used to resolve the problem of the sample thickness
holder is used in this technique, as shown in Fig. 11(a)              increasing at higher tilt angles.
[36,51í53]. It allows the sample to rotate over a tilt
range of 360°. The size of specimens manufactured by a                3.3 Conical tilt tomography
FIB technique is approximately 100í300 nm. The                             LANZAVECCHIA et al [54] have reported a
on-axis rotation holder is able to eliminate the missing              conical tilt tomography. Through this method, the
wedge, thus, the fanning and elongation effects will be               missing wedge in Fourier space will be reduced to the
minimized [36,51í53]. KE et al [53] have demonstrated                 missing cone. Figure 16 shows the schematic diagrams
slices of the 3D reconstruction from carbon nanotubes                 of a single-axis tilt tomography, a dual-axis tilt
(CNTs) inside semiconductor contact holes using two                   tomography and a conical tilt tomography in Fourier
different tilt series: ±90° (Figs. 15(a), (c)) and ±70°               space, from which it can be seen that the conical tilt
Fig. 15 Slices of CNTs grown inside semiconductor contact holes at two different positions using full ±90° tilt series (a) and (c), and
±70° tilt series (b) and (d), respectively) [53]
                           Mao-hua LI, et al/Trans. Nonferrous Met. Soc. China 24(2014) 3031í3050                     3043
                                                                high anisotropy. However, the drawbacks of a conical tilt
                                                                tomography are also present, such as increasing the
                                                                number of projections, only partly eliminating the
                                                                missing information and requiring a more elaborate
                                                                specimen stage.
Fig. 17 Orthoslices through 3D reconstruction of porous layer of La2Zr2O7 based on SIRT (top), and DART (bottom) from tilt series
with missing wedge of 0°, 20°, 40°, and 60° [58]
Fig. 18 3D reconstructions of unsupported Ga–Pd nanoparticles by HAADF-STEM tomography: (a) Orthoslices using CS-ET;
(b) Orthoslices using SIRT (b); (aƍ, bƍ) Corresponding 3D reconstructions [63]
                            Mao-hua LI, et al/Trans. Nonferrous Met. Soc. China 24(2014) 3031í3050                            3045
CS-ET       has    several    outstanding    advantages:
compensating for the missing wedge, being segmented                 5 Ongoing research and future perspectives
easily and requiring fewer projections. As a result, the
high fidelity of 3D reconstruction can be obtained based            5.1 New technique assisted ET
on CS-ET.                                                                A novel technology of electron beam precession is
                                                                    able to reduce multiple scattering effects by rocking the
4.3 EST                                                             beam around the optical axis [65,66]. REBLED et al [67]
     EST acquires a tilt series of 2D projections with              have carried out successfully 3D reconstruction of Sn
equal slope increments [55í57]. The technique makes                 precipitates embedded in Al matrix based on beam
use of pseudopolar fast Fourier transform and the                   precession assisted ET, as shown in Fig. 20, in which the
oversampling method with an iterative algorithm. LEE et             left column is without beam precession and the right
al [56] have demonstrated 3D reconstruction of single               column is with beam precession. Obviously, the images
keyhole limpet hemocyanin (KLH) particles by EST for                in right column present more homogenous background
the first time. Figure 19(a) shows iso-surface renderings           and more strong contrast features of the precipitate than
in different orientations of the model. Figures 19(b)í(g)           those in left column. Especially, in Fig. 20(c), the
show 3D reconstruction of KLH particles based on                    precipitate is invisible due to the combination of the
WBP-full, WBP-full-denoising, EST-full, WBP-2/3,                    precipitate inherent diffraction contrast condition and the
WBP-2/3-denoising and EST-2/3, respectively. The                    presence of thickness and bending artifacts, whereas in
WBP reconstructions have a few holes in this region                 Fig. 20(cƍ) the precipitate is clearly visible. These results
(green arrows) in Figs. 19(c) and (f), but the EST                  indicate that the combination of beam precession and
reconstructions are more continuous and smoother in                 BF-TEM can eliminate kinetic effect caused by multiple
Figs. 19(d) and (g). Experimental results indicate that the         scattering, artifacts of thickness and bending. Therefore,
reconstruction by EST has less noise, clear boundaries              high quality reconstruction can be obtained by this
and continuous density. In general, there are some                  technique.
advantages in EST. Firstly, the EST method can
alleviate the missing wedge artifacts. Secondly, it is              5.2 Multi-signal and multi-dimensional electron
suitable for beam-sensitive samples since a limited                      tomography
number of projections are enough to acquire precise                      Transmission electron microscope can generate
reconstruction.                                                     various signals such as elastically scattered electrons,
Fig. 19 Iso-surfaces of known higher resolution model viewed from different directions (a), 3D reconstructions of based on WBP-full
(b), WBP-full-denoising (c), EST-full (d), WBP-2/3 (e), WBP-2/3-denoising (f) and EST-2/3 (g) (2/3 denotes reduced-dose
reconstructions) [56]
3046                          Mao-hua LI, et al/Trans. Nonferrous Met. Soc. China 24(2014) 3031í3050
Fig. 20 BF-TEM images of Sn precipitate at 35° (a, aƍ), +9° (b, bƍ), í9° (c, cƍ) and í43° (d, dƍ) tilt angle [67]
                           Mao-hua LI, et al/Trans. Nonferrous Met. Soc. China 24(2014) 3031í3050                                    3047
inelastically scattered electrons, and X-ray. ET can
reconstruct materials at nanoscale by utilizing mass
thickness maps [4í8], diffraction maps [9í13], elemental
maps [13í28,32í37], thickness maps [17] and bulk
plasmon maps [17], etc. With the development of the
elemental mapping techniques in TEM, STEM and
UEM, ET can obtain 4D reconstruction through
space-energy maps [14í19, 32í37], energy-time maps
[68í70] and space-time maps [38,39].
           ܜ䖯ⱘ⬉ᄤᮁሖᠿᦣᡔᴃᴤ᭭⾥ᄺЁⱘথሩ
          ——Ѣ䗣ᇘ⬉ᄤᰒᖂ䬰ᠿᦣ䗣ᇘ⬉ᄤᰒᖂ䬰
               हਔ ܟ1ēཷཛྷூ 1ēܻ μ 1ē ບ 1ē ฟ 1ē ۂੜ 1ēథ ذޟ2
       ᨬ㽕˖ҟ㒡њѢ䗣ᇘ⬉ᄤᰒᖂ䬰ᠿᦣ䗣ᇘ⬉ᄤᰒᖂ䬰ⱘ⬉ᄤᮁሖᠿᦣᡔᴃᴤ᭭⾥ᄺ乚ඳⱘ᳔ᮄ䖯ሩDŽ䆺
       㒚ഄᦣ䗄⾡⬉ᄤᮁሖᠿᦣⱘ㉏ൟ˖Ѣ䗣ᇘ⬉ᄤᰒᖂ䬰ⱘᮁሖᠿᦣᡔᴃࣙᣀᯢഎᮁሖᠿᦣǃᱫഎᮁሖᠿᦣǃᔅ
       ᴳᱫഎᮁሖᠿᦣǃ⦃ᔶᱫഎᮁሖᠿᦣ㛑䞣䖛Ⓒᮁሖᠿᦣ˗Ѣᠿᦣ䗣ᇘ⬉ᄤᰒᖂ䬰ⱘᮁሖᠿᦣᡔᴃࣙᣀ催㾦⦃
       ᔶᱫഎᮁሖᠿᦣǃ⦃ᔶᱫഎᮁሖᠿᦣˈ䴲݅Ḑᯢഎᮁሖᠿᦣǃ⬉ᄤ㛑䞣ᤳ༅䈅ᮁሖᠿᦣ X ᇘ㒓㛑䈅ᮁሖᠿᦣDŽ
       䘧њӬ࣪ⱘؒ䕀㋏߫ˈ↨བঠ䕈ؒ䕀ǃৠ䕈ؒ䕀ǃ䫹ᔶؒ䕀ҹঞㄝ᭰⥛ؒ䕀ㄝDŽᘏ㒧њܜ䖯ⱘ䞡ᵘㅫ⊩ࣙᣀ⾏
       ᬷ䗁ҷ䞡ᵘᡔᴃǃय़㓽Ӵᛳㅫ⊩ҹঞㄝ᭰⥛ㅫ⊩DŽ᳔ৢˈᦤߎњ⬉ᄤᮁሖᠿᦣᡔᴃᴤ᭭⾥ᄺЁⱘথሩ䍟DŽ
       ݇䬂䆡˖⬉ᄤᮁሖᠿᦣᡔᴃ˗ᴤ᭭⾥ᄺ˗䗣ᇘ⬉ᄤᰒᖂ䬰˗ᠿᦣ䗣ᇘ⬉ᄤᰒᖂ䬰
                                                                                     (Edited by Xiang-qun LI)