Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Vacancy diffusion”

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 55 records · Page 3

Phase Transition in a Memristive Suspended MoS 2 Monolayer Probed by Opto- and Electro-Mechanics

Semiconducting monolayers of a 2D material are able to concatenate multiple interesting properties into a single component. Here, by combining opto-mechanical and electronic measurements, we demonstrate the presence of a partial 2H–1T' phase transition in a suspended 2D monolayer membrane of MoS 2 . Electronic transport shows unexpected memristive properties in the MoS 2 membrane, in the absence of any external dopants. A strong mechanical softening of the membrane is measured concurrently and may only be related to the 2H–1T' phase transition, which imposes a 3% directional elongation of the topological 1T' phase with respect to the semiconducting 2H. We note that only a few percent 2H–1T' phase switching is sufficient to observe measurable memristive effects. Our experimental results combined with first-principles total energy calculations indicate that sulfur vacancy diffusion plays a key role in the initial nucleation of the phase transition. Our study clearly shows that nanomechanics represents an ultrasensitive technique to probe the crystal phase transition in 2D materials or thin membranes. Lastly, a better control of the microscopic mechanisms responsible for the observed memristive effect in MoS 2 is important for the implementation of future devices.

2D materials↗

Interfacial Reactivity-Triggered Oscillatory Lattice Strains of Nanoalloys

Understanding the structure evolution of nanoalloys under reaction conditions is vital to the design of active and durable catalysts. Herein, we report an operando measurement of the dynamic lattice strains of dual-noble-metal alloyed with an earth-abundant metal as a model electrocatalyst in a working proton-exchange membrane fuel cell using synchrotron high-energy X-ray diffraction coupled with pair distribution function analysis. Additionally, the results reveal an interfacial reaction-triggered oscillatory lattice strain in the alloy nanoparticles upon surface dealloying. Analysis of the lattice strains with an apparent oscillatory irregularity in terms of frequency and amplitude using time-frequency domain transformation and theoretical calculation reveals its origin from a metal atom vacancy diffusion pathway to facilitate realloying upon dealloying. This process, coupled with surface metal partial oxidation, constitutes a key factor for the nanoalloy's durability under the electrocatalytic oxygen reduction reaction condition, which serves as a new guiding principle for engineering durable or self-healable electrocatalysts for sustainable fuel cell energy conversion.

