Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “INTERSTITIALS”

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 415 records · Page 23

Mechanism of vacancy formation induced by hydrogen in tungsten

We report a hydrogen induced vacancy formation mechanism in tungsten based on classical molecular dynamics simulations. We demonstrate the vacancy formation in tungsten due to the presence of hydrogen associated directly with a stable hexagonal self-interstitial cluster as well as a linear crowdion. The stability of different self-interstitial structures has been further studied and it is particularly shown that hydrogen plays a crucial role in determining the configuration of SIAs, in which the hexagonal cluster structure is preferred. Energetic analysis has been carried out to prove that the formation of SIA clusters facilitates the formation of vacancies. Such a mechanism contributes to the understanding of the early stage of the hydrogen blistering in tungsten under a fusion reactor environment.

36 MATERIALS SCIENCE↗

Structural stability and defect energetics of ZnO from diffusion quantum Monte Carlo

We have applied the many-body ab initio diffusion quantum Monte Carlo (DMC) method to study Zn and ZnO crystals under pressure and the energetics of the oxygen vacancy, zinc interstitial, and hydrogen impurities in ZnO. We show that DMC is an accurate and practical method that can be used to characterize multiple properties of materials that are challenging for density functional theory (DFT) approximations. DMC agrees with experimental measurements to within 0.3 eV, including the band-gap of ZnO, the ionization potential of O and Zn, and the atomization energy of O{sub 2}, ZnO dimer, and wurtzite ZnO. DMC predicts the oxygen vacancy as a deep donor with a formation energy of 5.0(2) eV under O-rich conditions and thermodynamic transition levels located between 1.8 and 2.5 eV from the valence band maximum. Our DMC results indicate that the concentration of zinc interstitial and hydrogen impurities in ZnO should be low under n-type and Zn- and H-rich conditions because these defects have formation energies above 1.4 eV under these conditions. Comparison of DMC and hybrid functionals shows that these DFT approximations can be parameterized to yield a general correct qualitative description of ZnO. However, the formation energy of defects in ZnO evaluated with DMC and hybrid functionals can differ by more than 0.5 eV.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dissociative diffusion mechanism in vacancy-rich materials according to mass action kinetics

We conducted two sets of diffusion-reaction numerical simulations using a finite difference method (FDM) in order to investigate fast impurity diffusion via interstitial sites in vacancy-rich materials such as Cu(In,Ga)Se 2 (CIGS) and Cu 2 ZnSn(S, Se) 4 (CZTSSe or CZTS) via the dissociative diffusion mechanism where the interstitial diffuser ultimately reacts with a vacancy to produce a substitutional. The first set of simulations extends the standard interstitial-limited dissociative diffusion theory to vacancy-rich material conditions where vacancies are annihilated in large amounts, introducing non-equilibrium vacancy concentration profiles. The second simulation set explores the vacancy-limited dissociative diffusion where impurity incorporation increases the equilibrium vacancy concentration. In addition to diffusion profiles of varying concentrations and shapes that were obtained in all simulations, some of the profiles can be fitted with the constant- and limited-source solutions of Fick’s second law despite the non-equilibrium condition induced by the interstitial-vacancy reaction. The first set of simulations reveals that the dissociative diffusion coefficient in vacancy-rich materials is inversely proportional to the initial vacancy concentration. In the second set of numerical simulations, impurity-induced changes in the vacancy concentration lead to distinctive diffusion profile shapes. The simulation results are also compared with published data of impurity diffusion in CIGS. And according to the characteristic properties of diffusion profiles from the two set of simulations, experimental detection of the dissociative diffusion mechanism in vacancy-rich materials may be possible.

36 MATERIALS SCIENCE↗

The role of water vapor during the synthesis of hydrogen doped In 2 O 3

