Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Liquid crystals”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Generalized representative structures for atomistic systems

A new method is presented to generate atomic structures that reproduce the essential characteristics of arbitrary material systems, phases, or ensembles. Previous methods allow one to reproduce the essential characteristics (e.g. the chemical disorder) of a large random alloy within a small crystal structure. The ability to generate small representations of random alloys, along with the restriction to crystal systems, results from using the fixed-lattice cluster correlations to describe structural characteristics. A more general description of the structural characteristics of atomic systems is obtained using complete sets of atomic environment descriptors. These are used within for generating representative atomic structures without restriction to fixed lattices. A general data-driven approach is provided here utilizing the atomic cluster expansion (ACE) basis. The N-body ACE descriptors are a complete set of atomic environment descriptors that span both chemical and spatial degrees of freedom and are used within for describing atomic structures. The generalized representative structure (GRS) method presented within generates small atomic structures that reproduce ACE descriptor distributions corresponding to arbitrary structural and chemical complexity. It is shown that systematically improvable representations of crystalline systems on fixed parent lattices, amorphous materials, liquids, and ensembles of atomic structures may be produced efficiently through optimization algorithms. With the GRS method, we highlight reduced representations of atomistic machine-learning training datasets that contain similar amounts of information and small 40–72 atom representations of liquid phases. The ability to use GRS methodology as a driver for informed novel structure generation is also demonstrated. The advantages over other data-driven methods and state-of-the-art methods restricted to high-symmetry systems are highlighted.

atomic cluster expansion↗
DOE OSTI2022

Thermoplasticity of metallic glasses: Processing and applications

Owing to their extremely low shear viscosity and fragility, bulk metallic glasses (BMGs) are envisioned as new net-shaping material replacing conventionally used crystalline metals for various applications. The first part of this review describes the general technique and kinetics of thermoplastic forming. The second part elaborates on the thermal processing capability of metallic glasses from atomic- to meter-scale. The micro-/nano-structure obtained after various quenching rates followed by temperature, time and heating rate dependent thermoplastic forming is described in the third section. The deformation behavior and flow kinetics of BMGs in terms of composition, mold features, and applied conditions are elaborated in the fourth section. The variation of volume and enthalpy at the supercooled liquid region, strength retention of metallic glasses compared to other conventional metals and alloys, and kinetics of liquid fragility are given in the fifth section. Finite element modeling and molecular dynamics simulations of high-temperature deformation in BMGs are presented in the sixth part. Finally, thermally-formed BMGs used for different applications, including energy, biomedical and micro-optics, are presented in the final part. Altogether, this review provides an overview of shaping capabilities and modifications in the macro-scale properties and short-to-medium range order of BMGs upon thermoplastic forming.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Reverse order-disorder transition of Janus particles confined in two dimensions

In this study Janus particles with different patch sizes, confined to two dimensions, generate a series of patterns of interest to the field of nanoscience. Here we observe reverse melting, where for some densities the system melts under cooling. For a broad range of hydrophobic patch sizes ( 60° < θ 0 < 90°), a reentrant transition from solid to liquid and then to an ordered phase emerges as temperature (T) decreases due to the formation of rhombus chains at low T. This reentrant phase has pseudo long-range orientational order but short-range translational order, similar to a hexatic phase. Our work provides guidelines to study the melting and assembly of Janus particles in two dimensions, as well as mechanisms to generate phases with specific symmetry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

NbZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys niobium-zirconium (Nb-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Nb and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on NERSC-Perlmutter and OLCF-Summit using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Cases 1-10 have been run without generating the file 'vaspout.h5'. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for binary alloys is AXBY, where A and B refer to the constituents, whereas X and Y are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: FILES contained in each subdirectory with name "case-N" where N ranges between 11 and 100, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. 0.POSCAR: input file that defines the atomic structure of a system 4. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 7. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. OUTCAR: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system FILES contained in each subdirectory with name "case-N" where N ranges between 1 and 10, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. {ID}.POSCAR: input file that defines the atomic structure of a system at the beginning of ID execution of geometry optimization with PREC=LOW 4. {ID}.CONTCAR: output file that provides the atomic positions and cell parameters at the end of ID execution of geometry optimization with PREC=LOW in the INCAR file 5. {ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of geometry optimization that has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.{ID}.out: file with diagnostic information about the execution of the ID execution of the geometry optimization with precision variable set to PREC=Low in the INCAR file 7. N{ID}.POSCAR: input file that defines the atomic structure of a system after the geometry optimization run at low precision. This represents the input for the ID execution of the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. N{ID}.CONTCAR: output file that provides the atomic positions and cell parameters after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 9. N{ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.{ID}.out: file with diagnostic information about the ID execution of geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216. REFERENCES (1) Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994) (3) Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. materials science 6, 15–50 (1996) (4) Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. review B 54, 11169 (1996) (5) Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2024

TaZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys tantalum-zirconium (Ta-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Ta and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on NERSC-Perlmutter and OLCF-Summit using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Cases 1-10 have been run without generating the file 'vaspout.h5'. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for binary alloys is AXBY, where A and B refer to the constituents, whereas X and Y are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: FILES contained in each subdirectory with name "case-N" where N ranges between 11 and 80, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. 0.POSCAR: input file that defines the atomic structure of a system 4. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 7. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. OUTCAR: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system FILES contained in each subdirectory with name "case-N" where N ranges between 1 and 10 and between 81 and 100, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. {ID}.POSCAR: input file that defines the atomic structure of a system at the beginning of ID execution of geometry optimization with PREC=LOW 4. {ID}.CONTCAR: output file that provides the atomic positions and cell parameters at the end of ID execution of geometry optimization with PREC=LOW in the INCAR file 5. {ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of geometry optimization that has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.{ID}.out: file with diagnostic information about the execution of the ID execution of the geometry optimization with precision variable set to PREC=Low in the INCAR file 7. N{ID}.POSCAR: input file that defines the atomic structure of a system after the geometry optimization run at low precision. This represents the input for the ID execution of the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. N{ID}.CONTCAR: output file that provides the atomic positions and cell parameters after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 9. N{ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.{ID}.out: file with diagnostic information about the ID execution of geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216. REFERENCES (1) Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994) (3) Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. materials science 6, 15–50 (1996) (4) Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. review B 54, 11169 (1996) (5) Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2020

Orientation-dependent stabilization of MgCr 2 O 4 spinel thin films

AB 2 O 4 normal spinels with a magnetic B site can host a variety of magnetic and orbital frustrations leading to spin-liquid phases and field-induced phase transitions. Here, we report the epitaxial growth of (111)-oriented MgCr 2 O 4 thin films. By characterizing the structural and electronic properties of films grown along the (001) and (111) directions, the influence of growth orientation has been studied. Despite distinctly different growth modes observed during deposition, the comprehensive characterization reveals no measurable disorder in the cation distribution nor multivalency issue for Cr ions in either orientation. Contrary to a naive expectation, the (111) stabilized films exhibit a smoother surface and a higher degree of crystallinity than (001)-oriented films. The preference in growth orientation is explained within the framework of heteroepitaxial stabilization in connection to a significantly lower (111) surface energy. Furthermore, these findings open broad opportunities in the fabrication of two-dimensional kagome-triangular heterostructures with emergent magnetic behavior inaccessible in bulk crystals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Relaxation dynamics of water in the vicinity of cellulose nanocrystals

Water molecules near cellulose nanocrystals (CNCs; produced via the sulfuric acid-catalyzed hydrolysis of wood pulp) are believed to relax slower than those in the bulk liquid, which may result in unique properties of CNC aqueous dispersions. Herein this study analyzed the polarization behavior of water molecules in CNC aqueous dispersions and other reference samples using a dielectric relaxation spectroscopy (DRS) technique in the microwave frequency range (0.2–20 GHz). As the CNC concentration increases, two slow relaxation components become prominent. The comparison with DRS data of aqueous dispersions of nanoporous silica, polyvinyl alcohol (PVA), and hairy CNC (HCNC) with amorphous chains protruding from both ends suggested that these slow relaxation modes of water near CNC surfaces cannot be attributed to direct hydrogen bonding interactions with the hydroxyl (OH) groups exposed and immobilized at the solid surface. Instead, they are similar to the water molecules interacting with OH groups attached to flexible polymer chains. Molecular dynamics (MD) simulations of the polarization behavior of water near the (110) facet of cellulose Iβ crystals confirmed that the interactions of water molecules with the cellulose crystal surface do not cause slower relaxations in the frequency range studied via the DRS. These results indicated that the CNC surface cannot be depicted with the crystallographic facets of cellulose Iβ; instead, it resembles a polymer-brush surface on which the short glucan residues or fragments of the strong acid-catalyzed hydrolysis process are swollen and extended into the aqueous phase.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Hydrazinoacetic acid is a biosynthetic precursor of the bacterially produced nitramine, N -nitroglycine

Nitramines [R(R′)N–NO 2 ; R,R′=H or alkyl] are valuable synthetic products, but knowledge of the biosynthetic processes that generate these compounds is limited. This work sought to elucidate the biosynthesis of a nitramine natural product, N-nitroglycine (NNG) by Streptomyces noursei . Stable isotope studies showed that S. noursei cells supplemented with L-(ε- 15 N)lysine, ( 15 N)glycine, or ( 13 C)hydrazinoacetic acid (HAA) incorporated 67%, 88%, and 67% of the isotope label into NNG, respectively, indicating that these compounds are biosynthetic precursors of NNG. Liquid chromatography coupled tandem mass spectrometry (LC-MS/MS) of 15 N-Lys-labeled NNG confirmed that the nitro nitrogen of NNG originates from Lys. Bioinformatics analysis of the S. noursei genome showed evidence for a biosynthetic gene cluster (BGC) that contained machinery for HAA biosynthesis ( nngKLM ), consistent with the results of the isotope labeling. In vitro reconstitution of the gene products produced HAA. The borders of this BGC were defined by cross-referencing the predicted BGC with previously published differential proteomics data. Furthermore, we show that azaserine is produced alongside NNG in S. noursei cultures, linking the two biosynthetic pathways via a proposed nitrosamine biosynthetic intermediate. Finally, the oxygen balance for NNG is −20.2% for the formation of carbon dioxide (CO 2 ), which is comparable to that of hexahydro-1,3,5- trinitro-1,3,5-triazene (common name: RDX; −21.6%). Crystal structure data of NNG indicate that the unit crystalizes as a pure material, not a hydrate, suggesting a favorable energetic crystallization phase. The combined results suggest a route that, with further development, could lead to sustainable production of energetic nitramines via synthetic biology or biocatalytic approaches.

biosynthesis↗
DOE OSTI2023

Universal and interpretable classification of atomistic structural transitions via unsupervised graph learning

Materials processing often occurs under extreme dynamic conditions leading to a multitude of unique structural environments. These structural environments generally occur at high temperatures and/or high pressures, often under non-equilibrium conditions, which results in drastic changes in the material's structure over time. Computational techniques, such as molecular dynamics simulations, can probe the atomic regime under these extreme conditions. However, characterizing the resulting diverse atomistic structures as a material undergoes extreme changes in its structure has proved challenging due to the inherently non-linear relationship between structures as large-scale changes occur. Here, we introduce SODAS++, a universal graph neural network framework, that can accurately and intuitively quantify the atomistic structural evolution corresponding to the transition between any two arbitrary phases. We showcase SODAS++ for both solid–solid and solid–liquid transitions for systems of increasing geometric and chemical complexity, such as colloidal systems, elemental Al, rutile and amorphous TiO 2 , and the non-stoichiometric ternary alloy Ag 26 Au 5 Cu 19 . Finally, we show that SODAS++ can accurately quantify all transitions in a physically interpretable manner, showcasing the power of unsupervised graph neural network encodings for capturing the complex and non-linear pathway, a material's structure takes as it evolves.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Strain release by 3D atomic misfit in fivefold twinned icosahedral nanoparticles with amorphization and dislocations

Multiple twinning to form fivefold twinned nanoparticles in crystal growth is common and has attracted broad attention ranging from crystallography research to physical chemistry and materials science. Lattice-misfit strain and defects in multiple twinned nanoparticles (MTP) are key to understand and tailor their electronic properties. However, the structural defects and related strain distributions in MTPs are poorly understood in three dimensions (3D). Here, we show the 3D atomic misfit and strain relief mechanism in fivefold twinned icosahedral nanoparticles with amorphization and dislocations by using atomic resolution electron tomography. We discover a two-sided heterogeneity in variety of structural characteristics. A nearly ideal crystallographic fivefold face is always found opposite to a less ordered face, forming Janus-like icosahedral nanoparticles with two distinct hemispheres. The disordered amorphous domains release a large amount of strain. Molecular dynamics simulations further reveal the Janus-like icosahedral nanoparticles are prevalent in the MTPs formed in liquid-solid phase transition. This work provides insights on the atomistic models for the modelling of formation mechanisms of fivefold twinned structures and computational simulations of lattice distortions and defects. We anticipate it will inspire future studies on fundamental problems such as twin boundary migration and kinetics of structures in 3D at atomic level.

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

Understanding supercooling mechanism in sodium sulfate decahydrate phase-change material

Salt hydrate-based phase-change materials are considered promising for future heat storage applications in residential heating/cooling systems. Smooth phase transition from the liquid to solid phase and vice versa is essential for effective heat exchanger; however, supercooling in salt hydrates delays the onset of liquid–solid phase transition. We investigate the molecular level mechanism of supercooling in sodium sulfate decahydrate (SSD). SSD is a complex salt hydrate whose properties are governed by electrostatic forces that include pure Coulombic interactions as well as hydrogen bonds. Experimentally, we examine the importance of a nucleator in reducing supercooling temperatures. We investigated the effect of various mass concentrations of a borax nucleator on a decrease of supercooling temperatures. Molecular dynamics simulation techniques are used to obtain a basic understanding of supercooling in SSD. We observe that by introducing borax as a nucleator, there is a decrease in the supercooling temperature before nucleation. Our molecular dynamics simulations show that long-range electrostatics between sodium and sulfate ion pairs and that with polar water molecules is responsible for delayed nucleation in SSD that results in supercooling, and also, dynamics of charged molecules slows down. The lack of crystallization leads to amorphous structures in supercooled SSD.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Data-Driven Kinetic Reaction Networks for Separation Chemistry

Understanding complex, multistep chemical reactions at the molecular level is a major challenge whose solution would greatly benefit the design and optimization of numerous chemical processes. The separation of rare-earth (4f) and actinide (5f) elements is an example where improving our chemical understanding is important for designing and optimizing new chemistries, even with a limited number of observations. Here, in this work, we leverage data-driven artificial intelligence and machine-learning approaches to develop kinetic reaction networks that describe the liquid–liquid extraction mechanism of uranium using N,N-di-2-ethylhexyl-isobutyramide (DEHiBA). Specifically, we compare and contrast the properties of two classes of models: (1) purely data-driven models that are regularized using chemistry-agnostic, L1 regression and (2) chemistry-informed models that are regularized using relative reaction energies provided by quantum mechanical calculations. We observe that purely data-driven models are unbiased, simple, and accurate in their predictions of experimental measurements when provided with sufficient data but are difficult to fully constrain and interpret. In contrast, chemistry-informed models exhibit significantly improved chemical interpretability and consistency, providing a detailed description of the separation process while achieving high accuracy through ensemble averaging. Overall, the dominant species predicted to be extracted into the organic phase is UO 2 (NO 3 ) 2 (DEHiBA) 2 , agreeing with experimental slope analysis, thermodynamic modeling, EXAFS, and crystal structures. This work demonstrates that leveraging the fundamental structure of the problem can lead to efficient learning schemes that provide both accurate predictions and chemical insights at a low computational cost.

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

High-entropy engineering of the crystal and electronic structures in a Dirac material

Dirac and Weyl semimetals are a central topic of contemporary condensed matter physics, and the discovery of new compounds with Dirac/Weyl electronic states is crucial to the advancement of topological materials and quantum technologies. Here we show a widely applicable strategy that uses high configuration entropy to engineer relativistic electronic states. We take the A MnSb 2 ( A = Ba, Sr, Ca, Eu, and Yb) Dirac material family as an example and demonstrate that mixing of Ba, Sr, Ca, Eu and Yb at the A site generates the compound (Ba 0.38 Sr 0.14 Ca 0.16 Eu 0.16 Yb 0.16 )MnSb 2 (denoted as A 5 MnSb 2 ), giving access to a polar structure with a space group that is not present in any of the parent compounds. A 5 MnSb 2 is an entropy-stabilized phase that preserves its linear band dispersion despite considerable lattice disorder. Although both A 5 MnSb 2 and A MnSb 2 have quasi-two-dimensional crystal structures, the two-dimensional Dirac states in the pristine A MnSb 2 evolve into a highly anisotropic quasi-three-dimensional Dirac state triggered by local structure distortions in the high-entropy phase, which is revealed by Shubnikov–de Haas oscillations measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Lowering the Activation Barriers for Lithium-Ion Conductivity through Orientational Disorder in the Cyanide Argyrodite Li6PS5CN

Rapid advancements in safe and high-energy-density energy storage are predicated on identifying new solid-state ion conductors with low activation energies and high ionic conductivities for all-solid-state battery technologies. Halide argyrodites are among some of the top candidates for solid-state electrolytes, as they can achieve ionic conductivities that approach liquid electrolytes. Incorporating dynamic pseudohalide species in argyrodite solid electrolytes presents an exciting opportunity to exploit lattice dynamics as a design principle to modulate the ion conduction properties of solid-state ion conductors. In the present study, we have prepared the new argyrodite Li 6 PS 5 CN containing orientationally disordered cyanide ions. The new cyanide argyrodite Li 6 PS 5 CN exhibits an activation barrier to Li-ion transport of 471 ± 25 meV and a room-temperature ionic conductivity of 6(2) × 10 –5 S cm –1 in comparison to the activation barrier of 502 ± 16 meV and an ionic conductivity of 2.3(1) × 10 –4 S cm –1 measured for the bromide analogue Li 6 PS 5 Br. Structural studies of both compounds by high-resolution X-ray diffraction indicate that Li 6 PS 5 CN and Li 6 PS 5 Br adopt nearly identical crystal structures with similar lattice parameters, which indicates that lower activation barriers in Li 6 PS 5 CN arise due to the cyanide ion itself rather than due to changes in the geometry of conduction pathways in the local lithium environment. The orientational disorder of the quadrupolar cyanide ion in Li 6 PS 5 CN points to a complex interplay of lattice polarizability and molecular dynamics that lower the activation barrier for lithium-ion conductivity in the cyanide argyrodite.

25 ENERGY STORAGE↗
DOE OSTI2022

The Jefferson Lab Eta Factory (JEF) experiment

The Jefferson Lab Eta Factory (JEF) experiment is an upcoming experiment designed to run in Hall D at Jefferson Lab using an upgraded GlueX spectrometer to study various decays of the ? meson. The GlueX spectrometer consists of a ?2 Tesla solenoid magnet housing a liquid hydrogen target and drift chambers used for tracking charged particles and an array of lead glass blocks (the Forward Calorimeter (FCAL)) downstream of the magnet for detecting neutral particles. The decays of the ? meson can be used to measure the light quark mass ratio via the ???+???0 channel and allow access to higher-order terms in Chiral Perturbation Theory via the ???0?? channel. The decays can also be used to constrain new charge-conjugation violating/parity conserving (CVPC) reactions and to search for signatures of dark matter. In particular, the ???0?? channel can be used to search for lepto-phobic dark vector (B) bosons in the reaction ??B? (B??0?) or dark scalar (S) bosons in the reaction ???0S (S???). Studying the rare radiative decay channel ???0?? requires replacing the 4×4×45 cm3 lead glass blocks in the inner region of the FCAL with 2×2×20 cm3 lead tungstate crystals, which will provide twice better position and energy resolution. This talk will describe the JEF physics program and the upgrade to the FCAL.

Taylor, Simon↗
DOE OSTI2022

Experimental observation of open structures in elemental magnesium at terapascal pressures

Investigating how solid matter behaves at enormous pressures, such as those found in the deep interiors of giant planets, is a great experimental challenge. Over the past decade, computational predictions have revealed that compression to terapascal pressures may bring about counter-intuitive changes in the structure and bonding of solids as quantum mechanical forces grow in influence. Although this behavior has been observed at modest pressures in the highly compressible light alkali metals, it has not been established whether it is commonplace among high-pressure solids more broadly. We used shaped laser pulses at the National Ignition Facility to compress elemental Mg up to 1.3 TPa, which is approximately four times the pressure at the Earth’s core. By directly probing the crystal structure using nanosecond-duration X-ray diffraction, we found that Mg changes its crystal structure several times with non-close-packed phases emerging at the highest pressures. Our results demonstrate that phase transformations of extremely condensed matter, previously only accessible through theoretical calculations, can now be experimentally explored.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Synthesis of Anhydrous Acetates for the Components of Nuclear Fuel Recycling in Dialkylimidazolium Acetate Ionic Liquids [plus supplemental information]

A series of anhydrous acetate salts with uranium {[C 2 C 1 im][UO 2 (OAc) 3 ] (1), [C 2 C 2 im][UO 2 (OAc) 3 ] (2), and [C 4 C 1 im][UO 2 (OAc) 3 ] (3)}, lanthanides {[C 2 C 2 im] 2 [La(OAc) 5 ] (4) and [C 2 C 1 im] 2 [Nd(OAc) 5 ] (5)}, and strontium {[C 2 C 1 im] n [Sr- (OAc) 3 ] n (6)} (where C 2 C 1 im = 1-ethyl-3-methylimidazolium, C 2 C 2 im = 1,3-diethylimidazolium, C 4 C 1 im = 1-butyl-3-methylimi-dazolium, and OAc = acetate) have been prepared and structurally characterized. Both lanthanides and strontium are common components of the nuclear fuel waste, and their separation from uranium is an important but still challenging task. A new synthetic method with dialkylimidazolium acetate ionic liquids (ILs) as the solvent has been developed for the direct synthesis of homoleptic acetates from the corresponding hydrates and, unexpectedly, hardly soluble f-element oxides. Although the group of characterized compounds shows perfect structural variability, all actinide and lanthanide metal ions form monomeric complex anions where the metal cation coordinates to five ligands including two oxygen atoms in the case of uranium, as is commonly observed for uranyl compounds. Crystallographic analyses revealed that the complex [UO 2 (OAc) 3 ] - anions possess rather standard D 3h symmetry featuring a hexagonal-bipyramidal coordination environment, while the lanthanide anions [Ln(OAc) 5 ] 2- are fully asymmetric and the Ln 3+ cations are 10-coordinated in the form of a distorted bicapped tetragonal antiprism. This is the first report of lanthanide ions coordinated in this fashion. For Sr 2+ , 9-fold coordination through oxygen atoms in the form of a strongly distorted tricapped trigonal prism is observed. Furthermore, the crystallization of anhydrous, homoleptic, anionic acetate complexes from such a large variety of different metal salts appears to be due to the properties of dialkylimidazolium acetate ILs themselves, including enhanced basicity from the high concentration of free anions and their greater affinity for hydrogen-bonding solutes relative to metal cations.

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

Electrochemical Oxidation in Garnet-Type Solid Electrolyte by Formation of Point Defects

All-solid-state batteries hold greater promise for improving safety and energy density over conventional battery technology employing organic liquid electrolytes. One of the required features of a Li + conducting solid electrolyte is electrochemical stability, attained thermodynamically or kinetically, within the targeted operating voltage and temperature ranges. Therefore, understanding of the oxidative or reductive degradation mechanism is important to allow the design of stable solid electrolyte materials. This work contributes to building an understanding of the oxidative degradation mechanism in lithium solid electrolytes at cell operating conditions. Here, we have focused on resolving the oxidative decomposition mechanism of Al-doped lithium garnet Li 6.28 Al 0.24 La 3 Zr 2 O 12 (LLZO) as a state-of-the-art inorganic ceramic electrolyte. By combining experimental and computational analyses, we show that oxidation of LLZO occurs by simultaneous loss of oxygen and lithium from the structure, resulting in substoichiometric LLZO, at a moderate temperature (80 °C) and a high electrode potential (4.3 V vs Li/Li + ). Based on X-ray absorption and diffraction analyses, we find that the zirconium coordination shells in LLZO contract while the crystal structure experiences positive chemical strain upon electrochemical oxidation. The results from ex situ structural characterization of both the local structure and crystal symmetry are supported by a substoichiometric LLZO with lithium and oxygen vacancies, modeled by density functional theory (DFT) calculations. These chemical and structural changes in LLZO suppress effective lithium-ion conductivity by an order of magnitude. Formation of lithium and oxygen vacancies in LLZO upon electrochemical oxidation is different from prior thermodynamic predictions of phase decomposition of LLZO. The difference here is that the experiments were conducted at near-room temperature, which can hinder the kinetics of phase separation, and thus, the resultant LLZO solid electrolyte is still single-phase but substoichiometric in Li and O. In conclusion, these findings contribute an important degradation mechanism of the electrolyte, relevant for practical operational conditions of solid-state batteries.

36 MATERIALS SCIENCE↗