Engineering Papers⌕ Search

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.

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

At least 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Synthesizing intermetallic thermoelectric materials using metal flux growth method

With zero energy poverty by 2050 goals, exploring innovative energy technologies is essential. Thermoelectric materials (TEMs) have the potential to harness waste heat and convert it to electricity (e.g., for spacecraft power, sensors, and cryogenic cooling). Two limitations of thermoelectric compounds are their low efficiency and limited pool of materials. Optimizing the efficiency of TEMs is challenging due to the interdependent relationship between electrical and thermal conductivity. A deeper understanding of electronic properties and lattice thermal conductivity is necessary to “understand, predict and ultimately control” the thermoelectric properties of TEMs. Understanding the influence of light-atom-doping and intrinsic defects on the electronic properties of the f-element containing TEMs are important questions that need to be answered to design and synthesize TEMs with optimal parameters. The purpose of this work is to understand how doping f-element based TEMs with light elements, such as carbon and boron, would impact the crystal structure, bond lengths, band gaps, and thermoelectric characteristics. High-quality f¬-element containing intermetallic single crystals are grown in metal flux media. In molten metal flux growth method, one or more low-melting metals such as aluminum, gallium, tin, and bismuth are used in excess as the reaction medium. As these are solution-state reactions, we can use lower temperatures than solid-state reactions and isolate the kinetic products instead of the thermodynamically stable compound. Then these are characterized using X-ray diffraction, scanning electron microscopy, energy dispersive spectroscopy, magnetization, electrical and thermal conductivities, and thermoelectric power.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2023

tinyIFD: A High-Throughput Binding Pose Refinement Workflow Through Induced-Fit Ligand Docking

A critical step in structure-based drug discovery is predicting whether and how a candidate molecule binds to a model of a therapeutic target. However, substantial protein side chain movements prevent current screening methods, such as docking, from accurately predicting the ligand conformations and require expensive refinements to produce viable candidates. Here, we present the development of a high-throughput and flexible ligand pose refinement workflow, called “tinyIFD”. The main features of the workflow include the use of specialized high-throughput, small-system MD simulation code mdgx.cuda and an actively learning model zoo approach. We show the application of this workflow on a large test set of diverse protein targets, achieving 66% and 76% success rates for finding a crystal-like pose within the top-2 and top-5 poses, respectively. We also applied this workflow to the SARS-CoV-2 main protease (M pro ) inhibitors, where we demonstrate the benefit of the active learning aspect in this workflow.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Workhorse minimally empirical dispersion-corrected density functional with tests for weakly bound systems: r 2 SCAN + rVV 10

SCAN+rVV10 has been demonstrated to be a versatile van der Waals (vdW) density functional that delivers good predictions of both energetic and structural properties for many types of bonding. Recently, the r 2 SCAN functional was devised as a revised form of SCAN with improved numerical stability. In this work, we refit the rVV10 functional to optimize the r 2 SCAN+rVV10 vdW density functional and test its performance for molecular interactions and layered materials. Our molecular tests demonstrate that r 2 SCAN+rVV10 outperforms its predecessor SCAN+rVV10 in both efficiency (numerical stability) and accuracy. This good performance is also found in lattice-constant predictions. In comparison with benchmark results from higher-level theories or experiments, r 2 SCAN+rVV10 yields excellent interlayer binding energies and phonon dispersions for layered materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Direct Prediction of Phonon Density of States With Euclidean Neural Networks

Abstract Machine learning has demonstrated great power in materials design, discovery, and property prediction. However, despite the success of machine learning in predicting discrete properties, challenges remain for continuous property prediction. The challenge is aggravated in crystalline solids due to crystallographic symmetry considerations and data scarcity. Here, the direct prediction of phonon density‐of‐states (DOS) is demonstrated using only atomic species and positions as input. Euclidean neural networks are applied, which by construction are equivariant to 3D rotations, translations, and inversion and thereby capture full crystal symmetry, and achieve high‐quality prediction using a small training set of examples with over 64 atom types. The predictive model reproduces key features of experimental data and even generalizes to materials with unseen elements, and is naturally suited to efficiently predict alloy systems without additional computational cost. The potential of the network is demonstrated by predicting a broad number of high phononic specific heat capacity materials. The work indicates an efficient approach to explore materials' phonon structure, and can further enable rapid screening for high‐performance thermal storage materials and phonon‐mediated superconductors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Large Destabilization of (TiVNb)-Based Hydrides via (Al, Mo) Addition: Insights from Experiments and Data-Driven Models

