Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic Monte Carlo”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Unraveling the Early-Stage Ordering of Krypton Solid Bubbles in Molybdenum: A Multimodal Study

Self-organization of defects such as fission gas bubbles in materials can lead to high inventory capacity for fission gas storage and help mitigate swelling caused by fission gases in nuclear fuel materials under radiation in nuclear reactors. Here, we report the physical mechanism of self-organization of krypton (Kr) gas bubbles in molybdenum (Mo) under ion implantation. The ion fluence and temperature-dependent formation of Kr solid bubble superlattice (SBS) in Mo were investigated by using both synchrotron-based small-angle X-ray scattering (SAXS) and transmission electron microscopy (TEM). Early stage self-organization of gas bubbles is observed at a fluence of 2.5 × 10 16 ions/cm 2 at temperatures of 300–400 °C. The bubble lattice constant increases with increasing implantation temperature from 300 to 400 °C. Both experiments and atomic kinetic Monte Carlo modeling indicate that the Kr solid bubbles are weakly ordered in comparison to previously studied helium (He) gas bubble superlattice (GBS) while the lattice constant are relatively smaller for Kr SBS compared to that of He GBS. The irradiation conditions suggest that spinodal decomposition, which is a form of phase separation, probably precedes gas bubble ordering in Mo. Altogether, our work sheds light on the formation mechanism of noble gas superlattice toward the development of radiation-tolerant materials which are important for the design of advanced nuclear reactors.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Revealing Intermetallic Active Sites of PtNi Nanocatalysts for Reverse Water Gas Shift Reaction

Bimetallic nanoparticles have gained great interest due to their potential applications in catalysis. By combining Lattice Boltzmann Modeling (LBM), Density Functional Theory (DFT) calculations, and kinetic Monte Carlo (KMC) simulation, here we rationalize the superior performance of supported Pt-Ni alloy nanostructures over Pt and Ni alone during the reverse Water Gas Shift (rWGS) reaction, which was observed previously. The LBM simulation including the diffusion and flow of the reacting gas molecules within the supporting nanoporous channels in SBA-15 indicates that the simple combination of contributions from Ni and Pt alone is not enough to result in the observed promotion in rWGS activity and selectivity of the Pt-Ni alloy. Instead, the synergy between Pt and Ni via a unique intermetallic motif is responsible according to the DFT and KMC simulation. In conclusion, compared to the monometallic counterparts, the intermetallic-introduced ensemble effect at the low-coordinated Pt-Ni hybrid site of Pt-Ni alloy nanocatalysts enables the enhanced activation of CO 2 and promoted removal of CO/H 2 O as gas phase product.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Evolution and Degradation Patterns of Electrochemical Cells Based on the Analysis of Interfacial Phenomena at Li Metal Anode/Electrolyte Interfaces

In this work, we report the results of a theoretical–computational analysis of the solid electrolyte interphase (SEI) growth and degradation dynamics occurring in lithium metal batteries during cycling. We use ab initio-kinetic Monte Carlo simulations to generate a synthetic data set, which is analyzed by machine learning methods. We aim to determine: (i) how modifications in interfacial interaction energies between solid electrolyte interphase (SEI) blocks and between Li ions and SEI facets impact the Coulombic efficiency (CE) of the battery and (ii) what factors, including reactions, microscopic transport, and other interfacial events, may lead to cell performance “failure” during prolonged charge and discharge cycles, signaled as a sharp decay in the CE over cycling. The demonstration of our approach is done on a cell including a Li metal surface interfacing with a previously introduced state-of-the-art electrolyte, and the idea can be applied to any electrochemical system. Outcomes include the identification of the leading chemical, physical, and structural variables causing cell failure and relating them to the electrolyte formulation, thus paving the way to future more refined analysis and electrolyte design.

batteries↗
DOE OSTI2025

First-Principles Dissociation Pathways of BCl 3 on the Si(100)-2 × 1 Surface

