Search NASA⌕ Search

SEARCH · Search NASA

Results for “density functional theory (DFT) calculations”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Towards Quantum-Chemical Level Calculations of SARS-CoV-2 Spike Protein Variants of Concern by First Principles Density Functional Theory

The spike protein (S-protein) is a crucial part of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), with its many domains responsible for binding, fusion, and host cell entry. In this review we use the density functional theory (DFT) calculations to analyze the atomic-scale interactions and investigate the consequences of mutations in S-protein domains. We specifically describe the key amino acids and functions of each domain, which are essential for structural stability as well as recognition and fusion processes with the host cell; in addition, we speculate on how mutations affect these properties. Such unprecedented large-scale ab initio calculations, with up to 5000 atoms in the system, are based on the novel concept of amino acid–amino acid-bond pair unit (AABPU) that allows for an alternative description of proteins, providing valuable information on partial charge, interatomic bonding and hydrogen bond (HB) formation. In general, our results show that the S-protein mutations for different variants foster an increased positive partial charge, alter the interatomic interactions, and disrupt the HB networks. We conclude by outlining a roadmap for future computational research of biomolecular virus-related systems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Informing air–carbon ablation modeling with theoretical calculations of atomic oxygen and nitrogen interacting with carbon surfaces

To understand the gas–surface chemistry above the thermal protection system of a hypersonic vehicle, it is necessary to map out the kinetics of key elementary reaction steps. In this work, extensive periodic density functional theory (DFT) calculations are performed to elucidate the interaction of atomic oxygen and nitrogen with both the basal plane and edge sites of highly oriented pyrolytic graphite (HOPG). Reaction energies and barriers are determined for adsorption, desorption, diffusion, recombination, and several reactions. These DFT results are compared with the most recent finite-rate model for air–carbon ablation. Our DFT results corroborated some of the parameters used in the model but suggest that further refinement may be necessary for others. The calculations reported here will help to establish a predictive kinetic model for the complex reaction network present under hypersonic flight conditions.

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

Origin of Stabilization of Ligand-Centered Mixed Valence Ruthenium Azopyridine Complexes: DFT Insights for Neuromorphic Applications

Redox-driven conductance changes are critical processes in molecular- and coordination-complex-based memristive thin films and devices that are envisioned for neuromorphic technologies, but fundamental mechanisms of conductance switching are not fully understood. Here, we explore charge disproportionation (CD) processes in [Ru II L 2 ](PF 6 ) 2 molecular systems that intrinsically involve interfragment charge transfer (IFCT). Using a combination of ab initio molecular dynamics simulation (AIMD), time-dependent density functional theory (TD-DFT), and density functional theory (DFT) calculations, we investigate the electron transfer mechanisms and the roles of temperature and cell volumetric expansion in facilitating the counterion movements and electronic transitions required for low-cost IFCT and charge redistribution. A detailed analysis of the density of states and TD-DFT calculations highlights that unpaired electrons play a crucial role in low-energy transitions, with the azo (N=N) groups of the ligand serving as the primary sites for electronic transport between molecular fragments, further stabilizing the asymmetric state. Localization of added electrons on azo ligands occurs with negligible change at the Ru centers, supported by atomic volume expansions up to +4.74 bohr 3 , and goes along with a progressive reduction of the HOMO−LUMO gap across redox states, suggesting enhanced conductivity. The TD-DFT analysis reveals a dominant IFCT excitation at 2082.76 nm in the doubly reduced (22) state, while a stabilization energy of 1.20 eV of the asymmetric (13) state relative to the symmetric (22) state is predicted by constrained DFT. Periodic DFT and AIMD simulations emulating a molecular film show that the stabilization of the asymmetric state, relative to a symmetric one, translates in net charge separation values (order of ∼0.33 e) that are strongly linked to increased counterion mobility (average counterion displacements exceeding 0.7 Å per atom during CD events) and the involvement of azo groups in electron redistribution. These findings, which align with previously reported experimental and computational data, provide key insights into the IFCT mechanisms and electronic transport facilitated by azo groups, with important implications for redox-driven memristive and neuromorphic technologies.

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

Electronic structure of compressively strained thin film La$_2$PrNi$_2$O$_7$