Hydrogen-doped In 2 O 3 synthesized using magnetron sputtering or atomic layer deposition in the presence of water vapor has high transparency in the near infra-red region and mobility values that are almost three to four times of those in commercially used In 2 O 3 :Sn at a much lower carrier concentration. However, simple questions like what happens when water molecules enter In 2 O 3 or why the dominant charge carrier changes from an oxygen vacancy in In 2 O 3 :Sn to H + in In 2 O 3 are not clear. Using hybrid functional based density functional theory calculations, in this study we show that water molecules spontaneously split into H + and OH – . The H + interstitial bonds with an oxygen in In 2 O 3 , and the OH – interstitial occupies a vacant lattice site or is trapped by an oxygen vacancy. Thus, in agreement with experimental findings, our results suggest that the concentration of oxygen vacancies decreases as the partial pressure of water increases. The charge transition level, ϵ ( – / + ) , of a hydrogen interstitial and the Fermi level at which the formation energies of OH – and H + are equal lie above the conduction band minimum, suggesting that the H + interstitial has a lower formation energy. In addition, water molecules do not exhibit any propensity to form clusters, meaning that H + ions can diffuse and remain spatially separated from each other.

36 MATERIALS SCIENCE↗

First-principles study of intrinsic point defects and Xe impurities in uranium monocarbide

Based on density functional theory (DFT) calculations, we perform an extensive investigation of intrinsic point defects and Xe impurities in uranium monocarbide (UC). The DFT calculations involve both the conventional generalized gradient approximation (GGA) and the GGA+U approach with the Hubbard parametric term (U), using up to 5×5×5 supercells. GGA calculations for the formation energy of intrinsic defects demonstrate the significant effect of using larger supercells than in previous studies. Results confirm that the ⟨111⟩ and ⟨100⟩ dumbbell interstitials are the most stable interstitial configurations for U and C, respectively. The interstitial mechanisms are favored for self-diffusion of both uranium and carbon and diffusion of Xe under equilibrium conditions. Calculations also reveal that the Xe substitutional defect at the C lattice site tends to adopt an off-site configuration, which can be interpreted as a Xe interstitial–C vacancy complex. We also utilize GGA+U to assess the impact of effective U parameter (Ueff) on the results. Moreover, we introduce a method to estimate the carbon chemical potential by fitting the phase diagram composition data and propose a selection Ueff=1.25 eV based on the experimental Xe diffusion activation energy. With this approach, GGA+U calculations reproduce the available experimental data for the formation energy of the carbon Frenkel pair and can explain the stepwise recovery of intrinsic properties and burst Xe release behavior in UC observed in annealing experiments.

Huang, Gui-Yang (ORCID:0000000301447376)↗

Phosphorus diffusion and deactivation during SiGe oxidation

Dopant profiles near the semiconductor–oxide interface are critical for microelectronic device performance. As the incorporation of Si 1−x Ge x into transistors continues to increase, it is necessary to understand the behavior of dopants in Si 1−x Ge x . In this paper, the diffusion and electrical activation of phosphorus within a strained, single-crystal Si 0.7 Ge 0.3 layer on Si during oxidation are reported. Both layers were uniformly doped, in situ, with an average phosphorus concentration of 4 × 10 19 atoms/cm 3 . After high-temperature oxidation, secondary ion mass spectrometry measurements revealed that the bulk of the phosphorus diffuses out of only the SiGe layer and segregates at the oxidizing SiGe–SiO 2 interface. Hall effect measurements corroborate the observed phosphorus loss and show that the phosphorus diffusing to the oxidizing interface is electrically inactive. Through density functional theory (DFT) calculations, it is shown that phosphorus interstitials prefer sites near the SiGe–SiO 2 interface. Finally, based on a combination of experimental data and DFT calculations, we propose that the phosphorus atoms are displaced from their lattice sites by Ge interstitials that are generated during SiGe oxidation. The phosphorus atoms then migrate toward the SiGe–SiO 2 interface through a novel mechanism of hopping between Ge sites as P–Ge split interstitials. Once they reach the interface, they are electrically inactive, potentially in the form of interstitial clusters or as part of the reconstructed interface or oxide.

Thornton, Chappel S.↗

Identifying native point defect configurations in α-alumina