BCl 3 is a promising acceptor precursor for atomic-precision δ-doping of silicon, as it has been observed to rapidly dissociate into boron doped into the silicon surface and surface chlorine, which can be removed upon annealing. The chemical pathway and the resulting kinetics, through which BCl 3 adsorbs and dissociates on silicon, however, have only been partially explained. Here, in this work, we use density functional theory to expand the dissociation reactions of BCl 3 to include reactions that take place across multiple silicon dimer rows and reactions which end in a bare B atom either at the surface, substituted for a surface silicon, or in a subsurface position. We further simulate the resulting scanning tunneling microscopy images for each of these BCl x dissociation fragments, demonstrating that they often display distinct features that may allow for relatively confident experimental identification. Finally, we input the full dissociation pathway for BCl 3 into a kinetic Monte Carlo model, which simulates realistic reaction pathways as a function of environmental conditions, such as the pressure and temperature of dosing. We find that BCl 2 is broadly dominant at low temperatures, while high temperatures and ample space on the silicon surface for dissociation encourage the formation of bridging BCl fragments and B substitutions on the surface. This work provides the chemical mechanisms for understanding atomic-precision doping of Si with B, enabling a number of relevant quantum applications, such as bipolar nanoelectronics, acceptor-based qubits, and superconducting Si.