The discovery of superconductivity in the bulk nickelates under high pressure is a major advance in physics. The recent observation of superconductivity at ambient pressure in compressively strained bilayer nickelate thin films has now enabled direct characterization of the superconducting phase through angle resolved photoemission spectroscopy (ARPES). Here we present an in-situ ARPES study of compressively strained La$_2$PrNi$_2$O$_7$ films grown by oxide molecular beam epitaxy, and the ozone treated counterparts with an onset T$_c$ of 40 K, supplemented with results from pulsed laser deposition films with similar T$_c$. We resolve a systematic strain-driven electronic band shift with respect to that of bulk crystals, in qualitative agreement with density functional theory (DFT) calculations. However, the strongly renormalized flat 3$d_{z2}$ band shifts a factor of 5-10 smaller than anticipated by DFT. Furthermore, it stays ~70 meV below the Fermi level, contradicting the expectation that superconductivity results from the high density of states of this band at the Fermi level. We also observed a non-trivial k$_z$ dispersion of the cuprate-like 3$d_{x2-y2}$ band. Combined with results from both X-ray diffraction and DFT, we suggest that the strained films are under ~5 GPa effective pressure, considerably larger than the naïve expectation from the DFT relaxed structure. Finally, the ~70 meV energy position is intriguingly close to the collective mode coupling more prominently seen in thin films, in the energy range of both oxygen related phonons and the maximum of the spin excitation spectrum.

FOS: Physical sciences↗
DOE OSTI2022

Defective Sn-Zn perovskites through bio-directed routes for modulating CO 2 RR

Sn and Zn oxides can catalyze electrochemical CO 2 reduction reactions (CO 2 RR) to produce useful chemicals. Herein, we report on perovskite-type tin-zinc oxides (TZO) synthesized by a biomineralization approach under benign aqueous conditions. The resulting TZO catalysts are disordered and possess a high quantity of defects, which serve as highly reactive active sites in CO 2 RR. Selectively toward formate or CO is tunable toward the relative amounts of Sn and Zn and their associated O defects. In-situ X-ray absorption spectroscopy measurements reveal electronic and structural transformations consistent with oxygen vacancy formation at Sn-O-Zn bridges as a function of cathodic potential that strongly correlates to CO2RR selectivity. Further, Density functional theory (DFT) calculations demonstrate that introducing an oxygen vacancy in TZO influences drastically influences formate/CO selectivity, consistent with our in-situ XAS findings and catalytic performance. Overall, these findings demonstrate the potential of biomineralization to fabricate perovskite-type metal oxides for energy conversion reactions.

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

Atomic-Level Structure of Mesoporous Hexagonal Boron Nitride Determined by High-Resolution Solid-State Multinuclear Magnetic Resonance Spectroscopy and Density Functional Theory Calculations

Mesoporous hexagonal boron nitride (p-BN) has received significant attention over the last decade as a promising candidate for water cleaning/pollutant removal and hydrogen storage applications. In this work, high-resolution solid-state NMR spectroscopy and plane-wave density-functional theory (DFT) calculations are used to obtain an atomic-level description of p-BN. 1 H– 15 N or 1 H– 14 N heteronuclear (HETCOR) correlation experiments recorded with either conventional NMR at room temperature or dynamic nuclear polarization surface-enhanced spectroscopy (DNP-SENS) at ca. 100 K reveal NB 2 H, NBH 2 , NBH 3 + species residing on the edges of BN sheets. Ultra-high field 35.2 T 11 B NMR spectroscopy was used to resolve 11 B NMR signals from BN 3 , BN 2 O x (OH) 1–x (x = 0–1), BNO x (OH) 2–x (x = 0–2), BO x (OH) 3–x (x = 0–3), and BO x (OH) 4–x – (x = 0–4). Importantly, 2D 11 B dipolar double-quantum–single-quantum homonuclear correlation spectra reveal that many pore/defect sites are composed of boron oxide/hydroxide clusters connected to the BN framework through BN 2 O units. 1D and 2D 11 B{ 15 N} HETCOR NMR experiments, in addition to plane-wave DFT calculations of nine different structural models, further confirm the assignment of all NMR signals. The detailed structure determination of the pore and edge/defect sites within p-BN should further enable the rational design and development of next-generation p-BN-based materials. In addition, the techniques outlined here should be applicable to determine structure within other porous and/or boron-based materials.

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

Novel ceramic capacitors with ultrahigh energy density and efficiency (Final Technical Report)