Intimately intertwined atomic and electronic structures of point defects govern diffusion-limited corrosion and underpin the operation of optoelectronic devices. For some materials, complex energy landscapes containing metastable defect configurations challenge first-principles modeling efforts. Here, we thoroughly reevaluate native point defect geometries for the illustrative case of α-Al 2 O 3 by comparing three methods for sampling candidate geometries in density functional theory calculations: displacing atoms near a naively placed defect, initializing interstitials at high-symmetry points of a Voronoi decomposition, and Bayesian optimization. We find symmetry-breaking distortions for oxygen vacancies in some charge states, and we identify several distinct oxygen split-interstitial geometries that help explain literature discrepancies involving this defect. We also report a surprising and, to our knowledge, previously unknown trigonal geometry favored by aluminum interstitials in some charge states. Importantly, these new configurations may have transformative impacts on our understanding of defect migration pathways in aluminum-oxide scales protecting metal alloys from corrosion. Overall, the Voronoi scheme appears most effective for sampling candidate interstitial sites because it always succeeded in finding the lowest-energy geometry identified in this study, although no approach found every metastable configuration. Finally, we show that the position of defect levels within the band gap can depend strongly on the defect geometry, underscoring the need to conduct careful searches for ground-state geometries in defect calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Thermodynamic Modeling of Complex Solid Solutions in the Lu-H-N System via Graph Neural Network Accelerated Monte Carlo Simulations

Metal hydrides are important across diverse applications, such as hydrogen storage, batteries, gas sensors, nuclear reactions, and high-temperature superconductivity. Previous computational studies of metal hydrides under extreme pressures, e.g., 𝑂⁡(10 2 ) ⁢GPa, usually treat them as stoichiometric compounds without considering interstitial lattice disorder. As pressures become more moderate in the 𝑂⁡(10 0 ) ⁢GPa and below range, hydrogen disorder at interstitial lattice sites becomes prominent, e.g., in the N-doped Lu hydride that was recently claimed superconducting near 1 GPa. Further adding compositional complexity from alloying and/or multielement interstitial occupation makes elucidating pressure- and temperature-dependent observables intractable by first-principles calculations alone. We therefore propose a lattice graph neural-network surrogate modeling approach to predict configuration- and pressure-dependent equation-of-state properties. Their efficiency permits Monte Carlo simulations to calculate Gibbs energies and pressure-dependent phase diagrams, thereby revealing insights into the synthesis conditions required for achieving desired phase equilibria. We demonstrate this concept for the compositionally complex cubic Lu(H,N,Va) 3 system where three constituents (hydrogen, nitrogen and vacancy) have disordered multielement interstitial occupancies and insights into pressure-dependent phase equilibria are critically needed, e.g., N-doping levels can significantly lower dehydrogenation temperatures and provide a new strategy to optimize hydrogen-storage alloys. This work can improve the thermodynamic understanding of the Lu-H-N system and help rational synthesis of N-doped Lu hydrides, but more generally demonstrates an efficient approach to model pressure-dependent thermodynamics of multicomponent solid solutions.

Monte Carlo methods↗

Natural off-stoichiometry causes carrier doping in half-Heusler filled tetrahedral structures

