Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “vacancy binding”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Atomistic modeling of lanthanide diffusion in refractory body-centered cubic molybdenum

Lanthanide fission products can strongly interact with candidate cladding alloys, but their transport properties in refractory metals remain poorly understood. Here, in this work, we investigate the atomic-scale diffusion behavior of La, Ce, Pr, and Nd in body-centered cubic (bcc) molybdenum, a potential candidate for advanced nuclear cladding. Self-consistent mean-field transport modeling is performed to evaluate the fission product transport and vacancy mobility, informed by first-principles and nudged elastic band calculations of vacancy formation energies, migration barriers, and solute–vacancy binding characteristics. Compared with bcc Fe, lanthanide solutes in bcc Mo exhibit slower tracer diffusion due to higher vacancy formation and migration energies. Furthermore, the calculations reveal that the influence of fission products on migration barriers in bcc Mo are not as extensive in range compared to bcc Fe. Among the studied lanthanides, La exhibits the strongest vacancy binding while also being the fastest diffuser in Mo. These findings highlight how refractory bcc alloys can reduce fission product infiltration, offering valuable insight into the development of durable cladding systems for advanced reactors.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Diffusion Behavior of Oversized Fission Products in bcc Fe Cladding: A First-Principles Study

Fuel-Cladding Chemical Interaction (FCCI) poses significant challenges in nuclear reactors, where fission products from nuclear fuel interact with Fe-based cladding materials, potentially compromising their structural integrity. This study investigates the diffusion behavior of oversized fission products, Pr, Nd, Ce, and La, within bcc Fe cladding using density functional theory (DFT), nudged elastic band (NEB) method, and self-consistent mean field (SCMF) theory. Our results reveal significant long-range vacancy binding energies, particularly up to the 6th nearest neighbor, with La exhibiting the strongest binding affinity, followed by Nd, Ce, and Pr. The NEB calculations indicate significant high barriers for the dissociation of 1nn vacancy-solute pairs for all fission products. The tracer diffusion coefficients of these fission products was derived in Arrhenius form. The significant trapping effect of vacancies by a very dilute amount of fission products reduces vacancy mobility, leading to an oversaturation of point defects, void nucleation, and swelling. These are critical issues for irradiated cladding materials. The tracer diffusion coefficients indicate that Nd diffuses the fastest, followed by La, Ce, and Pr. This study provides essential insights for developing advanced cladding materials and design strategies to mitigate FCCI, ultimately enhancing nuclear reactor safety and performance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Diffusion behavior of lanthanide fission products in bcc Fe cladding: A first-principles study

Fuel-cladding chemical interaction poses significant challenges in nuclear reactors, where fission products generated from nuclear fuel interact with Fe-based cladding materials, potentially compromising their structural integrity. This study investigates the diffusion behavior of lanthanide fission products, Lanthanum (La), Cerium (Ce), Praseodymium (Pr), and Neodymium (Nd), within body-centered cubic (bcc) Fe cladding using the density functional theory, nudged elastic band method, and self-consistent mean field theory. Our results reveal significant vacancy binding energies, particularly with the 1st and 2nd nearest neighbors, which diminish beyond the 5th nearest neighbor, with La exhibiting the strongest binding affinity, followed by Nd, Ce, and Pr. The nudged elastic band calculations indicate significant high barriers for the dissociation of 1st nearest neighbor vacancy-solute pairs for all fission products. The tracer diffusion coefficients of these fission products were derived in an Arrhenius form, with a magnetic correction that accounts for the high-temperature paramagnetic state. The significant trapping effect of vacancies caused by a very dilute concentration of fission products reduces vacancy mobility, potentially leading to modifications in point defect supersaturation, void nucleation, and swelling under irradiation. These represent critical challenges for irradiated cladding materials. The tracer diffusion coefficients indicate that Nd diffuses the fastest, followed by La, Ce, and Pr. Furthermore, this study provides essential insights for understanding fission product transport in cladding materials and informs future design strategies to mitigate fuel-cladding chemical interaction, ultimately enhancing nuclear reactor safety and performance.

Diffusion↗
DOE OSTI2025

