SEARCH · Search NASA

Results for “metastable defect configurations”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

11 records

All sources NASA NTRS DOE OSTI

DOE OSTI2023

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↗

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↗

DOE OSTI2024

Chiral, Topological, and Knotted Colloids in Liquid Crystals

The geometric shape, symmetry, and topology of colloidal particles often allow for controlling colloidal phase behavior and physical properties of these soft matter systems. In liquid crystalline dispersions, colloidal particles with low symmetry and nontrivial topology of surface confinement are of particular interest, including surfaces shaped as handlebodies, spirals, knots, multi-component links, and so on. These types of colloidal surfaces induce topologically nontrivial three-dimensional director field configurations and topological defects. Director switching by electric fields, laser tweezing of defects, and local photo-thermal melting of the liquid crystal host medium promote transformations among many stable and metastable particle-induced director configurations that can be revealed by means of direct label-free three-dimensional nonlinear optical imaging. The interplay between topologies of colloidal surfaces, director fields, and defects is found to show a number of unexpected features, such as knotting and linking of line defects, often uniquely arising from the nonpolar nature of the nematic director field. This review article highlights fascinating examples of new physical behavior arising from the interplay of nematic molecular order and both chiral symmetry and topology of colloidal inclusions within the nematic host. Furthermore, the article concludes with a brief discussion of how these findings may lay the groundwork for new types of topology-dictated self-assembly in soft condensed matter leading to novel mesostructured composite materials, as well as for experimental insights into the pure-math aspects of low-dimensional topology.

chirality↗

NASA NTRS1995

Analysis of the Defect Structure of B2 Feal Alloys

The Bozzolo, Ferrante and Smith (BFS) method for alloys is applied to the study of the defect structure of B2 FeAI alloys. First-principles Linear Muffin Tin Orbital calculations are used to determine the input parameters to the BFS method used in this work. The calculations successfully determine the phase field of the B2 structure, as well as the dependence with composition of the lattice parameter. Finally, the method is used to perform 'static' simulations where instead of determining the ground state configuration of the alloy with a certain concentration of vacancies, a large number of candidate ordered structures are studied and compared, in order to determine not only the lowest energy configurations but other possible metastable states as well. The results provide a description of the defect structure consistent with available experimental data. The simplicity of the BFS method also allows for a simple explanation of some of the essential features found in the concentration dependence of the heat of formation, lattice parameter and the defect structure.

Bozzolo, Guillermo↗

DOE OSTI2024

Structure and migration of heavily irradiated grain boundaries and dislocations in Ni in the athermal limit

The microstructural evolution at and near preexisting grain boundaries (GBs) and dislocations in materials under high radiation doses is still poorly understood. In this work, we use the creation relaxation algorithm (CRA) developed for atomistic modeling of high-dose irradiation in bulk materials to probe the athermal limit of saturation of GB and dislocation core regions under irradiation in fcc Ni. We find that, upon continuously subjecting a single dislocation or GB to Frenkel pair creation in the athermal limit, a local steady-state disordered defect structure is reached with excess properties that fluctuate around constant values. Case studies are given for a straight screw dislocation which elongates into a helix under irradiation and several types of low- and high-angle GBs, which exhibit coupled responses such as absorption of extrinsic dislocations, roughening, and migration. A positive correlation is found between the initial GB energy and the local steady-state GB energy under irradiation across a wide variety of GB types. Metastable GB structures with similar density in the defect core region but different initial configurations are found to converge to the same limiting structure under CRA. The mechanical responses of pristine and irradiated dislocations and GB structures are compared under an applied shear stress. Irradiated screw and edge dislocations are found to exhibit a hardening response, migrating at larger flow stresses than their pristine counterparts. Mobile GBs are found to exhibit softening or hardening responses depending on GB character. Although some GBs recover their initial pristine structures upon migration outside of the radiation zone, many GBs sustain different flow stresses corresponding to altered mobile core structures.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Native point defects in HgCdTe infrared detector material: Identifying deep centers from first principles

