Advanced Electron Microscopy

Jill Banfield

Background

As one of the most powerful atomic-scale imaging methods, high-resolution (HR) transition electron microscopy (TEM) was employed early for mineral analysis and rapidly extended to understand mineral reactions and complex defects (Veblen et al. 1993). In recent years, advances in electron optics, detectors, sample environments as well as data acquisition and analysis methods have enabled the extraction of more information and the study of challenging materials. In particular, the use of cryogenic sample preparation and imaging methods, as well of the development of low-dose imaging techniques, provided opportunities for structural analyses of fully hydrated fine-grained minerals, particularly clays. Our group is collaborating closely with scientists at the National Center for Electron Microscopy in the Molecular Foundry to develop and apply these methods, as well as new directions, to the systems of study in this proposal.

Cryogenic Transmission Electron Microscopy and Tomography

CryoTEM is routinely used to elucidate molecular structures in the life sciences and was the basis of the 2017 Nobel Prize in chemistry. Very fast cryogenic freezing preserves hydrated samples in a near-native state and reduces irradiation damage by decreasing the diffusion length of reactive species generated by inelastic electron scattering (Hattne et al. 2018) making otherwise sensitive samples amenable to real space atomic-scale characterization. The prevailing paradigm in the life sciences is to average hundreds to thousands of images of nominally identical individual molecules as a means to achieve high resolution while circumventing the deleterious effects of electron-induced beam damage on any individual molecule. This paradigm is not applicable for systems in which the objects of interest are unique, which includes most non-biological samples, and is especially true for clay minerals. We have developed an imaging paradigm for nanominerals that takes advantage of the increased dose robustness and native-state preservation of cryo-frozen samples, while enabling the structural resolution to better than 2 Å and tomographic reconstruction in 3D.

Figure 1 Average structure unit from high-resolution cryo-TEM image of an individual 2:1 clay layer in aqueous solution, with atomic model overlain (Whittaker, Banfield et al., in preparation). Seven adjacent unit cells shown for clarity. Teal = octahedral sheet; blue = tetrahedral sheet; purple = sodium; interlayer H2O not shown. Scale bar = 5 Å.

High electron density contrast between clay mineral nanoparticles and the surrounding aqueous medium produces strong amplitude contrast that can yield near atomic-resolution images of the basal planes of single smectite layers and stacking geometries of water-separated layers (Fig. 1). Moreover, cryo-TEM data increasingly appear able to directly locate the positions of cations within hydrated smectite crystalline hydrates (Whittaker, Lammers, Banfield et al. 2019) and even at the outer clay surfaces. Average structural units constructed from hundreds to thousands of unit cells within a single mineral layer reveal details of the electric double layer (EDL) structure with nearly-atomic resolution (Fig. 1).

Nanomineral structures, even those of layered minerals like clays, are inherently three-dimensional. We have classified defect structures from individual clay particles using cryo electron tomography (CET) to show that there are a hierarchy of pore geometries that are not currently accounted for in conceptual models of clay rich media (Whittaker, Banfield et al. 2020). In particular, topological defects such as dislocations, in which a layer terminates within a particle forcing adjacent layers to bend around it, give rise to pores with diameters that are larger than the equilibrium basal spacing in neighboring interlayers, but smaller than what are conventionally considered to be ‘macropores’, i.e., spaces in between particles. The diffusivities of charged species are highly dependent on the pore diameter, which suggests that these pore structures may play a heretofore unanticipated role in solute transport in clay-rich media.

Low-dose HAADF-STEM

High-angle annular dark-field (HAADF) scanning transmission electron microscopy (STEM) (Fig. 2) provides high-resolution imaging with intensity sensitive to atomic number, Zn, where 1<n<2 depending on imaging-condition and material (Hartel et al. 1996). HAADF-STEM has been used to image high-Z element distributions in Cs-bearing phlogopite, Fe-rich celadonite, Fe-rich interstratified minerals of 7-Å berthierine and 14-Å chlorite. Clay minerals are generally very electron-beam sensitive, however. For example, the damage of halloysite starts and completes at an electron dose of ~1950 and ~8800 e2, respectively (Kogure 2016).

