Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Crystal structure prediction”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Determination of Vibrational Motions Driving Photoinduced Electron Transfer Reactions in Molecular Crystals and Organic Thin Films

The objective of this research is to develop an understanding of how nuclear motions can be used to drive and control photo-induced charge transfer and singlet fission reactions in molecular crystals and thin films. By following the structural evolution of reacting molecules and nanostructures on the timescale of their nuclear motion, we will determine what interplay of nuclear coordinates is most efficient in driving electron transfers, a necessary first step in most photovoltaic and photocatalytic processes. We will identify which nuclear coordinates effectively promote desirable electron transfer processes, or ineffectively dissipate photon energy into unwanted thermal energy, thus guiding the design and implementation of molecular systems tailored for high efficiency processes. Specifically, we aim to identify which nuclear coordinates are responsible for driving electron transfer and singlet fission reactions at different points along the reaction coordinate, particularly in regimes which deviate from Marcus theory predictions. We will do so by investigating systems with electron transfer and singlet fission occurring from non-thermalized states due to strong donor-acceptor electronic coupling. First, we will prove that specific vibrational modes are responsible for charge transfer processes in molecular crystals. Secondly, we will investigate how vibrational coherences and specific nuclear motions can be used to drive high yield singlet fission processes. Finally, we will determine how complex film morphology affects long-range charge transport, a key limiting factor in our ability to make gains in efficiency.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

De novo design of drug-binding proteins with predictable binding energy and specificity

The de novo design of small molecule–binding proteins has seen exciting recent progress; however, high-affinity binding and tunable specificity typically require laborious screening and optimization after computational design. We developed a computational procedure to design a protein that recognizes a common pharmacophore in a series of poly(ADP-ribose) polymerase–1 inhibitors. One of three designed proteins bound different inhibitors with affinities ranging from <5 nM to low micromolar. X-ray crystal structures confirmed the accuracy of the designed protein-drug interactions. Molecular dynamics simulations informed the role of water in binding. Binding free energy calculations performed directly on the designed models were in excellent agreement with the experimentally measured affinities. We conclude that de novo design of high-affinity small molecule–binding proteins with tuned interaction energies is feasible entirely from computation.

Science & Technology - Other Topics↗
DOE OSTI2020

Modeling the role of local crystallographic correlations in microstructures of Ti-6Al-4V using a correlated structure visco-plastic self-consistent polycrystal plasticity formulation

This paper presents a multi-level crystal plasticity-based simulation framework for modeling mechanical response and microstructure evolution of Ti-6Al-4V with α-lath/lamellar microstructures. Here, the model is a correlated structure visco-plastic self-consistent (CS-VPSC) formulation linking three scales: a single crystals micro-scale, a lath/lamellar colony meso-scale, and a lath/lamellar aggregate macro-scale. A selected hardening law for the evolution of critical resolved shear stress per slip system used in CS-VPSC is phenomenological. However, it adjusts the resistances of basal and prismatic slip systems based on the geometry of slip transfer between adjacent lamellae. Consistent with experimental evidences, the resolved shear stress on the pyramidal slip planes is dependent not only on the stress in the direction of slip but also on the two orthogonal shear stress components and the three normal stress components (non-Schmid effects). Electron backscatter diffraction (EBSD) data in conjunction with a procedure relying on α→β phase transformation is used to construct paired variants of α-laths/lamellae satisfying their local crystallographic correlations. The procedure fits volume fractions of individual laths/lamellae with the experimental EBSD data and selects a distribution of habit planes between adjacent variants with respect to the loading direction. The simulation framework is applied to interpret the deformation behavior in tension and compression along two sample directions of Ti-6Al-4V fabricated via laser powder bed fusion. Moreover, the model is used to simulate texture evolution during rolling of the material to large strains. It is demonstrated that the model is capable of predicting plastic anisotropy/asymmetry and the concomitant texture evolution. While the model reveals a significant effect of habit plane inclination with respect to the loading direction on yield stress, the comparison of the data and model predictions shows that a random distribution of habit planes fits the flow response. It is further inferred that the tension-compression asymmetry arises from the non-Schmid effects.