The half-Heusler filled tetrahedral structures (FTSs) are zinc-blende-like compounds, where an additional atom is filling its previously empty interstitial site. The FTSs having 18 valence electrons per formula unit are an emerging family of functional materials, whose intrinsic doping trends underlying a wide range of electronic functionalities are yet to be understood. Interestingly, even pristine compounds without any attempt at impurity/chemical doping exhibit intriguing trends in the free carriers they exhibit. Applying the first principles theory of doping to a few prototype compounds in the A IV B X C IV and A IV B IX C V groups, we describe the key ingredients controlling the materials’ propensity for both intrinsic and extrinsic doping: (a) The spontaneous deviations from 1:1:1 stoichiometry reflect predictable thermodynamic stability of specific competing phases. (b) Bulk ABC compounds containing 3d elements in the B position (ZrNiSn and ZrCoSb) are predicted to be naturally 3 d rich. The B = 3d interstitials are the prevailing shallow donors, whereas the potential acceptors (e.g., Zr vacancy and Sn-on-Zr antisite) are ineffective electron killers, resulting in an overall uncompensated n-type character, even without any chemical doping. In these materials, the band edges are “natural impurity bands” due to non-Daltonian off-stoichiometry, such as B interstitials, not intrinsic bulk controlled states as in a perfect crystal. (c) Bulk ABC compounds containing 5d elements in the B position (ZrPtSn, ZrIrSb, and TaIrGe) are predicted to be naturally C rich and A poor. This promotes the hole-producing C-on- A antisite defects rather than B-interstitial donors. The resultant p-type character (without chemical doping) therein is “latent” for C = Sn and Sb; however, as the C-on-A hole-producing acceptors are rather deep and p typeness is manifest only at high temperature or via impurity doping. In contrast, in TaIrGe (B = Ir, 5d), the prevailing hole-producing Ge-on-Ta antisite (C-on- A) is shallow, making it a real p-type compound. This general physical picture establishes the basic trends of carriers in this group of materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Dislocation loop evolution in Kr-irradiated ThO 2

Here, the early stage of microstructural evolution of ThO 2 , under krypton irradiation at 600, 800, and 1000°C, was investigated using in situ transmission electron microscopy (TEM). Dislocation loops grew faster, whereas their number density decreased with increasing irradiation temperature. Loop density was found to decrease with ion dose. Interstitial dislocation loops, including Frank loops with Burgers vector of $\textit{a}/3\langle111\rangle$ and perfect loops with Burgers vector of $\textit{a}/2\langle110\rangle$, were determined by traditional TEM and atomic resolution–scanning TEM techniques. Atomistic and mesoscale level modeling are performed to interpret experimental observations. The migration energy barriers of defects in ThO 2 were calculated using density-functional theory. The energetics of different dislocation loop types were studied using molecular dynamics simulations. Loop density and diameter were analyzed using a kinetic rate theory model that considers stoichiometric loop evolution. This analysis reveals that loop growth is governed by the mobility of cation interstitials, whereas loop nucleation is determined by the mobility of anion defects. Lastly, a rate theory model was used to extract the diffusion coefficients of thorium interstitials, oxygen interstitials, and vacancies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Ab initio prediction of rapid kinetics of Fe impurities in δ -Pu

Here, we study the formation energies of iron impurities in δ-Pu within spin–orbital-polarized density functional theory (SOP-DFT). The thermodynamic solubility limit of iron in δ-Pu is calculated, indicating low miscibility. We show that surprisingly, Fe impurities at equilibrium are almost equally likely to occupy octahedral interstitial sites or substitutional sites, with slight preference for the former. In contrast, we find the energy of the tetrahedral interstitial Fe to be nearly 1 eV higher than the octahedral one. We explore the energy landscape for Fe impurity hopping diffusion and conclude that Fe impurities in δ-Pu are divided into two populations: (i) Immobile substitutional Fe impurities and (ii) highly mobile interstitial Fe impurities. The latter, (ii), migrate between octahedral interstitial sites with an energy barrier of around 0.2 eV. The energy barrier for exchange between the two populations is calculated to exceed 0.7 eV. Finally, we discuss the role of magnetic order on the impurity energetics.

Atomic structure↗

Migration of defect clusters and xenon-vacancy clusters in uranium dioxide

The possible transition states, minimum energy paths and migration mechanisms of defect clusters and xenon-vacancy defect clusters in uranium dioxide have been investigated using the dimer and the nudged elastic-band methods. The nearby O atom can easily hop into the oxygen vacancy position by overcoming a small energy barrier, which is much lower than that for the migration of a uranium vacancy. A simulation for a vacancy cluster consisting of two oxygen vacancies reveals that the energy barrier of the divacancy migration tends to decrease with increasing the separation distance of divacancy. For an oxygen interstitial, the migration barrier for the hopping mechanism is almost three times larger than that for the exchange mechanism. Xe moving between two interstitial sites is unlikely a dominant migration mechanism considering the higher energy barrier. A net migration process of a Xe-vacancy pair containing an oxygen vacancy and a xenon interstitial is identified by the NEB method. We expect the oxygen vacancy-assisted migration mechanism to possibly lead to a long distance migration of the Xe interstitials in UO2. The migration of defect clusters involving Xe substitution indicates that Xe atom migrating away from the uranium vacancy site is difficult.