High-entropy alloys (HEAs) represent an interesting alloying strategy that can yield exceptional performance properties needed across a variety of technology applications, including hydrogen storage. Examples include ultrahigh volumetric capacity materials (BCC alloys → FCC dihydrides) with improved thermodynamics relative to conventional high-capacity metal hydrides (like MgH 2 ), but still further destabilization is needed to reduce operating temperature and increase system-level capacity. Here, in this work, we demonstrate efficient hydride destabilization strategies by synthesizing two new Al 0.05 (TiVNb) 0.95–x Mo x (x = 0.05, 0.10) compositions. We specifically evaluate the effect of molybdenum (Mo) addition on the phase structure, microstructure, hydrogen absorption, and desorption properties. Both alloys crystallize in a bcc structure with decreasing lattice parameters as the Mo content increases. The alloys can rapidly absorb hydrogen at 25 °C with capacities of 1.78 H/M (2.79 wt %) and 1.79 H/M (2.75 wt %) with increasing Mo content. Pressure-composition isotherms suggest a two-step reaction for hydrogen absorption to a final fcc dihydride phase. The experiments demonstrate that increasing Mo content results in a significant hydride destabilization, which is consistent with predictions from a gradient boosting tree data-driven model for metal hydride thermodynamics. Furthermore, improved desorption properties with increasing Mo content and reversibility were observed by in situ synchrotron X-ray diffraction, in situ neutron diffraction, and thermal desorption spectroscopy.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Elucidating and contrasting the mechanisms for Mg and Ca sulfate ion-pair formation with multi-level embedded quantum mechanics/molecular dynamics simulations

Here, solutions and minerals containing sulfate (SO 4 2- ), and Ca 2+ and Mg 2+ cations, are ubiquitous throughout the lithosphere and are significant components of seawater, thus presenting a prototypical system for the study of strong electrolytes and crystal nucleation mechanisms. However, despite their relative abundance, key questions remain unanswered about the most fundamental atomic-level steps of their mineralization pathways and aqueous dynamics. Here, we carry out enhanced sampling multi-level molecular dynamics (MD) embedded correlated wavefunction theory simulations to elucidate ion-pairing mechanisms for Mg–SO 4 and Ca–SO 4 in concentrated aqueous solution, accurately capturing effects arising from both structural dynamics and electron exchange–correlation. We predict contact-ion-pair formation to be barrierless and highly exoergic for Ca–SO 4 , in agreement with its minimal solubility, whereas for Mg–SO 4 , solvent-shared and contact ion pairs have similar free energies, qualitatively consistent with its higher solubility. Finally, we demonstrate that brief high-temperature pre-equilibration may be utilized to accelerate convergence of free energies in blue-moon-ensemble enhanced-sampling MD.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Too Many Materials and Too Many Applications: An Experimental Problem Waiting for a Computational Solution

Finding the best material for a specific application is the ultimate goal of materials discovery. However, there is also the reverse problem: when experimental groups discover a new material, they would like to know all the possible applications this material would be promising for. Computational modeling can aim to fulfill this expectation, thanks to the sustained growth of computing power and the collective engagement of the scientific community in developing more efficient and accurate workflows for predicting materials' performances. We discuss the impact that reproducibility and automation of the modeling protocols have on the field of gas adsorption in nanoporous crystals. We envision a platform that combines these tools and enables effective matching between promising materials and industrial applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Unconventional pressure-induced magnetic transitions and mechanical properties in non-magnetic LaFeSi