Figure 2 Illustration of HAADF- and 4D-STEM imaging mode. ⍺ is the convergence angle of electron (E) probe. β is the inner collection angle of the annular detector. The figure is adapted from Fang et al., 2019.
Figure 3 HAADF-STEM image of amesite along [010]. T and O stand for the tetrahedral and octahedral layer, respectively. The green arrows in mark the oxygen atomic columns. Zhang, Banfield et al., in preparation

Recently, our group imaged montmorillonite along [001] with a total electron dose less than ~100 e2 and achieved a resolution of 1.2 Å in HR-TEM (Whittaker, Banfield et al. 2018) and we applied the low-dose technique in HAADF‒STEM achieving 1-Å resolution with an electron dose of ~6,000 e2.

We used low-dose HAADF‒STEM to study stacking disorder in amesite (a Cr-rich 1:1 layer silicate that is isostructural with lizardite) and for the first time were able not only to distinguish the cations but also resolve the oxygen atoms in the octahedral and tetrahedral sheets. In Fig. 3, the bright dots between the two darker dots marked by green arrows are octahedral cation columns (O), and those darker dots are oxygen columns.

The image quality combined with sensitivity to atomic number has revealed the transformation from amesite to chlorite to be associated with the repartitioning of Cr among octahedral sites.

Capability Development

Electron Multiple Scattering Correction for Tomography Reconstruction

Tomographic methods for the reconstruction of three dimensional objects like cells from two dimensional images are well-established in the life sciences. However, the high amplitude contrast that makes cryoTEM a powerful technique for visualizing hydrated cations in the interlayer of clay particles also limits the applicability of traditional tomographic reconstruction techniques that do not explicitly treat multiple scattering that is common in ‘thick and heavy’ environmental samples. In collaboration with Dr. Colin Ophus at the National Center for Electron Microscopy (NCEM) and Dr. Laura Waller's computational imaging group at UC-Berkeley, we are applying a multiple-scattering tomographic reconstruction algorithm (Ren et al. 2020) to reveal the 3D structure of hydrated clay particles at dimensions smaller than the thickness of an individual layer. With improvements in data processing undertaken in the upcoming proposal period, we expect to achieve lattice resolution (<5 Å) of each of the tens to hundreds of clay mineral layers within a volume of approximately 1 µm x 700 nm x 250 nm. Geometry-dependent microscopic theories of clay layer interactions will then be applied directly to authentic clay mineral structures to determine intermolecular forces that underlie these structures. This will allow for the explicit determination of properties such as permeability, swelling pressure, and shear modulus for which the microscopic origins are currently simulated or inferred from ensemble measurements. This imaging paradigm holds great promise for the atomic scale characterization of particles and interfaces of broad geological significance, including defects, cracks, interfaces, and hierarchical structures.

Information Extraction from Large Tomographic Datasets

The goal of the tomographic imaging project over the three year proposal period is to accurately predict macroscopic properties from microscopic tomographic models. In principle, tomographic data sets contain sufficient multiscale information to accomplish this goal. However, the wealth of information contained within a CET dataset, which is typically tens to hundreds of gigabytes, exceeds current capabilities to analyze it in full. Two new directions arise for expanding the scope of CET data analysis. Initial steps will focus on the automated segmentation of clay mineral layers from CET using deep learning algorithms implemented in Weka (open source) and Dragonfly™ to create multi-scale structural models that can integrate with complimentary characterization techniques such as XRD. Structural models of segmented layers will be used generate three-dimensional pair distance distribution functions that will inform new theories for X-ray scattering and diffraction based on observed microstructures. Establishing a direct link between microscopic features and macroscopic measurements will allow for detailed analysis of clay mineral microstructures using readily available techniques such as XRD without requiring CET to be acquired for each condition of interest. In an orthogonal effort, nanoscale information about the location of isomorphic substitutions, structure of the EDL, effect of curvature and strain on EDL structure, and relative orientation of interacting layers will be determined by in-painting atomic-resolution information using learned features from the data. Average structural units will be calculated using sub-tomogram averaging implemented in IMOD (Mastronarde & Held 2017) and RELION 3.0 (Scheres 2012). These average features will be classified into distinguishable structural motifs and refined iteratively to generate a basis set of atomic-resolution clay mineral unit structures.