Defect clusters↗

Development of a robust modeling tool for radiation-induced segregation in austenitic stainless steels

Irradiation-assisted stress corrosion cracking (IASCC) of austenitic stainless steels in Light Water Reactor (LWR) components has been linked to changes in grain boundary composition due to irradiation induced segregation (RIS). This work developed a robust RIS modeling tool to account for thermodynamics and kinetics of the atom and defect transportation under combined thermal and radiation conditions. The diffusion flux equations were based on the Perks model formulated through the linear theory of the thermodynamics of irreversible processes. Both cross and non-cross phenomenological diffusion coefficients in the flux equations were considered and correlated to tracer diffusion coefficients through Manning’s relation. The preferential atom-vacancy coupling was described by the mobility model, whereas the preferential atom-interstitial coupling was described by the interstitial binding model. The composition dependence of the thermodynamic factor was modeled using the CALPHAD approach. Detailed analysis on the diffusion fluxes near and at grain boundaries of irradiated austenitic stainless steels suggested the dominant diffusion mechanism for chromium and iron is via vacancy, while that for nickel can swing from the vacancy to the interstitial dominant mechanism. The diffusion flux in the vicinity of a grain boundary was found to be greatly influenced by the composition gradient formed from the transient state, leading to the oscillatory behavior of alloy compositions in this region. This work confirms that both vacancy and interstitial diffusion, and segregation itself, have important roles in determining the microchemistry of Fe, Cr, and Ni at irradiated grain boundaries in austenitic stainless steels.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Investigating the Role of Defects in α-Pu Phase Transitions using Density Functional Theory

Before reaching its relatively low melting point of 913 K, plutonium undergoes six solid state phases: simple monoclinic alpha, body-centered monoclinic beta, face-centered orthorhombic gamma, face centered cubic delta, body-centered tetragonal delta-prime, and body-centered cubic epsilon phases. Due to their low symmetry, the monoclinic alpha and beta phases have not been studied extensively, however it has been previously hypothesized that the transition between alpha and beta is caused by defects in the lattices due to the incommensurate number of atoms. Additionally, the martensitic delta to alpha phase transformation in plutonium-gallium alloys induces an expanded alpha lattice known as alpha prime, most likely due to the insoluble gallium impurity atom. This work explores how defects such as plutonium interstitials and gallium interstitials and substitutionals may affect the alpha plutonium lattice, providing a better understanding of the mechanisms for the alpha-beta and martensitic phase transformations through density functional theory (DFT). Results indicate that plutonium interstitials within the alpha lattice induce beta features structurally and electronically, such as volume expansion and bonding changes, and gallium interstitials and substitutionals within the alpha lattice induce volume changes as observed experimentally in alpha prime.

36 MATERIALS SCIENCE↗

Modeling of radiation damage in silicon solar cells

One MeV electron irradiation produces preponderantly isolated vacancy interstitial pairs. If neither of these defects is mobile, the concentration of each grows linearly with fluence. Annealing of damage depends on the nature of the damage. Vacancy interstitial pairs which are bound by an interaction such that they mutually annihilate rather than dissociate are termed close pairs; close pair recovery usually occurs at a lower temperature than the temperature at which long distance defect migration occurs. Annealing of the remaining frozen in damage occurs when a temperature is reached where the vacancy or interstitial is mobile; usually the interstitial is more mobile than the vacancy. The recovery occurs in two regimes which may be resoluable.

Oehrlein, G.↗

Radiation damage and defect behavior in proton irradiated lithium-counterdoped n+p silicon solar cells