Campbell, Quinn T. [Sandia National Laboratories (↗
DOE OSTI2023

Theoretical Perspective of Promoting Direct Methane-to-Methanol Conversion at Complex Metal Oxide–Metal Interfaces

Direct methane conversion to methanol has been considered as an effective and economic way to release greenhouse effects and current high demand on methanol in industry. However, the process has long been challenging due to lack of viable catalysts to compromise the activation of methane that typically occurs high temperatures and retaining of produced methanol that requires mild conditions. Here, this perspective demonstrates an effective strategy to promote direct methane to methanol conversion by engineering the active sites and chemical environments at complex metal oxide - copper oxide - copper interfaces. Such effort strongly depends on the extensive theoretical studies by combining density functional theory (DFT) calculations and kinetic Monte Carlo (KMC) simulations to provide in-depth understanding of reaction mechanism and active sites, which built the strong basis to enable the identification of design principles and advance the catalyst optimization for selective CH 4 -to-CH 3 OH conversion.

03 NATURAL GAS↗
DOE OSTI2023

Aggregation of Size-Selected Oxide Clusters Deposited onto Au(111)

Kinetic Monte Carlo (kMC) simulations along with density functional theory (DFT) calculations were used to investigate the aggregation of size-selected Nb 3 O y (y = 5, 6, 7) clusters deposited onto the Au(111) surface. Recent STM experiments showed that the cluster binding sites and sizes of the cluster assemblies on the Nb 3 O y /Au(111) surfaces strongly depend on the stoichiometry of the clusters, i.e., the oxygen-to-niobium ratio. To better understand the origins of these differences, kMC simulations of the nucleation and growth of cluster assemblies were performed using energy barriers for diffusion and intercluster interactions estimated from DFT calculations of cluster binding and dimerization energies, respectively. Comparisons of the kMC simulations with STM images of the as-deposited Nb 3 O y /Au(111) surfaces at RT and after high temperature annealing were used to further optimize the energetics and gauge the importance of nearest neighbor interactions. The kMC simulations demonstrate that the assembly of Nb 3 O y clusters on Au(111) are largely controlled by the magnitude of the barriers for diffusion and interparticle-bond formation, while changes at higher temperatures are sensitive to the binding energies between nearest neighbors. Simulations for the Nb 3 O 5 and Nb 3 O 6 clusters, which exhibit smaller cluster assembly sizes in STM, required larger diffusion barriers as well as different barriers for interparticle binding, which reflected differences in DFT calculated dimerization energies. Furthermore, the results demonstrate the effectiveness of combined DFT and kMC calculations for understanding how the stoichiometry affects the aggregation of small oxide clusters on a metal surface.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Critical Hydrogen Coverage Effect on the Hydrogenation of Ethylene Catalyzed by δ-MoC(001): An Ab Initio Thermodynamic and Kinetic Study

Highly selective oxidation of methane to methanol has long been challenging in catalysis. Here, we reveal key steps for the promotion by water when tuning the selectivity of a well-defined CeO 2 /Cu 2 O/Cu(111) catalyst from CO and CO 2 to methanol under a reaction environment with methane, oxygen, and water. Ambient-pressure x-ray photoelectron spectroscopy showed that water added to CH 4 and O 2 led to surface CH 3 O and accelerated methanol production. These results were consistent with density functional theory calculations and kinetic Monte Carlo simulations which showed that water preferentially dissociates over the active Ce ions at the CeO 2 -Cu 2 O/Cu(111) interface. The adsorbed OH blocked O-O bond cleavage that would dehydrogenate CH 3 O to CO and CO 2 , and it directly converted this species to methanol, while O 2 now reoxidized the reduced surface. Water adsorption also displaced the produced methanol into the gas phase.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Highly Selective Methane to Methanol Conversion on Inverse SnO 2 /Cu 2 O/Cu(111) Catalysts: Unique Properties of SnO 2 Nanostructures and the Inhibition of the Direct Oxidative Combustion of Methane

Direct methane to methanol (CH 4 → CH 3 OH) conversion in heterogeneous catalysis has been a long-standing challenge due to the difficulties in equalizing the activation of methane and protection of the methanol product at the same reaction conditions. Here, we report an inverse catalyst, consisting of small structures of SnO 2 (0.5-1 nm in size) dispersed on Cu 2 O/Cu(111), for highly selective CH 3 OH production from CH 4 . This system was investigated by combining theoretical [density functional theory calculations (DFT), kinetic Monte Carlo simulations (KMC)] and experimental methods [scanning tunneling microscopy (STM), ambient-pressure X-ray photoelectron spectroscopy (AP-XPS)]. The DFT and AP-XPS studies showed that on SnO 2 /Cu 2 O/Cu(111) the conversion of CH 4 by oxygen (O 2 ) preferred complete combustion to carbon dioxide (CO 2 ). The addition of water (H 2 O) enhanced the production of CH 3 OH to nearly 100% selectivity in KMC simulations. This trend was consistent with results of AP-XPS. The presence of water in the reaction environment rendered an extremely high amount of methoxy species (*CH 3 O), a precursor for CH 3 OH production. Further, the high CH 3 OH selectivity of SnO 2 /Cu 2 O/Cu(111) reflected the unique atomic and electronic structure of the supported SnO 2 nanoparticles. As a result, the O 2 adsorption and dissociation, and thus the full combustion of CH 4 to CO 2 , was completely suppressed; while the H 2 O dissociative adsorption was still feasible, providing active hydroxyl species for a truly selective CH 4 to CH 3 OH conversion.

03 NATURAL GAS↗
DOE OSTI2023

Correlating Experimentally Determined Hydrogen Binding Energy with Hydrogen Evolution Activity over Metal Monolayers on Molybdenum Nitride

It is well established that hydrogen binding energy (HBE) is a key descriptor for hydrogen evolution reaction (HER) activity, and such a relationship is a useful tool for searching efficient and cost-effective HER catalysts. However, in almost all cases, the HBE values are obtained from density functional theory (DFT) calculations. Here in this study, temperature programmed desorption (TPD) was used to experimentally determine the HBE values of metal monolayers supported on molybdenum nitride (Mo 2 N), and electrochemical measurements were performed on the same surfaces. Combined DFT and kinetic Monte Carlo (kMC) simulations were used to validate the trend observed with TPD and the electrochemical HER activity. Depositing one monolayer of Pt on Mo 2 N led to similar HBE values seen for bulk Pt, and electrochemical measurements showed that monolayer Pt on Mo 2 N had similar HER activity to that of bulk Pt. Similar studies were also performed for monolayer Pd on Mo 2 N to confirm the correlation. This work also demonstrates that Mo 2 N is a promising support to reduce precious metal loading in HER catalysts.

08 HYDROGEN↗
DOE OSTI2025

From Micro-environments to Macroscopic Effects: How the Alkaline Hydrogen Evolution Reaction Drives Cu Cathodic Corrosion

Cathodic corrosion of copper (Cu) has posed a significant challenge for over a century, impeding various technological progresses such as electrochemical conversion of CO 2 (eCO 2 RR) into fuels and other value-added carbon products. Here, in this study, employing a combined Density Functional Theory (DFT) and kinetic Monte Carlo (kMC) simulation approach, we delve into the atomistic level mechanism driving this phenomenon in Cu. Our hypothesis posits the pivotal role of alkaline hydrogen evolution reaction (HER) in facilitating cathodic corrosion in Cu. We rigorously develop a pH-dependent hydroxide (OH) adsorption mechanism and calculate the equilibrium OH coverage (𝜃 OH ) at varying pH levels, the thermodynamic stability of subsurface oxygen (O sub ), as well as the Cu-vacancy mediated diffusion of subsurface oxygens (O sub ). Through comprehensive analysis, we establish correlation among various microenvironments, including oxygen diffusion in subsurface layers, pH-dependent OH adsorption, and Cu dissolution into the electrolyte as (Cu-OH) complexes. Furthermore, our investigation explores the correlation between surface coordination environment of active sites and cathodic corrosion of Cu. Finally, by integrating DFT-derived thermodynamic data into a kMC model, we successfully predict the formation of experimentally observed corrosion pits on Cu-surfaces. This combined approach not only advances our fundamental understanding of Cu cathodic corrosion but also offers insights crucial for developing effective corrosion mitigation strategies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Methanol Synthesis from CO 2 Hydrogenation over a Potassium-Promoted Cu x O/Cu(111) ( x ≤ 2) Model Surface: Rationalizing the Potential of Potassium in Catalysis

Alkalis have been reported as a promotor in the heterogeneous catalysis, being able to enhance the activity and selectivity of catalysts. The effective utilization of alkalis in catalyst optimization requires the fundamental understanding of the underlying mechanism. In this work, we take a potassium (K)-modified Cu x O/Cu(111) (x ≤ 2) model surface as a case study to rationalize the nature of K during the carbon dioxide hydrogenation using combined density functional theory (DFT) calculation and the kinetic Monte Carlo (KMC) simulation. Our result demonstrates the significant tuning of selectivity from carbon monoxide to methanol on going from Cu(111) to K-modified Cu x O/Cu(111). The deposited K + stabilizes the Cu x O thin film under the reducing condition of carbon dioxide hydrogenation. More importantly, our study reveals that K + acts as an active center for selective tuning in the binding, an accelerator for charge transfer, and a mediator for the electron tunneling. As a result, the K-modified Cu x O/Cu(111) opens a methanediol [H 2 C(OH) 2 ]-mediated formate pathway to facilitate the selective conversion of carbon dioxide to methanol. Our study develops the intrinsic rules of design to tune the catalytic performance using alkali metals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Modeling the Photoelectrochemical Evolution of Lead-Based, Mixed-Halide Perovskites Due to Photosegregation

Lead-based, mixed-halide perovskites such as methylammonium lead iodide-bromide [MAPb(I 1-x Br x ) 3 ] undergo anion photosegregation under illumination. This is observed as low band gap photoluminescence from photogenerated iodine-rich domains due to favorable band offsets that induce carrier funneling into them. Unfortunately, theoretical rationalizations of mixed-halide photosegregation are complicated by biases inherent to photoluminescence-based observations. Recent compositionally-weighted X-ray diffraction (XRD) measurements now reveal broad distributions of photosegregated stoichiometries not captured by existing photosegregation models. To better bridge experiment and theory, we perform kinetic Monte Carlo (KMC) simulations of photosegregation within the context of a band gap-based thermodynamic model, which has previously accounted for numerous experimental observations. Our KMC simulations are modified to consider high carrier density Fermi-Dirac statistics that result from carrier funneling and accumulation within photosegregated I-rich domains. Obtained KMC results reproduce broad x terminal distributions seen experimentally and illustrate how they are characterized by a central, heavily I-enriched stoichiometry. I-rich domain “drifting” during photosegregation rationalizes the long photosegregation timescales seen experimentally with drifting simultaneously producing a wake of variable stoichiometry I-rich inclusions that form the lion’s share of stoichiometric heterogeneities seen in compositionally-weighted XRD measurements. These simulations and accompanying rationalizations further reveal a general criterion for realizing favorable free energies to induce demixing. Central to the criterion is the statistical occupation of low gap inclusions in the parent alloy by excitations. Furthermore, the resulting model thus provides a general framework for conceptualizing mixed-halide perovskite light and temperature sensitivities, mediated by photocarriers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Cesium-Induced Active Sites for C–C Coupling and Ethanol Synthesis from CO 2 Hydrogenation on Cu/ZnO(000$\bar{1}$) Surfaces

The efficient conversion of carbon dioxide, a major air pollutant, into ethanol or higher alcohols is a big challenge in heterogeneous catalysis, generating great interest in both basic scientific research and commercial applications. Here, we report the facilitated methanol synthesis and the enabled ethanol synthesis from carbon dioxide hydrogenation on a catalyst generated by codepositing Cs and Cu on a ZnO(000$\bar{1}$) substrate. A combination of catalytic testing, X-ray photoelectron spectroscopy (XPS) measurements, and calculations based on density functional theory (DFT) and kinetic Monte Carlo (KMC) simulation was used. The results of XPS showed a clear change in the reaction mechanism when going from Cs/Cu(111) to a Cs/Cu/ZnO(000$\bar{1}$) catalyst. The Cs-promoting effect on C–C coupling is a result of a synergy among Cs, Cu, and ZnO components that leads to the presence of CH x and CH y O species on the surface. Furthermore, according to the DFT-based KMC simulations, the deposition of Cs introduces multifunctional sites with a unique structure at the Cu–Cs–ZnO interface, particularly being able to promote the interaction with CO 2 and thus the methanol synthesis predominantly via the formate pathway. More importantly, it tunes the CHO binding strongly enough to facilitate the HCOOH decomposition to CHO via the formate pathway, but weakly enough to allow further hydrogenation to methanol. The fine-tuning of CHO binding also enables a close alignment of a CHO pair to facilitate the C–C coupling and eventually ethanol synthesis. Our study opens new possibilities to allow the highly active and selective conversion of carbon dioxide to higher alcohols on widely used and low-cost Cu-based catalysts.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Selective Methane Oxidation to Methanol on ZnO/Cu 2 O/Cu(111) Catalysts: Multiple Site-Dependent Behaviors

Due to the abundance of natural gas in our planet, a major goal is to achieve a direct methane to methanol conversion at medium to low temperatures using mixtures of methane and oxygen. Here, we report an efficient catalyst, ZnO/Cu 2 O/Cu(111), for this process investigated using a combination of reactor testing, scanning tunneling microscopy, ambient-pressure x-ray photoemission spectroscopy, density functional calculations, and kinetic Monte Carlo simulations. The catalyst is capable of methane activation at room temperature and transforms mixtures of methane and oxygen to methanol at 450 K with a selectivity of ~30%. This performance is not seen for other heterogenous catalysts which usually require the addition of water to enable a significant conversion of methane to methanol. The unique coarse structure of the ZnO islands supported on a Cu 2 O/Cu(111) substrate provides a collection of multiple centers that display different catalytic activity during the reaction. ZnO-Cu 2 O step sites are active centers for methanol synthesis when exposed to CH 4 and O 2 due to an effective O-O bond dissociation, which enables a methane-to-methanol conversion with a reasonable selectivity. Upon addition of water, the defected O-rich ZnO sites, introduced by Zn vacancies, show superior behavior towards methane conversion and enhance the overall methanol selectivity to over 80%. Thus, in this case, the surface sites involved in a direct CH 4 → to CH 3 OH conversion are different from those engaged in methanol formation without water. Furthermore, the identification of the site-dependent behavior of ZnO/Cu 2 O/Cu(111) opens a design strategy for guiding efficient methane reforming with high methanol selectivity.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Activation and Conversion of Methane to Syngas over ZrO 2 /Cu(111) Catalysts near Room Temperature

Enzymatic systems achieve the catalytic conversion of methane at room temperature under mild conditions. Here, in this study, varying thermodynamic and kinetic parameters, we show that the reforming of methane by water (MWR, CH 4 + H 2 O → CO + 3H 2 ) and the water-gas shift reaction (WGS, CO + H 2 O → H 2 + CO 2 ), two essential processes to integrate fossil fuels toward a H2 energy loop, can be achieved on ZrO 2 /Cu(111) catalysts near room temperature. Measurements of ambient-pressure X-ray photoelectron spectroscopy and mass spectrometry, combined with density functional calculations and kinetic Monte Carlo simulations, were used to study the behavior of the inverse oxide/metal catalysts. The superior performance is associated with a unique zirconia-copper interface, where multifunctional sites involving zirconium, oxygen, and copper work coordinatively to dissociate methane and water at 300 K and move forward the MWR and WGS processes.

03 NATURAL GAS↗
DOE OSTI2024

Tailoring Hierarchical Structure and Rare Earth Affinity of Compositionally Identical Polymers via Sequence Control

Macromolecule sequence, structure, and function are inherently intertwined. While well-established relationships exist in proteins, they are more challenging to define for synthetic polymer nanoparticles due to their molecular weight, sequence, and conformational dispersities. Furthermore, to explore the impact of sequence on nanoparticle structure, we synthesized a set of 16 compositionally identical, sequence-controlled polymers with distinct monomer patterning of dimethyl acrylamide and a bioinspired, structure-driving di(phenylalanine) acrylamide (FF). Sequence control was achieved through multiblock polymerizations, yielding unique ensembles of polymer sequences which were simulated by kinetic Monte Carlo simulations. Systematic analysis of the global (tertiary- and quaternary-like) structure in this amphiphilic copolymer series revealed the effect of multiple sequence descriptors: the number of domains, the hydropathy of terminal domains, and the patchiness (density) of FF within a domain, each of which impacted both chain collapse and the distribution of single- and multichain assemblies. Furthermore, both the conformational freedom of chain segments and local-scale, β-sheet-like interactions were sensitive to the patchiness of FF. To connect sequence, structure, and target function, we evaluated an additional series of nine sequence-controlled copolymers as sequestrants for rare earth elements (REEs) by incorporating a functional acrylic acid monomer into select polymer scaffolds. We identified key sequence variables that influence the binding affinity, capacity, and selectivity of the polymers for REEs. Collectively, these results highlight the potential of and boundaries of sequence control via multiblock polymerizations to drive primary sequence ensembles hierarchical structures, and ultimately the functionality of compositionally identical polymeric materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Unveiling the Origin of Morphological Instability in Topologically Complex Electrocatalytic Nanostructures

Coarsening and degradation phenomena in metals have largely focused on thermally driven processes, such as bulk and surface diffusion. However, dramatic coarsening has been reported in high-surface-area, nanometer-sized Pt-based catalysts during potential cycling in an electrolyte at room temperature─a temperature too low for the process to be explained purely by surface mobility values measured in both vacuum and electrolytes (∼10 –22 and ∼10 –18 cm 2 /s, respectively). This morphological evolution must be due to a different mechanism for mass transport that is sensitive to electrochemical conditions (e.g., electrolyte composition, potential limits, and scan rate). However, there have been no notable studies of electrochemically induced coarsening in nanometer-sized electrocatalysts. Here, we unveil the origins of coarsening in an electrolyte through coupled in situ experiments and atomistic kinetic Monte Carlo (kMC) simulations. Our work demonstrates electrochemical coarsening is driven by two concurrent mechanisms that can be explained at the atomistic level: (i) dissolution/redeposition during the reduction of an oxidized species and (ii) rapid surface diffusion of undercoordinated atoms.

gold↗
DOE OSTI2026

Atomic Evolution of Hydrogen Intercalation Wave Dynamics in Palladium Nanocrystals Revealed by Liquid-Phase Transmission Electron Microscopy

Solute-intercalation-induced phase separation creates spatial heterogeneities in host materials, a phenomenon ubiquitous in batteries, hydrogen storage, and other energy devices. Despite many efforts, probing intercalation processes at the atomic scale has been a significant challenge. By utilizing liquid-phase transmission electron microscopy (TEM), we study hydrogen (de)intercalation in palladium nanocrystals as a model system and have achieved unprecedented atomic-resolution imaging of hydrogen intercalation wave dynamics. Our observations reveal that intercalation wave mechanisms, instead of shrinking-core mechanisms, prevail at ambient temperature for palladium nanocubes ranging from ∼60 nm down to ∼10 nm. Systematic image analysis uncovers the atomic evolution of the hydrogen intercalation wave, transitioning from nonplanar and inclined boundaries to those closely aligned with {100} planes. Our kinetic Monte Carlo simulations demonstrate that the observed intercalation wave dynamics correspond to sorption pathways minimizing the lattice mismatch strain at the phase boundary. In conclusion, unveiling the atomic intercalation pathways holds profound implications for engineering intercalation-mediated devices and advancements in energy sciences.

Lee, Daewon [Lawrence Berkeley National Laboratory↗