Search NASA⌕ Search

SEARCH · Search NASA

Results for “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 37 records · Page 2

The XtalOpt Evolutionary Algorithm for Crystal Structure Prediction

Significant progress has been made in the field of a priori crystal structure prediction, with a number of recent remarkable success stories. Herein, we briefly outline the methods that have been developed for finding the global minimum structure and interesting local minima without the need for experimental information. Focus is placed on describing the XTALOPT evolutionary algorithm (EA) developed in our group towards this end. XTALOPT is published under well-known open-source licenses, and the EA searches can be analyzed via the AVOGADRO chemical editor and visualizer. We describe new algorithmic developments that have made it possible to predict the structures of ever-more complex crystalline lattices. Benchmark tests, which clearly illustrate how the new developments improve the success rate and accelerate the discovery of the global minimum structure, are carried out. Lastly, we describe how XTALOPT has been employed to predict novel ternary hydrides that have the propensity for high-temperature superconductivity under pressure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

How Accurate Can Crystal Structure Predictions Be for High-Energy Molecular Crystals?

Molecular crystals have shallow potential energy landscapes, with multiple local minima separated by very small differences in total energy. Predicting molecular packing and molecular conformation in the crystal generally requires ab initio methods of high accuracy, especially when polymorphs are involved. We used dispersion-corrected density functional theory (DFT-D) to assess the capabilities of an evolutionary algorithm (EA) for the crystal structure prediction (CSP) of well-known but challenging high-energy molecular crystals (HMX, RDX, CL-20, and FOX-7). While providing the EA with the experimental conformation of the molecule quickly re-discovers the experimental packing, it is more realistic to start instead from a naïve, flat, or neutral initial conformation, which reflects the limited experimental knowledge we generally have in the computational design of molecular crystals. By doing so, and using fully flexible molecules in fully variable unit cells, we show that the experimental structures can be predicted in fewer than 20 generations. Nonetheless, one must be aware that some molecular crystals have naturally hindered evolutions, requiring as many attempts as there are space groups of interest to predict their structures, and some may require the accuracy of all-electron calculations to discriminate between closely ranked structures. To save resources in this computationally demanding process, we showed that a hybrid xTB/DFT-D approach could be considered in a subsequent study to push the limits of CSP beyond 200+ atoms and for cocrystals.

42 ENGINEERING↗

Validation of finite element and boundary element methods for predicting structural vibration and radiated noise

Analytical and experimental validation of methods to predict structural vibration and radiated noise are presented. A rectangular box excited by a mechanical shaker was used as a vibrating structure. Combined finite element method (FEM) and boundary element method (BEM) models of the apparatus were used to predict the noise radiated from the box. The FEM was used to predict the vibration, and the surface vibration was used as input to the BEM to predict the sound intensity and sound power. Vibration predicted by the FEM model was validated by experimental modal analysis. Noise predicted by the BEM was validated by sound intensity measurements. Three types of results are presented for the total radiated sound power: (1) sound power predicted by the BEM modeling using vibration data measured on the surface of the box; (2) sound power predicted by the FEM/BEM model; and (3) sound power measured by a sound intensity scan. The sound power predicted from the BEM model using measured vibration data yields an excellent prediction of radiated noise. The sound power predicted by the combined FEM/BEM model also gives a good prediction of radiated noise except for a shift of the natural frequencies that are due to limitations in the FEM model.

Seybert, A. F.↗

Predicting structures of cross-linked condensation polymers

Mathematical procedure is used to predict structure of cross-linked condensation polymer differentiated from an additional polymer resulting from specific reaction. Procedure will greatly reduce amount of empirical formulation and testing needed to produce desired product.

Marsh, H. E.↗

Experimental validation of finite element and boundary element methods for predicting structural vibration and radiated noise

This research report is presented in three parts. In the first part, acoustical analyses were performed on modes of vibration of the housing of a transmission of a gear test rig developed by NASA. The modes of vibration of the transmission housing were measured using experimental modal analysis. The boundary element method (BEM) was used to calculate the sound pressure and sound intensity on the surface of the housing and the radiation efficiency of each mode. The radiation efficiency of each of the transmission housing modes was then compared to theoretical results for a finite baffled plate. In the second part, analytical and experimental validation of methods to predict structural vibration and radiated noise are presented. A rectangular box excited by a mechanical shaker was used as a vibrating structure. Combined finite element method (FEM) and boundary element method (BEM) models of the apparatus were used to predict the noise level radiated from the box. The FEM was used to predict the vibration, while the BEM was used to predict the sound intensity and total radiated sound power using surface vibration as the input data. Vibration predicted by the FEM model was validated by experimental modal analysis; noise predicted by the BEM was validated by measurements of sound intensity. Three types of results are presented for the total radiated sound power: sound power predicted by the BEM model using vibration data measured on the surface of the box; sound power predicted by the FEM/BEM model; and sound power measured by an acoustic intensity scan. In the third part, the structure used in part two was modified. A rib was attached to the top plate of the structure. The FEM and BEM were then used to predict structural vibration and radiated noise respectively. The predicted vibration and radiated noise were then validated through experimentation.