Two lithium-counterdoped n+p silicon solar cells with different lithium concentrations were irradiated by 10-MeV protons. Cell performance was measured as a function of fluence, and it was found that the cell with the highest concentration of lithium had the highest radiation resistance. Deep level transient spectroscopy which showed two deep level defects that were lithium related. Relating the defect energy levels obtained from this study with those from earlier work using 1-MeV electron irradiation shows no correlation of the defect energy levels. There is one marked similarity: the absence of the boron-interstitial-oxygen-interstitial defect. This consistency strengthens the belief that lithium interacts with oxygen to prevent the formation of the boron interstitial-oxygen interstitial defect. The results indicate that, in general, addition of lithium in small amounts to the p-base of a boron doped silicon solar cell such that the base remains p-type, tends to increase the radiation resistance of the cell.

Stupica, John↗

Physiological water model development

The water of the human body can be categorized as existing in two main compartments: intracellular water and extracellular water. The intracellular water consists of all the water within the cells and constitutes over half of the total body water. Since red blood cells are surrounded by plasma, and all other cells are surrounded by interstitial fluid, the intracellular compartment has been subdivided to represent these two cell types. The extracellular water, which includes all of the fluid outside of the cells, can be further subdivided into compartments which represent the interstitial fluid, circulating blood plasma, lymph, and transcellular water. The interstitial fluid surrounds cells outside of the vascular system whereas plasma is contained within the blood vessels. Avascular tissues such as dense connective tissue and cartilage contain interstitial water which slowly equilibrates with tracers used to determine extracellular fluid volume. For this reason, additional compartments are sometimes used to represent these avascular tissues. The average size of each compartment, in terms of percent body weight, has been determined for adult males and females. These compartments and the forces which cause flow between them are presented. The kidneys, a main compartment, receive about 25 percent of the cardiac output and filters out a fluid similar to plasma. The composition of this filtered fluid changes as it flows through the kidney tubules since compounds are continually being secreted and reabsorbed. Through this mechanism, the kidneys eliminate wastes while conserving body water, electrolytes, and metabolites. Since sodium accounts for over 90 percent of the cations in the extracellular fluid, and the number of cations is balanced by the number of anions, considering the renal handling sodium and water only should sufficiently describe the relationship between the plasma compartment and kidneys. A kidney function model is presented which has been adapted from a previous model of normal renal function in man. To test the validity of the proposed kidney model, results predicted by the model will be compared to actual data involving injected or ingested fluids and subsequent urine flow rates. Comparison of the model simulation to actual data following the ingestion of 1 liter of water is shown. The model simulation is also shown with actual data following the intravenous infusion of hypertonic saline.

Doty, Susan↗

The Influence of C and Si on the Flow Behavior of NiAl Single Crystals

Alloys based on the intermetallic compound NiAl are considered potential replacements for Ni and Co-based superalloys in high temperature structural applications due to their excellent oxidation resistance, low densities, high thermal conductivities, and increased melting points. Unfortunately, NiAl exhibits low tensile ductility at room temperature and low strengths at elevated temperatures which have combined to hinder its development. Recent efforts, have revealed that NiAl in the presence of sufficient solute levels, is subject to the phenomenon of strain aging which manifests itself as: sharp yield points, abnormally low strain rate sensitivities (SRS), plateaus or peaks in yield stress and work hardening rate as a function of temperature, flow stress transients upon an upward change in strain rate, reduced tensile elongations at elevated temperatures, and serrated stress-strain curves. Though recent efforts via either alloying or the removal of interstitial impurities, have resulted in consistent room-temperature tensile elongations exceeding 5% and the elimination of serrated flow, the effects of particular substitutional and interstitial elements and the mechanisms by which they might enhance or hinder the mechanical properties remain unknown. Consequently, the purpose of the present paper is to provide a preliminary assessment of the influence of common substitutional and interstitial impurities on the deformation behavior of NiAl. To accomplish this goal a series of NiAl single crystal alloys containing various interstitial solutes were prepared and their mechanical properties were evaluated between 77 and 1100 K. Because Si is a common impurity in conventional purity single crystals grown by the Bridgman method, Si concentrations were also varied in order to determine the influence of this element.

Weaver, M. L.↗