Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Vacancies”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Point defect evolution under irradiation: Finite size effects and spatio-temporal correlations

The evolution of point defect concentrations under irradiation is controlled by their diffusion properties, and by their formation and elimination mechanisms. The latter includes the mutual recombination of vacancies and interstitials, and the elimination of point defects at sinks. Two models are traditionally used to predict the defect concentration evolution, the standard rate theory (SRT) and the atomistic kinetic Monte Carlo (AKMC). We show in this work a large discrepancy in the defect concentrations between both methods when the average number of defects in the AKMC simulation box is close or lower than one. The reason is that AKMC naturally captures strong space and time correlations between vacancies and interstitials generated by the finite size of the periodic simulation box. These correlations strongly affect the recombination rate and the point defect concentrations. SRT fails to predict such correlations, and the corresponding solution deviates from the more accurate solution given by the AKMC simulation under similar conditions. These finite size correlation effects are strong when the elimination of point defects occurs by recombination only, but can still be significant in the presence of sinks. In order to account for the spacio-temporal correlation in a continuum framework, we introduce a Correlated Pair Theory (CPT). This theory fully takes into account the correlations between vacancy and interstitial pairs and predicts point defect concentrations in good agreement with AKMC simulations. Inversely, if the goal is for the AKMC to reproduce defect concentrations in bulk using small simulation boxes, the naturally occurring correlations need to be corrected. We show here that the CPT can be used to modify the elimination rates in the AKMC simulations, so as to yield point defect concentrations in agreement with SRT.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Cluster dynamics modeling of irradiation growth in single crystal Zr

Irradiation growth of single-crystal Zr is modeled with a reduced set of cluster dynamics equations. Nucleation and growth of basal and prismatic dislocation loops are both accounted for in the model. Equations are developed for the time evolution of single point defects, a limited number of point defect clusters, interstitial and vacancy loop nucleation rates, as well as for the growth of vacancy loops on basal planes and interstitial loops on prismatic planes. Reduction of the usual infinite hierarchy of cluster dynamics equations to the simple set studied here is justified on the physical basis of the stability of small loops once nucleated. This simplified cluster dynamics model with a small number of adjustable parameters avoids the complexity of explicit representation of higher order point defect clusters. The model shows consistency with experimental observations of the following aspects: (1) the growth rates of Zr crystals along the a- and c-axes; (2) the onset dose for breakaway irradiation growth; (3) the saturation dislocation loop densities of vacancy and interstitial loops; (4) the effects of cold work and temperature on irradiation growth.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Correlations for the specific heat capacity of ( U x Pu 1 - x ) 1 - y Gd y O 2 - z derived from molecular dynamics

We report UO 2 is the primary conventional fuel used in most nuclear reactors with Gd 2 O 3 commonly added as a burnable absorber to produce a more level power distribution in the reactor core at the beginning of operation. It can also be mixed with other actinide oxides to produce mixed oxide (MOx) fuel. In this study, molecular dynamics simulations were used to predict the specific heat capacity of Gd-doped PuO 2 , UO 2 and (U, Pu)O 2 MOx accommodating Gd 3+ substituted at cation sites via two charge compensation mechanisms - oxygen vacancy formation and the oxidation of U 4+ to U 5+ . The specific heat capacity values for PuO 2 and UO 2 are in good agreement with other studies showing a distinct peak at high temperatures - above 1800 K. As Gd 3+ is added, the peak height reduces for each composition considered. An analytical fit was applied to the data where Gd 3+ was fully charge compensated by either oxygen vacancies or U 5+ . The expression was then validated by predicting the specific heat capacity for three compositions of (Ux Pu 1-x ) 1-y Gd y O 2-z containing both oxygen vacancies and U 5+ , and compared to molecular dynamics data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Microstructural evolution of compositionally complex solid-solution alloys under in-situ dual-beam irradiation

