Efficient and Versatile Model for Vibrational STEM-EELS

Systematic Absences of Optical Phonon Modes in Phonon Dispersion Measured by Electron Microscopy

Phonon dispersion relations are widely used to elucidate the vibrational properties of materials. As an emerging technique, momentum-resolved vibrational spectroscopy in scanning transmission electron microscopy offers an unparalleled approach to explore q-dependent phonon behavior at local structures. In this study, we systematically investigate the phonon dispersion of monolayer graphene across several Brillouin zones (BZs) using momentum-resolved vibrational spectroscopy and find that the optical phonon signals vanish at the Γ points with indices (hk0) satisfying h+2k=3n (n denoted integers). Theoretical analysis reveals that the observed phenomena arise from the complete destructive interference of the scattered waves from different basis atoms. This observation, corroborated by the study of diamond, should be a general characteristic of materials composed of symmetrically equivalent pairs of the same elements. Moreover, our results emphasize the importance of multiple scattering in interpreting the vibrational signals in bulk materials. We demonstrate that the systematic absences and dynamic effects, which have not been much appreciated before, offer new insights into the experimental assessment of local vibrational properties of materials.

Footprints of atomic-scale features in plasmonic nanoparticles as revealed by electron energy loss spectroscopy

We present a first-principles theoretical study of the atomistic footprints in the valence electron energy loss spectroscopy (EELS) of nanometer-size metallic particles. Charge density maps of excited plasmons and EEL spectra for specific electron paths through a nanoparticle (Na380 atom cluster) are modeled using ab initio calculations within time-dependent density functional theory. Our findings unveil the atomic-scale sensitivity of EELS within this low-energy spectral range. Whereas localized surface plasmons (LSPs) are particularly sensitive to the atomistic structure of the surface probed by the electron beam, confined bulk plasmons (CBPs) reveal quantum size effects within the nanoparticle's volume. Moreover, we prove that classical local dielectric theories mimicking the atomistic structure of the nanoparticles reproduce the LSP trends observed in quantum calculations, but fall short in describing the CBP behavior observed under different electron trajectories.

Nanoscale Localized Phonons at Al2O3 Grain Boundaries

Nanoscale defects like grain boundaries (GBs) would introduce local phonon modes and affect the bulk materials' thermal, electrical, optical, and mechanical properties. It is highly desirable to correlate the phonon modes and atomic arrangements for individual defects to precisely understand the structure-property relation. Here we investigated the localized phonon modes of Al2O3 GBs by combination of the vibrational electron energy loss spectroscopy (EELS) in scanning transmission electron microscope and density functional perturbation theory (DFPT). The differences between GB and bulk obtained from the vibrational EELS show that the GB exhibited more active vibration at the energy range of <50 meV and >80 meV, and further DFPT results proved the wide distribution of bond lengths at GB are the main factor for the emergence of local phonon modes. This research provides insights into the phonon-defect relation and would be of importance in the design and application of polycrystalline materials.

Direct Visualization of Localized Vibrations at Complex Grain Boundaries

Grain boundaries (GBs) are a prolific microstructural feature that dominates the functionality of a wide class of materials. The functionality at a GB results from the unique atomic arrangements, different from those in the grain, that have driven extensive experimental and theoretical studies correlating atomic-scale GB structures to macroscopic electronic, infrared optical, and thermal properties. In this work, a SrTiO₃ GB is examined using atomic-resolution aberration-corrected scanning transmission electron microscopy and ultrahigh-energy-resolution monochromated electron energy-loss spectroscopy, in conjunction with density functional theory. This combination enables the correlation of the GB structure, nonstoichiometry, and chemical bonding with a redistribution of vibrational states within the GB dislocation cores. The new experimental access to localized GB vibrations provides a direct route to quantifying the impact of individual boundaries on macroscopic properties.

A comparison of molecular dynamics potentials used to account for thermal diffuse scattering in multislice simulations

Here we investigate electron scattering simulations with thermal displacements incorporated using molecular dynamics potentials. Specifically, we explore the sensitivity of electron scattering to the phonon band structure, or more explicitly interatomic forces. Silicon serves as the model material where we introduce thermal atomic displacements via empirical and machine-learned molecular dynamics interatomic potentials and compare them to finite-temperature density functional theory interatomic forces. We demonstrate that when molecular dynamics potentials do not sufficiently reproduce the correct phonon band structure, significant errors in the simulated diffraction and image intensities can occur. Moreover, for Si, we find that multislice simulations using machine-learned interatomic potentials are more accurate than empirical ones. In addition to the selected atomic potential, we demonstrate that the sensitivity to the phonon band structure also depends on the crystal zone axis, which can be used to enhance sensitivity to thermal displacements. Finally, we provide a sensitivity analysis with angle-resolved scanning transmission electron microscopy (STEM) to enhance image sensitivity to the details of the phonon band structure.