Antiferroelectric ceramics are a special class of material that have shown great potential as the dielectric in electrical capacitors due to their high energy- and power-density. During each charge-discharge cycle, the ceramic undergoes transformation to a ferroelectric phase and resumes its antiferroelectric phase. The hysteresis associated with the transitions leads to a mediocre energy efficiency and service lifetime of antiferroelectric capacitors and, hence, their almost absence in commercial products. Under the support of this research project, we first formulated a universal lattice-compatibility theory that included electrostatic polarization energy along with elastic energy and thermal energy to understand the origin of the hysteresis in antiferroelectric oxides. Guided by this compatibility theory, we conducted high-throughput density functional theory (DFT) calculations to assess chemical modifiers and their effect on crystal structures of 400+ PbZrO 3 -based compositions. Down-selected compositions were experimentally validated for their suppressed hysteresis and higher energy efficiency. The verified low-hysteresis compositions were then expanded to an antiferroelectric ceramic library with nearly 500 new compositions (more than 1,500 samples) using high-throughput experiments involving ceramic synthesis and property screening. The large quantity of data generated (theory and experimental) in these tasks were processed by machine-learning techniques and identified trends were fed to the next iteration. In the end, we successfully discovered four compositions with near-zero hysteresis, yielding a world-record energy efficiency of 98.2% at an energy density of 3.0 J/cm 3 . Furthermore, our antiferroelectric ceramic capacitor reaches 79.5 million charge-discharge cycles lifetime, a factor of 80 enhancement over previous antiferroelectric ceramics with large hysteresis. These research accomplishments have not only met the milestones set in the SOPO, but also led to two patent filings, three journal publications (one of them was in Advanced Materials, impact factor 29.4), and nine oral presentations at various venues. Through the course of the project, three postdocs, four Ph.D. students, and one M.S. student were trained. In short, our project established a new methodology in searching next-generation functional ceramics on the fundamental side and discovered several high-efficiency antiferroelectric compositions for capacitors on the applied side. Once fabricated into the multilayer form for commercial applications, these ceramic capacitors can potentially enable the high temperature high power density DC-link capacitors that are critical for the next generation inverters in electric vehicles. The project also significantly contributed to the nation’s workforce development in the STEM fields.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Theoretical Prediction and Interpretation of 237 Np Mössbauer Isomer Shifts

Different theoretical approaches for the calculation of 237 Np Mössbauer isomer shifts are investigated. The traditional contact density route is compared to a previously proposed alternative approach that uses energy derivatives with respect to the nuclear radius. Both approaches yield similar results as long as suitable basis sets augmented with large exponents and relativistic methods are used. Density functional theory (DFT) calculations do not show a strong dependency of the 237 Np isomer shift on the chosen functional. Wavefunction calculations show that dynamic electron correlation can be important when covalent bonding influences the isomer shift. Effects from spin–orbit coupling are small. The isomer shifts of ionic solids and Np(III) organometallic complexes are largely governed by the oxidation state of Np. Isomer shifts of organometallic Np(IV) complexes are strongly affected by donation bonding. Furthermore, detailed analysis of the wavefunction results with different active spaces demonstrates that correlation among the outer core Np and occupied ligand frontier orbitals contributes significantly to isomer shifts of Np(IV) compounds.

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

Bulk Fermi surfaces of the Dirac type-II semimetallic candidate NiTe 2

Here, we present a study on the Fermi surface of the Dirac type-II semimetallic candidate NiTe 2 via the temperature and angular dependence of the de Haas–van Alphen (dHvA) effect measured in single crystals grown through Te flux. In contrast to its isostructural compounds like PtSe 2 , band-structure calculations predict NiTe 2 to display a tilted Dirac node very close to its Fermi level that is located along the Γ to A high-symmetry direction within its first Brillouin zone. The angular dependence of the dHvA frequencies is found to be in agreement with the first-principles calculations when the electronic bands are slightly shifted with respect to the Fermi level (ε F ), and therefore provide support for the existence of a Dirac type-II node in NiTe 2 . Nevertheless, we observed mild disagreements between experimental observations and density functional theory (DFT) calculations as, for example, nearly isotropic and light experimental effective masses. This indicates that the dispersion of the bands is not well captured by DFT. Despite the coexistence of Dirac-like fermions with topologically trivial carriers, samples of the highest quality display an anomalous and large either linear or sublinear magnetoresistivity. Finally, this suggests that Lorentz invariance breaking Dirac-like quasiparticles dominate the carrier transport in this compound.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Can a deep-learning model make fast predictions of vacancy formation in diverse materials?