Seybert, A. F.↗

Quadrupolar NMR crystallography guided crystal structure prediction (QNMRX-CSP) of zwitterionic organic HCl salts

In this work, we benchmark quadrupolar NMR crystallography guided crystal structure prediction (QNMRX-CSP) for determining the crystal structures of two zwitterionic organic HCl salts, L-ornithine HCl ( Orn ) and L-histidine HCl·H 2 O ( Hist ). These salts present an interesting challenge for QNMRX-CSP, as gas-phase geometry optimizations used to generate starting structures for the organic zwitterionic fragments fail to capture their correct solid-state geometries. To overcome this limitation, geometry optimizations using the COSMO water-solvation model are employed to generate initial structural models. Using this approach, QNMRX-CSP yields structural models of the two zwitterionic organic HCl salts that closely match experimentally determined crystal structures. In addition, the application of QNMRX-CSP to Hist represents a further step toward the de novo structural determination of solvated organic HCl salts, as Hist is the first benchmark system of this type to include a water molecule as a component of its crystal structure. This work is significant for its potential application to the structural determination of active pharmaceutical ingredients, which often feature complex organic components and solvated solid forms.

Fleischer, Carl H. [Florida State Univ., Tallahass↗

Structure prediction of epitaxial inorganic interfaces by lattice and surface matching with Ogre

We present a new version of the Ogre open source Python package with the capability to perform structure prediction of epitaxial inorganic interfaces by lattice and surface matching. In the lattice matching step, a scan over combinations of substrate and film Miller indices is performed to identify the domain-matched interfaces with the lowest mismatch. Subsequently, surface matching is conducted by Bayesian optimization to find the optimal interfacial distance and in-plane registry between the substrate and the film. For the objective function, a geometric score function is proposed based on the overlap and empty space between atomic spheres at the interface. The score function reproduces the results of density functional theory (DFT) at a fraction of the computational cost. The optimized interfaces are pre-ranked using a score function based on the similarity of the atomic environment at the interface to the bulk environment. Final ranking of the top candidate structures is performed with DFT. Ogre streamlines DFT calculations of interface energies and electronic properties by automating the construction of interface models. The application of Ogre is demonstrated for two interfaces of interest for quantum computing and spintronics, Al/InAs and Fe/InSb.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Decoy selection for protein structure prediction via extreme gradient boosting and ranking

Background: Identifying one or more biologically-active/native decoys from millions of non-native decoys is one of the major challenges in computational structural biology. The extreme lack of balance in positive and negative samples (native and non-native decoys) in a decoy set makes the problem even more complicated. Consensus methods show varied success in handling the challenge of decoy selection despite some issues associated with clustering large decoy sets and decoy sets that do not show much structural similarity. Recent investigations into energy landscape-based decoy selection approaches show promises. However, lack of generalization over varied test cases remains a bottleneck for these methods. Results: We propose a novel decoy selection method, ML-Select, a machine learning framework that exploits the energy landscape associated with the structure space probed through a template-free decoy generation. The proposed method outperforms both clustering and energy ranking-based methods, all the while consistently offering better performance on varied test-cases. Moreover, ML-Select shows promising results even for the decoy sets consisting of mostly low-quality decoys. Conclusions: ML-Select is a useful method for decoy selection. This work suggests further research in finding more effective ways to adopt machine learning frameworks in achieving robust performance for decoy selection in template-free protein structure prediction.

59 BASIC BIOLOGICAL SCIENCES↗

CHARGED POINT DEFECTS AND DEFECT COMPLEXES IN PHOSPHORENE AND STRUCTURAL PREDICTION OF TAU-11, A DENSITY FUNCTIONAL THEORY APPROACH

With the advent of powerful computers, first principle calculations have become a robust and reliable tool to complement experiments in the study of materials. First principle calculations such as density functional theory can be used to guide experiments to save time and effort when it comes to new and advanced materials discovery. Study of charge defects in materials and crystal structure prediction are such areas in which density functional theory calculations facilitate the prediction of defect formation energy, charge transition levels and elemental site occupations. Experimentally, studying charge defects requires inferring results from various techniques. Likewise, determining site occupancies for similar elements using experimental techniques is challenging and sometimes impossible. In this work, we study the charged intrinsic and extrinsic defects and defect complexes in phosphorene and determine the site occupancies of Al, Fe and Si in the τ11 phase using density functional theory. In the defect work, we calculate the defect formation energy, charge transition level, binding energies and Stokes shift for vacancy, dopant substitution and dopant-vacancy defect complexes in phosphorene. We found that vacancy defect in phosphorene becomes negatively charged in n-type doping and may passivate the dopants and reduce carrier concentration and mobility. For non-metal dopants in phosphorene, we predict that O, S and Si prefer to form dopant-vacancy complexes removing the vacancy defect states from the band gap. Mn dopant-vacancy defect complex exhibits possibility of switching between two magnetic spin states. Lastly, using density functional theory, we complement the results from neutron powder diffraction to determine the site occupancies of Al and Si in τ11 phase which was used to determine the alloying element to develop a high temperature and low density alloy.

36 MATERIALS SCIENCE↗

On the numerical sensitivity of cellular automata grain structure predictions to large thermal gradients and cooling rates

Cellular automata (CA) models of as-solidified grain structure, originally developed and applied to casting, have become a common means of predicting grain structure resulting from Additive Manufacturing (AM) processes. The majority of these models are based on the decentered octahedron approach, which attempts to correct for the effect of grid anisotropy on the prediction of competitive solidification of dendritic grains. However, AM solidification occurs under cooling rates ($\dot{T}$) and thermal gradients (G) that are orders of magnitude larger than those encountered in casting, and no systematic investigation on the effect of the CA model cell size (Δx) and time step (Δt) on AM microstructure predictions has been performed. Here, in this study, such an investigation is first performed via simulation of individual grains of various crystallographic orientations with a fixed, unidirectional G, showing that CA prediction of the steady-state undercooling matched the expected values based on the interfacial response function at small G and deviated from the expected values at large G. Simulation of competitive growth of multiple grains showed a weakening of the predicted texture as G and Δx became large. Simulation of solidification under AM conditions, where G and $\dot{T}$ vary spatially across the melt pools, showed that not only does grain selection weaken and deviate from expectations at large Δx, but grains with crystallographic $\langle$100$\rangle$ aligned with the grid directions are more adversely affected by the temperature field discontinuities than grains with other crystallographic orientations. Despite the fact that the exact grain competition results depended on Δt, the overall texture development was notably less sensitive to Δt than Δx, provided that a reasonable value of Δt is selected based on the ratio of Δx to the maximum local solidification velocity in the simulation domain. Finally, from the directional solidification and AM simulation results, an analysis of computational cost compared to simulation resolution is performed based on an equation derived to quantify the relatively inaccuracy in grain selection based on the model and temperature field inputs. From this analysis, it is concluded that there is a need for algorithmic improvements to improve CA grain competition accuracy for large G processing conditions as sufficiently small Δx to resolve the necessary competition is intractable for many AM processing conditions.

36 MATERIALS SCIENCE↗

Elucidation of Bottom-up Growth of CaCO3 Involving Prenucleation Clusters from Structure Predictions and Decomposition of Globally Optimized (CaCO3)n Nanoclusters

Low-energy minima structures for (CaCO3)n, n = 28, are predicted using bottom-up genetic algorithms in conjunction with density functional theory electronic structure calculations, in comparison with the frozen and relaxed top-down clusters generated by cuts from the calcite, vaterite, and aragonite crystal structures. Similarities in structural motifs for the bottom-up and relaxed top-down are revealed using a fragment recognition technique. Fragment energy decomposition analysis shows that the bottom-up and relaxed top-down clusters belong to two classes of amorphous clusters with distinct intracluster energy distributions, despite their structural similarity. The bottom-up clusters with >20 formula units are surface stabilized with negative surface energy densities. In contrast, the top-down clusters are interior stabilized with positive surface energy densities. We prove that the sign of the surface energy density determines whether the nucleation reaction energy as a function of nuclear size has a maximum or a minimum. The surface-stabilized bottom-up clusters are proposed to be a type of prenucleation cluster at the minimum of the nucleation reaction energy. A mechanism for mineralization of CaCO3 involving prenucleation clusters and nonclassical growth pathway is proposed on the basis of our theoretical findings, which is consistent with previous titration experiments.

calcium carbonate, prenucleation cluster, biominer↗

Sensitivity of cellular automata grain structure predictions for high solidification rates

Understanding as-solidified grain structure development during constrained (G 0) alloy solidification conditions is a necessary component of understanding the larger process-microstructure relationships during additive manufacturing (AM) processing. This study applies a cellular automata (CA) model, for static solidification velocities and thermal gradients on the order of those expected in AM melt pools, to investigate grain shape development across ranges of input parameters governing temperature field evolution, nucleation, and growth. Variation in nucleated grain shape from columnar to equiaxed as a function of thermal gradient and nucleation density is modeled with CA and verified against the established transition model of Hunt and Gaumann. Furthermore, while thermal conditions and nucleation density are shown to have the largest roles on microstructure development, the sensitivity of grain structure was very unevenly distributed across parameter space. In particular, it was found that microstructure is very sensitive to mean nucleation undercooling and the interfacial response function in process space regimes near grain structure transitions, but entirely controlled by a combination of thermal gradient, solidification velocity, and nucleation density away from such transitions. The understanding of grain shape sensitivity obtained through CA modeling, particularly at large thermal gradient and nucleation density, will be necessary for accurate microstructure prediction during the extreme processing conditions and large input parameter uncertainty common to AM.

36 MATERIALS SCIENCE↗

ExaCA grain structure predictions for laser powder bed fusion processing

This dataset contains simulated cross-sections of laser powder bed fusion grain structures produced using the microstructure model ExaCA, in turn using time-temperature history data produced using the heat transport model AdditiveFOAM. These predictions show the grain structure using various permutations of hatch spacing and nucleation density. Also contained in the dataset are the input files necessary to reproduce the results. AdditiveFOAM (https://github.com/ORNL/AdditiveFOAM) and ExaCA (https://github.com/LLNL/ExaCA), both open source software, are required to reproduce the results in this dataset using the given input files. This DOI was updated 12/09/2025.

36 MATERIALS SCIENCE↗

Crystal structure prediction with host-guided inpainting generation and foundation potentials

Unconditional crystal structure generation with diffusion models faces challenges in identifying symmetric crystals as the unit cell size increases. Here, we present the crystal host-guided generation (CHGGen) framework to address this challenge through conditional generation using an inpainting method, which optimizes a fraction of atomic positions within a predefined and symmetrized host structure to improve the success rate for symmetric structure generation. By integrating inpainting structure generation with a foundation potential for structure optimization, we demonstrate the method on the ZnS–P 2 S 5 and Li–Si chemical systems, where the inpainting method generates a higher fraction of symmetric structures than unconditional generation. The practical significance of CHGGen extends to enabling the structural modification of crystal structures, particularly for systems with partial occupancy or intercalation chemistry. The inpainting method also allows for seamless integration with other generative models, providing a versatile framework for accelerating materials discovery.

Zhong, Peichen [University of California, Berkeley↗

Structure Prediction of Ionic Epitaxial Interfaces with Ogre Demonstrated for Colloidal Heterostructures of Lead Halide Perovskites

Colloidal epitaxial heterostructures are nanoparticles composed of two different materials connected at an interface, which can exhibit properties different from those of their individual components. Combining dissimilar materials offers exciting opportunities to create a wide variety of functional heterostructures. However, assessing structural compatibility–the main prerequisite for epitaxial growth–is challenging when pairing complex materials with different lattice parameters and crystal structures. This complicates both the selection of target heterostructures for synthesis and the assignment of interface models when new heterostructures are obtained. Here, we demonstrate Ogre as a powerful tool to accelerate the design and characterization of colloidal heterostructures. To this end, we implemented developments tailored for the high-efficiency prediction of epitaxial interfaces between ionic/polar materials, which encompass most colloidal semiconductors. These include the use of pre-screening candidate models based on charge balance at the interface and the use of a classical potential for fast energy evaluations, with parameters automatically calculated based on the input bulk structures. These developments are validated for perovskite-based CsPbBr 3 /Pb 4 S 3 Br 2 heterostructures, where Ogre produces interface models in excellent agreement with density functional theory and experiments. Furthermore, we use Ogre to rationalize the templating effect of CsPbCl 3 on the growth of lead sulfochlorides, where perovskite seeds induce the formation of Pb 4 S 3 Cl 2 rather than Pb 3 S 2 Cl 2 due to better epitaxial compatibility. Finally, combining Ogre simulations with experimental data enables us to unravel the structure and composition of the hitherto unsolved CsPbBr 3 /Bi x Pb y S z interface, and to assign a structure to several other reported metal halide- and oxide-based interfaces. The Ogre package is available on GitHub or via the OgreInterface desktop application, available for Windows, Linux, and Mac.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