Localized Phonon Densities of States at Grain Boundaries in Silicon

Since it is now possible to record vibrational spectra at nanometer scales in the electron microscope, it is of interest to explore whether extended defects in crystals such as dislocations or grain boundaries will result in measurable changes of the phonon densities of states (dos) that are reflected in the spectra. Phonon densities of states were calculated for a set of high angle grain boundaries in silicon. The boundaries are modeled by supercells with up to 160 atoms, and the vibrational densities of states were calculated by taking the Fourier transform of the velocity–velocity autocorrelation function from molecular dynamics simulations with larger supercells doubled in all three directions. In selected cases, the results were checked on the original supercells by comparison with the densities of states obtained by diagonalizing the dynamical matrix calculated using density functional theory. Near the core of the grain boundary, the height of the optic phonon peak in the dos at 60 meV was suppressed relative to features due to acoustic phonons that are largely unchanged relative to their bulk values. This can be attributed to the variation in the strength of bonds in grain boundary core regions where there is a range of bond lengths.

Theory of Atomic-Scale Vibrational Mapping and Isotope Identification with Electron Beams

Transmission electron microscopy and spectroscopy currently enable the acquisition of spatially resolved spectral information from a specimen by focusing electron beams down to a sub-angstrom spot and then analyzing the energy of the inelastically scattered electrons with few-meV energy resolution. This technique has recently been used to experimentally resolve vibrational modes in 2D materials emerging at mid-infrared frequencies. Here, on the basis of first-principles theory, we demonstrate the possibility of identifying single isotope atom impurities in a nanostructure through the trace that they leave in the spectral and spatial characteristics of the vibrational modes. Specifically, we examine a hexagonal boron nitride molecule as an example of application, in which the presence of a single isotope impurity is revealed through changes in the electron spectra, as well as in the space-, energy-, and momentum-resolved inelastic electron signal. We compare these results with conventional far-field spectroscopy, showing that electron beams offer superior spatial resolution combined with the ability to probe the complete set of vibrational modes, including those that are optically dark. Our study is relevant for the atomic-scale characterization of vibrational modes in materials of interest, including a detailed mapping of isotope distributions.

The abTEM code: transmission electron microscopy from first principles

Simulation of transmission electron microscopy (TEM) images or diffraction patterns is often required to interpret experimental data. Since nuclear cores dominate electron scattering, the scattering potential is typically described using the independent atom model, which completely neglects valence bonding and its effect on the transmitting electrons. As instrumentation has advanced, new measurements have revealed subtle details of the scattering potential that were previously not accessible to experiment. We have created an open-source simulation code designed to meet these demands by integrating the ability to calculate the potential via density functional theory (DFT) with a flexible modular software design. abTEM can simulate most standard imaging modes and incorporates the latest algorithmic developments. The development of new techniques requires a program that is accessible to domain experts without extensive programming experience. abTEM is written purely in Python and designed for easy modification and extension. The effective use of modern open-source libraries makes the performance of abTEM highly competitive with existing optimized codes on both CPUs and GPUs and allows us to leverage an extensive ecosystem of libraries, such as the Atomic Simulation Environment and the DFT code GPAW. abTEM is designed to work in an interactive Python notebook, creating a seamless and reproducible workflow from defining an atomic structure, calculating molecular dynamics (MD) and electrostatic potentials, to the analysis of results, all in a single, easy-to-read document.  This article provides ongoing documentation of abTEM development. In this first version, we show use cases for hexagonal boron nitride, where valence bonding can be detected, a 4D-STEM simulation of molybdenum disulfide including ptychographic phase reconstruction, a comparison of MD and frozen phonon modeling for convergent-beam electron diffraction of a 2.6-million-atom silicon system, and a performance comparison of our fast implementation of the PRISM algorithm for a decahedral 20000-atom gold nanoparticle.

The abTEM code: transmission electron microscopy from first principles

