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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
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 OSTI2021

The Li–F–H ternary system at high pressures

Evolutionary crystal structure prediction searches have been employed to explore the ternary Li-F-H system at 300 GPa. Metastable phases were uncovered within the static lattice approximation, with LiF 3 H 2 , LiF 2 H, Li 3 F 4 H, LiF 4 H 4 , Li 2 F 3 H and LiF 3 H lying within 50 meV/atom of the 0 K convex hull. All of these phases contain H n F¯ n+1 (n = 1; 2) anions, and Li + cations. Other structural motifs such as LiF slabs, $H$ $^{+}_{3}$ molecules and F δ- ions are present in some of the low enthalpy Li-F-H structures. The bonding within the H n F¯ n+1 molecules, which may be bent or linear, symmetric or asymmetric, is analyzed. The five phases closest to the hull are insulators, while LiF 3 H is metallic and predicted to have a vanishingly small superconducting critical temperature. Li 3 F 4 H is predicted to be stable at zero pressure. Furthermore, this study lays the foundation for future investigations of the role of temperature and anharmonicity on the stability and properties of compounds and alloys in the Li-F-H ternary system.

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 OSTI2021

Structural, Vibrational, and Electronic Properties of 1D-TlInTe 2 under High Pressure: A Combined Experimental and Theoretical Study

Analogous to 2D layered transition-metal dichalcogenides, the TlSe family of quasi-one dimensional chain materials with the Zintl-type structure exhibits novel phenomena under high pressure. In the present work, we have systematically investigated the high-pressure behavior of TlInTe 2 using Raman spectroscopy, synchrotron X-ray diffraction (XRD), and transport measurements, in combination with first principles crystal structure prediction (CSP) based on evolutionary approach. We found that TlInTe2 undergoes a pressure-induced semiconductor-to-semimetal transition at 4 GPa, followed by a superconducting transition at 5.7 GPa (with T c = 3.8 K). An unusual giant phonon mode (A g ) softening appears at ~10–12 GPa as a result of the interaction of optical phonons with the conduction electrons. The high-pressure XRD and Raman spectroscopy studies reveal that there is no structural phase transitions observed up to the maximum pressure achieved (33.5 GPa), which is in agreement with our CSP calculations. In addition, our calculations predict two high-pressure phases above 35 GPa following the phase transition sequence as I4/mcm (B37) → Pbcm → Pm$\bar{3}$m (B2). Electronic structure calculations suggest Lifshitz (L1 & L2-type) transitions near the superconducting transition pressure. Furthermore, our findings on TlInTe 2 open up a new avenue to study unexplored high-pressure novel phenomena in TlSe family induced by Lifshitz transition (electronic driven), giant phonon softening, and electron–phonon coupling.

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

Tuning chemical precompression: Theoretical design and crystal chemistry of novel hydrides in the quest for warm and light superconductivity at ambient pressures

Over the past decade, a combination of crystal structure prediction techniques and experimental synthetic work has thoroughly explored the phase diagrams of binary hydrides under pressure. The fruitfulness of this dual approach is demonstrated in the recent identification of several superconducting hydrides with T c s approaching room temperature. We start with an overview of the computational procedures for predicting stable structures and estimating their propensity for superconductivity. A survey of phases with high T c reveals some common structural features that appear conducive to the strong coupling of the electronic structure with atomic vibrations that leads to superconductivity. We discuss the stability and superconducting properties of phases containing two of these—molecular H 2 units mixed with atomic H and hydrogenic clathrate-like cages—as well as more unique motifs. Lastly, we argue that ternary hydride phases, whose exploration is still in its infancy, are a promising route to achieve simultaneous superconductivity at high temperatures and stability at low pressures. Several ternary hydrides arise from the addition of a third element to a known binary hydride structure through site mixing or onto a new site, and several more are based on altogether new structural motifs.

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 OSTI2020

Disorderly Conduct of Benzamide IV: Crystallographic and Computational Analysis of High Entropy Polymorphs of Small Molecules