Density Functional Theory Study of Iron–Oxygen Divacancies in Magnetite (Fe 3 O 4 ) and Hematite (Fe 2 O 3 )

Density functional theory (DFT) calculations are employed to investigate the formation energies, charge redistribution, and binding energies of iron–oxygen divacancies in magnetite (Fe 3 O 4 ) and hematite (Fe 2 O 3 ). For magnetite, we focus on the low-temperature phase to explore variations with local environments. Building on previous DFT calculations of the variations in formation energies for oxygen vacancies with local charge and spin order in magnetite, we extend this analysis to include octahedral iron vacancies before analyzing the iron–oxygen divacancies. We also assessed the relative stability of iron–oxygen divacancies by comparing their formation energies with those of individual vacancies. Our findings reveal a significant energetic driving force for the formation of divacancy clusters, particularly in magnetite, where divacancies in the +1 charge state exhibit formation energies comparable to those of neutral iron vacancies under oxidizing conditions. In hematite, the results indicate a strong tendency for oxygen vacancies to bind to iron vacancies. These results highlight the significance of iron–oxygen vacancy complexes in the transport properties of iron oxides, with particular relevance to diffusion mechanisms under irradiation conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Robustness of Vacancy-Bound Non-Abelian Anyons in the Kitaev Model in a Magnetic Field

Non-Abelian anyons in quantum spin liquids (QSLs) provide a promising route to fault-tolerant topological quantum computation. In the exactly solvable Kitaev honeycomb model, such anyons of the QSL state can be bound to nonmagnetic spin vacancies and endowed with non-Abelian statistics by an infinitesimal magnetic field. Here, we investigate how this approach for stabilizing non-Abelian anyons extends to a finite magnetic field represented by a proper Zeeman term. Through large-scale density-matrix renormalization group simulations, we compute the vacancy-anyon binding energy as a function of magnetic field for both the ferromagnetic and antiferromagnetic Kitaev models. Here, we find that anyon binding remains robust within the entire QSL phase for the ferromagnetic Kitaev model but breaks down already inside this phase for the antiferromagnetic Kitaev model. To compute a binding energy several orders of magnitude below the magnetic energy scale, we introduce both a refined definition and an extrapolation scheme based on carefully tailored perturbations.