Here, this work attempts to link the microstructural evolution of single-phase compositionally complex alloy (CCA) compositions under dual-beam irradiation to their Mn-content via the stacking-fault energy (SFE) and vacancy migration energies. Two alloys, Cr 18 Fe 27 Mn 27 Ni 28 and Cr 15 Fe 35 Mn 15 Ni 35 , along with less compositionally complex pure Ni and Fe 56 Ni 44 binary were irradiated at 500 and 600 °C under dual-beam 1 MeV Kr 2+ and 16 keV He + heavy-ions up to 7 displacements per atom (dpa) with a He/dpa ratio of 0.75 %/dpa using in situ transmission electron microscopy (TEM). Due to the bubble-stabilizing effect of implanted He, bubbles were observed in all irradiations, and populations of faulted interstitial loops were characterized in Cr 18 Fe 27 Mn 27 Ni 28 and Cr 15 Fe 35 Mn 15 Ni 35 . A reduction in swelling was observed in the two CCAs compared to pure Ni and Fe 56 Ni 44 . Although swelling increased from 500 to 600 °C in Fe 56 Ni 44 , Cr 15 Fe 35 Mn 15 Ni 35 and Cr 18 Fe 27 Mn 27 Ni 28 both swelled slightly more at 500 °C. This was attributed to the difference in vacancy mobility, stronger pinning effect of vacancies on He, and the sink strength of faulted dislocation loops. Faulted interstitial loops nucleated with a higher number density and dislocation line density in Cr 15 Fe 35 Mn 15 Ni 35 at both temperatures, and at 600 °C in both materials. The differences in faulted loop population and temperature effect on swelling are correlated to the Mn-content and the measured SFE (20.2 ± 6.7 mJ/m 2 for Cr 18 Fe 27 Mn 27 Ni 28 and 9.2 ± 3.4 mJ/m 2 for Cr 15 Fe 35 Mn 15 Ni 35 ).

36 MATERIALS SCIENCE↗
DOE OSTI2025

Back to functional hydrides: Effects of neutron-irradiated microstructure on hydrogen retention in yttrium hydride

Functional hydrides are promising candidates for advanced nuclear reactors, particularly in portable or transportable applications, due to their high hydrogen-retention capabilities, enabling efficient neutron moderation, and compact reactor design. However, hydrogen mobility in hydrides at elevated irradiation temperatures poses significant technological challenges, necessitating a comprehensive understanding of their irradiation behavior. Furthermore, this study investigated the microstructural and chemical stability of neutron-irradiated yttrium hydrides to assess their hydrogen-retention capacity. A targeted literature review was also conducted to contextualize neutron-irradiation effects on functional hydrides with regards to structural stability and hydrogen retention. Experimental characterizations revealed that, at high temperatures, irradiated hydrides retained their phase stability, which was likely enhanced by irradiation-induced microstructure evolution. Notably, an amorphous yttrium and oxygen -rich surface layer was present at the free surface of the hydride. Its thickness decreased while a continuous crystalline Y-O-rich layer was formed with increasing neutron damage. Additionally, the number density of dislocation loops and cavities generally increased as a function of neutron dose. First-principles calculations of hydrogen behavior within yttrium vacancy clusters in yttrium hydrides demonstrated vacancy-size-dependent hydrogen stability and configuration, highlighting the role of vacancy geometry in regulating hydrogen retention. Thus, the presence of irradiation-induced dislocation loops and cavities were hypothesized to improve hydrogen retention. Collectively, these findings advance the understanding of hydride behavior under neutron irradiation as well as their technological readiness for portable or transportable nuclear reactors.

Defect clusters↗
DOE OSTI2026

The impact of chemistry on anion migration in bixbyite-structured lanthanide oxides