36 MATERIALS SCIENCE↗
DOE OSTI2022

High Sensitivity of Non–Fullerene Organic Solar Cells Morphology and Performance to a Processing Additive

Although solvent additives are used to optimize device performance in many novel non-fullerene acceptor (NFA) organic solar cells (OSCs), the effect of processing additives on OSC structures and functionalities can be difficult to predict. Here, two polymer-NFA OSCs with highly sensitive device performance and morphology to the most prevalent solvent additive chloronaphthalene (CN) are presented. Devices with 1% CN additive are found to nearly double device efficiencies to 10%. However, additive concentrations even slightly above optimum significantly hinder device performance due to formation of undesirable morphologies. A comprehensive analysis of device nanostructure shows that CN is critical to increasing crystallinity and optimizing phase separation up to the optimal concentration for suppressing charge recombination and maximizing performance. Here, domain purity and crystallinity are highly correlated with photocurrent and fill factors. However, this effect is in competition with uncontrolled crystallization of NFAs that occur at CN concentrations slightly above optimal. This study highlights how slight variations of solvent additives can impart detrimental effects to morphology and device performance of NFA OSCs. Furthermore, successful scale-up processing of NFA-based OSCs will require extreme formulation control, a tuned NFA structure that resists runaway crystallization, or alternative methods such as additive-free fabrication.

14 SOLAR ENERGY↗
DOE OSTI2020

A Lanthanum‐Filled Carbon–Boron Clathrate

Abstract We report a carbon‐boron clathrate with composition 2 La@B 6 C 6 (LaB 3 C 3 ). Like recently reported SrB 3 C 3 , [1] single‐crystal X‐ray diffraction and computational modelling indicate that the isostructural La member crystallizes in the cubic bipartite sodalite structure (Type‐VII clathrate) with La atoms encapsulated within truncated octahedral cages composed of alternating carbon and boron atoms. The covalent nature of the B−C bonding results in a hard, incompressible framework, and owing to the balanced electron count, La 3+ [B 3 C 3 ] 3− exhibits markedly improved pressure stability and is a semiconductor with an indirect band gap predicted near 1.3 eV. A variety of different guest atoms may potentially be substituted within Type‐VII clathrate cages, presenting opportunities for a large family of boron‐stabilized, carbon‐based clathrates with ranging physical properties.

Strobel, Timothy A.↗
DOE OSTI2022

Numerical Study of Mineral Nucleation and Growth on a Substrate

We report nucleation and growth processes of minerals and other crystals can significantly affect one another due to the transport limitations and local depletion of reactive ions in the solution. Most numerical models and experimental measurements have typically focused on either growth or nucleation, but not both. In this work, we incorporate a heterogeneous nucleation process based on classical nucleation theory into a microcontinuum model that implements the Darcy–Brinkman–Stokes approach to study the interplay between nucleation and crystal growth on a substrate in diffusive systems. We demonstrate how the Damköhler number (reaction rate) and nucleation rate prefactor change the effective nucleation rate on a substrate. Higher surface growth rates deplete the solute concentration around the nuclei that appear initially on the substrate, creating islands that screen against further nucleation. The model predicts that measured nucleation rates may be affected by the history of crystal nucleation on the substrate. In the extreme case of high growth rates relative to diffusion, it predicts that the rate of subsequent nucleation is limited by reactant depletion. We introduce a nondimensional number α to represent the relation between surface propagation rate during growth and the heterogeneous nucleation rate. We show that it is important to control Damköhler number and α to achieve similar precipitation regimes at different reaction and nucleation rates. We suggest that the observed universality can guide the interpretation of experimental results on nucleation rates, since matching experiment can be achieved by tuning transport, reaction, and nucleation parameters simultaneously. In addition, we show how the bulk solution concentration affects the structure and topology of precipitation on a substrate.

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

DFT-guided flux synthesis of a family of layered titanates crystallizing in the lepidocrocite structure type