Xiao, Bo [Oak Ridge National Laboratory (ORNL), Oa↗
DOE OSTI2025

Effects of Minor Alloying Elements and Impurities on Long-Range Ordering in Ni-Cr-Based Alloys

There is growing concern about long-range ordered Ni2Cr phase forming during thermal aging of Ni-Cr-based alloys (e.g., 690, 625), which leads to lattice contraction, hardening, and embrittlement, possibly limiting the service life of light water reactor components. The timespan when this phase transformation occurs can vary in orders of magnitude between model and commercial alloys and is likely related to the presence of major (e.g., Fe, Mo) and minor alloying elements (e.g., Ti, Si, Nb). These minor alloying elements have been hypothesized to bind with vacancies and therefore limit diffusion. Contrarily, impurities (e.g., P, S) are believed to enhance vacancy diffusion and accelerate the Ni2Cr phase transformation because of a vacancy drag mechanism. This study implements experimental isothermal aging of model Ni-Cr binary and Ni-Cr-X ternary alloys, where X is Ti, Si, Nb, and P to help understand how each minor alloying element affects the Ni2Cr phase transformation timeline. Synchrotron-based X-ray diffraction and hardness measurements are used to quantify the evolution of the Ni2Cr phase transformation through a variety of short aging times to identify the onset of long-range order. These experimental results are compared to first principles simulations from literature to understand the mechanism by which each element impacts vacancy diffusion and the formation of the long-range ordered phase, Ni2Cr.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Understanding the effect of minor alloying elements on helium bubble formation in ferritic-martensitic steels

Ferritic-martensitic steels are promising structural materials for advanced nuclear reactors. To minimize long-term radioactivity, reduced-activation ferritic-martensitic steels have been developed by substituting high-activation elements like Ni and Mo with low-activation elements such as W. However, the impact of these alloying modifications on helium bubble formation, which plays a key role in material swelling, remains unclear. Here, in this study, we compared helium bubble formation in ferritic-martensitic steel T91 and reduced-activation ferritic-martensitic steel F82H. Both materials were irradiated with sequential 100 keV, 150 keV, and 200 keV helium ions to a dose of 0.5 dpa and a helium concentration of 9,000 appm at 500°C. The helium bubbles in F82H exhibited a larger average size and a lower density than those in T91, suggesting differences in minor alloying elements may influence the bubble growth. Here, to investigate the effects of these alloying elements, we characterized radiation-induced segregation near bubbles and grain boundaries. Prominent Ni-Mn-Si enriched clusters were found near bubbles in T91, while only Mn-Si enriched clusters were found near bubbles in F82H. In addition, the obvious Cr enrichment near grain boundaries was absent around bubbles in both steels. The different segregation trends among elements revealed the variations in element diffusion mechanisms and the different sink biases between bubbles and grain boundaries. Cr enrichment near grain boundaries is mostly driven by interstitial-mediated diffusion. However, since bubble growth relies on net vacancy flux, vacancy-mediated diffusion plays a dominant role in controlling element segregation near bubbles. Therefore, Cr enrichment was not found near bubbles. Because of preferential vacancy-drag diffusion for Ni, Si and Mn, these elements were enriched near bubbles. Due to the strong binding energies of vacancies with these solute atoms, the vacancy diffusivity can be reduced near these solutes. Therefore, the more prominent Ni-Si-Mn clustered near helium bubbles in T91 lead to stronger suppression of helium bubble growth compared to F82H.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Effects of heteroatom doping on hydrogen uptake in tungsten oxide

Redox-active transition metal oxides (TMOs) that can undergo proton-insertion coupled electron transfer (PICET) are promising candidates for catalyzing molecular conversion reactions which require the transfer of hydrogen atoms (or the thermochemical equivalent, H + , e − ). Herein, we studied the effects of isovalent (Mo 6+ ) and aliovalent (V 5+ and Nb 5+ ) heteroatom doping on the electrochemical PICET behavior of monoclinic tungsten oxide (WO 3 ). Cyclic voltammetry in aqueous acidic electrolytes shows that the addition of redox-active heteroatoms (Mo 6+ and V 5+ ) leads to systematic shifts in redox couple half-wave potentials (E 1/2 ), broadening, and an overall decrease in the current response. Conversely, the non-redox active heteroatom (Nb 5+ ) only reduces the current response with no observable peak-current broadening. This broadening is attributed to changes in the proton binding affinities of oxygen in different chemical environments, i.e., bridging different pairs of redox-active transition metal cations. We determined the hydrogen bond dissociation free energy (H BDFE) values to elucidate the thermodynamic effect of heteroatom substitution. Density functional theory calculations reveal a differentiation in the hydrogen binding and oxygen vacancy formation energies between heteroatom doped structures. The PICET-induced structural phase transitions of the pristine and doped samples were further probed with operando electrochemical X-ray diffraction (EC-XRD) and with ex situ chemical reduction. The broadening of the potential-dependent current response with increased heteroatom doping manifests in the operando EC-XRD results as prolonged structural regions where multiple hydrogen bronze phases exist and the appearance of cubic bronze phases at lower degrees of reduction compared to pristine WO 3 .

Holzapfel, Noah P. [North Carolina State Universit↗
DOE OSTI2025

Effect of Mg and Ni impurities on tritium diffusion in lithium ceramics through cluster dynamics simulations

This study investigates the impact of Mg and Ni doping on tritium diffusion in LiAlO 2 and LiAl 5 O 8 ceramics, that are used in tritium-producing burnable absorber rods (TPBARs). Utilizing Centipede simulations across a broad temperature range (500 K to 1250 K), we explore the interplay between defect dynamics, cluster formation, and tritium mobility. In LiAlO 2 , Mg doping significantly enhances tritium diffusivity by increasing tritium interstitial concentrations and diffusion coefficients of key species, thereby doubling the overall tritium diffusivity. Ni doping, while shifting the dominant defect to Li vacancies, maintains high tritium mobility due to the low binding energy of Li vacancy-tritium complexes, which ensures effective tritium migration. In LiAl 5 O 8 , Mg and Ni doping results in a slight reduction in the diffusion coefficients of key species, yet the dramatic increase in tritium interstitial concentrations compensates, leading to a net small increase in tritium diffusivity. In conclusion, the findings highlight the critical role of defects in tritium transport and the effect of Mg and Ni defects on the performance of these ceramics in demanding nuclear environments.

Cluster dynamics↗
DOE OSTI2025

Radiation-reduced diffusion of Nd in bcc Fe

Here, this study combines ab initio calculations and kinetic Monte Carlo simulations to investigate the bulk diffusion of Nd, an exemplar lanthanide element, in bcc Fe. Contrary to the conventional belief that irradiation accelerates diffusion processes, our study reveals that irradiation can unexpectedly reduce the diffusion of Nd due to the strong binding between substitutional Nd atoms and vacancies in bcc Fe. Combined with the high concentration of nonequilibrium vacancies, this promotes the abundant formation of Nd+3Va defect complexes containing a single Nd substitutional atom and three vacancies, accompanied by the depletion of Nd-vacancy pairs that mediate the diffusion of Nd.

Density functional theory (DFT)↗
DOE OSTI2025

Effect of Sn microalloying on the nucleation of L1 2 Al 3 Zr precipitates in a dilute aluminum-zirconium alloy

While L1 2 -Al 3 Zr nanoprecipitates provide a balance between strengthening and good electrical conductivity, the precipitation of L1 2 -Al 3 Zr in aluminum requires aggressive heat treatments. An improved age-hardening response was observed during isochronal aging of an Al-0.24Zr (wt%) alloy when microalloyed with Sn. A new mechanism termed Low melting point Element-Assisted Nucleation (LEAN) is proposed to explain the lower temperature nucleation of L1 2 -Al 3 Zr precipitates observed in this alloy based on the addition of a low melting point element, such as Sn. Characterization verified the first-principles density functional theory prediction that Zr and Sn atoms cluster during homogenization owing to the favorable binding energy of Zr-Sn-vacancy triplets. Direct microstructural observations revealed these clusters form Sn nanoprecipitates that assist the nucleation of L1 2 -Al 3 Zr at 200°C, where L1 2 -Al 3 Zr precipitation is not expected due to the low diffusivity of Zr atoms in Al. At higher temperatures (≳350°C), the acceleration of L1 2 -Al 3 Zr precipitation is driven by faster Zr diffusion in Al with Sn microalloying and the nuclei formed via the LEAN mechanism. In conclusion, this combination of mechanisms explains the improvement in age hardening through L1 2 -Al 3 Zr precipitation with Sn microalloying.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Probing the Effects of Dimethyl Aluminum Isopropoxide Surface Reaction Byproducts on Atomic Layer Deposition Nucleation

Atomic layer deposition (ALD) processes that leverage a myriad of metal–organic and complementary reactant combinations have been identified to realize precise and conformal thin film growth. However, the effects of the ALD reaction byproducts on nucleation and growth mechanisms are rarely considered. Site-selective ALD processes provide an opportunity for the detailed investigation of uniform surface sites with atomistic accuracy. Intentional pretreatment with a known ALD reaction byproduct – isopropanol – enables a significant improvement in the nucleation rate reproducibility of dimethylaluminum isopropoxide and water ALD on rutile TiO 2 (110). In situ spectroscopic ellipsometry reveals a partially reversible byproduct binding that is site-selective for TiO 2 (110) surface oxygen vacancies. First-principles calculations reveal surface site-specific thermodynamics for adsorption of isopropanol and water that may influence ALD nucleation. In conclusion, the sensitivity of site-selective ALD motivates consideration of secondary surface reactions when designing precision deposition processes, including area- or site-selective ALD reactions.

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

The effect of gamma ray irradiation on few layered MoSe2: A material for nuclear and space applications

In recent years, emerging two-dimensional (2D) materials, such as molybdenum diselenide (MoSe2), have been at the center of attention for many researchers. This is due to their unique and fascinating physicochemical properties that make them attractive in space and defense applications that include shielding harsh irradiation environments. In this study, we examined the effects of gamma (γ) rays at various doses on the structural, chemical, and optical properties of MoSe2 layers. After the samples were exposed to intense gamma radiation (from a 60Co source) with various exposure times to vary the total accumulated dosage (up to 100 kGy), Raman and photoluminescence spectroscopies were used to study and probe radiation-induced changes to the samples. When compared to pristine materials, very few changes in optical properties were typically observed, indicating good robustness with little sensitivity, even at relatively high doses of gamma radiation. The imaging using scanning electron microscopy revealed a number of nano-hillocks that were connected to substrate alterations. X-ray photoelectron spectroscopies revealed that Mo’s binding energies remained the same, but Se’s binding energies blueshifted. We associated this shift with the decrease in Se vacancies that occurred after irradiation as a result of Mo atoms creating adatoms next to Se atoms. When compared to pristine materials, very few changes in optical, chemical, and structural properties were typically observed. These findings highlight the inherent resilience of MoSe2 in hostile radioactive conditions, which spurs additional research into their optical, electrical, and structural characteristics as well as exploration for potential space, energy, and defense applications.

Materials Science↗
DOE OSTI2025

Energy Barriers for H 2 Uptake and Release by Defects in Boron Nitride

Density functional theory simulations (DFT) have been used in combination with the climbing image nudged elastic band method (CI-NEB) to determine the energy barriers required to dissociate H 2 at four types of h-BN defects and evaluate the potential for reacted H species to hop to adjacent available sites. The defects were selected based on previous work in which ab initio thermodynamics was used to identify specific defects with optimal thermodynamics for reversible H 2 storage, that is, the defects would dissociate molecular H 2 while not binding H species too strongly so that hydrogen release is favorable. The four defects investigated are the N and B monovacancies, which only offer either B or N defective sites respectively, as well as larger cluster vacancy defects named 3V(1B2N) and hexagonal 6V(3B3N), which present a mixture of B and N defect sites available for H binding. While homogeneous H 2 dissociation by B and N monovacancies have energy barriers larger than 0.4 eV, we found that mixed BN terminated defects have energy barrier lower than 0.4 eV for H 2 dissociation. In particular, H 2 dissociation by the hexagonal 6V(3B3N) defect, exclusively made of Frustrated Lewis Pairs (FLPs), is almost barrierless. The energy barrier for subsequent H hopping to available adjacent B or N sites of the defect were also calculated. It was found that the hopping of H species generally involves high energy barriers ranging from 0.26 eV to 2.57 eV. In conclusion, this work completes previous thermodynamic theoretical investigations and provides a more comprehensive picture of H 2 reaction with defective h-BN for potential H 2 storage and recovery.

Sassi, Michel J. [Pacific Northwest National Labor↗
DOE OSTI2026

The presence of Ru metal modifies the behavior of deuterium bound to pyridinic nitrogen in doped graphene-like materials

We used a combination of x-ray photoelectron spectroscopy (XPS) and density functional theory (DFT) to investigate how the presence of a ruthenium (Ru) metal substrate modifies the N-D bond strength of the pyridinic nitrogen species in nitrogen (N)-doped graphene (Gr) on Ru(0001) compared to metal-free highly-oriented pyrolytic graphite (HOPG) substrate. The N-dopants were introduced through low-energy N+/N2+ irradiation. N is embedded in the carbon materials in two predominant configurations: as graphitic N (GN: N substituted in the hexagonal C lattice) and pyridinic N (PN: substitutional N adjacent to a C vacancy). XPS data showed that upon atomic deuterium (D) exposure at 220 K, the PN peak shifted to a higher binding energy by +1.2 eV for HOPG and +1.0 eV for N-doped graphene on Ru(0001), while the GN peak remained unchanged, indicating that the D atoms bound solely to pyridinic N. 85% of PN sites on HOPG can be saturated with D atoms, whereas only ~30% of pyridinic N sites are able to bind D atoms in N-doped Gr/Ru(0001). Our DFT analysis shows that this difference is due to the coordination of PN to Ru atoms, which necessitates bond cleavage of the N-Ru interaction prior to D atom adsorption. D begins desorbing from N-HOPG at ~573 K and is fully desorbed by ~973 K, whereas desorption from N-doped graphene on Ru(0001) begins at ~290 K, with complete desorption observed at ~700 K, indicating that the Ru metal weakens the N-D bond strength. We also studied a high-surface-area, layered, and porous N-doped carbon material, internally labeled NC900, which was synthesized by pyrolysis of glucose and graphitic carbon nitride (g-C3N4) at 1173 K. D exposure caused a +1.1 eV PN peak shift in NC900, with no change in the GN peak. Moreover, the NC900 exhibited the same desorption behavior as HOPG, demonstrating that well-defined model systems can effectively capture the behavior of more complex N-doped carbon materials.

Alupothe Gedara, Buddhika S.↗
DOE OSTI2025

Defect Complexes in CrSBr Revealed Through Electron Microscopy and Deep Learning

Atomic defects underpin the properties of van der Waals materials, and their understanding is essential for advancing quantum and energy technologies. Scanning transmission electron microscopy is a powerful tool for defect identification in atomically thin materials, and extending it to multilayer and beam-sensitive materials would accelerate their exploration. Here, we establish a comprehensive defect library in a bilayer of the magnetic quasi-1D semiconductor CrSBr by combining atomic-resolution imaging, deep learning, and calculations. We apply a custom-developed machine learning work flow to detect, classify, and average point vacancy defects. This classification enables us to uncover several distinct Cr interstitial defect complexes, combined Cr and Br vacancy defect complexes, and lines of vacancy defects that extend over many unit cells. We show that their occurrence is in agreement with our computed structures and binding energy densities, reflecting the intriguing layer interlocked crystal structure of CrSBr. Our ab initio calculations show that the interstitial defect complexes give rise to highly localized electronic states. These states are of particular interest due to the reduced electronic dimensionality and magnetic properties of CrSBr and are, furthermore, predicted to be optically active. Our results broaden the scope of defect studies in challenging materials and reveal new defect types in bilayer CrSBr that can be extrapolated to the bulk and to over 20 materials belonging to the same FeOCl structural family.

deep learning↗
DOE OSTI2025

Nitrate-To-Ammonia Electroconversion at Neutral pH on Polycrystalline Vanadium Sulfide Derived from Vanadium Disulfide

The electrochemical nitrate reduction reaction (NO3RR) offers a pathway to produce NH3 for fuel and fertilizer from waste NO3-. In this work, a polycrystalline vanadium sulfide (VSx), which is derived from solvothermally grown and annealed VS2, is shown to exhibit excellent NO3RR activity (2.3 +- 0.6 mg.cm-2 geo..h-1 @ -0.92 VRHE) and Faradaic efficiency to NH4+ (69 +- 6% at -0.69 VRHE) in buffered neutral pH electrolyte containing 0.1 M NO3-. A variety of characterization techniques are leveraged to support the VSx assignment, including X-ray photoelectron spectroscopy, near-edge X-ray absorption fine structure spectroscopy, selected area electron diffraction, and X-ray diffraction measurements. The VS2 annealing step reduces the oxide character and generates VSx, which, based on the improved NO3RR activity, results in the creation of active sites for NO3- binding. To help shed light on NO3RR on VSx, VS2 is used as a model system, and a grand-canonical density functional theory (GC-DFT) investigation of VS2 shows strong evidence that S vacancies are active sites for NO3RR, where NO3- outcompetes H+ for adsorption at the S-vacancy sites. Moreover, GC-DFT results highlight a thermodynamically favorable reaction to generate NH4+ in an aqueous electrolyte at relevant cathodic potentials. As an annealed material, VSx may contain undersaturated V sites, which show an electronic structure similar to the theoretically calculated S-vacancy site of VS2, and these sites may contribute to the observed increase in NO3RR activity and selectivity for NH4+ on VSx versus unannealed VS2. Finally, kinetic isotope effect measurements suggest that the kinetic rate-limiting step of the NO3RR on VSx is not proton-coupled, indicating it may be the first electron transfer to adsorbed NO3*.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