Modeling Temperature-Dependent Electron Thermal Diffuse Scattering via Machine-Learned Interatomic Potentials and Path-Integral Molecular Dynamics.

Electron thermal diffuse scattering is shown to be sensitive to subtle changes in atomic vibrations and shows promise in assessing lattice dynamics at nanometer resolution. Here, we demonstrate that machine-learned interatomic potentials (MLIPs) and path-integral molecular dynamics can accurately capture the potential energy landscape and lattice dynamics needed to describe electron thermal diffuse scattering. Using SrTiO_{3} as a test bed at cryogenic and room temperatures, we compare electron thermal diffuse scattering simulations using different approximations to incorporate thermal motion. Only when the simulations are based on quantum mechanically accurate MLIPs in combination with path-integral molecular dynamics that include nuclear quantum effects is there excellent agreement with experiments.

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.

Tuning Surface Acidity of Mixed Conducting Electrodes: Recovery of Si-Induced Degradation of Oxygen Exchange Rate and Area Specific Resistance.

Metal oxides are an important class of functional materials, and for many applications, ranging from solid oxide fuel/electrolysis cells, oxygen permeation membranes, and oxygen storage materials to gas sensors (semiconducting and electrolytic) and catalysts, the interaction between the surface and oxygen in the gas phase is central. Ubiquitous Si-impurities are known to impede this interaction, commonly attributed to the formation of glassy blocking layers on the surface. Here, the surface oxygen exchange coefficient (kchem ) is examined for Pr0.1 Ce0.9 O2-δ (PCO), a model mixed ionic electronic conductor, via electrical conductivity relaxation measurements, and the area-specific resistance (ASR) by electrochemical impedance spectroscopy. It is demonstrated that even low silica levels, introduced by infiltration, depress kchem by a factor 4000, while the ASR increases 40-fold and we attribute this to its acidity relative to that of PCO. The ability to fully regenerate the poisoned surface by the subsequent addition of basic Ca- or Li-species is further shown. This ability to not only recover Si-poisoned surfaces by tuning the relative surface acidity of an oxide surface, but subsequently outperform the pre-poisoned response, promises to extend the operating life of materials and devices for which the catalytic oxygen/solid interface reaction is central.

Three-dimensional atomic packing in amorphous solids with liquid-like structure.

Liquids and solids are two fundamental states of matter. However, our understanding of their three-dimensional atomic structure is mostly based on physical models. Here we use atomic electron tomography to experimentally determine the three-dimensional atomic positions of monatomic amorphous solids, namely a Ta thin film and two Pd nanoparticles. We observe that pentagonal bipyramids are the most abundant atomic motifs in these amorphous materials. Instead of forming icosahedra, the majority of pentagonal bipyramids arrange into pentagonal bipyramid networks with medium-range order. Molecular dynamics simulations further reveal that pentagonal bipyramid networks are prevalent in monatomic metallic liquids, which rapidly grow in size and form more icosahedra during the quench from the liquid to the glass state. These results expand our understanding of the atomic structures of amorphous solids and will encourage future studies on amorphous-crystalline phase and glass transitions in non-crystalline materials with three-dimensional atomic resolution.

Capturing 3D atomic defects and phonon localization at the 2D heterostructure interface.

The three-dimensional (3D) local atomic structures and crystal defects at the interfaces of heterostructures control their electronic, magnetic, optical, catalytic, and topological quantum properties but have thus far eluded any direct experimental determination. Here, we use atomic electron tomography to determine the 3D local atomic positions at the interface of a MoS2-WSe2 heterojunction with picometer precision and correlate 3D atomic defects with localized vibrational properties at the epitaxial interface. We observe point defects, bond distortion, and atomic-scale ripples and measure the full 3D strain tensor at the heterointerface. By using the experimental 3D atomic coordinates as direct input to first-principles calculations, we reveal new phonon modes localized at the interface, which are corroborated by spatially resolved electron energy-loss spectroscopy. We expect that this work will pave the way for correlating structure-property relationships of a wide range of heterostructure interfaces at the single-atom level.