This manuscript describes atomistic calculations of oxygen vacancy and interstitial migration in bixbyite structured lanthanide oxides. We examine two types of compounds, one in which only one type of lanthanide cation is present and a second class in which two lanthanides are present in a 3:1 ratio as dictated by the symmetry of the bixbyite crystal structure. Using temperature accelerated dynamics and the nudged elastic band method, we quantify the role of chemistry on the energy barriers for the most important pathways for both vacancy and interstitial migration. We then analyze the impact of these variations on the overall diffusivity of each defect, quantifying the contribution of each pathway using the theory of kinosons. We find that vacancy mobility can vary by as much as three orders of magnitude through changes in chemistry at 500 K. Changes in interstitial mobility are more modest but can still vary by an order of magnitude. This points to the ability to tune the mass transport characteristics of these compounds through appropriate choices in chemistry. We have also included supplementary information containing the atomic structures of the relevant pathways.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Tuning Point Defects by Elastic Strain Modulates Nanoparticle Exsolution on Perovskite Oxides

Exsolution generates stable and catalytically active metal nanoparticles via phase precipitation out of a host oxide. An ability to control the size and dispersion of the exsolution particles is desirable for design of nanostructured (electro)catalysts. In this work, we demonstrate that tuning point defects by lattice strain affects both the thermodynamics and the kinetics of iron (Fe 0 ) exsolution on La 0.6 Sr 0.4 FeO 3 (LSF) thin film model. By combining in situ surface characterization and ab initio defect modeling, we show oxygen vacancy and Schottky defects to be the primary point defects formed upon Fe0 exsolution. Lattice strain tunes the formation energy, and thus the abundance of these defects, and alters the amount and size of the resulting exsolution particles. In addition, we find that the density of exsolved nanoparticles matches the concentration of oxygen vacancy pairs, thus pointing to the surface oxygen vacancy pairs as preferential nucleation sites for exsolution. The tensile-strained LSF with a facile formation of these critical point defects results in a higher Fe 0 metal concentration, a larger density of nanoparticles, and a reduced particle size at its surfaces. These results provide important mechanistic insights and highlight the role of point-defect engineering in designing nanostructured catalysts in energy and fuel conversion technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Investigating the Electronic Structure of Prospective Water-Splitting Oxide BaCe 0.25 Mn 0.75 O 3-δ before and after Thermal Reduction

BaCe 0.25 Mn 0.75 O 3-δ (BCM), a non-stoichiometric oxide with a layered perovskite-like crystal structure, has recently emerged as a prospective contender for application in renewable energy harvesting by solar thermochemical hydrogen generation. Using solar-thermal energy and a reducing environment, oxygen vacancies can be created in high-temperature BCM, and the reduced crystal so obtained can, in turn, produce H 2 by stripping oxygen from H 2 O. Therefore, a first step toward understanding the working mechanism and optimizing the performance of BCM is a thorough and comparative analysis of the electronic structure of the pristine and the reduced material. In this paper, we probe the electronic structure of BCM using the combined effort of first-principles calculations and experimental O K-edge X-ray absorption spectroscopy (XAS). The computed projected density of states (PDOS) and orbital plots are used to propose a simplified model for orbital mixing between the oxygen and metal atoms. With the help of state-of-the-art simulations, we are able to find the origins of the XAS peaks and categorize them on the basis of contribution from Ce and Mn. For the reduced crystal, the calculations show that the change in electron density resulting from the reduction is strongly localized around the oxygen vacancy. Experimental measurements reveal a marked lowering of the first O K-edge peak in the reduced crystal. Using theoretical analysis, this is shown to result from lifting of spin degeneracy in the absorption peaks as well as from a diminished O 2p contribution to the frontier unoccupied orbitals, in accordance with the tight binding scheme. The simulated results serve as a reference for the extent of spectral change as a function of the percentage of oxygen vacancies in the reduced crystal. Here, our study paves the way for the investigation of the working mechanism of BCM and for computational and experimental efforts aimed at design and discovery of efficient water-splitting oxides.

08 HYDROGEN↗
DOE OSTI2024

Correlated Displacement of Dynamic Elastic Dipoles Produces Nonclassical Electrostriction in Zr-Doped Ceria