Here, this work demonstrates that two consecutive magnetic transitions occur in non-magnetic LaFeSi compound under negative hydrostatic pressure (small volume expansion). Electronic structure properties of LaFeSi were calculated using density functional theory (DFT) to understand the origin of these magnetic transitions. Mechanical properties of LaFeSi, especially elastic properties (shear and Young’s modulus, and Poisson’s ratio), were determined by DFT distortion calculations to elucidate lattice anisotropy and stability. The antiferromagnetic transition at –7.52 GPa, predicted in our work for LaFeSi, is consistent with the reported experimentally observed antiferromagnetic ground states in CeMnSi and LaMnSi compounds having similar crystal geometry. This finding of unexpected magnetic ordering predicts that layered non-magnetic materials, like LaFeSi, may become magnetically active with lattice expansion.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Modeling Single-Crystal Battery Materials: From Fundamental Understanding to Performance Evaluation

The performance of rechargeable batteries is fundamentally influenced by the physicochemical properties and microstructural features of their key material components. Recent experimental advancements have highlighted the potential of single-crystal (SC) morphologies to address inherent limitations of polycrystalline (PC) electrodes and solid-state electrolytes, offering tunable charge transport kinetics and improved cell cycling performance. Here, this review examines how state-of-the-art computational modeling, from atomistic and mesoscale to continuum-level approaches, including machine learning methodologies, has been utilized to investigate the critical factors governing the electrochemical behavior of SC battery materials. We explore how predictive modeling can elucidate the processing–structure–property–performance relationships of SC cathodes, anodes, and solid-state electrolytes, with a focus on unique SC characteristics such as crystallographic anisotropy, size effects, and facet-dependent properties. Additionally, we identify limitations in commonly used modeling techniques and discuss strategies to address these challenges. By integrating high-fidelity simulations with experimental insights, this review aims to outline a clear path for the rational design and optimization of SC battery components, paving the way for accelerated advancements in energy storage technologies.

Materials science↗
DOE OSTI2024

Solvation and oxidation effects on the crystal structure and morphology of tetraoxolene-based materials

Owing to their rich redox behavior and strong metal chelating ability, tetraoxolene ligands have been established as one of the foundational building blocks for multifunctional metal–organic materials. Here, we show how simple and often overlooked synthetic parameters can be used to control the structures of transition metal and lanthanide-based metal–tetraoxolene materials across multiple length scales. Through the synthesis of twelve new compounds, we provide a comprehensive survey detailing how the choice of solvent, initial ligand redox state, and in situ oxidant impact the local coordination geometry and chain architecture, as well as the crystal size and shape (e.g., rods vs. platelets). This work represents an important step towards the synthesis of new metal–tetraoxolene materials with predictable architectures and, therefore, targeted functionality.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Anomalous elastic softening in ferroelectric hafnia under pressure

his study employs first-principles density-functional theory (DFT) calculations to explore the elastic and mechanical properties of ferroelectric hafnia (HfO 2 ) in its polar orthorhombic 𝑃⁢𝑐⁢𝑎⁢2 1 phase under varying hydrostatic pressure conditions up to 30 GPa. Utilizing a plane-wave basis set and Perdew-Burke-Ernzerhof generalized-gradient approximation for solids in our DFT calculations, we investigate both pure and yttrium-substituted HfO 2 . Our findings reveal an anomalous reduction in the 𝐶 33 component of the elastic tensor with increasing pressure, which becomes significant above 15 GPa and signals a potential pressure-driven structural phase transition at higher pressure. The analysis of atomic displacements under pressure sheds light on the unusual mechanical behavior and phase stability of this material. Additionally, we observe a transition from an indirect band gap to a direct band gap with increasing pressure, which could have significant implications for optical applications. Here, the effects of yttrium substitution on the mechanical and electronic properties are further examined, revealing that yttrium substitution softens the elastic response of this material and reduces the electronic band gap. These results enhance our understanding of elastic and mechanical responses of ferroelectric hafnia and its potential for applications in microelectronics, piezoelectric devices, and nonvolatile ferroelectric random-access memories. Further experimental validation is recommended to confirm our predictions and explore the practical implications of the observed phase transitions and electronic behavior of the ferroelectric hafnia under high-pressure conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