The presence of point defects, such as vacancies, plays an important role in materials design. Here, we explore the extrapolative power of a graph neural network (GNN) to predict vacancy formation energies. We show that a model trained only on perfect materials can also be used to predict vacancy formation energies (E vac ) of defect structures without the need for additional training data. Such GNN-based predictions are considerably faster than density functional theory (DFT) calculations and show potential as a quick pre-screening tool for defect systems. To test this strategy, we developed a DFT dataset of 530 E vac consisting of 3D elemental solids, alloys, oxides, semiconductors, and 2D monolayer materials. We analyzed and discussed the applicability of such direct and fast predictions. We applied the model to predict 192 494 E vac for 55 723 materials in the JARVIS-DFT database. Our work demonstrates how a GNN-model performs on unseen data.

2D materials↗
DOE OSTI2021

Ab Initio Study of Atomic and Electronic Structure of Promising Ba4XMn3O12 (X = Nb, Ce, Pr) Oxides for Solar Thermochemical Hydrogen Production

The two-step metal oxide water-splitting cycle is one of the most viable approach for Solar Thermochemical Hydrogen (STCH) production. Challenges exist in finding suitable oxides that can satisfy thermodynamics of the STCH redox cycle under viable range of temperatures and partial pressures. Recently, BXM family members, i.e., Ba4NbMn3O12 (BNM), Ba4CeMn3O12 (BCM), Ba4PrMn3O12 (BPM) are discovered to exhibit promising STCH performance. However, the magnetic degrees of freedom of their experimental crystal structures is not characterized, and only hypothetical models exist for their electronic structure. We performed Monte-Carlo sampling for the magnetic spin arrangement of Mn atoms to provide explicit atomic structure models for their ground state 12R polytype and determined the most stable spin configuration for each member after structure relaxation using density functional theory (DFT) calculations. We also elucidate the most stable charge configuration (Mn3+/Mn4+ arrangement) for BNM among the existing hypothesis in literature. We investigate these strongly correlated oxides using different levels of theory (DFT+U, Hybrid-DFT) to develop better understanding of their electronic structure. These calculations allowed us to develop a model based on molecular orbital interactions that help explain the observed changes in their experimental x-ray spectroscopy measurements. Thereby, help link differences in electronic structure of BXM family members to their varying water-splitting behavior. These structure models will further facilitate future in-depth defect studies to model STCH redox cycle of these oxides.

atomic and magnetic structure↗
DOE OSTI2021

Suppression of antiferromagnetic order and strong ferromagnetic spin fluctuations in Ca(Co 1-x Ni x ) 2-y As 2 single crystals

CaCo 2–y As 2 is a unique itinerant system having strong magnetic frustration. Here, we report the effect of electron doping on the physical properties resulting from Ni substitutions for Co. The single crystals of Ca(Co 1–x Ni x ) 2–y As 2 were characterized by single-crystal x-ray diffraction, energy-dispersive x-ray spectroscopy, magnetization M versus temperature T, magnetic field H, time t, and heat capacity C p (H, T) measurements. The A-type antiferromagnetic (AFM) transition temperature T N = 52 K for x = 0 decreases to 22 K with only 3% Ni substitution and is completely suppressed for x > 0.16. For 0.11 ≤ x ≤ 0.52 strong ferromagnetic (FM) fluctuations develop as revealed by magnetic susceptibility χ(T) = M(T) / H measurements. For x = 0.11 and 0.16 competing AFM and FM interactions result in a reentrant spin-glass behavior below T N , as evidenced by the observations of thermomagnetic hysteresis and magnetic relaxation. Enhanced FM fluctuations are also found for the x = 0.21 and 0.31 crystals, where χ c increases significantly at low T. A large χ anisotropy in these compositions where χ c is up to a factor of two larger than χ ab suggests that the FM spin fluctuations are quasi-1D in nature. Weak ferromagnetic contributions to the magnetization are found at T = 2 K for x = 0.11 –0.31. Heat-capacity Cp(T) measurements reveal the presence of FM quantum spin fluctuations for 0.11 ≤ x ≤ 0.52, where a logarithmic T dependence of Cp(T) / T is observed at low T. Here, the suppression of AFM order by the development of strong FM fluctuations in Ca(Co 1–x Ni x ) 2–y As 2 crystals suggests the presence of a FM quantum-critical point at x ≈ 0.20. Our density-functional theory (DFT) calculations confirm that FM fluctuations are enhanced by Ni substitutions for Co in CaCo 2–y As 2 . The Sommerfeld electronic heat-capacity coefficient is enhanced for x = 0, 0.21, and 0.42 by about a factor of two compared to DFT calculations of the density of states (DOS) at the Fermi energy, suggesting an enhancement of the DOS from electron-phonon and/or electron-electron interactions. The crystals with x > 0.52 do not exhibit FM spin fluctuations or magnetic order at T ≥ 1.8 K, which was found from the DFT calculations to arise from a Stoner transition. Superconductivity is not observed above 1.8 K for any of the compositions. Neutron-diffraction studies of crystals with x = 0.11 and 0.16 in the crossover regime (0.1 ≲ x ≲ 0.2) show no evidence of A-type ordering as observed in the parent compound with x = 0. Furthermore, no other common magnetic structures, such as ferromagnetic (FM), helical stacking of in-plane FM layers, or in-plane AFM structure, are found with an ordered moment greater than the uncertainty of 0.05 μ B per transition-metal atom.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Magnetic states of the quasi-one-dimensional iron chalcogenide Ba 2 FeS 3

