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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The effect of composition on helium bubble nucleation in complex concentrated alloys

Molecular dynamics simulations of helium implantation have been performed in the MoNbTaTi complex concentrated alloy (CCA) to determine the impact of composition on helium clustering and bubble nucleation. It was found that compositions with lower interstitial atomic volume resulted in smaller cluster sizes, higher helium-to-vacancy ratios, and larger migration barriers. As the atomic volume decreases, it is harder to accommodate additional helium in the lattice and for helium to diffuse between interstitial sites, resulting in smaller cluster sizes. In particular, increasing molybdenum directly correlated with a decreased interstitial atomic volume and increased migration barrier. Niobium exhibited the opposite trend, where more niobium in the material resulted in larger interstitial atomic volumes and cluster sizes but lower migration barriers. The atomic volume was determined to be an indicator of the susceptibility of the material to helium bubble growth and could potentially be used as a metric to screen compositions for resistance to helium damage. In this way, properties such as available interstitial volumes could be used to screen materials more rapidly given the vast compositional space of CCAs.

McCarthy, Megan Jeanne [Sandia National Laboratori↗
NASA NTRS1978

Investigation of the fracture mechanism of Ti-5Al-2.5Sn at cryogenic temperatures

Fractography and metallographic sectioning were used to investigate the influence of microstructure on the fracture mechanism and fracture toughness (KIC) of normal interstitial and extra low interstitial (ELI) Ti-5Al-25Sn at 20 K (-423 F) and 77 K (-320 F). Plates of each grade were mill annealed at 815 C followed by either air or furnace cooling. These variations in composition and cooling rate resulted in differences in the volume fraction and internal structure of the dispersed beta phase and in the ordering of the alpha matrix. The ELI alloys were tougher than the normal interstitial plates. KIC of the furnace-cooled ELI plate was 25% lower than that of the air-cooled ELI material. Variations in cooling rate had no influence of KIC of the normal interstitial alloys. Fractography showed that a large portion of the fracture surfaces were covered with elongated dimples. Metallographic sections of specimens deformed at 77 K showed that these features form at the intersections of slip bands or deformation twins with grain or twin boundaries. Ordering and higher interstitial levels increase the local strain in slip bands resulting in void nucleation at lower macroscopic strains and lower KIC values.

Van Stone, R. H.↗
NASA NTRS1983

Lithification opf gas-rich chondrite regolith breccias by grain boundary and localized shock melting

The fine-grained matrices (less than 150 microns) of 14 gas-rich ordinary chondrile regolith breccias were studied in an attempt to decipher the nature of the lithification process that converted loose regolith material into consolidated breccias. It is found that there is a continuouos gradation in matrix textures from nearly completely clastic (class A) to highly cemented (class C) breccias in which the remining clasts are completely surrounded by interstitial, shock-melted material. It is concluded that this interstitial material is formed by shock melting in the porous regolith. In general, the abundances of solar-wind-implanted He-4 and Ne-20 are inversely correlated with the abundance of intenstitial, shock-melted, feldspathic material. Chondrites with the highest abundance of interstitial, melted material (class C) experienced the highest shock pressures and temperatures and suffered the most extensive degassing. It is this interstitial, feldspathic melt that lithifies and cements the breccias together; those breccias with very little interstitial melt (class A) are the most porous and least consolidated.

Bischoff, A.↗
NASA NTRS2015

Kinetic Monte Carlo Simulation of Oxygen Diffusion in Ytterbium Disilicate

Ytterbium disilicate is of interest as a potential environmental barrier coating for aerospace applications, notably for use in next generation jet turbine engines. In such applications, the transport of oxygen and water vapor through these coatings to the ceramic substrate is undesirable if high temperature oxidation is to be avoided. In an effort to understand the diffusion process in these materials, we have performed kinetic Monte Carlo simulations of vacancy-mediated and interstitial oxygen diffusion in Ytterbium disilicate. Oxygen vacancy and interstitial site energies, vacancy and interstitial formation energies, and migration barrier energies were computed using Density Functional Theory. We have found that, in the case of vacancy-mediated diffusion, many potential diffusion paths involve large barrier energies, but some paths have barrier energies smaller than one electron volt. However, computed vacancy formation energies suggest that the intrinsic vacancy concentration is small. In the case of interstitial diffusion, migration barrier energies are typically around one electron volt, but the interstitial defect formation energies are positive, with the result that the disilicate is unlikely to exhibit experience significant oxygen permeability except at very high temperature.

Coatings↗
DOE OSTI2025

Impact of amorphous pockets on displacement damage evolution in silicon

Silicon has long been known to exhibit amorphization in response to heavy particle bombardment. For doses below the total amorphization threshold, partial amorphization is observed in the form of scattered amorphous pockets. While extensive research has gone into modeling the formation and evolution of amorphous pockets in response to irradiation, no studies yet investigate their impact on the evolution of other damage such as interstitial supersaturation and clustering. In this study, we survey the impact of amorphous pockets on defect evolution in silicon when treated as static sinks. MD is first used to show that amorphous pockets provide energetically favorable sites for point defects relative to the crystalline bulk, supporting the hypothesis that they act as sinks. A 0-D cluster dynamics model is then constructed, taking an interstitial clustering model from the literature and including amorphous pockets as a sink species. We conduct our survey for temperatures between 30 and 400 °C and sink strengths between 1 to 6 x 10 10 cm −2 . Both implantation- and radiation-induced damage states are investigated using interstitial and vacancy concentrations as initial condition variables. We find that, due to the differing migration rates of the interstitial and the vacancy, amorphous pockets have a non-monotonic impact on the final damage state depending on the effective sink strength of the amorphous pockets, resulting in increased damage formation in regimes of intermediate amorphization. In conclusion, this result emphasizes the important role of amorphous pockets in governing the evolution of damage in partially amorphized crystalline materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Coupled cluster and dislocation dynamics modeling of microstructure evolution in irradiated materials

We develop here a coupled cluster and dislocation dynamics framework to study the microstructure evolution of irradiated materials. The framework not only accounts for the three dimensional diffusion of radiation-generated clusters, but also their interaction with dislocation networks and the resultant climb motion of discrete dislocations within finite crystals. The framework is solved with a superposition solution scheme, and is applied to investigate the evolution of the irradiation-induced dislocation loops in zirconium (Zr), considering the effects of various bias factors including the diffusion anisotropy difference (DAD) of interstitials and interstitial clusters, the dislocation bias of defects to discrete dislocation segments, and the production bias of defects from the radiation cascade. We find that the DAD is the most critical factor influencing the kinetics of the loop evolution in Zr, while the recombination/interaction of mobile defects can induce a strong spatial dependence of the loop evolution together with the DAD. Here, the method is also adopted to study the evolution of interstitial $\langle$a$\rangle$ and vacancy $\langle$c$\rangle$ dislocation loop ensembles consistent with the microstructure observed during irradiation-induced growth of Zr. Our findings not only reveal the spatial dependence of the size and ellipticity of the dislocation loops, but also suggest a limit on the anisotropy factor of interstitials to reproduce the co-growth of $\langle$a$\rangle$ and $\langle$c$\rangle$ loops in zirconium, in good agreement with experimental observations and other simulation results.

Bias factors↗
DOE OSTI2024

Parameterization of vacancy production rate in phase-field models of fission gas bubble evolution in nuclear fuel

Phase-field modeling has increasingly been used to study microstructural evolution in fission gas bubbles in nuclear fuel to improve understanding of fission gas release. To improve computational efficiency, often only vacancies and gas atoms are included as defect species. In this case, the net effects of vacancy and interstitial production, recombination, and biased sink absorption are included as a net vacancy source, or net vacancy source combined with an effective sink. However, there has been a lack of clarity on what parameter values should be used for these approaches to best match the more complete physical picture that includes interstitials and vacancies. Here, we compare a phase-field model of void growth to analytical models for the source-only and source plus sink approach to gain insight into how the phase-field models can be parameterized effectively. The source-only approach provides greater flexibility to match growth rates determined from the full vacancy-interstitial picture. A strategy was developed for determining the value of the net vacancy source term by comparing to an analytical model that includes vacancy and interstitial production, recombination, and biased sink absorption. Finally, this strategy can be used to parameterize phase-field models of fission gas bubble growth.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Directing Polymorphism of Colloid Crystals Using Conformational Entropy of Polymer Chains

Controlling the polytypes of close-packed structures of spherical colloids is still a challenging problem despite their wide occurrence. Here, in this work, we show that systematic engineering of the polytype structures of close-packed colloids is possible by using the conformational entropy of polymer chains confined in the interstitial space of colloid crystals. Our interstitial space analysis shows that the hexagonal close-packed (HCP) structures offer larger local interstitial space domains, and the structure director chains in favor of HCP counteract the entropic advantages of the face-centered cubic (FCC) lattices. Using model block copolymer colloids and the known lattice entropy of FCC, a proportionality parameter, β CP = 1.91 × 10 –3 ± 3.67 × 10 –4 , for quantifying the conformational entropy contribution toward HCP structures is extracted. This work demonstrates that the interstitial space of colloid crystals serves as a new structure engineering tool for the self-assembly of colloids.

36 MATERIALS SCIENCE↗
DOE OSTI2025

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↗
DOE OSTI2025

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↗
NASA NTRS1980

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.↗
NASA NTRS1987

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↗
NASA NTRS1993

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↗
NASA NTRS1996

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.↗
NASA NTRS1992

Ventilation-perfusion matching during exercise

In normal subjects, exercise widens the alveolar-arterial PO2 difference (P[A-a]O2) despite a more uniform topographic distribution of ventilation-perfusion (VA/Q) ratios. While part of the increase in P(A-a)O2 (especially during heavy exercise) is due to diffusion limitation, a considerable amount is caused by an increase in VA/Q mismatch as detected by the multiple inert gas elimination technique. Why this occurs is unknown, but circumstantial evidence suggests it may be related to interstitial pulmonary edema rather than to factors dependent on ventilation, airway gas mixing, airway muscle tone, or pulmonary vascular tone. In patients with lung disease, the gas exchange consequences of exercise are variable. Thus, arterial PO2 may increase, remain the same, or fall. In general, patients with advanced chronic obstructive pulmonary disease (COPD) or interstitial fibrosis who exercise show a fall in PO2. This is usually not due to worsening VA/Q relationships but mostly to the well-known fall in mixed venous PO2, which itself results from a relatively smaller increase in cardiac output than VO2. However, in interstitial fibrosis (but not COPD), there is good evidence that a part of the fall in PO2 on exercise is caused by alveolar-capillary diffusion limitation of O2 transport; in COPD (but not interstitial fibrosis), a frequent additional contributing factor to the hypoxemia of exercise is an inadequate ventilatory response, such that minute ventilation does not rise as much as does CO2 production or O2 uptake, causing arterial PCO2 to increase and PO2 to fall.

Non-NASA Center↗
NASA NTRSDate not supplied

Tin-Essako 001: A Metal-Rich Ureilite?

Introduction: Metal-rich achondrites include a variety of types, and likely have a variety of origins. Models range from gravitational mixing at the core-mantle boundaries of differentatiated asteroids, to complex impact mixing scenarios. We describe a new type of metal-rich achondrite that might be the first metal-rich ureilite. Sample: Tin-Essako (TE) 001 (~4.3 g) was found in Mali in 2020 and purchased by Jay Piatek in 2021. It was classified as a metal-rich ungrouped achondrite, with olivine and oxygen isotope (d18O=8.2810‰, d17O=3.718‰, avg. 3) compositions suggesting affinity to ureilites [1]. We studied one polished section (~187.7 mm2) of TE 001. Petrography: TE 001 consists of ~60% metal and 40% silicates, heterogeneously distributed. The metal is largely fresh, but iron oxides (presumably terrestrial) occur along one edge and in some patches and veins in the interior. The silicates are dominantly olivine (≥90%), with melt-textured areas of plagioclase + Si-rich glass. Minor phases include chromite and carbon. Olivine occurs as rounded grains (up to ~2.5 mm) in metal, commonly with rims of melt-textured plagioclase + glass. Olivine also occurs in more massive areas having a “honeycomb” texture, with rounded olivine “cells” surrounded by an interstitial network of reduced olivine riddled with tiny metal grains, plus melt-textured plagioclase + Si-rich glass. Chromite occurs as subhedral to rounded grains; smaller grains (30-250 mm) are included in olivine and a few larger grains (400-500 mm) are isolated within metal. A carbon phase occurs as lacy-textured rims around olivine grains in metal, or small patches within metal. Mineral Compositions: The olivine (excluding interstitial areas) is Fo 73.9±0.5, with 0.25±0.02 wt.% CaO, 0.31±0.01 wt% Cr2O3, 0.01 wt.% NiO, and molar Fe/Mn=49.8±2.2 (53 analyses). Olivine in interstitial areas has Fo up to at least 91. Smaller chromite grains have Fe# (molar Fe/[Fe+Mg]) = 0.54±0.01, Cr# (molar Cr/[Cr+Al]) = 0.52±0.01, 0.52±0.02 wt.% V2O3 and 0.21±0.04 wt% ZnO (28 analyses). One larger grain is zoned from Fe# = 0.45, Cr# =0.52 to Fe# = 0.40, Cr# =0.57, and contains thin Al-rich lamellae not resolved by EMPA. One irregularly shaped patch of chromite included in olivine has Fe# =0.26 and contains no ZnO. Plagioclase laths are An ~53-60, with ≤0.01 wt.% K2O. Glass contains 75-76 wt% SiO2 and ~16 wt% Al2O3. The metal contains 5.2±0.2 wt% Ni, 0.46±0.02 wt.% Co, and 0.01±0.01 wt.% Cr, with Si and P below detection (168 analyses). Discussion: The olivine + chromite assemblage in TE 001 is similar to the most ferroan ureilites (Fo ~75-79 [2]), as are the oxygen isotopes [1]. The presence of a carbon phase supports this, although the identity of this phase (graphite as in ureilites?) remains to be determined. The “honeycomb” textured areas, in particular, the presence of olivine “reduction rims,” resemble shock-smelted olivine areas in ureilites [3], but interstitial melt-textured plagioclase laths + glass like those in TE 001 have not been reported in such (or any) ureilites. Olivine in TE 001 is marginally more ferroan than in the most FeO-rich ureilite, with Fe-Mg-Mn composition offset from the trend of olivine + low-Ca pyroxene ureilites similar to augite-bearing ureilites [4]. CaO and Cr2O3 contents are in the range of those in ureilite olivine [ ] , though Cr2O3 is at the extreme low end of the range [5]. Chromites (except the unusual one) have similar Fe# to the most ferroan primary chromites in ureilites [2], but distinctly lower Cr# (0.52 vs. 0.71). Metal compositions are with the range for metal in ureilites [7]. The absence of pyroxene and sulfide, and the high abundance of metal in TE 001, are unlike ureilites. One possibility is that a ferroan, chromite-bearing, pyroxene-poor ureilite was invaded (possibly due to impact) by a metallic liquid (low S content suggests very high temperature), resulting in complete melting of pyroxene and smelting of olivine, with rapid recrystallization of melted silicate as plagioclase + glass. Alternatively, a pre-existing metal-rich ureilite assemblage may have been impact melted. Additional types of data will be obtained to evaluate these hypotheses (i.e., is the metal indigenous?) and assess affinity to ureilites.

olivine↗
DOE OSTI2024

Negative thermal expansion coefficient and amorphization in defective 4H-SiC

This paper presents thermal expansion coefficient (TEC) and amorphization in 4H-SiC containing point defects at different concentrations. We considered vacancy defects, interstitial defects, and Frenkel pair defects and investigated the thermomechanical response of the lattice over a wide range of temperatures using classical molecular dynamics simulations. The results show that 4H-SiC with vacancy defects exhibits a negative TEC above a critical defect density of around 9% (irrespective of the temperature). With interstitial defects, it exhibits a positive TEC (regardless of the defect density), and with Frenkel pair defects it shows a transition from positive TEC to negative TEC for a defect density greater than 8%. The coupling between temperature-induced expansion and defect-introduced stress in the lattice forms the mechanistic basis for the observed variation in TEC. Furthermore, the specific heat decreases rapidly with an increase in defect density at room temperature, with the highest sensitivity of the lattice observed for the Frenkel pair defects followed by interstitial defects and then by vacancy defects. Finally, these findings highlight the critical implications of defects on thermal expansion behavior of 4H-SiC with applications in radiation environments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Be₂C as a neutron moderator for molten salt reactors: A first-principles study of structural, electronic, and defect properties

Beryllium carbide (Be 2 C), valued for its high neutron moderation efficiency and low absorption cross section, is a promising high-temperature neutron moderator for molten salt reactors. Its practical adoption, however, demands significant technological maturation, requiring comprehensive theoretical and experimental studies of its response to different conditions, including high temperature and irradiation. Here, we report initial results on the fundamental properties and radiation-induced defects of Be 2 C, focusing on antisites, vacancies, interstitial atoms, and Frenkel pairs in the Be and C sublattices. Using density functional theory (DFT) and ab initio molecular dynamics (AIMD), we calculate the defects formation and binding energies, evaluating their dependence on the supercell size, charge states, and chemical environment. In general, carbon defects exhibit higher formation energies, greater sensitivity to cell size, and stronger impacts on the density of states compared to beryllium defects, with charged state the effects being more pronounced. Static DFT reveals multiple metastable interstitial configurations, while AIMD identifies ground states as C-C <100> dumbbells and octahedral Be interstitials. In conclusion, the diversity of metastable configurations and defect states complicates the diffusion mechanisms, requiring further molecular dynamics analysis to elucidate the mechanisms and rates of radiation-induced atomic transport, as well as the structural stability of Be2C.

Ab initio modeling↗