K 0.72 Fe 0.72 Ti 1.28 O 4 , a layered mixed-metal oxide belonging to the lepidocrocite structure type, was grown out of a molten KCl-KF mixture at 800 °C. The compound crystallizes in orthorhombic space group Cmcm and exhibits two-dimensional layers, composed of mixed Fe/Ti octahedra, separated by K ions that occupy the interlayer layer space and maintain charge balance. The targeted flux synthesis of the Ni, Cu, and Zn analogs was carried out after density functional theory (DFT) calculations predicted their formation to be energetically favorable, yielding K 0.77 Ni 0.38 Ti 1.62 O 4 , K 0.77 Cu 0.38 Ti 1.62 O 4 , and K 0.73 Zn 0.36 Ti 1.64 O 4 . Finally, all four compositions can be prepared via conventional solid-state technique by mixing and heating stoichiometric amounts of KNO 3 , appropriate transition metal oxide precursors, and TiO 2 .

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

High-Pressure Equation of State of 1,3,5-triamino-2,4,6-trinitrobenzene: Insights into the Monoclinic Phase Transition, Hydrogen Bonding, and Anharmonicity

The high-pressure equation of state (EOS) of energetic materials (EMs) is important for continuum and mesoscale models of detonation performance and initiation safety. Additionally, obtaining a high-fidelity EOS of the insensitive EM 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) has proven to be difficult because of challenges in experimental characterization at high pressures (HPs). In this work, powder X-ray diffraction patterns were fitted using the recently discovered monoclinic I2/a phase above 4 GPa, which shows that TATB is less compressible than when indexed with the triclinic $P\bar{1}$ phase. First-principles calculations were performed with Perdew–Burke–Ernzerhof (PBE) and PBE0 functionals including thermal effects using the $P\bar{1}$ phase. PBE0 improves the description of hydrogen bonding and thus predicts accurate planar a and b lattice parameters under ambient conditions. However, discrepancies in the predicted lattice parameters above 4–10 GPa compared with experimental measurements indexed with $P\bar{1}$ are further evidence of a structural modification at high pressure. Layer sliding defects are formed during molecular dynamics simulations, which induces an anharmonic effect on the thermal expansion of the c lattice parameter. In short, the results provide several insights into determining high-fidelity EOS parameters for TATB and other molecular crystals.

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

Machine learning of 27Al NMR electric field gradient tensors for crystalline structures from DFT

NMR crystallography has emerged as a promising technique for the determination and refinement of atomic coordinates in crystal structures. The crystal structure of compounds containing quadrupolar nuclei, such as 27Al, can be improved by directly comparing solid-state NMR measurements to DFT computations of the electric field gradient (EFG) tensor. The non-negligible computational cost of these first-principles calculations limits the applicability of this method to all but the most well-defined structures. We developed a fast, low-cost machine learning model to predict EFG parameters based on local structural motifs and elemental parameters. We computed 8081 EFG tensors from 1681 27Al crystalline solids using DFT and benchmarked them against 105 experimentally measured 27Al sites. Surprisingly, simple local geometric features dominate the predictive performance of the resulting random-forest model, yielding an R2 value of 0.98 and an RMSE of 0.61 MHz for CQ, the quadrupolar coupling constant. This model accuracy should enable pre-refining future structural assignments before finally validating with first-principles calculations. Such a catalogue of 27Al NMR tensors can serve as a tool for researchers assigning complex NMR spectra influenced by the nuclear electric quadrupole interaction.

Sun, He↗
DOE OSTI2021

Monte carlo simulations and experiments of all-silica zeolite LTA assembly combining structure directing agents that match cage sizes

In this work, we investigated the influence of organic structure-directing agents (OSDAs) on the formation rates of all-silica zeolite LTA using both simulations and experiments, to shed light on the crystallization process. We compared syntheses using one OSDA with a diameter close to the size of the large cavity in LTA, and two OSDAs of diameters matching the sizes of both the small and large LTA cavities. Reaction-ensemble Monte Carlo (RxMC) simulations predict a speed up of LTA formation using two OSDAs matching the LTA pore sizes; this qualitative result is confirmed by experimental studies of crystallization kinetics, which find a speedup in all-silica LTA crystallization of a factor of 3. Analyses of simulated rings and their Si–O–Si angular energies during RxMC crystallizations show that all ring sizes in the faster crystallization exhibit lower angular energies, on average, than in the slower crystallization, explaining the origin of the speedup through packing effects.

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