Quasi-one-dimensional iron-based ladders and chains, with the 3d iron electronic density n=6, are attracting considerable attention. Recently, a new iron chain system Ba 2 FeS 3 , also with n=6, was prepared under high-pressure and high-temperature conditions. Here the magnetic and electronic phase diagrams are theoretically studied for this quasi-one-dimensional compound. Based on first-principles calculations, a strongly anisotropic one-dimensional electronic band behavior near the Fermi level was observed. In addition, a three-orbital electronic Hubbard model for this chain was constructed. Introducing the Hubbard and Hund couplings and studying the model via the density matrix renormalization group (DMRG) method, we studied the ground-state phase diagram. A robust staggered ↑-↓-↑-↓ AFM region was unveiled in the chain direction, consistent with our density functional theory (DFT) calculations. Furthermore, at intermediate Hubbard U coupling strengths, this system was found to display an orbital selective Mott phase (OSMP) with one localized orbital and two itinerant metallic orbitals. At very large U/W (W=bandwidth), the system displays Mott insulator characteristics, with two orbitals half-filled and one doubly occupied. Our results for high pressure Ba 2 FeS 3 provide guidance to experimentalists and theorists working on this one-dimensional iron chalcogenide chain material.

1-dimensional systems↗
DOE OSTI2025

Mechanistic Insights into Acetate Selectivity on Intermetallic CuPd(110) in CO Reduction

Experimental studies demonstrate that CuPd(110) uniquely favors acetate formation during CO reduction (CORR), contrasting with the preference for ethylene on Cu surfaces. To elucidate this selectivity, we employed explicit solvation density functional theory (DFT) calculations to investigate the reaction mechanism from both thermodynamic and kinetic angles. Here, our findings reveal that on CuPd(110), the acetate-pathway intermediate H 2 CCO is thermodynamically favored at experimental conditions, while CHCHO─a precursor to ethylene─is preferred on Cu(111). Beyond thermodynamics, we find that H 2 CCO is kinetically accessible under the experimental conditions on CuPd(110), aiding acetate formation. Electron density difference analyses further corroborate distinct protonation preferences supporting this mechanism. We propose a thermodynamic screening parameter based on the Gibbs free energy, G H 2 CCO < G CHCHO , as a guide for designing Cu-based catalysts with enhanced acetate selectivity. These results offer critical mechanistic insights into the CORR product distribution and a rational framework for future catalyst design.

Acetate↗
DOE OSTI2023

Understanding the influence of solvents on the Pt-catalyzed hydrodeoxygenation of guaiacol