Determining the three-dimensional atomic structure of an amorphous solid.

Amorphous solids such as glass, plastics and amorphous thin films are ubiquitous in our daily life and have broad applications ranging from telecommunications to electronics and solar cells1-4. However, owing to the lack of long-range order, the three-dimensional (3D) atomic structure of amorphous solids has so far eluded direct experimental determination5-15. Here we develop an atomic electron tomography reconstruction method to experimentally determine the 3D atomic positions of an amorphous solid. Using a multi-component glass-forming alloy as proof of principle, we quantitatively characterize the short- and medium-range order of the 3D atomic arrangement. We observe that, although the 3D atomic packing of the short-range order is geometrically disordered, some short-range-order structures connect with each other to form crystal-like superclusters and give rise to medium-range order. We identify four types of crystal-like medium-range order-face-centred cubic, hexagonal close-packed, body-centred cubic and simple cubic-coexisting in the amorphous sample, showing translational but not orientational order. These observations provide direct experimental evidence to support the general framework of the efficient cluster packing model for metallic glasses10,12-14,16. We expect that this work will pave the way for the determination of the 3D structure of a wide range of amorphous solids, which could transform our fundamental understanding of non-crystalline materials and related phenomena.

Correlating the three-dimensional atomic defects and electronic properties of two-dimensional transition metal dichalcogenides.

The electronic, optical and chemical properties of two-dimensional transition metal dichalcogenides strongly depend on their three-dimensional atomic structure and crystal defects. Using Re-doped MoS2 as a model system, here we present scanning atomic electron tomography as a method to determine three-dimensional atomic positions as well as positions of crystal defects such as dopants, vacancies and ripples with a precision down to 4 pm. We measure the three-dimensional bond distortion and local strain tensor induced by single dopants. By directly providing these experimental three-dimensional atomic coordinates to density functional theory, we obtain more accurate electronic band structures than derived from conventional density functional theory calculations that relies on relaxed three-dimensional atomic coordinates. We anticipate that scanning atomic electron tomography not only will be generally applicable to determine the three-dimensional atomic coordinates of two-dimensional materials, but also will enable ab initio calculations to better predict the physical, chemical and electronic properties of these materials.

Observing crystal nucleation in four dimensions using atomic electron tomography.

Nucleation plays a critical role in many physical and biological phenomena that range from crystallization, melting and evaporation to the formation of clouds and the initiation of neurodegenerative diseases1-3. However, nucleation is a challenging process to study experimentally, especially in its early stages, when several atoms or molecules start to form a new phase from a parent phase. A number of experimental and computational methods have been used to investigate nucleation processes4-17, but experimental determination of the three-dimensional atomic structure and the dynamics of early-stage nuclei has been unachievable. Here we use atomic electron tomography to study early-stage nucleation in four dimensions (that is, including time) at atomic resolution. Using FePt nanoparticles as a model system, we find that early-stage nuclei are irregularly shaped, each has a core of one to a few atoms with the maximum order parameter, and the order parameter gradient points from the core to the boundary of the nucleus. We capture the structure and dynamics of the same nuclei undergoing growth, fluctuation, dissolution, merging and/or division, which are regulated by the order parameter distribution and its gradient. These experimental observations are corroborated by molecular dynamics simulations of heterogeneous and homogeneous nucleation in liquid-solid phase transitions of Pt. Our experimental and molecular dynamics results indicate that a theory beyond classical nucleation theory1,2,18 is needed to describe early-stage nucleation at the atomic scale. We anticipate that the reported approach will open the door to the study of many fundamental problems in materials science, nanoscience, condensed matter physics and chemistry, such as phase transition, atomic diffusion, grain boundary dynamics, interface motion, defect dynamics and surface reconstruction with four-dimensional atomic resolution.