SEARCH · Engineering Papers

Results for “crystal structure prediction”

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.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2026

AI-assisted rapid crystal structure generation towards a target local environment

In material design, traditional crystal structure prediction approaches are expensive as they require extensive structural sampling through expensive energy minimization methods. Emerging artificial intelligence (AI) generative models have shown great promise in rapidly generating realistic crystals, but they typically handle only a few tens of atoms per unit cell. To overcome this limitation, we introduce a symmetry-informed approach, the Local Environment Geometry-Oriented Crystal Generator (LEGO-xtal). Our method generates initial structures using AI models trained on an augmented dataset, and then optimizes them using structure descriptors rather than energy-based optimization. We demonstrate its effectiveness by expanding from 25 known low-energy sp2 carbon allotropes to over 1700, all within 0.5 eV/atom of the ground-state energy of graphite. This framework offers a generalizable strategy for the targeted design of materials with modular building blocks, such as metal-organic frameworks and battery materials.

Ridwan, Osman Goni [University of North Carolina a↗

DOE OSTI2024

Suppression of Disorder in Benzamide and Thiobenzamide Crystals by Fluorine Substitution

Disorder is a common feature of molecular crystals that complicates determination of structures and can potentially affect electric and mechanical properties. Suppression of disorder is observed in otherwise severely disordered benzamide and thiobenzamide crystals by substituting hydrogen with fluorine in the ortho-position of the phenyl ring. Fluorine occupancies of 20–30% are sufficient to suppress disorder without changing the packing motif. Crystal structure prediction calculations reveal a much denser lattice energy landscape for benzamide compared to 2-fluorobenzamide, suggesting that fluorine substitution makes disorder less likely.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Structurally Constrained Evolutionary Algorithm for the Discovery and Design of Metastable Phases

Metastable materials are abundant in nature and technology, showcasing remarkable properties that inspire innovative materials design. However, traditional crystal structure prediction methods, which rely solely on energetic factors to determine a structure’s fitness, are not suitable for predicting the vast number of potentially synthesizable phases that represent a local minimum corresponding to a state in thermodynamic equilibrium. Here, we present a new approach for the prediction of metastable phases with specific structural features, and interface this method with the XTALOPT evolutionary algorithm. Our method relies on structural features that include the local crystalline order (e.g., the coordination number or chemical environment), and symmetry (e.g., Bravais lattice and space group) to filter the breeding pool of an evolutionary crystal structure search. The effectiveness of this approach is benchmarked on three known metastable systems: XeN 8 , with a two-dimensional polymeric nitrogen sublattice, brookite TiO 2 , and a high pressure BaH 4 phase that was recently characterized. Additionally, a newly predicted metastable melaminate salt, P1¯WC 3 N 6 , was found to possess an energy that is lower than two phases proposed in a recent computational study. Here, the method presented here could help in identifying the structures of compounds that have already been synthesized, and developing new synthesis targets with desired properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Force Field X: A computational microscope to study genetic variation and organic crystals using theory and experiment

Force Field X (FFX) is an open-source software package for atomic resolution modeling of genetic variants and organic crystals that leverages advanced potential energy functions and experimental data. FFX currently consists of nine modular packages with novel algorithms that include global optimization via a many-body expansion, acid–base chemistry using polarizable constant-pH molecular dynamics, estimation of free energy differences, generalized Kirkwood implicit solvent models, and many more. Applications of FFX focus on the use and development of a crystal structure prediction pipeline, biomolecular structure refinement against experimental datasets, and estimation of the thermodynamic effects of genetic variants on both proteins and nucleic acids. The use of Parallel Java and OpenMM combines to offer shared memory, message passing, and graphics processing unit parallelization for high performance simulations. Overall, the FFX platform serves as a computational microscope to study systems ranging from organic crystals to solvated biomolecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Combined First-Principles and Experimental Investigation into the Reactivity of Codeposited Chromium–Carbon under Pressure

High-pressure synthesis in the diamond anvil cell suffers from the lack of a general approach for the control of precursor stoichiometry and homogeneity. Here, we present results from a new method we have developed that uses magnetron cosputtering to prepare stoichiometrically precise and atomically mixed amorphous films of Cr:C. Laser-heated diamond anvil cell experiments carried out on a flake of this sample at pressures between 13.5 and 24.3 GPa lead to the observation of Cr 3 C (Pnma) over the entire pressure range–in good agreement with our in-house theoretical predictions–but also reveal two other metastable phases that were not expected: a novel monoclinic chromium carbide phase and the NaCl-type CrC (Fm3̅m) phase. The unexpected stability of CrC is investigated by using first-principles methods, revealing a large stabilizing effect tied to substoichiometry at the carbon site. These results offer an important case study into the current limitations of crystal structure prediction methods with regard to phase complexity and bolster the growing need for advanced theoretical approaches that can more completely survey experimentally unexplored phase space.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Superconductivity in the Li-B-C system at 100 GPa

Layer Li-B-C compounds have been shown to have feasible superconductivity. Using the adaptive genetic algorithm, we predict the structures of the Li-B-C system at 100 GPa. Here we identify several low-enthalpy metallic phases with stoichiometries of Li B 2 C, Li B 3 C, Li 2 BC 2 , Li 3 B 2 C 3 , Li 3 BC, and Li 5 BC. Using a fast screening method of electron-phonon interaction, we find that Li B 3 C is a promising candidate for superconductivity. The consecutive calculations using the full Brillouin zone confirm the existence of the strong electron-phonon coupling (EPC) in this system. The anharmonic B-C phonon modes near the zone center provide the major contribution to the EPC. The EPC constant is 1.40, and the estimated critical temperature is 22 K. In this paper, we indicate that superconductivity can also happen without a layered structural motif in the Li-B-C system. We also demonstrate an effective strategy for crystal structure prediction of superconducting materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Solid State Reduction Driven Synthesis of Mn Containing Multi-principal Component Alloys

In recent years, high entropy alloys (HEAs), also known as multi-principal component alloys (MPCAs) have emerged as a new and exciting class of materials. This paper reports on the solid state reduction synthesis of a series of CoFeNiMn-based MPCA compositions, starting from a mixture of the corresponding oxides. One of the aims of the study was to test whether the degree of reduction of MnO, a highly stable oxide, could be enhanced by tailoring the alloy composition. Specifically, the influence of Ni content was studied because Ni exhibits a significant negative enthalpy of mixing with Mn. High purity precursor powders of Co(OH) 2 , Fe 2 O 3 , MnO 2 , and NiO were milled and mixed using standard ceramic processing methods. Here, the nominal sample compositions (assuming complete oxide reduction) were (CoFeMn) x Ni (1-x) , for x = 0, 0.083, 0.166, and 0.25. The oxide samples were subjected to a series of isothermal reduction anneals in flowing 3 pct H 2 –Ar at 1100 °C. The resulting microstructures were characterized using scanning electron microscopy (SEM), X-ray energy dispersive spectroscopy (EDS) and X-ray diffraction (XRD). The composition of the resulting MPCAs was determined quantitatively using wavelength dispersive spectroscopy (WDS) in the electron microprobe. The study revealed that for each of the initial oxide compositions studied, it was possible to achieve an MPCA with ~ 25 at. pct Mn. These results were found to be consistent with the predictions of a thermodynamic model whereby a negative enthalpy of mixing (ΔH mix ), combined with a contribution from configurational entropy, can offset a positive free energy of reduction (ΔG red ). The incorporation of vibrational entropy into first principles calculations was found to have a significant effect on the predicted crystal structure of the MPCAs.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Electronic structure and energetics of a heterodimeric BChl g ′/Chl a ′ special pair generated by exposure of Heliomicrobium modesticaldum to dioxygen

Heliobacteria are anoxygenic phototrophs that have a Type I homodimeric reaction center containing bacteriochlorophyll g (BChl g). Previous experimental studies have shown that in the presence of light and dioxygen, BChl g is converted into 8 1 -OH-chlorophyll a F (hereafter Chl a F ), with an accompanying loss of light-driven charge separation. These studies suggest that the reaction center only loses the ability to transfer electrons once both BChl g′ molecules of the P 800 special pair have been converted to Chl a F ′. Here, the present work confirms that the partially converted BChl g′/Chl a F ′ special pair remains functional in samples exposed to dioxygen by demonstrating its presence using hyperfine couplings obtained from Q-band 1 H ENDOR, 2D 14 N HYSCORE and DFT methods. The DFT calculations of the BChl g′/BChl g′ homodimeric primary donor, which are based on the recently published X-ray crystal structure, predict that the unpaired electron spin is equally delocalized over both BChl g′ molecules and provide an excellent match to the experimental hyperfine couplings of the anaerobic samples. Exposure to dioxygen leads to substantial changes in the hyperfine interactions, indicative of greater localization of the unpaired electron spin. The measured hyperfine couplings are reproduced in the DFT calculations by replacing one of the BChl g′ molecules of the primary donor with a Chl a F ′ molecule. The calculations reveal that the spin density becomes localized on BChl g′ in the heterodimeric primary donor. Time-dependent DFT calculations demonstrate that conversion of either or both of the accessory BChl g molecules and/or one of the BChl g′ molecules of P800 to Chl a F ′ results in minor effects on the energy of the charge-separated states. In contrast, if both of the BChl g′ molecules of P 800 are converted a large increase in the energy of the charge-separated state occurs. This suggests that the reaction center remains functional when only one half of the dimer is converted, however, conversion of both halves of the P 800 dimer leads to loss of function

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Computational discovery of two-dimensional rare-earth iodides: promising ferrovalley materials for valleytronics

Two-dimensional Ferrovalley materials with intrinsic valley polarization are rare but highly promising for valley-based nonvolatile random access memory and valley filter devices. These ferromagnetic materials exhibit valleys at or near the Fermi level with intrinsic magnetism. The strong coupling between magnetism and spin–orbit coupling induces intrinsic valley polarization. Using Kinetically Limited Minimization, an unconstrained crystal structure prediction algorithm, and prototype sampling based on first-principles calculations, we have discovered new Ferrovalley materials, rare-earth iodides RI 2 , where R is a rare-earth element belonging to Sc, Y, or La-Lu, and I is Iodine. The rare-earth iodides are layered and demonstrate either 2H, 1T, or 1T d phase as the ground state in bulk, analogous to transition metal dichalcogenides (TMDCs). The calculated exfoliation energy of monolayers (MLs) is comparable to that of graphene and TMDCs, suggesting possible experimental synthesis. The MLs in the 2H phase exhibit ferromagnetism due to unpaired electrons in d and f orbitals. Throughout the rare-earth series, d bands have valley polarization at K and $\overline{K}$ points in the Brillouin zone in the vicinity of the Fermi level. Large intrinsic valley polarization in the range of 15–143 meV without external stimuli is observed in these Ferrovalley materials, which can be enhanced further by applying an in-plane bi-axial strain. These valleys can selectively be probed and manipulated for information storage and processing, potentially offering superior performance beyond conventional electronics and spintronics. Here we further show that the 2H ferromagnetic phase of RI 2 MLs possesses non-zero Berry curvature and exhibits anomalous valley Hall effect with considerable anomalous Hall conductivity. Our work will incite exploratory synthesis of the predicted Ferrovalley materials and their application in valleytronics and beyond.

2D materials↗

DOE OSTI2023

Structure, stability, and superconductivity of N-doped lutetium hydrides at kbar pressures

Here, the structure of the material responsible for the room temperature and near ambient pressure superconductivity reported in an N-doped lutetium hydride has not been conclusively determined. Herein, density functional theory calculations are performed in an attempt to uncover what it might be. Guided by a range of strategies including crystal structure prediction and modifications of existing structure types, we present an array of Lu-N-H phases that are dynamically stable at experimentally relevant pressures. Although none of the structures found are thermodynamically stable, and none are expected to remain superconducting above ~17 K at 10 kbar, a number of metallic compounds with fcc Lu lattices – as suggested by the experimental X-ray diffraction measurements of the majority phase – are identified. The system whose calculated equation of states matches best with that measured for the majority phase is fluorite-type LuH 2 , whose 10 kbar superconducting critical temperature was estimated to be 0.1 K using the Allen-Dynes modified McMillan equation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Constraints on 5 f -electron magnetism in Ga-stabilized $δ$-Pu from x-ray magnetic circular dichroism

Density-functional theory models of δ-Pu accurately predict crystal structure, phonon density of states, and unit cell volumes but require magnetic degrees of freedom which have never been experimentally verified. Some models invoke an on-site cancellation of spin and orbital moments, engendering a near-zero bulk magnetization, undetectable by most probes. Here, we employ x-ray magnetic circular dichroism at the Pu M 4,5 edges to directly probe spin and orbital moments using the magneto-optical sum rules. The data show no dichroism within experimental error, constraining polarized moments at 6 T and 3 K to μ L,S <0.1⁢μ B . Finally, these experiments point to the absence of even unconventional spin-orbit compensated order in δ-Pu.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Structural analysis and functional evaluation of the disordered ß–hexosyltransferase region from Hamamotoa (Sporobolomyces) singularis

Hamamotoa (Sporobolomyces) singularis codes for an industrially important membrane bound ß-hexosyltransferase (BHT), (BglA, UniprotKB: Q564N5) that has applications in the production of natural fibers such as galacto-oligosaccharides (GOS) and natural sugars found in human milk. When heterologously expressed by Komagataella phaffii GS115, BHT is found both membrane bound and soluble secreted into the culture medium. In silico structural predictions and crystal structures support a glycosylated homodimeric enzyme and the presence of an intrinsically disordered region (IDR) with membrane binding potential within its novel N-terminal region (1–110 amino acids). Additional in silico analysis showed that the IDR may not be essential for stable homodimerization. Thus, we performed progressive deletion analyses targeting segments within the suspected disordered region, to determine the N-terminal disorder region’s impact on the ratio of membrane-bound to secreted soluble enzyme and its contribution to enzyme activity. The ratio of the soluble secreted to membrane-bound enzyme shifted from 40% to 53% after the disordered N-terminal region was completely removed, while the specific activity was unaffected. Furthermore, functional analysis of each glycosylation site found within the C-terminal domain revealed reduced total secreted protein activity by 58%–97% in both the presence and absence of the IDR, indicating that glycosylation at all four locations is required by the host for the secretion of active enzyme and independent of the removed disordered N-terminal region. Overall, the data provides evidence that the disordered region only partially influences the secretion and membrane localization of BHT.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Superconductivity in CH 4 and BH – 4 containing compounds derived from the high-pressure superhydrides

Inspired by the synthesis of the high-pressure Fm3m LaH 10 superconducting superhydride, systematic density functional theory (DFT) calculations are performed to study ternaries that could be derived from it by replacing two of the hydrogen atoms with boron or carbon and varying the identity of the electropositive element. Though many of the resulting alkali-metal and alkaline-earth MC 2 H 8 phases are predicted to be dynamically stable at mild pressures, their superconducting critical temperatures (T c s) are low because their metallicity results from the filling of an electride-like band. Substitution with a trivalent element leads to phases with substantial metal d- character at the Fermi level whose T c s are typically above 40 K. Here, among the MB 2 H 8 phases examined, KB 2 H 8 , RbB 2 H 8 and CsB 2 H 8 are predicted to be dynamically stable at very mild pressures, and their stability is rationalized by a DFT-Chemical Pressure analysis that elucidates the role of the M atom size. Quantum anharmonic effects strongly affect the properties of KB 2 H 8 , the highest predicted T c compound, near 10 GPa, but molecular dynamics simulations reveal it would decompose below its T c at this pressure. Nonetheless, at ca. 50 GPa KB 2 H 8 is predicted to be thermally stable with a superconducting figure of merit surpassing that of the recently synthesized LaBeH 8 .

36 MATERIALS SCIENCE↗

DOE OSTI2023

Accelerating crystal structure determination with iterative AlphaFold prediction

Experimental structure determination can be accelerated with artificial intelligence (AI)-based structure-prediction methods such as AlphaFold . Here, an automatic procedure requiring only sequence information and crystallographic data is presented that uses AlphaFold predictions to produce an electron-density map and a structural model. Iterating through cycles of structure prediction is a key element of this procedure: a predicted model rebuilt in one cycle is used as a template for prediction in the next cycle. This procedure was applied to X-ray data for 215 structures released by the Protein Data Bank in a recent six-month period. In 87% of cases our procedure yielded a model with at least 50% of C α atoms matching those in the deposited models within 2 Å. Predictions from the iterative template-guided prediction procedure were more accurate than those obtained without templates. It is concluded that AlphaFold predictions obtained based on sequence information alone are usually accurate enough to solve the crystallographic phase problem with molecular replacement, and a general strategy for macromolecular structure determination that includes AI-based prediction both as a starting point and as a method of model optimization is suggested.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2009

Computational Study on Full-length Human Ku70 with Double Stranded DNA: Dynamics, Interactions and Functional Implications

The Ku70/80 heterodimer is the first repair protein in the initial binding of double-strand break (DSB) ends following DNA damage, and is a component of nonhomologous end joining repair, the primary pathway for DSB repair in mammalian cells. In this study we constructed a full-length human Ku70 structure based on its crystal structure, and performed 20 ns conventional molecular dynamic (CMD) simulations on this protein and several other complexes with short DNA duplexes of different sequences. The trajectories of these simulations indicated that, without the topological support of Ku80, the residues in the bridge and C-terminal arm of Ku70 are more flexible than other experimentally identified domains. We studied the two missing loops in the crystal structure and predicted that they are also very flexible. Simulations revealed that they make an important contribution to the Ku70 interaction with DNA. Dislocation of the previously studied SAP domain was observed in several systems, implying its role in DNA binding. Targeted molecular dynamic (TMD) simulation was also performed for one system with a far-away 14bp DNA duplex. The TMD trajectory and energetic analysis disclosed detailed interactions of the DNA-binding residues during the DNA dislocation, and revealed a possible conformational transition for a DSB end when encountering Ku70 in solution. Compared to experimentally based analysis, this study identified more detailed interactions between DNA and Ku70. Free energy analysis indicated Ku70 alone is able to bind DNA with relatively high affinity, with consistent contributions from various domains of Ku70 in different systems. The functional implications of these domains in the processes of Ku heterodimerization and DNA damage recognition and repair can be characterized in detail based upon this analysis.

Hu, Shaowen↗

DOE OSTI2025

Explainable Synthesizability Prediction of Inorganic Crystal Polymorphs Using Large Language Models

Abstract We evaluate the ability of machine learning to predict whether a hypothetical crystal structure can be synthesized and explain those predictions to scientists. Fine‐tuned large language models (LLMs) trained on a human‐readable text description of the target crystal structure perform comparably to previous bespoke convolutional graph neural network methods, but better prediction quality can be achieved by training a positive‐unlabeled learning model on a text‐embedding representation of the structure. An LLM‐based workflow can then be used to generate human‐readable explanations for the types of factors governing synthesizability, extract the underlying physical rules, and assess the veracity of those rules. These explanations can guide chemists in modifying or optimizing non‐synthesizable hypothetical structures to make them more feasible for materials design.

Kim, Seongmin [Department of Chemical and Biologic↗

DOE OSTI2025

Explainable Synthesizability Prediction of Inorganic Crystal Polymorphs Using Large Language Models

Kim, Seongmin [Department of Chemical and Biologic↗

DOE OSTI2024

Ultrafast Yttrium Hydride Chemistry at High Pressures via Non-equilibrium States Induced by an X-ray Free Electron Laser

Controlling the formation and stoichiometric content of the desired phases of materials has become of central interest for a variety of fields. The possibility of accessing metastable states by initiating reactions by X-ray-triggered mechanisms over ultrashort time scales has been enabled by the development of X-ray free electron lasers (XFELs). Utilizing the exceptionally high-brilliance X-ray pulses from the EuXFEL, we report the synthesis of a previously unobserved yttrium hydride under high pressure, along with nonstoichiometric changes in hydrogen content as probed at a repetition rate of 4.5 MHz using time-resolved X-ray diffraction. Furthermore, exploiting non-equilibrium pathways, we synthesize and characterize a hydride in a Weaire–Phelan structure type at pressures as low as 125 GPa, predicted using a crystal structure search, with a hydrogen content of 4.0–5.75 hydrogens per cation, that is enthalpically metastable on the convex hull.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