SEARCH · Search NASA
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.
Abstract not provided
Finding an optimal match between two different crystal structures underpins many important materials science problems, including describing solid-solid phase transitions and developing models for interface and grain boundary structures. Here, we formulate the matching of crystals as an optimization problem where the goal is to find the alignment and the atom-to-atom map that minimize a given cost function such as the Euclidean distance between the atoms. We construct an algorithm that directly solves this problem for large finite portions of the crystals and retrieves the periodicity of the match subsequently. We demonstrate its capacity to describe transformation pathways between known polymorphs and to reproduce experimentally realized structures of semi-coherent interfaces. Additionally, from our findings, we define a rigorous metric for measuring distances between crystal structures that can be used to properly quantify their geometric (Euclidean) closeness.
Thermoelectric materials enable the energy conversion of waste heat into electricity, helpful to relieve global energy crisis. Here, we report a systematic investigation on high-temperature thermoelectric gadolinium selenides, cubic Gd 3-x Se 4 (x = 0.16, 0.21 and 0.25) and orthorhombic Gd 2 Se 3-y (y = 0.02, 0.06 and 0.08). High energy synchrotron x-ray diffraction and total scattering have been used to investigate the crystallographic and local structures. Atomic-scale clusters of Gd vacancy in the cubic phase are observed by employing the reverse Monte Carlo simulation. For cubic Gd 3-x Se 4 , adjusting Gd vacancy triggers the effect of multiple conduction bands, confirmed by the increase in effective masses and theoretical calculations. A reasonable peak zT value of 0.27 is achieved at 850 K for Gd 3-x Se 4 (x = 0.16). On the other hand, tuning Se vacancy enables the optimization of electron concentration for the orthorhombic Gd 2 Se 3-y . More significantly, its low deformation potential (Ξ = 12eV) gives rise to enhanced electron mobility and higher peak zT value of 0.54 at 850 K for Gd 2 Se 3-y (y = 0.02). Intriguingly, a higher zT of 1.2 at 1200 K is reasonably predicted by quality factor analysis. Finally, this work extends the scope of high-temperature thermoelectric materials and facilitates the exploration of novel high-temperature thermoelectric materials
Abstract Bacteriophage T4 is a classic model system for studying the mechanisms of DNA processing. A key protein in T4 DNA processing is the gp32 single-stranded DNA-binding protein. gp32 has two key functions: it binds cooperatively to single-stranded DNA (ssDNA) to protect it from nucleases and remove regions of secondary structure, and it recruits proteins to initiate DNA processes including replication and repair. Dda is a T4 helicase recruited by gp32, and we purified and crystallized a gp32–Dda–ssDNA complex. The low-resolution structure revealed how the C-terminus of gp32 engages Dda. Analytical ultracentrifugation analyses were consistent with the crystal structure. An optimal Dda binding peptide from the gp32 C-terminus was identified using surface plasmon resonance. The crystal structure of the Dda–peptide complex was consistent with the corresponding interaction in the gp32–Dda–ssDNA structure. A Dda-dependent DNA unwinding assay supported the structural conclusions and confirmed that the bound gp32 sequesters the ssDNA generated by Dda. The structure of the gp32–Dda–ssDNA complex, together with the known structure of the gp32 body, reveals the entire ssDNA binding surface of gp32. gp32–Dda–ssDNA complexes in the crystal are connected by the N-terminal region of one gp32 binding to an adjacent gp32, and this provides key insights into this interaction.
We present a materials generation framework that couples a symmetry-conditioned variational autoencoder with a differentiable SO(3) power spectrum objective to steer candidates toward a specified local environment under the crystallographic constraints. In particular, we implement a fully differentiable pipeline that performs batch-wise optimization on both direct and latent crystallographic representations. Using the GPU acceleration, the implementation achieves about fivefold speed compared to our previous CPU workflow, while yielding comparable outcomes. In addition, we introduce the optimization strategy that alternatively performs optimization on the direct and latent crystal representations. This dual-level relaxation approach can effectively escape local minima defined by different objective gradients, thus increasing the success rate of generating complex structures satisfying the target local environments. This framework can be extended to systems consisting of multi-components and multi-environments, providing a scalable route to generate material structures with the target local environment.
Atomic structures and adatom geometries of surfaces encode information about the thermodynamics and kinetics of the processes that lead to their formation, and which can be captured by a generative physical model. In this work, we develop a workflow based on a machine-learning-based analysis of scanning tunneling microscopy images to reconstruct the atomic and adatom positions, and a Bayesian optimization procedure to minimize statistical distance between the chosen physical models and experimental observations. We optimize the parameters of a 2- and 3-parameter Ising model describing surface ordering and use the derived generative model to make predictions across the parameter space. For concentration dependence, we compare the predicted morphologies at different adatom concentrations with the dissimilar regions on the sample surfaces that serendipitously had different adatom concentrations. The proposed workflow can be used to reconstruct the thermodynamic models and associated uncertainties from the experimental observations of materials microstructures. The code used in the manuscript is available at https://github.com/saimani5/Adatom_interactions.
Packing motifs—patterns in how molecules orient relative to one another in a crystal structure—are an important concept in many subdisciplines of materials science because of correlations observed between specific packing motifs and properties of interest. That said, packing motif data sets have remained small and noisy due to intensive manual labeling processes and insufficient labeling schemes. The most prominent labeling algorithms calculate relative interplanar angles of nearest neighbor molecules to determine the packing motif of a molecular crystal, but this simple approach can fail when neighbors are naively sampled isotropically around the crystal structure. To remedy this issue, here we propose an optimization algorithm, which rotates the molecular crystal structure to find representative molecules that inform the packing motif. We package this algorithm into an automated framework—Autopack—which both optimally rotates the crystal structure and labels the packing motif based on the appropriate neighboring molecules. In this work, we detail the Autopack framework and its performance, which shows improvements compared to previous state-of-the-art labeling methods, providing the first quantitative point of comparison for packing motif labeling algorithms. Furthermore, using Autopack (available at https://ipo.llnl.gov/technologies/software/autopack), we perform the first large-scale study of potential relationships between chemicals’ compositions and packing motifs, which shows that these relationships are more complex than previously hypothesized from studies that used only tens of polycyclic aromatic hydrocarbon molecules. Autopack’s capabilities help pose next steps for crystal engineering research focusing not only on a molecule’s adoption of a specific packing motif but also on new structure–property relationships.
Topology optimization (TO) is commonly applied to design the unit cells of periodic structures. For example, metamaterials, lattice structures, phononic crystals (PhC), and photonic crystals (PC) have all been previously designed via TO. Unfortunately, the optimal structures for certain design objectives, e.g., bandgaps, are often impossible to manufacture as they have disconnected regions or “islands” of solid material (ISM) that are not self-supporting. Additionally, designs with enclosed void space (EVS) are problematic for additive manufacturing (AM) since support material or pre-sintered powder cannot be removed after manufacturing. We present a series of constraints that may be incorporated into any TO framework to ensure structures are self-supporting without enclosed voids. Additionally, we employ homogenization-based constraints that allow the designer to tune the elastic stiffness and isotropy of the optimized design. The proposed constraints are evaluated on example microstructures and utilized in a simple optimization test problem to highlight their abilities and limitations so that guidelines for appropriate combinations of constraints may be proposed. Effective constraint combinations are demonstrated on the design of 3D photonic crystals for maximum bandgap subject to manufacturing and stiffness constraints.
We present a workflow that traces the path from the bulk structure of a crystalline material to assessing its performance in carbon capture from coal's postcombustion flue gases. This workflow is applied to a database of 324 covalent-organic frameworks (COFs) reported in the literature, to characterize their CO 2 adsorption properties using the following steps: (1) optimization of the crystal structure (atomic positions and unit cell) using density functional theory, (2) fitting atomic point charges based on the electron density, (3) characterizing the pore geometry of the structures before and after optimization, (4) computing carbon dioxide and nitrogen isotherms using grand canonical Monte Carlo simulations with an empirical interaction potential, and finally, (5) assessing the CO 2 parasitic energy via process modeling. The full workflow has been encoded in the Automated Interactive Infrastructure and Database for Computational Science (AiiDA). Both the workflow and the automatically generated provenance graph of our calculations are made available on the Materials Cloud, allowing peers to inspect every input parameter and result along the workflow, download structures and files at intermediate stages, and start their research right from where this work has left off. In particular, our set of CURATED (Clean, Uniform, and Refined with Automatic Tracking from Experimental Database) COFs, having optimized geometry and high-quality DFT-derived point charges, are available for further investigations of gas adsorption properties. We plan to update the database as new COFs are being reported.
Explore the source record for details and available documents.
Here, we developed a PyTorch-based architecture called HydraGNN that implements graph convolutional neural networks (GCNNs) to predict the formation energy and the bulk modulus for models of solid solution alloys for various atomic crystal structures and relaxed volumes. We trained the GCNN surrogate model on a dataset for nickel–niobium (NiNb) generated by the embedded atom model (EAM) empirical interatomic potential for demonstration purposes. The dataset was generated by calculating the formation energy and the bulk modulus as a prototypical elastic property for optimized geometries starting from initial body-centered cubic (BCC), face-centered cubic (FCC), and hexagonal compact packed (HCP) crystal structures, with configurations spanning the possible compositional range for each of the three types of initial crystal structures. Numerical results show that the GCNN model effectively predicts both the formation energy and the bulk modulus as function of the optimized crystal structure, relaxed volume, and configurational entropy of the model structures for solid solution alloys.
Metal organic frameworks (MOFs) are porous crystalline materials that display a wide variety of physical and chemical properties. Their single crystal structure determination is often challenging because in most cases micro- or nano-sized crystals spontaneously form upon MOF synthesis, which cannot be recrystallized. The production of larger single crystals for structure determination involves optimizing, and thus modifying, the conditions of synthesis, in which success cannot be guaranteed. Failure to produce crystals suitable for single-crystal X-ray diffraction leaves the 3D structure of the MOF compound unknown, and scientists must resort to more challenging structure solution methods based on X-ray powder or electron diffraction data. These laborious tasks can be avoided by using serial crystallography techniques which merge data collected on many micro-crystals. Here, we report the application of three synchrotron serial crystallography methods. We call these “mesh”, “grid” and “mesh&collect” scans. “Still” images (no rotation) are collected in the mesh scan approach, whereas small rotational wedges are collected in the grid scan method. The third protocol, mesh&collect, combines the acquisition of still images and rotational wedges. Using these means, we determine the ab initio structure of benchmark MOFs, MIL-100(Fe) and ZIF-8, that differ largely in unit cell size. These methods are expected to be widely applicable and facilitate structure determination of many MOF microcrystalline systems.
Nitrogen is a primary constituent of Earth’s atmosphere and an essential component of many biological processes. Recent developments in understanding N incorporation into silicate minerals suggest that the solid earth and lithosphere could provide a substantial reservoir for the planetary N budget. Furthermore, determining the concentration and distribution of low abundances of N in silicates presents a significant challenge in microanalysis. We tested capabilities for high spatial resolution quantification of low levels of N using EPMA with the goal of developing standards to be used in coupled N and noble gas analyses of silicate minerals. We investigated the concentration and distribution of N in a suite of ammonium-bearing silicate minerals that reflect the dominant N-bearing phases of the lithosphere: tobelite, buddingtonite, ammonioleucite, tsaregorodtsevite, and hyalophane. These minerals form during NH4-rich alteration of K-bearing silicate minerals, and the bonding environment of N in each mineral is likely to differ based on crystal structure. Using an optimized analytical routine for N quantification enables sensitivity approaching 0.02 weight% N and high-resolution quantitative mapping of sub-micron variations in N contents.
Chimeric antigen receptor (CAR) T cells have emerged as a promising class of therapeutic agents, generating remarkable responses in the clinic for a subset of human cancers. One major challenge precluding the wider implementation of CAR therapy is the paucity of tumor-specific antigens. Here, we describe the development of a CAR targeting the tumor-specific isocitrate dehydrogenase 2 (IDH2) with R140Q mutation presented on the cell surface in complex with a common human leukocyte antigen allele, HLA-B*07:02. Engineering of the hinge domain of the CAR, as well as crystal structure-guided optimization of the IDH2 R140Q- HLA-B*07:02-targeting moiety, enhances the sensitivity and specificity of CARs to enable targeting of this HLA-restricted neoantigen. This approach thus holds promise for the development and optimization of immunotherapies specific to other cancer driver mutations that are difficult to target by conventional means.
Nitrogen is a primary constituent of Earth’s atmosphere and an essential component of many biological processes. Recent developments in understanding N incorporation into silicate minerals suggest that the solid earth and lithosphere could provide a substantial reservoir for the planetary N budget. Furthermore, determining the concentration and distribution of low abundances of N in silicates presents a significant challenge in microanalysis. We tested capabilities for high spatial resolution quantification of low levels of N using EPMA with the goal of developing standards to be used in coupled N and noble gas analyses of silicate minerals. We investigated the concentration and distribution of N in a suite of ammonium-bearing silicate minerals that reflect the dominant N-bearing phases of the lithosphere: tobelite, buddingtonite, ammonioleucite, tsaregorodtsevite, and hyalophane. These minerals form during NH 4 -rich alteration of K-bearing silicate minerals, and the bonding environment of N in each mineral is likely to differ based on crystal structure. Using an optimized analytical routine for N quantification enables sensitivity approaching 0.02 weight% N and high-resolution quantitative mapping of sub-micron variations in N contents.
Abstract Measuring the similarity between two arbitrary crystal structures is a common challenge in crystallography and materials science. Although there are an infinite number of ways to mathematically relate two crystal structures, only a few are physically meaningful. Here we introduce both a geometry-based and a symmetry-adapted similarity metric to compare crystal structures. Using crystal symmetry and combinatorial optimization we describe an algorithm to arrive at the structural relationship that minimizes these similarity metrics across all possible maps between any pair of crystal structures. The approach makes it possible to (i) identify pairs of crystal structures that are identical, (ii) quantitatively measure the similarity between crystal structures, and (iii) find and rank structural transformation pathways between any pair of crystal structures. We discuss the advantages of using the symmetry-adapted cost metric over the geometric cost. Finally, we show that all known structural transformation pathways between common crystal structures are recovered with the mapping algorithm. The methodology presented in this study will be of value to efforts that seek to catalogue crystal structures, identify structural transformation pathways or prune large first-principles datasets used to parameterize on-lattice Hamiltonians.
Not Available
This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called • BCC_opt • FCC_opt • HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string “A_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters “A_B_Câ€, three files are provided: • A_B_C_opt.xyz: The optimized geometries in xyz format • A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. • A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. • A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.