Molecular underpinnings of ssDNA specificity by Rep HUH-endonucleases and implications for HUH-tag multiplexing and engineering

Abstract Replication initiator proteins (Reps) from the HUH-endonuclease superfamily process specific single-stranded DNA (ssDNA) sequences to initiate rolling circle/hairpin replication in viruses, such as crop ravaging geminiviruses and human disease causing parvoviruses. In biotechnology contexts, Reps are the basis for HUH-tag bioconjugation and a critical adeno-associated virus genome integration tool. We solved the first co-crystal structures of Reps complexed to ssDNA, revealing a key motif for conferring sequence specificity and for anchoring a bent DNA architecture. In combination, we developed a deep sequencing cleavage assay, termed HUH-seq, to interrogate subtleties in Rep specificity and demonstrate how differences can be exploited for multiplexed HUH-tagging. Together, our insights allowed engineering of only four amino acids in a Rep chimera to predictably alter sequence specificity. These results have important implications for modulating viral infections, developing Rep-based genomic integration tools, and enabling massively parallel HUH-tag barcoding and bioconjugation applications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Helicity-dependent photocurrents in the chiral Weyl semimetal RhSi

Weyl semimetals are crystals in which electron bands cross at isolated points in momentum space. Associated with each crossing point (or Weyl node) is a topological invariant known as the Berry monopole charge. The circular photogalvanic effect (CPGE), whereby circular polarized light generates a helicity-dependent photocurrent, is a notable example of a macroscopic property that emerges directly from the topology of the Weyl semimetal band structure. Recently, it was predicted that the amplitude of the CPGE associated with optical transitions near a Weyl node is proportional to its monopole charge. In chiral Weyl systems, nodes of opposite charge are nondegenerate, opening a window of wavelengths where the CPGE resulting from uncompensated Berry charge can emerge. Here, we report measurements of CPGE in the chiral Weyl semimetal RhSi, revealing a CPGE response in an energy window that closes at 0.65 eV, in agreement with the predictions of density functional theory.

36 MATERIALS SCIENCE↗
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 OSTI2021

Dynamic Barriers to Crystallization of Calcium Barium Carbonates

Metastable carbonates play important roles in geochemistry, biomineralization and serve as model systems for nonclassical theories of nucleation and growth. Balcite (Ca 0.5 Ba 0.5 CO 3 ) is a remarkable carbonate phase that is isostructural with a high-temperature modification of calcite (CaCO 3 ), yet can be synthesized at ambient conditions. Here, we investigate crystallization pathways in the Ba-Ca-CO 3 -H 2 O system, with a focus on the transformation of amorphous calcium barium carbonate (ACBC) to balcite and subsequent decomposition into the equilibrium calcite (CaCO3) and witherite (BaCO 3 ) phases. Density functional theory calculations show that balcite is an unstable solid solution (Ca 1-x BaxCO 3 , R 3¯ m) in the range 0.17 < x < 0.5, but is accessible through the amorphous ACBC precursor for x ≲ 0.5, and predict its decomposition into calcite and witherite. We confirm this pathway experimentally but found demixing to proceed slowly and remain incomplete even after 9 months. Nucleation kinetics of balcite from ACBC was assessed using a microfluidic assay, where increasing barium content led to a surprising increase in the balcite nucleation rate, despite decreasing thermodynamic driving force. We attribute crystallization rates that dramatically accelerate with time to changes in interfacial structure and composition during coarsening of the amorphous precipitate. By carefully quantifying the thermodynamic and kinetic contributions in the multistep crystallization of a metastable carbonate, we produce insights that allow us to better interpret the formation and persistence of metastable minerals in natural and synthetic environments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Lattice Structure and Dynamics of Sparse Molecular Crystals: OsO 4 and RuO 4