Benzamide, a simple derivative of benzoic acid and a common intermediate of pharmaceutical compounds, was reported to form two polymorphs in 1832, but the single crystal structure of the more stable form was not solved until 1959. Nearly 50 years later, the second form was characterized by powder diffraction, followed shortly thereafter by characterization of a third form, a polytype of the most thermodynamically stable Form I. These two new forms, Forms II and III, are metastable. Herein, we describe a fourth polymorph, Form IV, discovered by melt crystallization concurrently with its crystallization under confinement at small length scales (<10 nm), where it is stable indefinitely. Form III exists under confinement in larger pores, and melting point data for different pore sizes corroborate the existence of Form IV below 10 nm. Form IV is highly disordered, precluding indexing of powder diffraction data other than hk0 reflections. Nonetheless, a combination of powder X-ray diffraction and computational crystal structure prediction reveals that Form IV contains a 2D motif resembling that of Form II, but with longrange order in the third dimension masked by ubiquitous stacking faults. This approach relies on distilling a large number of candidate structures to a few possible disorder models based on benzamide tetrads that organize in 2D parquet-like tiles, with organization along the third dimension, that can be modeled with various stacking fault configurations having distinct intermolecular interactions and translations in the dimension orthogonal to the tiling planes. These observations reveal a bewildering crystallographic complexity for such a simple molecule. Nonetheless, the approach described herein demonstrates that challenging structures that may be abandoned prematurely because of poor crystallinity, twinning, or disorder can be solved.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2020

Route to high-$T_c$ superconductivity via CH 4 -intercalated H 3 S hydride perovskites

While exploring potential superconductors in the C-S-H ternary system using first-principles crystal structure prediction methods, we uncovered a class of hydride perovskites based on the intercalation of methane into an H 3 S framework. These intriguing H 3 S–CH 4 structures emerge as metastable at ~100 GPa. Further, electron-phonon coupling calculations indicate that phases with CSH 7 stoichiometry are potential superconductors with $T_c$ values ranging from 100 K to 190 K at megabar pressures. The results are expected to guide the experimental search for new high-$T_c$ superconductors, including those stable at lower pressures than previously documented superconducting hydrides such as H 3 S and LaH 10 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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 OSTI2021

Benchmarking Coordination Number Prediction Algorithms on Inorganic Crystal Structures

Coordination numbers and geometries form a theoretical framework for understanding and predicting materials properties. Algorithms to determine coordination numbers automatically are increasingly used for machine learning and automatic structural analysis. In this work, we introduce MaterialsCoord, a benchmark suite containing 56 experimentally-derived crystal structures (spanning elements, binaries, and ternary compounds) and their corresponding coordination environments as described in the research literature. We also describe CrystalNN, a novel algorithm for determining near neighbors. We compare CrystalNN against 7 existing near-neighbor algorithms on the MaterialsCoord benchmark, finding CrystalNN to perform similarly to several well-established algorithms. For each algorithm, we also assess computational demand and sensitivity towards small perturbations that mimic thermal motion. Finally, we investigate the similarity between bonding algorithms when applied to the Materials Project database. Finally, we expect that this work will aid the development of coordination prediction algorithms as well as improve structural descriptors for machine learning and other applications.

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

Prediction of chlorine and fluorine crystal structures at high pressure using symmetry driven structure search with geometric constraints

The high-pressure properties of fluorine and chlorine are not yet well understood because both are highly reactive and volatile elements, which have made conducting diamond anvil cell and x-ray diffraction experiments a challenge. Here, we use ab initio methods to search for stable crystal structures of both elements at megabar pressures. We demonstrate how symmetry and geometric constraints can be combined to efficiently generate crystal structures that are composed of diatomic molecules. Furthermore, our algorithm extends the symmetry driven structure search method by adding constraints for the bond length and the number of atoms in a molecule while still maintaining generality. As a method of validation, we have tested our approach for dense hydrogen and reproduced the known molecular structures of Cmca-12 and Cmca-4. We apply our algorithm to study chlorine and fluorine in the pressure range of 10 GPa–4000 GPa while considering crystal structures with up to 40 atoms per unit cell. We predict chlorine to follow the same series of phase transformations as elemental iodine from Cmca to Immm to Fm3¯m, but at substantially higher pressures. We predict fluorine to transition from a C2/c to Cmca structure at 70 GPa, to a novel orthorhombic and metallic structure with P42/mmc symmetry at 2500 GPa, and finally to its cubic analog form with Pm3¯n symmetry at 3000 GPa.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