Wu, Zhi-Peng [State University of New York (SUNY),↗

On the role of the zirconia/ceria interface in the degradation of solid oxide electrolysis cells

Gd-doped Ceria (GDC) barrier layers are widely used in solid oxide electrolysis cells to prevent reactions between the Y-stabilized Zirconia (YSZ) electrolyte and oxygen electrode. In some cases, electrolysis operation causes fracture near the GDC/YSZ interface. Here, we explain this effect using a numerical calculation, showing that a peak oxygen pressure is observed near the GDC/YSZ interface during electrolysis operation. Moreover, GDC/YSZ interdiffusion and possible zirconate phase formation occurring during electrode firing yield a low oxygen ion conductivity region that introduces a much higher peak oxygen partial pressure at the interdiffusion layer/YSZ interface, thereby exacerbating degradation. The results indicate that the peak pressure can be reduced and, thereby, damage can be mitigated by altering processing in order to minimize interdiffusion, decreasing the GDC layer thickness, and reducing oxygen electrode polarization resistance or altering electrolyte transport properties, e.g., changing the Gd doping level of GDC to increase the oxygen vacancy diffusivity.

Physics↗

Enhancing the metal–insulator transition in VO 2 heterostructures with graphene interlayers

In this work, studies were conducted on epitaxial VO 2 thin films to assess to the effect of remote epitaxy on the metal–insulator transition (MIT). The epitaxial VO 2 heterostructures were synthesized on both bare Al 2 O 3 (0001) substrates and Al 2 O 3 substrates coated with two monolayer-thick graphene. While both systems exhibit the MIT, the film grown by remote epitaxy on graphene demonstrates improved transport properties. Electrical transport measurements show that the on/off ratio is enhanced by a factor of ~7.5 and the switching temperature window is narrower for VO 2 thin films grown on graphene. By characterizing the heterostructures with a suite of structural, chemical, and spectroscopic tools, we find that the graphene interlayer inhibits oxygen vacancy diffusion from Al 2 O 3 (0001) during the VO 2 growth, resulting in improved electrical behavior at the MIT.

36 MATERIALS SCIENCE↗

Development of interatomic potential and effect of ordering on defect properties in CrMnV

Developing materials that can withstand extreme environments, such as high radiation doses and elevated temperatures, is crucial for next-generation particle accelerators, including the 2.4 MW Long-Baseline Neutrino Facility. High-Entropy Alloys have emerged as promising candidates for beam window materials due to their superior mechanical strength, corrosion resistance, and radiation tolerance. In this study, we focus on the Cr–Mn–V alloy system, developing and employing machine-learning interatomic potentials (MLIPs) to investigate the formation of an ordered phase and its influence on defect properties. Using hybrid Monte Carlo-Molecular Dynamics simulations, we observe the formation of a B2-ordered phase at lower temperatures, consistent with Density Functional Theory (DFT) predictions. Ordered structures display a bimodal distribution of migration energies and reduced mean square displacement values, indicating suppressed vacancy diffusion. Our results also show that the migration energy barrier varies based on the atomic species, with Mn and V exhibiting the highest and lowest average barriers, respectively. These findings suggest that atomic ordering inhibits defect mobility, potentially enhancing the radiation resistance of CrMnV alloys. The validated MLIP provides a reliable framework for simulations that are faster than traditional DFT while maintaining the accuracy required to study defect and ordering properties.

36 MATERIALS SCIENCE↗

Measurement of Ehrlich-Schwoebel barrier contribution to the self-organized formation of ordered surface patterns on Ge(001)

Normal-incidence 1-keV Ar + ion bombardment leads to amorphization and ultrasmoothing of Ge at room temperature, but at elevated temperatures the Ge surface remains crystalline and is unstable to the formation of self-organized nanoscale patterns of ordered pyramid-shaped pits. The physical phenomenon distinguishing the high-temperature patterning from room-temperature ultrasmoothing is believed to be a surface instability due to the Ehrlich-Schwoebel barrier for diffusing vacancies and adatoms, which is not present on the amorphous material. This real-time grazing-incidence small-angle x-ray scattering study compares smoothing of a prepatterned Ge sample at room temperature with patterning of an initially flat Ge sample at an elevated temperature. In both experiments, when the nanoscale structures are relatively small in height, the average kinetics can be explained by a linear theory. Furthermore, the linear theory coefficients, indicating surface stability or instability, were extracted for both experiments. A comparison between the two measurements allows estimation of the contribution of the Ehrlich-Schwoebel barrier to the self-organized formation of ordered nanoscale patterns on crystalline Ge surfaces.

36 MATERIALS SCIENCE↗

Atomistic Mechanism Underlying the Si ( 111 ) - ( 7 × 7 ) Surface Reconstruction Revealed by Artificial Neural-Network Potential

The 7×7 reconstruction of the Si(111) surface represents arguably the most fascinating surface reconstruction so far observed in nature. Yet, the atomistic mechanism underpinning its formation remains unclear after it was discovered sixty years ago. Experimentally, it is observed post priori so that analysis of its formation mechanism can only be carried out in analogy with archaeology. Theoretically, density-functional theory (DFT) correctly predicts the Si(111)-(7×7) ground state but is impractical to simulate its formation process; while empirical potentials failed to produce it as the ground state. Developing an artificial neural-network potential of DFT quality, we carried out accurate large-scale simulations to unravel the formation of the Si(111)-(7×7) surface. We reveal a possible step-mediated atom-pop rate-limiting process that triggers massive nonconserved atomic rearrangements, most remarkably, a critical process of collective vacancy diffusion that mediates a sequence of selective dimer, corner-hole, stacking-fault, and dimer-line pattern formation, to fulfill the 7×7 reconstruction. Our findings may not only solve the long-standing mystery of this famous surface reconstruction but they also illustrate the power of machine learning in studying complex structures.

36 MATERIALS SCIENCE↗

Final-stage densification kinetics of direct current–sinteredZrB2

Final-stage sintering was analyzed for nominally phase pure zirconium diboride synthesized by borothermal reduction of high-purity ZrO2. Analysis was conducted on ZrB2 ceramics with relative densities greater than 90% using the Nabarro–Herring stress–directed vacancy diffusion model. Temperatures of 1900°C or above and an applied uniaxial pressure of 50 MPa were required to fully densify ZrB2 ceramics by direct current sintering. Ram travel data were collected and used to determine the relative density of the specimens during sintering. Specimens sintered between 1900 and 2100°C achieved relative densities greater than 97%, whereas specimens sintered below 1900°C failed to reach the final stage of sintering. The average grain size ranged from 1.0 to 14.7 μm. The activation energy was calculated from the slope of an Arrhenius plot that used the Kalish equation. The activation energy was 162 ± 34 kJ/mol, which is consistent with the activation energy for dislocation movement in ZrB2. The diffusion coefficients for dislocation motion that controls densification were 5.1 × 10−6 cm2/s at 1900°C and 5.1 × 10−5 cm2/s at 2100°C, as calculated from activation energy and average grain sizes. This study provides evidence that the dominant mechanism for final-stage sintering of ZrB2 ceramics is dislocation motion.

36 MATERIALS SCIENCE↗

The Critical Stack Pressure to Alter Void Generation at Li/Solid-Electrolyte Interfaces during Stripping

The lithium stripping process generates vacancies, which may accumulate as voids and lead to uneven current distribution and dendrite growth in the following plating cycles. A stack pressure is typically required during stripping, but how to optimize the stack pressure is not clear. In this work, extremely lithiophilic Li/Li 2 O and lithiophobic Li/LiF interfaces were used to reveal the combining effect of interface interaction and stack pressure induced lithium creep on the stripping critical current density (CCD). A multiscale simulation scheme with Density Functional Theory (DFT), kinetic Monte Carlo (KMC) simulations, and an analytical model was developed. The analytical model predicted lithiophobic interfaces require a higher stack pressure than lithiophilic interfaces to reach the same CCD. The KMC simulations also showed higher stack pressure is needed at lithiophobic interfaces to accelerate Li vacancy diffusion into the bulk and maintain a flat surface. This stack pressure needs to be high enough to alter the Li forward-and-backward hopping barriers at the interface. This multiscale simulation scheme illustrates the importance to include the chemical-mechanical effects during Li stripping morphology evolution. It can be used to design ideal interlayer coating materials to maintain a flat Li surface during cycling.

25 ENERGY STORAGE↗

Interpretation of Ion Irradiation and Neutron Irradiation Damage in Additively Manufactured 316 Stainless Steel using Multiscale Modeling

The accelerated adoption of nuclear energy necessitates advanced manufacturing technologies, such as additive manufacturing, to meet heightened supply chain requirements and support innovative reactor technologies. Due to the unique microstructural characteristics of additively manufactured materials under distinct solidification conditions, comprehensive evaluation of their performance in reactor environments is essential. The Advanced Materials and Manufacturing Technologies program under the Department of Energy's Office of Nuclear Energy focuses on understanding the irradiation performance and damage evolution of laser powder bed fusion 316 stainless steel, with an emphasis on integrating ion and neutron irradiation data to accelerate the development and qualification of materials for advanced nuclear reactor applications. While ion irradiation is a cost- and time-effective method, modeling and simulation are required to interpret the data for the broader range of irradiation conditions encountered in advanced reactors. In fiscal year 2025, integrated multiscale modeling and simulations were conducted to assess irradiation damage in additively manufactured 316 stainless steel. Key outcomes include predictions of chromium enrichment at grain boundaries, nickel enrichment at dislocation cell walls and void surfaces, and heterogeneous void evolution under ion and neutron irradiation conditions. Cluster dynamics simulations revealed the coarsening of voids at high irradiation temperatures and the suppression of void growth by high network dislocation density, while also demonstrating significant growth and coarsening of voids and self-interstitial atom loops at low dose rates. Machine learning-accelerated atomistic simulations highlighted the impact of the local environment and chromium concentration on vacancy diffusivity, providing key insights on the influence of composition on void swelling and radiation-induced segregation. Additionally, molecular dynamics simulations demonstrated the presence of defect production bias and a significant effect of carbon content on defect cluster behavior. These combined efforts aim to predict the performance of additively manufactured materials under various reactor conditions, supporting their qualification for nuclear reactor applications by interpreting ion irradiation data. This report underscores the potential of integrated multiscale modeling to analyze ion irradiation data in the effort to accelerate the qualification of additively manufactured materials for nuclear reactor components.

316 stainless steel↗

Sluggish, chemical bias and percolation phenomena in atomic transport by vacancy and interstitial diffusion in Ni-Fe alloys

In this paper, we report the results of extended atomistic modeling of intrinsic mobility of point defects and associated atomic transport in Ni-Fe model binary alloys. We consider the effects of composition and temperature and present evidence of the sluggish and chemically biased diffusion, and percolation effects occurring in atomic transport via the vacancy and interstitial migration mechanisms. The results are analyzed and discussed in the light of previous studies and some experimental observations. It is demonstrated that the sluggish diffusion, the chemically biased diffusion, and the percolation are interlinked phenomena that are defined by the chemical complexity of particular alloys. Methods for predicting these phenomena in multicomponent alloys are discussed.We report a fundamental understanding of sluggish diffusion, chemically-biased diffusion, as well as percolation phenomena, in Ni-Fe random alloys for vacancy and interstitial atom migration mechanisms.

36 MATERIALS SCIENCE↗

First-principles investigation of lanthanides diffusion in HCP zirconium via vacancy-mediated transport

The diffusion of lanthanide fission products plays an important role in the growth of the fuel-cladding chemical interaction (FCCI) region in metallic fuels. The use of a Zr interdiffusion barrier may mitigate the transport of lanthanides from the fuel to the cladding, but the efficacy of such a liner is not yet known. In this paper, the stability and vacancy-mediated diffusion of La, Ce, Pr, and Nd in hexagonal close-packed (HCP) Zr is investigated via density functional theory (DFT) calculations and self-consistent mean field (SCMF) analysis. DFT is used to calculate the formation, binding, and migration energies of vacancies and vacancy-solute pairs. The DFT energetics are used in the KineCluE code to calculate the Onsager transport coefficients. La is found to be the fastest diffusing species in HCP Zr and experiences an almost isotropic diffusion behavior. The other three species (Ce, Pr, and Nd) demonstrate anisotropic diffusion where the diffusion in the basal planes is significantly faster than that along the c-axis. Further, the calculated lanthanide diffusivities in HCP Zr are fitted to an Arrhenius relation and the activation energies and prefactors are reported for the first time. Furthermore, the vacancy drag and the segregation tendencies were analyzed using the calculated off-diagonal transport coefficients. According to our vacancy-mediated diffusion model, lanthanides are expected to be enriched at vacancy sinks at low temperatures, while at high temperatures, lanthanides are depleted at sinks and will preferably diffuse into the bulk. The enrichment/depletion transition temperature depends on the diffusion direction (basal or axial) and hence will be controlled by the grain texture and orientation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Influence of Hydrogen Isotopes on Vacancy Formation and Antisite Defect Diffusion in Palladium and Vanadium Metals

Density functional theory calculations are performed to study the interaction between hydrogen isotopes, vacancy, and antisite defects in Pd and V. Various defect configurations and defect-defect distances are systematically explored. Binding energies and migration energy barriers are calculated and tabulated. The study provides atomistic data for subsequent mesoscale simulations of tritium, vacancy, and antisite diffusion. In Pd, a strong binding is found between a vacancy and tritium (0.16 eV). This tritium binding increases the thermal concentration of vacancies by a factor of ~10 at 500 ?C. The increase in vacancy concentration enhances V antisite diffusion in Pd by a factor of ~6 at 500 ?C. The influence of tritium is even stronger in V, with a tritium and vacancy binding energy of 0.38 eV. Such a strong binding increases the thermal concentration of vacancies by a factor of ~300 at 500 ?C. The increase in vacancy concentration enhances Pd antisite diffusion in V by a factor of ~640 at 500 ?C. Vanadium and Pd exhibit a strong driving force to intermix with a formation energy of -1.57 eV for V antisite in Pd and -1.05 eV for Pd antisite in V. The results suggest vanadium diffusion into Pd is energetically stronger than the reverse. Zero-point-energy corrections are taken into account and calculations for hydrogen and tritium are presented.

Setyawan, Wahyu↗

Accelerated kinetic Monte Carlo: A case study; vacancy and dumbbell interstitial diffusion traps in concentrated solid solution alloys

Vacancy and self-interstitial atomic diffusion coefficients in concentrated solid solution alloys can have a non-monotonic concentration dependence. Here, the kinetics of monovacancies and $\langle$100$\rangle$ dumbbell interstitials in Ni–Fe alloys are assessed using lattice kinetic Monte Carlo (kMC). The non-monotonicity is associated with superbasins, which impels using accelerated kMC methods. Detailed implementation prescriptions for first passage time analysis kMC (FPTA-kMC), mean rate method kMC (MRM-kMC), and accelerated superbasin kMC (AS-kMC) are given. The accelerated methods are benchmarked in the context of diffusion coefficient calculations. The benchmarks indicate that MRM-kMC underestimates diffusion coefficients, while AS-kMC overestimates them. In this application, MRM-kMC and AS-kMC are computationally more efficient than the more accurate FPTA-kMC. Our calculations indicate that composition dependence of migration energies is at the origin of the vacancy’s non-monotonic behavior. In contrast, the difference between formation energies of Ni–Ni, Ni–Fe, and Fe–Fe dumbbell interstitials is at the origin of their non-monotonic diffusion behavior. Additionally, the migration barrier crossover composition—based on the situation where Ni or Fe atom jumps have lower energy barrier than the other one—is introduced. KMC simulations indicate that the interplay between composition dependent crossover of migration energy and geometrical site percolation explains the non-monotonic concentration-dependence of atomic diffusion coefficients.

36 MATERIALS SCIENCE↗

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↗

Bulk Diffusion of Cl through O Vacancies in α-Cr 2 O 3 : A Density Functional Theory Study

Many metals form passive oxides which are broken or weakened by anionic-induced degradation leading to material failure and anionic diffusion is an important step of this degradation process. The diffusion of anionic species through oxides involves a combination of diffusion along grain boundaries, cracks, and channels, as well as bulk diffusion via point defects which is the focus of this study. Using density functional theory, we study bulk diffusion of Cl through O vacancies in α -Cr 2 O 3 as a model system for passive metal oxides. Little is known about Cl diffusion and bonding characteristics, so we benchmark our work through comparison to numerous studies on O diffusion in α -Cr 2 O 3 to analyze similarities and differences between the O and Cl diffusion in the passivation and degradation process of α -Cr 2 O 3 respectively. Unlike O diffusion, the lowest diffusion barrier for Cl is cross-plane diffusion between two (0001) planes through a vacant cation site but the much shorter in-plane diffusion path within the same coordination polyhedron has 35% higher barrier. This work provides the basis for considering the contributions of Cl bulk-diffusion in the overall diffusion kinetics of Cl through α -Cr 2 O 3 .

Oware Sarfo, Kofi↗