Simulation of transmission electron microscopy (TEM) images or diffraction patterns is often required to interpret experimental data. Since nuclear cores dominate electron scattering, the scattering potential is typically described using the independent atom model, which completely neglects valence bonding and its effect on the transmitting electrons. As instrumentation has advanced, new measurements have revealed subtle details of the scattering potential that were previously not accessible to experiment. We have created an open-source simulation code designed to meet these demands by integrating the ability to calculate the potential via density functional theory (DFT) with a flexible modular software design. abTEM can simulate most standard imaging modes and incorporates the latest algorithmic developments. The development of new techniques requires a program that is accessible to domain experts without extensive programming experience. abTEM is written purely in Python and designed for easy modification and extension. The effective use of modern open-source libraries makes the performance of abTEM highly competitive with existing optimized codes on both CPUs and GPUs and allows us to leverage an extensive ecosystem of libraries, such as the Atomic Simulation Environment and the DFT code GPAW. abTEM is designed to work in an interactive Python notebook, creating a seamless and reproducible workflow from defining an atomic structure, calculating molecular dynamics (MD) and electrostatic potentials, to the analysis of results, all in a single, easy-to-read document.  This article provides ongoing documentation of abTEM development. In this first version, we show use cases for hexagonal boron nitride, where valence bonding can be detected, a 4D-STEM simulation of molybdenum disulfide including ptychographic phase reconstruction, a comparison of MD and frozen phonon modeling for convergent-beam electron diffraction of a 2.6-million-atom silicon system, and a performance comparison of our fast implementation of the PRISM algorithm for a decahedral 20000-atom gold nanoparticle.

ab initio description of bonding for transmission electron microscopy

The simulation of transmission electron microscopy (TEM) images or diffraction patterns is often required to interpret their contrast and extract specimen features. This is especially true for high-resolution phase-contrast imaging of materials, but electron scattering simulations based on atomistic models are widely used in materials science and structural biology. Since electron scattering is dominated by the nuclear cores, the scattering potential is typically described by the widely applied independent atom model. This approximation is fast and fairly accurate, especially for scanning TEM (STEM) annular dark-field contrast, but it completely neglects valence bonding and its effect on the transmitting electrons. However, an emerging trend in electron microscopy is to use new instrumentation and methods to extract the maximum amount of information from each electron. This is evident in the increasing popularity of techniques such as 4D-STEM combined with ptychography in materials science, and cryogenic microcrystal electron diffraction in structural biology, where subtle differences in the scattering potential may be both measurable and contain additional insights. Thus, there is increasing interest in electron scattering simulations based on electrostatic potentials obtained from first principles, mainly via density functional theory, which was previously mainly required for holography. In this Review, we discuss the motivation and basis for these developments, survey the pioneering work that has been published thus far, and give our outlook for the future. We argue that a physically better justified ab initio description of the scattering potential is both useful and viable for an increasing number of systems, and we expect such simulations to steadily gain in popularity and importance.

Four-dimensional vibrational spectroscopy for nanoscale mapping of phonon dispersion in BN nanotubes

Directly mapping local phonon dispersion in individual nanostructures can advance our understanding of their thermal, optical, and mechanical properties. However, this requires high detection sensitivity and combined spatial, energy and momentum resolutions, thus has been elusive. Here, we demonstrate a four-dimensional electron energy loss spectroscopy technique, and present position-dependent phonon dispersion measurements in individual boron nitride nanotubes. By scanning the electron beam in real space while monitoring both the energy loss and the momentum transfer, we are able to reveal position- and momentum-dependent lattice vibrations at nanometer scale. Our measurements show that the phonon dispersion of multi-walled nanotubes is locally close to hexagonal-boron nitride crystals. Interestingly, acoustic phonons are sensitive to defect scattering, while optical modes are insensitive to small voids. This work not only provides insights into vibrational properties of boron nitride nanotubes, but also demonstrates potential of the developed technique in nanoscale phonon dispersion measurements.

Lattice resolution of vibrational modes in the electron microscope

The combination of aberration correction and ultra high energy resolution with monochromators has made it possible to record images showing lattice resolution in phonon modes, both with a displaced collection aperture and more recently with an on -axis collection aperture. In practice the objective aperture has to include Bragg reflections that correspond to the observed lattice image spacings, and the specimen has to be sufficiently thick for adequate phonon scattered intensity. There has been controversy as to whether the images with the on axis detector are really a consequence of lattice resolution in a phonon mode or just a transfer of information from an image that was formed by elastically scattered electrons. We present results of calculations based on a theory that includes the possibility of dynamical electron diffraction for both incident and scattered electrons and the full phonon dispersion relation. We show that Umklapp scattering from the second Brillouin Zone back to the first Brillouin Zone is necessary for lattice resolution with the on axis detector and that it is therefore reasonable to attribute the lattice resolution to the phonon scattering.