Bio-oil derived from fast pyrolysis of lignocellulosic biomass needs to be deoxygenated to become a substitute for petroleum fuels. Here, we study the hydrodeoxygenation mechanism of guaiacol, a bio-oil model compound derived from the lignin fraction of biomass, on Pt(1 1 1) terrace sites in the presence of water, diethyl ether, 1-butanol, and n-hexane as solvent. Using first-principles periodic density functional theory (DFT) calculations and mean-field microkinetic reactor modeling, a detailed reaction mechanism is investigated targeting various products such as catechol, phenol, anisole, benzene, cyclohexanone, and cyclohexanol. Solvent phase DFT outcomes are mostly similar to that of the vapor phase; however, microkinetic modeling results suggest that rate controlling species and transition states differ somewhat in the various reaction environments. Catechol was found to be the major aromatic product across all reaction environments. Over Pt(1 1 1), unsaturated monooxygenate production from catechol is unlikely.

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

Prussian blue analogues as platform materials for understanding and developing oxygen evolution reaction electrocatalysts

Transition metal based materials containing Fe have drawn great attention as oxygen evolution reaction (OER) catalysts. The nature of the electrocatalytic active species remains under debate due to the ambiguous physicochemical properties of the catalyst materials, such as the oxidation states and crystal structures. Here, in order to address this issue, transition metal Prussian blue analogues (TM-PBA, Na(TM)(Fe)(CN) 6 , TM = V, Fe, Co, and Ni) with an isomorphous structure are investigated for OER catalysis. Our combined experimental measurements and density functional theory (DFT) calculations reveal that TM-PBAs exhibit volcano-like OER activity with Ni-PBA located near the top of the volcano. Such a volcano-like activity profile can be attributed to the distinctive binding energy difference between *O and *OH on different TM-PBAs surfaces. This research demonstrates that TM-PBAs can be used as platform materials for understanding structure-property-activity relationships in OER catalysts.

25 ENERGY STORAGE↗
DOE OSTI2024

Exploring Spin‐Orbit Effects in a [Cu 6 Tl] + Nanocluster Featuring an Uncommon Tl−H Interaction

Reaction of [CuH(PPh 3 )] 6 with 1 equiv. of Tl(OTf) results in formation of [Cu 6 TlH 6 (PPh 3 ) 6 ][OTf] ([1]OTf]), which can be isolated in good yields. Variable-temperature 1 H NMR spectroscopy, in combination with density functional theory (DFT) calculations, confirms the presence of a rare Tl−H orbital interaction. According to DFT, the 1 H chemical shift of the Tl-adjacent hydride ligands of [1] + includes 7.7 ppm of deshielding due to spin-orbit effects from the heavy Tl atom. In conclusion, this study provides valuable new insights into a rare class of metal hydrides, given that [1][OTf] is only the third isolable species reported to contain a Tl−H interaction.

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

Quantitative nanoscale tracking of oxygen vacancy diffusion inside single ceria grains by in situ transmission electron microscopy

Oxygen vacancy formation and migration in ceria is critical to its electrochemical and catalytic properties in systems for chemical and energy transformation, but its quantification is rather challenging especially at atomic-scale because of disordered distribution. Here we report a rational approach to track oxygen vacancy diffusion in single grains of pure and Sm-doped ceria at -20 °C to 160 °C using in situ (scanning) transmission electron microscopy ((S)TEM). To create a gradient in oxygen vacancy concentration, a small region (~30 nm in diameter) inside a ceria grain is reduced to the C-type CeO 1.68 phase by the ionization or radiolysis effect of a high-energy electron beam. The evolution in oxygen vacancy concentration is then mapped through lattice expansion measurement using scanning nano-beam diffraction or 4D STEM at a spatial resolution better than 2 nm; this allows direct determination of local oxygen vacancy diffusion coefficients in a very small domain inside pure and Sm-doped ceria at different temperatures. Further, the activation energies for oxygen transport are determined to be 0.59, 0.66, 1.12, and 1.27 eV for pure CeO 2 , Ce 0.94 Sm 0.06 O 1.97 , Ce 0.89 Sm 0.11 O 1.945 , and Ce 0.8 Sm 0.2 O 1.9 , respectively, implying that activation energy increases due to impurity scattering. The results are qualitatively supported by density functional theory (DFT) calculations. In addition, our in situ TEM investigation reveals that dislocations impede oxygen vacancy diffusion by absorbing oxygen vacancies from the surrounding areas and pinning them locally. With more oxygen vacancies absorbed, dislocations show extended strain fields with local tensile zone sandwiched between the compressed ones. Therefore, dislocation density should be reduced in order to minimize the resistance to oxygen vacancy diffusion at low temperatures.

36 MATERIALS SCIENCE↗