OsO 4 and RuO 4 are molecular oxides with unique tetrameric structures and rare +8 oxidation states. Accurately modeling their properties remains challenging for density functional theory (DFT) due to weak intertetramer interactions, which standard functionals fail to capture. Here, in this work, we show that the van der Waals (vdW)-corrected density functional (vdW-DF-optB86b) provides structural parameters that are much closer to experimental values than the standard generalized gradient approximation, with volume predictions that fall within the experimentally observed range. Phonon band structure analysis shows that the inclusion of vdW interactions stabilizes soft phonon modes, highlighting the importance of dispersion corrections for accurate predictions of lattice dynamics. Experimental measurements of the phonon density of states for OsO 4 , obtained via inelastic neutron scattering, demonstrate good agreement with our vdW-DF-optB86b calculations. These results validate OsO 4 and RuO 4 as valuable benchmarks for structural and vibrational calculations via vdW-corrected DFT methods and offer insights for studying the broader class of sparse molecular materials.

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

Modeling the Environment-Dependent Kinetics of Oxygen Reduction Reaction – a Continuum Model for Electric Double Layer

Here, for proton-exchange-membrane fuel cells (PEMFCs) to achieve broad commercialization, improved energy-conversion efficiency with minimal Pt-based electrocatalyst is required. Because the sluggish rate of oxygen reduction reaction (ORR) limits the efficiency of PEMFCs, the efficiency improvement requires a better understanding of ORR kinetics and mechanism to design better catalyst. To understand the ORR mechanism, theoretical and experimental analyses have been conducted. While previous studies reasonably explained the catalyst-dependent activity on single crystal catalysts in 0.1 M perchloric acid solution, the explicit effect of electrolyte and related microenvironments is not thoroughly understood. The change in the electrolyte alters the electric-double-layer (EDL) structure and thus the local microenvironment at the electrode/electrolyte interface. Thus, the structure of the EDL should be carefully analyzed to uncover the electrolyte-dependent reaction kinetics. In this talk, we propose a multiscale continuum model to predict the EDL structure and examine the effect of perchloric acid concentration on ORR activity on Pt (111). The model includes Density Potential Functional Theory (DPFT) for electron density and Modified Poisson Boltzmann equation for species’ density and electric potential. Also, the interaction between adsorbents and electric field is taken into account by minimizing the grand potential. After model validation with experimentally measured double-layer capacity data as a function of applied potential and concentration, the effect of the perchloric acid concentration (0.02 M – 0.2 M) on ORR activity is analyzed and discussed. It is shown that the model reproduces the specific activity obtained in the experiments when assuming the oxygen adsorption is limiting the rate, which can be attributed to the large energetic barrier for solvent reorganization. Then, extension of the model to PEMFC ionomer electrolytes will be introduced. Overall, the model framework and findings provide insights into the ORR mechanism and guidance on how to tailor catalyst materials for increased PEMFC performance.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2022

Competition of magnetocrystalline anisotropy of uranium layers and zigzag chains in UNi 0.34 Ge 2 single crystals

Structural and thermodynamic properties of single-crystalline UNi 1-x Ge 2 with x = 0.66 have been investigated by measuring magnetization, specific heat, and thermal expansion over a wide range of temperatures and magnetic fields. The measurements revealed the emergence of a long-range antiferromagnetic ordering of uranium magnetic moments below the Néel temperature T N = 45.5 (1) K and the existence of two easy axes in the studied compound, namely b and c, which correspond to the planes of the uranium zig-zag chains. A magnetic field applied along these two crystallographic directions induces in the system a first-order metamagnetic phase transition (from antiferromagnetic to field-polarized paramagnetism), and the width of the magnetic hysteresis associated with that transition reaches as much as 40 kOe at the lowest temperatures. A magnetic phase diagram developed from the experimental data showed that the metastable region associated with that magnetic hysteresis forms a funnel that narrows toward the N´eel point in a zero magnetic field. The four-layer Ising model has successfully predicted the colinear antiferromagnetic structure in UNi 0.34 Ge 2 (known from earlier reports), its magnetic phase diagram, and temperature and field variations of its magnetization. Moreover, it suggests that the first-order phase transition extends down to zero magnetic field, although it is barely detectable in the experiments performed in low magnetic fields. According to this model, the second-order phase transition occurs in the compound only in zero field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