Four-Dimensional STEM Ptychography

The dose-efficient technique 4D–STEM ptychography could help to beat the trade-off between SNR and beam damage (Fang et al. 2019). In 4D–STEM, a 2D pixelated detector is used to record the diffraction pattern (Fig. 2), rather than an ADF detector integrating parts of the diffraction pattern and giving one signal at each probe position. One problem of 4D-STEM is that the commercialized 2D pixelated detector has a readout time of 1 ms (1 kHz), which significantly limits its application for clay minerals at the atomic scale because of sample drift caused by sample charging and long readout time. Recently Lawrence Berkeley National Laboratory has developed a 4D-camera working at a very fast readout speed (Ciston et al. 2019), 87 kHz. Using this new 4D-camera, we can record the diffraction pattern at each probe position with a dwell time as short as 12 μs. From the recorded 4D dataset, 2D for scanning probe plus 2D for diffraction pattern, the projected potential of the structure could be reconstructed using the ptychography method, which may be the most effective way to observe light elements like oxygen (Yang et al. 2017). Using this cutting-edge STEM technique, we anticipate resolving all the atomic columns along both [010] and especially [010]

Applications to Proposal Objectives

Electron microscopy methods will have numerous applications through all the proposal topics. Key examples include:

Topic 1 High-resolution Tomography of Nanopores in Biogenic Carbonates


Topic 2 Hydrated Structure of Fresh, Aged and Healed Single Cracks


Topic 3 Defect Stacking Structures of 2:1 and 1:1 Clays


References


Ciston, J., Johnson, I. J., Draney, B. R., Ercius, P., Fong, E., Goldschmidt, A., . . . Denes, P. (2019). The 4D Camera: Very High Speed Electron Counting for 4D-STEM. Microscopy and Microanalysis, 25(S2), 1930-1931.
Fang, S., Wen, Y., Allen, C. S., Ophus, C., Han, G. G. D., Kirkland, A. I., . . . Warner, J. H. (2019). Atomic electrostatic maps of 1D channels in 2D semiconductors using 4D scanning transmission electron microscopy. Nature Communications, 10(1), 1127.
Hartel, P., Rose, H., & Dinges, C. (1996). Conditions and reasons for incoherent imaging in STEM. Ultramicroscopy, 63(2), 93-114
Hattne, J., Shi, D., Glynn, C., Zee, C.-T., Gallagher-Jones, M., Martynowycz, M. W., . . . Gonen, T. (2018). Analysis of Global and Site-Specific Radiation Damage in Cryo-EM. Structure, 26(5), 759-766.e754.
Mastronarde, D. N., & Held, S. R. (2017). Automated tilt series alignment and tomographic reconstruction in IMOD. Journal of Structural Biology, 197(2), 102-113.
Ren, D., Ophus, C., Chen, M., & Waller, L. (2020). A multiple scattering algorithm for three dimensional phase contrast atomic electron tomography. Ultramicroscopy, 208, 112860.
Scheres, S. H. W. (2012). RELION: Implementation of a Bayesian approach to cryo-EM structure determination. Journal of Structural Biology, 180(3), 519-530.
Veblen, D. R., Banfield, J. F., Guthrie, G. D., Heaney, P. J., Ilton, E. S., Livi, K. J. T., & Smelik, E. A. (1993). High-Resolution and Analytical Transmission Electron Microscopy of Mineral Disorder and Reactions. science, 260(5113), 1465.
Whittaker, M. L., Comolli, L. R., Gilbert, B., & Banfield, J. F. (2020). Layer size polydispersity in hydrated montmorillonite creates multiscale porosity networks. Applied Clay Science, 190, 105548.
Whittaker, M. L., Lammers, L. N., Carrero, S., Gilbert, B., & Banfield, J. F. (2019). Ion exchange selectivity in clay is controlled by nanoscale chemical–mechanical coupling. Proceedings of the National Academy of Sciences, 116(44), 22052.
Yang, H., MacLaren, I., Jones, L., Martinez, G. T., Simson, M., Huth, M., . . . Nellist, P. D. (2017). Electron ptychographic phase imaging of light elements in crystalline materials using Wigner distribution deconvolution. Ultramicroscopy, 180, 173-179.