By combining experimental data with density functional theory-based ab initio molecular dynamics modeling, this work provides evidence that nonclassical electrostriction in isovalent Zr-doped ceria is due to the correlated anharmonic motion of dynamic elastic dipoles associated with multiple [ZrO 8 ]-local bonding units with a high Zr concentration (Zr 0.1 Ce 0.9 O 2 ). Introduction of 0.5 mol % trivalent or divalent codopants (Sc, Yb, La, or Ca) reduces the longitudinal electrostriction strain coefficient by more than a factor of 10, produces a 3-fold decrease in the relative dielectric permittivity, and increases the elastic modulus. Since these changes depend neither on the radius nor on the valency of the codopant, we conclude that the responsible species are charge-compensating oxygen vacancies (VO). For trivalent dopants (Do 0.005 Zr 0.1 Ce 0.895 O 1.9975 ), oxygen vacancies are present at a concentration ratio 1:40 with respect to Zr, giving, for random distribution, a characteristic interaction distance of ≤2.3 unit cells (1.2 nm). Oxygen vacancies participate in [ZrO 7 -V O ] local bonding units, disrupting the correlated dynamic displacements of the connected [ZrO 8 ]-local bonding units. Finally, such correlated motion of dynamic elastic dipoles may also explain the exponential increase in the longitudinal electrostriction strain coefficient with an increase in Zr concentration to <0.2 mole fraction and must be taken into account for further development of nonclassical electrostrictors based on Zr-doped ceria.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Fundamental Thermodynamic, Kinetic, and Mechanical Properties of Lithium and Its Alloys

Lithium alloying reactions are beneficial in promoting uniform plating and stripping of lithium metal in all-solid-state batteries. First-principles calculations are performed to predict thermodynamic, kinetic, and mechanical properties of lithium and several important Li–M alloys (M = Mg, Ag, Zn, Al, Ga, In, Sn, Sb, and Bi). While the Li–Mg binary system forms a solid solution, most other lithium–metal alloys prefer stoichiometric intermetallic compounds with common local motifs that enable fast Li diffusion. Lithium and Li-rich alloys exhibit an unusually flat energy landscape along paths that connect BCC to close-packed structures like FCC and HCP, with important implications for mechanical properties. Very low migration barriers for Li diffusion that rival those of superion conductors are predicted, both in pure Li and in Li–M intermetallics. However, vacancy concentration, which is crucial for substitutional diffusion, is predicted to be low in metallic Li and most Li–M intermetallics. Compounds such as B32 LiAl and LiGa as well as D0 3 Li 3 Sb and Li 3 Bi exhibit structural vacancies at higher ends of their voltage windows, which together with low migration barriers leads to exceptionally high Li mobilities. In the Li–Mg solid solution, the addition of Mg is found to decrease the vacancy tracer diffusion coefficient by an order of magnitude.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Role of Short-Range Order on Diffusion Coefficients in the Li–Mg Alloy

Li–Mg alloys are important because of their beneficial role in fostering uniform plating and stripping of lithium in all-solid-state batteries. The alloy Li x Mg 1–x forms a solid solution on the BCC crystal structure when the lithium content is greater than x ≈ 0.3. The activation barriers of lithium and magnesium exchanges with a vacancy, crucial for substitutional diffusion, are predicted to be exceptionally low and almost identical, with negligible dependence on the alloy composition. The equilibrium vacancy concentration at room temperature is predicted to be very low, and it also remains almost constant with no dependence on Mg content in the alloy (for x Li ≥ 0.5). Nevertheless, both experiments and kinetic Monte Carlo simulations indicate that the tracer diffusion coefficients decrease by almost an order of magnitude with the addition of Mg to the alloy. In this contribution, the crucial role that chemical short-range order plays in affecting the diffusion coefficients is studied. Chemical short-range order is found to increase the effective activation barrier for lithium and magnesium diffusion by making successive atomic hops with vacancies correlated.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Defect Diffusion Graph Neural Networks for Materials Discovery in High-Temperature Energy Applications