ReaxFF Parameter Set for Boron Clusters and Icosahedral Boron Crystals: Comparison with Density Functional Theory and Machine-Learning Potentials

Icosahedral boron materials, which include regular icosahedra of 12 boron atoms have gained increasing attention due to their potential applications as superhard materials, semiconductors, and energy storage media. However, the synthesis of high quality crystals of these materials has been a major barrier to the development of these applications. To enable computational prediction of synthesis conditions yielding high-quality icosahedral boron crystals, herein we tested and refined a set of ReaxFF parameters for the nucleation and growth of such crystals. We focused on matching the relative energies of small boron clusters obtained by density functional theory since such small clusters and similar motifs are likely present in crystal nuclei and at the interface of growing crystals. Using a training set of B 80 clusters, including a low-energy core–shell structure containing a B 12 icosahedron core and a high-energy single-shell structure produced in preliminary ReaxFF simulations, the ReaxFF parameter set was refined to better reproduce energies calculated by density functional theory (DFT). Among existing ReaxFF parameter sets and the machine-learning interatomic potentials MACE-MP-0, MACE-MP-0b3, MACE-MPA-0, PFP v7.0.0, and SevenNet-MF-ompa, only our new parameter set and PFP v7.0.0 correctly ranked these B 80 clusters. This refinement led to improved agreement with DFT for a test set of 58 clusters consisting of 8–103 boron atoms. Furthermore, our refined parameter set yielded greater local icosahedral structure than the previously existing ReaxFF parameter set for larger scale simulations of crystallization from supercooled liquid boron. Additionally, simulations of solid boron in contact with molten nickel using our refined ReaxFF parameters yielded a boron solubility value that agrees moderately well with experimental expectations, while the previous boron parameters gave a value that was much too low.

boron↗
DOE OSTI2020

Mie-Resonant Three-Dimensional Metacrystals

Optical metamaterials, engineered to exhibit electromagnetic properties not found in natural materials, may enable new light-based applications including cloaking and optical computing. While there have been significant advances in the fabrication of two-dimensional metasurfaces, planar structures create nontrivial angular and polarization sensitivities, making omnidirectional operation impossible. Although three-dimensional (3D) metamaterials have been proposed, their fabrication remains challenging. In this work, we use colloidal crystal engineering with DNA to prepare isotropic 3D metacrystals from Au nanocubes. We show that such structures can exhibit refractive indices as large as ~8 in the mid-infrared, far greater than that of common high-index dielectrics. Additionally, we report the first observation of multipolar Mie resonances in metacrystals with well-formed habits, occurring in the mid-infrared for sub-micron metacrystals, which we measured using synchrotron infrared microspectroscopy. Finally, we predict that arrays of metacrystals could exhibit negative refraction. The results present a promising platform for engineering devices with unnatural optical properties.

3D metamaterials↗
DOE OSTI2021

X-ray diffraction and equation of state of the C–S–H room-temperature superconductor

X-ray diffraction indicates that the structure of the recently discovered carbonaceous sulfur hydride (C–S–H) room-temperature superconductor is derived from previously established van der Waals compounds found in the H 2 S–H 2 and CH 4 –H 2 systems. Crystals of the superconducting phase were produced by a photochemical synthesis technique, leading to the superconducting critical temperature T c of 288 K at 267 GPa. X-ray diffraction patterns measured from 124 to 178 GPa, within the pressure range of the superconducting phase, are consistent with an orthorhombic structure derived from the Al 2 Cu-type determined for (H 2 S) 2 H 2 and (CH 4 ) 2 H 2 that differs from those predicted and observed for the S–H system at these pressures. The formation and stability of the C–S–H compound can be understood in terms of the close similarity in effective volumes of the H 2 S and CH 4 components, and denser carbon-bearing S–H phases may form at higher pressures. Furthermore, the results are crucial for understanding the very high superconducting T c found in the C–S–H system at megabar pressures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Ferromagnetic MnBi 4 Te 7 obtained with low-concentration Sb doping: A promising platform for exploring topological quantum states