We investigate the native point defects in the long-wavelength infrared (LWIR) detector material Hg0.75Cd0.25Te using a dielectric-dependent hybrid density functional combined with spin–orbit coupling. Characterizing these point defects is essential as they are responsible for intrinsic doping and nonradiative recombination centers in the detector material. The dielectric-dependent hybrid functional allows for an accurate description of the bandgap (Eg) for Hg1−xCdxTe (MCT) over the entire compositional range, a level of accuracy challenging with standard hybrid functionals. Our comprehensive examination of the native point defects confirms that cation vacancies VHg(Cd) are the primary sources of p-type conductivity in the LWIR material given their low defect formation energies and the presence of a shallow acceptor level (−/0) near the valence-band maximum. In addition to the shallow acceptor level, the cation vacancies exhibit a deep charge transition level (2−/−) situated near the midgap, characteristic of nonradiative recombination centers. Our results indicate that Hg interstitial could also be a deep center in the LWIR MCT through a metastable configuration under the Hg-rich growth conditions. While an isolated Te antisite does not show deep levels, the formation of VHg–TeHg defect complex introduces a deep acceptor level within the bandgap.

Physics↗

NASA NTRS2003

A Dynamic Pathway for Stone-Wales Bond Rotation on Carbon Nanotubes through Diamond-Like Bonds

A new lower energy barrier with a two-step pathway of Stone-Wales (SW) ,ond rotation on carbon nanotubes (CNTs) is found through molecular dynamics (MD) simulations of CNTs under tension. The first step involves going over to a stable sp3-like metastable configuration with half rotated and partially tilted C-C bond. The second step involves going over to the fully rotated C-C bond with the formation of a SW defect in the nanotube. The energy barrier for this two-step dynamic pathway is significantly lower than the previously known static barrier for in-plane rotation of the C-C bond on a tensile strained (> 4%) CNT.

Wei, Chen-Yu↗

DOE OSTI2025

Atomic-scale structure of ZrO 2 : Formation of metastable polymorphs

Metastable phases can exist within local minima in the potential energy landscape when they are kinetically “trapped” by various processing routes, such as thermal treatment, grain size reduction, chemical doping, interfacial stress, or irradiation. Despite the importance of metastable materials for many technological applications, little is known about the underlying structural mechanisms of the stabilization process and atomic-scale nature of the resulting defective metastable phase. Investigating ion-irradiated and nanocrystalline zirconia with neutron total scattering experiments, we show that metastable tetragonal ZrO 2 consists of an underlying structure of ferroelastic, orthorhombic nanoscale domains stabilized by a network of domain walls. The apparent long-range tetragonal structure that can be recovered to ambient conditions is only the configurational ensemble average of the underlying orthorhombic domains. This structural heterogeneity with a distinct short-range order is more broadly applicable to other nonequilibrium materials and provides insight into the synthesis and recovery of functional metastable phases with unique physical and chemical properties.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2022

Stabilized Synthesis of 2D Verbeekite: Monoclinic PdSe 2 Crystals with High Mobility and In-Plane Optical and Electrical Anisotropy

PdSe 2 has a layered structure with an unusual, puckered Cairo pentagonal tiling. Its atomic bond configuration features planar 4-fold-coordinated Pd atoms and intralayer Se–Se bonds that enable polymorphic phases with distinct electronic and quantum properties, especially when atomically thin. PdSe 2 is conventionally orthorhombic, and direct synthesis of its metastable polymorphic phases is still a challenge. Here, we report an ambient-pressure chemical vapor deposition approach to synthesize metastable monoclinic PdSe 2 . Monoclinic PdSe 2 is shown to be synthesized selectively under Se-deficient conditions that induce Se vacancies. These defects are shown by first-principles density functional theory calculations to reduce the free energy of the metastable monoclinic phase, thereby stabilizing it during synthesis. Further, the structure and composition of the monoclinic PdSe 2 crystals are identified and characterized by scanning transmission electron microscopy imaging, convergent beam electron diffraction, and electron energy loss spectroscopy. Polarized Raman spectroscopy of the monoclinic PdSe 2 flakes reveals their strong in-plane optical anisotropy. Electrical transport measurements show that the monoclinic PdSe 2 exhibits n-type charge carrier conduction with electron mobilities up to ~298 cm 2 V –1 s –1 and a strong in-plane electron mobility anisotropy of ~1.9. The defect-mediated growth pathway identified in this work is promising for phase-selective direct synthesis of other 2D transition metal dichalcogenides.