Here, the migration of crystallographic defects dictates material properties and performance for a plethora of technological applications. Density functional theory (DFT)-based nudged elastic band (NEB) calculations are a powerful computational technique for predicting defect migration activation energy barriers, yet they become prohibitively expensive for high-throughput screening of defect diffusivities. Without introducing hand-crafted (i.e., chemistry- or structure-specific) descriptors, we propose a generalized deep learning approach to train surrogate models for NEB energies of vacancy migration by hybridizing graph neural networks with transformer encoders and simply using pristine host structures as input. With sufficient training data, computationally efficient and simultaneous inference of vacancy defect thermodynamics and migration activation energies can be obtained to compute temperature-dependent vacancy diffusivities and to down-select candidates for more thorough DFT analysis or experiments. Thus, as we specifically demonstrate for potential water-splitting materials, candidates with desired defect thermodynamics, kinetics, and host stability properties can be more rapidly targeted from open-source databases of experimentally validated or hypothetical materials.

14 SOLAR ENERGY↗
DOE OSTI2023

Quantum Simulations of Fermionic Hamiltonians with Efficient Encoding and Ansatz Schemes

Here, we propose a computational protocol for quantum simulations of fermionic Hamiltonians on a quantum computer, enabling calculations on spin defect systems which were previously not feasible using conventional encodings and a unitary coupled-cluster ansatz of variational quantum eigensolvers. We combine a qubit-efficient encoding scheme mapping Slater determinants onto qubits with a modified qubitcoupled cluster ansatz and noise-mitigation techniques. Our strategy leads to a substantial improvement in the scaling of circuit gate counts and in the number of required qubits, and to a decrease in the number of required variational parameters, thus increasing the resilience to noise. We present results for spin defects of interest for quantum technologies, going beyond minimum models for the negatively charged nitrogen vacancy center in diamonds and the double vacancy in 4H silicon carbide (4H-SiC) and tackling a defect as complex as negatively charged silicon vacancy in 4H-SiC for the first time.

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

Excited State Properties of Point Defects in Semiconductors and Insulators Investigated with Time-Dependent Density Functional Theory

Here, we present a formulation of spin-conserving and spin-flip hybrid time-dependent density functional theory (TDDFT), including the calculation of analytical forces, which allows for efficient calculations of excited state properties of solid-state systems with hundreds to thousands of atoms. We discuss an implementation on both GPU- and CPU-based architectures along with several acceleration techniques. We then apply our formulation to the study of several point defects in semiconductors and insulators, specifically the negatively charged nitrogen-vacancy and neutral silicon-vacancy centers in diamond, the neutral divacancy center in 4H silicon carbide, and the neutral oxygen-vacancy center in magnesium oxide. Our results highlight the importance of taking into account structural relaxations in excited states in order to interpret and predict optical absorption and emission mechanisms in spin defects.

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

Structure and Reactivity of Pristine and Reduced Spinel CoFe 2 O 4 (001)/(100) Surfaces

Cobalt ferrite, CoFe 2 O 4 (CFO), nanocrystals are efficient and competitive anode materials in the field of electrochemical water splitting. Using density functional theory with on-site Hubbard U corrections (DFT+U), we have investigated the structural, electronic and magnetic properties of CFO (001)/(100) surfaces, as well as their reactivities towards water adsorption. Special attention has been focused on the formation of oxygen vacancies (V O ), due to their key role in the oxidation activity of metal oxides, often based on the Mars-van-Krevelen mechanism. Our results show that vacancy formation is easiest at oxygen sites that are not bound to tetrahedrally-coordinated Fe. Water adsorbs mainly in molecular form on the Co/Fe metal cations, whereas it dissociates at defects. In comparison to other spinels, CFO is similar to NiFe 2 O 4 , exhibiting relatively low energy cost of V O formation and a strong affinity of the vacancies towards water. Furthermore, these findings suggest that CFO may be a more promising oxidation catalyst than NiCo 2 O 4 and Co 3 O 4 .

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

NO 2 Interactions with MoO 3 and CuO at Atmospherically Relevant Pressures