The tuning of the magnetic phase, chemical potential, and structure is crucial to observe diverse exotic topological quantum states in Mn Bi 2 Te 4 (Bi 2 Te 3 ) m (m = 0–3). Here we show a ferromagnetic (FM) phase with a chiral crystal structure in Mn (Bi 1–x Sb x ) 4 Te 7 , obtained via tuning the growth conditions and Sb concentration. Unlike previously reported Mn (Bi 1–x Sb x ) 4 Te 7 , which exhibits FM transitions only at high Sb doping levels, our samples show FM transitions (T C = 13.5 K) at 15%–27% doping levels. Furthermore, our single-crystal x-ray-diffraction structure refinements find Sb doping leads to a chiral structure with the space group of P3, contrasted with the centrosymmetric $P\bar{3}m1$ crystal structure of the parent compound MnBi 4 Te 7 . Through angle-resolved photoemission spectroscopy measurements, we also demonstrated that the nontrivial band topology is preserved in the Sb-doped FM samples. Given that the nontrivial band topology of this system remains robust for low Sb doping levels, our success in making FM Mn(Bi 1–x Sb x ) 4 Te 7 with x = 0.15, 0.175, 0.2, and 0.27 paves the way for realizing the predicted topological quantum states, such as the axion insulator and Weyl semimetals. Additionally, we also observed magnetic glassy behavior in both antiferromagnetic MnBi 4 Te 7 and FM Mn (Bi 1–x Sb x )4 Te 7 samples, which we believe originates from cluster spin-glass phases coexisting with long-range antiferromagnetic/FM orders. Further, we have also discussed how the antisite Mn ions impact the interlayer magnetic coupling and how FM interlayer coupling is stabilized in this system.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A constitutive model incorporating grain refinement strengthening on metallic alloys

Surface nano-crystallization techniques have been recently developed as one of the most efficient ways to optimize materials’ structure, and therefore develop the local and global mechanical behavior as to increase strength without compromising ductility. Here, in this work, we present a constitutive model incorporating grain refinement hardening to simulate the nano-crystallization technique, specifically, surface mechanical attrition treatment. The computation is implemented using user-defined VUMAT subroutines. As an example of its application, a geometry model with full coverage of random impacts are employed. The results show that the model has rather precise predictability of grain size evolution during plastic deformation. The readily embedded with a computational code of material dynamics enables this novel model to be a promising tool to study the dynamic evolution of microstructures under plastic deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Regioselectivity mechanism of the Thunbergia alata Δ 6 -16:0-acyl carrier protein desaturase

Plant plastidial acyl–acyl carrier protein (ACP) desaturases are a soluble class of diiron-containing enzymes that are distinct from the diiron-containing integral membrane desaturases found in plants and other organisms. The archetype of this class is the stearoyl-ACP desaturase which converts stearoyl-ACP into oleoyl (18:1Δ 9 cis )-ACP. Several variants expressing distinct regioselectivity have been described including a Δ 6 -16:0-ACP desaturase from black-eyed Susan vine ( Thunbergia alata ). We solved a crystal structure of the T. alata desaturase at 2.05 Å resolution. Using molecular dynamics (MD) simulations, we identified a low-energy complex between 16:0-ACP and the desaturase that would position C6 and C7 of the acyl chain adjacent to the diiron active site. The model complex was used to identify mutant variants that could convert the T. alata Δ 6 desaturase to Δ 9 regioselectivity. Additional modeling between ACP and the mutant variants confirmed the predicted regioselectivity. Furthermore, to validate the in-silico predictions, we synthesized two variants of the T. alata desaturase and analyzed their reaction products using gas chromatography-coupled mass spectrometry. Assay results confirmed that mutants designed to convert T. alata Δ 6 to Δ 9 selectivity exhibited the predicted changes. In complementary experiments, variants of the castor desaturase designed to convert Δ 9 to Δ 6 selectivity lost some of their Δ9 desaturation ability and gained the ability to desaturate at the Δ 6 position. The computational workflow for revealing the mechanistic understanding of regioselectivity presented herein lays a foundation for designing acyl-ACP desaturases with novel selectivities to increase the diversity of monoenes available for bioproduct applications.

59 BASIC BIOLOGICAL SCIENCES↗