2D material↗

DOE OSTI2025

Unraveling the Atomic Mechanism of the Crystalline Phase‐Dependent Structural Features and Special Spectral Design of α‐, β‐, and Ɛ‐Ga₂O₃

Atomic‐scale phase transformations profoundly influence the functional properties of Ga₂O₃ polymorphs. By combining irradiation experiments with microstructure characterization and theoretical approaches, phase‐specific energy‐dissipation pathways in α‐, β‐, and ε‐Ga₂O₃ are uncovered and strategies for targeted property design are outlined. Competing antiphase boundaries (APBs) and twin domain boundaries (TDBs) promote irreversible α→ε interconversion through domain fragmentation. In β‐Ga₂O₃, defect‐induced stress gradients drive two distinct local transformations: surface Ga‐aggregated β→δ that stabilizes transient states, and latent‐track‐confined β→κ phase transition with recoverable distortions via cation reordering. Under electronic excitation, β‐Ga₂O₃ forms nanohillocks via robust GaO₆ octahedra (high density/strong Ga─O bonds), while α/ε‐Ga₂O₃ generates nanopores from tetrahedral Ga looseness (low bonding energy), highlighting phase‐dependent surface dynamics shaped by atomic packing and bonding anisotropy. Defect‐regulated recombination suppresses visible photoluminescence in α/β‐Ga₂O₃, whereas in ε‐Ga₂O₃ bandgap narrowing of ΔE: 0.30 eV is observed, enhancing emission. Linking phase‐dependent defect‐carrier interactions and metastable‐phase engineering in Ga₂O₃ enables property optimization for power‐electronics and optoelectronics devices.

electronic state configuration↗

DOE OSTI2024

Accelerating defect predictions in semiconductors using graph neural networks

First-principles computations reliably predict the energetics of point defects in semiconductors but are constrained by the expense of using large supercells and advanced levels of theory. Machine learning models trained on computational data, especially ones that sufficiently encode defect coordination environments, can be used to accelerate defect predictions. Here, we develop a framework for the prediction and screening of native defects and functional impurities in a chemical space of group IV, III–V, and II–VI zinc blende semiconductors, powered by crystal Graph-based Neural Networks (GNNs) trained on high-throughput density functional theory (DFT) data. Using an innovative approach of sampling partially optimized defect configurations from DFT calculations, we generate one of the largest computational defect datasets to date, containing many types of vacancies, self-interstitials, anti-site substitutions, impurity interstitials and substitutions, as well as some defect complexes. We applied three types of established GNN techniques, namely crystal graph convolutional neural network, materials graph network, and Atomistic Line Graph Neural Network (ALIGNN), to rigorously train models for predicting defect formation energy (DFE) in multiple charge states and chemical potential conditions. We find that ALIGNN yields the best DFE predictions with root mean square errors around 0.3 eV, which represents a prediction accuracy of 98% given the range of values within the dataset, improving significantly on the state-of-the-art. We further show that GNN-based defective structure optimization can take us close to DFT-optimized geometries at a fraction of the cost of full DFT. The current models are based on the semi-local generalized gradient approximation-Perdew–Burke–Ernzerhof (PBE) functional but are highly promising because of the correlation of computed energetics and defect levels with higher levels of theory and experimental data, the accuracy and necessity of discovering novel metastable and low energy defect structures at the PBE level of theory before advanced methods could be applied, and the ability to train multi-fidelity models in the future with new data from non-local functionals. The DFT-GNN models enable prediction and screening across thousands of hypothetical defects based on both unoptimized and partially optimized defective structures, helping identify electronically active defects in technologically important semiconductors.

Rahman, Md Habibur (ORCID:000000027705984X)↗