NO x concentrations in some geographic regions are harmful to human health. Gas filters to trap NO x and other toxic chemicals contain metal oxides, including MoO 3 and CuO. These materials are also being investigated for NO x gas sensors. In a step to understand the fundamental adsorption mechanism in sensors and the effect on binding site availability in gas filters, ambient-pressure X-ray photoelectron spectroscopy (APXPS) was used to study the interaction of NO 2 with polycrystalline MoO 3 and CuO surfaces under pressures up to 0.01 Torr (14 parts per million volume (ppmv)). Density functional theory-based computational modeling was performed to reveal the mechanisms of NO 2 interactions with the MoO 3 (010) and CuO(111) surfaces to aid interpretation of the experimental results. With pressure dependence, NO 2 interacts with reduced Mo 5+ atoms generated by oxygen vacancies and abstracts hydrogen atoms from hydroxyl groups on MoO 3 without accumulating N-containing species on the surface; vacancy-induced electronic states in the band gap are also removed, hinting toward an increase in the resistivity of the material. N-containing species begin accumulating on the CuO surface at atmospherically relevant pressures of 140 ppbv. NO 2 only decomposes at oxygen vacancy sites of CuO. The nitrogen species leave the CuO surface upon evacuation, highlighting the importance of in situ surface characterization when studying gas sensing and adsorption mechanisms. Finally, these results imply that NO 2 removes hydroxyl and O vac binding sties on these materials when used in gas filtration and sensing applications. Furthermore, the results show the key role of O vac sites in the gas sensing mechanism of MoO 3 and highlight the potential of APXPS for further studies of gas sensors.

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

LaMnO 3 Dopants for Efficient Thermochemical Water Splitting Identified by Density Functional Theory Calculations

While ceria is the standard material for two-step water splitting, perovskites are emerging as viable alternatives. In this work, based on the orthorhombic LaMnO 3 supercell, we substitute Li Na K Rb Mg Ca Sr Ba on the A-sites (La sites) and Al Ga In Mg Zn on the B-sites (Mn sites) at a concentration of 37.5%. The range of temperature and oxygen partial pressure at which each composition is stable is predicted. For compositions that are stable in relevant temperature and pressure ranges, the oxygen vacancy formation energies are determined for all of the oxygen vacancy site positions available in the computational supercell. Mg, Ca, Sr, and Ba A-site-substituted LaMnO 3 and Al and In B-site-substituted LaMnO 3 meet these two criteria for candidates in solar-thermal water splitting applications. Finally, oxygen vacancy formation energy can also be controlled by adjusting the doping strategy.

08 HYDROGEN↗
DOE OSTI2024

Variation in Cation Adsorption Mechanism Controlled by Chemical and Structural Heterogeneities at the Quartz (101)–Water Interface

Mineral–water interfacial reactions are central to chemical processes that control the fate of nutrients and contaminants in natural environments. Mineral surfaces commonly have complex structures and compositions whose impact on interfacial reactivity is poorly understood. Here, in this work, we investigated the effects of surface heterogeneities on Rb + sorption at the quartz (101)–10 mM RbCl solution interface at pH 9.8 using in situ high-resolution X-ray reflectivity. Two surface locales (i.e., Spots A and B) having distinct interfacial structures were chosen: Spot A was characterized by its low defect density (≤20% topmost Si vacancies) and Rb + adsorption occurred predominantly as an inner-sphere complex. In comparison, Spot B had a higher defect density (~50% vacancies) and was covered with poorly crystalline SiO 2 . A substantially larger Rb + uptake (i.e., 7-times higher coverage) was observed on this defective surface where Rb + incorporated in the vacancy sites (confirmed by density functional tight binding-based molecular dynamics simulations) or adsorbed directly on the disordered film. These results provide a direct quantification of how surface heterogeneity influences the geochemical behavior of mineral–water interfaces, in particular highlighting the important role of chemical and structural defects on the sorbate speciation and coverage at silicate mineral surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