SEARCH · Engineering Papers
Results for “electronic structure methods”
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.
Spectroscopic characterization of Mn 2+ and Cd 2+ coordination to phosphorothioates in the conserved A9 metal site of the hammerhead ribozyme
Phosphorothioate modifications have widespread use in the field of nucleic acids. As substitution of sulfur for oxygen can alter metal coordination preferences, the phosphorothioate metal-rescue experiment is a powerful method for identifying metal coordination sites that influence specific properties in a large RNAs. The A9/G10.1 metal binding site of the hammerhead ribozyme (HHRz) has previously been shown to be functionally important through phosphorothioate rescue experiments. While an A9-S Rp substitution is inhibitory in Mg 2+ , thiophilic Cd 2+ rescues HHRz activity. Mn 2+ is also often used in phosphorothioate metal-rescue studies but does not support activity for the A9-S Rp HHRz. Here, we use EPR, electron spin-echo envelope modulation (ESEEM), and X-ray absorption spectroscopic methods to directly probe the structural consequences of Mn 2+ and Cd 2+ coordination to R p and S p phosphorothioate modifications at the A9/G10.1 site in the truncated hammerhead ribozyme (tHHRz). The results demonstrate that while Cd 2+ does indeed bind to S in the thio-substituted ligand, Mn 2+ coordinates to the non–sulfur oxo group of this phosphorothioate, regardless of isomer. Computational models demonstrate the energetic preference of Mn—O over Mn—S coordination in metal-dimethylthiophosphate models. In the case of the tHHRz, the resulting Mn 2+ coordination preference of oxygen in either R p or S p A9 phosphorothioates differentially tunes catalytic activity, with Mn—O coordination in the A9-S Rp phosphorothioate enzyme being inhibitory.
In situ nanoscale imaging of moiré superlattices in twisted van der Waals heterostructures
Abstract Direct visualization of nanometer-scale properties of moiré superlattices in van der Waals heterostructure devices is a critically needed diagnostic tool for study of the electronic and optical phenomena induced by the periodic variation of atomic structure in these complex systems. Conventional imaging methods are destructive and insensitive to the buried device geometries, preventing practical inspection. Here we report a versatile scanning probe microscopy employing infrared light for imaging moiré superlattices of twisted bilayers graphene encapsulated by hexagonal boron nitride. We map the pattern using the scattering dynamics of phonon polaritons launched in hexagonal boron nitride capping layers via its interaction with the buried moiré superlattices. We explore the origin of the double-line features imaged and show the mechanism of the underlying effective phase change of the phonon polariton reflectance at domain walls. The nano-imaging tool developed provides a non-destructive analytical approach to elucidate the complex physics of moiré engineered heterostructures.
Structure, folding and flexibility of co-transcriptional RNA origami
RNA origami is a method for designing RNA nanostructures that can self-assemble through co-transcriptional folding with applications in nanomedicine and synthetic biology. However, to advance the method further, an improved understanding of RNA structural properties and folding principles is required. Here we use cryogenic electron microscopy to study RNA origami sheets and bundles at sub-nanometre resolution revealing structural parameters of kissing-loop and crossover motifs, which are used to improve designs. In RNA bundle designs, we discover a kinetic folding trap that forms during folding and is only released after 10 h. Furthermore, exploration of the conformational landscape of several RNA designs reveal the flexibility of helices and structural motifs. Finally, sheets and bundles are combined to construct a multidomain satellite shape, which is characterized by individual-particle cryo-electron tomography to reveal the domain flexibility. Together, the study provides a structural basis for future improvements to the design cycle of genetically encoded RNA nanodevices.
The influence of lattice defects, recombination, and clustering on thermal transport in single crystal thorium dioxide
Thermal transport is a key controlling property for thorium dioxide in many applications where defect-generating radiation fields are present. An understanding of the effect of nanoscale lattice defects on thermal transport in this material is currently unavailable due to a lack of high-quality, single crystal material from which unit processes may be investigated. In this work, a series of thorium dioxide single crystals are exposed to 2 MeV proton irradiation at room temperature and 600°C to create microscale regions with varying densities and types of point and extended defects. Defected regions are investigated using spatial domain thermoreflectance to quantify the change in thermal conductivity as a function of ion fluence as well as transmission electron microscopy and Raman spectroscopy to interrogate the structure of the generated defects. Together, this com- bination of methods provides important initial insight into defect formation, recombination, and clustering in thorium dioxide and the effect of those defects on thermal transport. These methods also provide a promising pathway for the quantification of the smallest-scale defects which cannot be captured using traditional microscopy techniques and play and outsized role in degrading thermal performance.
Kinetic modeling of hot tail runaway electron generation during plasma disruptions using the JOREK code
The generation of runaway electrons (REs) during disruptions poses a significant challenge for the operation of tokamaks. The production of these high-energy electrons can cause substantial damage, particularly when the plasma current is high, making it a critical concern for ITER. For the high-temperature plasmas anticipated in ITER, the primary generation of REs may be dominated by the hot tail mechanism, which consists of the acceleration of hot electrons from the pre-disruption population which have not yet thermalized with the bulk following the rapid cooling of the plasma. To account for the significant 3D effects on RE production, a hot tail modeling framework has been developed within the non-linear 3D extended MHD code JOREK. This paper presents the structure of this framework, which is based on test electrons evolving in MHD fields. The verification of the method shows good agreement with the reference DREAM code for 0D test cases, as well as for axisymmetric simulations of 15 MA ITER H-mode disruption scenarios. Furthermore, a proof-of-principle application to a DIII-D case demonstrates the framework’s capability to capture for the first time the hot tail generation in 3D MHD simulations in realistic geometry. Preliminary results suggest that the production of REs is significantly reduced by stochastic losses.
Machine learning and TDDFT software for stopping power computation
(SF-24-012) Stopping power describes the rate that a material slows radiation particles passing through it and is useful in designing many technologies. Few organizations can perform new measurements, which require significant resources and rare equipment, and all others rely on coarse approximations rendered from pre-existing data. Methods for computing stopping power in new materials, such as time-dependent density functional theory (TD-DFT), have only recently (circa-2015) become available but are too computationally costly to use frequently enough to have a pronounced impact. We have created a method that opens a pathway to computing stopping power without any need for experimental data by combining electronic structure computations and machine learning.
Fragment-based determination of a proteinase K structure from MicroED data using ARCIMBOLDO_SHREDDER
Structure determination of novel biological macromolecules by X-ray crystallography can be facilitated by the use of small structural fragments, some of only a few residues in length, as effective search models for molecular replacement to overcome the phase problem. Independence from the need for a complete pre-existing model with sequence similarity to the crystallized molecule is the primary appeal of ARCIMBOLDO , a suite of programs which employs this ab initio algorithm for phase determination. Here, the use of ARCIMBOLDO is investigated to overcome the phase problem with the electron cryomicroscopy (cryoEM) method known as microcrystal electron diffraction (MicroED). The results support the use of the ARCIMBOLDO_SHREDDER pipeline to provide phasing solutions for a structure of proteinase K from 1.6 Å resolution data using model fragments derived from the structures of proteins sharing a sequence identity of as low as 20%. ARCIMBOLDO_SHREDDER identified the most accurate polyalanine fragments from a set of distantly related sequence homologues. Alternatively, such templates were extracted in spherical volumes and given internal degrees of freedom to refine towards the target structure. Both modes relied on the rotation function in Phaser to identify or refine fragment models and its translation function to place them. Model completion from the placed fragments proceeded through phase combination of partial solutions and/or density modification and main-chain autotracing using SHELXE . The combined set of fragments was sufficient to arrive at a solution that resembled that determined by conventional molecular replacement using the known target structure as a search model. This approach obviates the need for a single, complete and highly accurate search model when phasing MicroED data, and permits the evaluation of large fragment libraries for this purpose.
Role of shell composition and morphology in achieving single-emitter photostability for green-emitting “giant” quantum dots
The use of the varied chemical reactivity of precursors to drive the production of a desired nanocrystal architecture has become a common method to grow thick-shell graded alloy quantum dots (QDs) with robust optical properties. Conclusions on their behavior assume the ideal chemical gradation and uniform particle composition. In this work, advanced analytical electron microscopy (high-resolution scanning transmission electron microscopy coupled with energy dispersive spectroscopy) is used to confirm the nature and extent of compositional gradation and these data are compared with performance behavior obtained from single-nanocrystal spectroscopy to elucidate structure, chemical-composition, and optical-property correlations. In particular, the evolution of the chemical structure and single-nanocrystal luminescence was determined for a time-series of graded-alloy “CdZnSSe/ZnS” core/shell QDs prepared in a single-pot reaction. In a separate step, thick (~6 monolayers) to giant (>14 monolayers) shells of ZnS were added to the alloyed QDs via a successive ionic layer adsorption and reaction (SILAR) process, and the impact of this shell on the optical performance was also assessed. By determining the degree of alloying for each component element on a per-particle basis, we observe that the actual product from the single-pot reaction is less “graded” in Cd and more so in Se than anticipated, with Se extending throughout the structure. The latter suggests much slower Se reaction kinetics than expected or an ability of Se to diffuse away from the initially nucleated core. It was also found that the subsequent growth of thick phase-pure ZnS shells by the SILAR method was required to significantly reduce blinking and photobleaching. However, correlated single-nanocrystal optical characterization and electron microscopy further revealed that these beneficial properties are only achieved if the thick ZnS shell is complete and without large lattice discontinuities. Thus, we identify the necessary structural design features that are required for ideal light emission properties in these green-visible emitting QDs.
De novo synthesis and near atomic resolution imaging of host immune receptors critical for pathogen recognition (Abbreviated Final Report)
The innate immune system serves as the body’s first line of defense against invading pathogens, responding rapidly through the deployment of immune cells at common sites of infection, such as the skin and airways. These immune-cell sentinels express a range of highly conserved receptors, including Toll-like receptors (TLRs), which recognize and bind to pathogen-derived components. This recognition event triggers intracellular signaling cascades that initiate and coordinate immune responses. Despite their significance, the complete structural characterization of full-length TLRs, including their extracellular domain (ECD), transmembrane domain (TMD), and Toll/interleukin-1 receptor (TIR) domain, remains incomplete. In this study, we examined the expression and isolation of human TLR4 incorporated into nanodiscs using two approaches: a cell-free synthesis system and a cell-based transfection strategy employing Expi293F cells derived from the human embryonic kidney lineage. Our findings demonstrate that NLP-bound human TLR4 produced via the cell-based method yielded functional protein suitable for time-resolved single-particle cryo-electron microscopy (cryo-EM). This advancement enables structural characterization of full-length TLR4, maintaining the integrity of its extracellular domain, transmembrane domain, and Toll/Interleukin-1 receptor (TIR) domain.
Ground State Properties of the Wide Band Gap Semiconductor Beryllium Sulfide (BeS)
We report the results from self-consistent calculations of electronic, transport, and bulk properties of beryllium sulfide (BeS) in the zinc-blende phase, and employed an ab-initio local density approximation (LDA) potential and the linear combination of atomic orbitals (LCAO). We obtained the ground state properties of zb-BeS with the Bagayoko, Zhao, and Williams (BZW) computational method, as enhanced by Ekuma and Franklin (BZW-EF). Our findings include the electronic energy bands, the total (DOS) and partial (pDOS) densities of states, electron and hole effective masses, the equilibrium lattice constant, and the bulk modulus. The calculated band structure clearly shows that zb-BeS has an indirect energy band gap of 5.436 eV, from Γ to a point between Γ and X, for an experimental lattice constant of 4.863 Å. This is in excellent agreement with the experiment, unlike the findings of more than 15 previous density functional theory (DFT) calculations that did not perform the generalized minimization of the energy functional, required by the second DFT theorem, which is inherent to the implementation of our BZW-EF method.
LoTToR: An Algorithm for Missing-Wedge Correction of the Low-Tilt Tomographic 3D Reconstruction of a Single-Molecule Structure
A single-molecule three-dimensional (3D) structure is essential for understanding the thermal vibrations and dynamics as well as the conformational changes during the chemical reaction of macromolecules. Individual-particle electron tomography (IPET) is an approach for obtaining a snap-shot 3D structure of an individual macromolecule particle by aligning the tilt series of electron tomographic (ET) images of a targeted particle through a focused iterative 3D reconstruction method. The method can reduce the influence on the 3D reconstruction from large-scale image distortion and deformation. Due to the mechanical tilt limitation, 3D reconstruction often contains missing-wedge artifacts, presented as elongation and an anisotropic resolution. Here, we report a post-processing method to correct the missing-wedge artifact. This low-tilt tomographic reconstruction (LoTToR) method contains a model-free iteration process under a set of constraints in real and reciprocal spaces. A proof of concept is conducted by using the LoTToR on a phantom, i.e., a simulated 3D reconstruction from a low-tilt series of images, including that within a tilt range of ±15°. The method is validated by using both negative-staining (NS) and cryo-electron tomography (cryo-ET) experimental data. A significantly reduced missing-wedge artifact verifies the capability of LoTToR, suggesting a new tool to support the future study of macromolecular dynamics, fluctuation and chemical activity from the viewpoint of single-molecule 3D structure determination.
Fabrication, thermal analysis, and heavy ion irradiation resistance of epoxy matrix nanocomposites loaded with silane-functionalized ceria nanoparticles
This paper describes a detailed understanding of how nanofillers function as radiation barriers within the polymer matrix, and how their effectiveness is impacted by factors such as composition, size, loading, surface chemistry, and dispersion. Here, we designed a comprehensive investigation of heavy ion irradiation resistance in epoxy matrix composites loaded with surface-modified ceria nanofillers, utilizing tandem computational and experimental methods to elucidate radiolytic damage processes and relate them to chemical and structural changes observed through thermal analysis, vibrational spectroscopy, and electron microscopy. A detailed mechanistic examination supported by FTIR spectroscopy data identified the bisphenol A moiety as a primary target for degradation reactions. Results of computational modeling by the Stopping Range of Ions in Matter (SRIM) Monte Carlo simulation were in good agreement with damage analysis from surface and cross-sectional SEM imaging. All metrics indicated that ceria nanofillers reduce the damage area in polymer nanocomposites, and that nanofiller loading and homogeneity of dispersion are key to effective damage prevention. The results of this study represent a significant pathway for engineered irradiation tolerance in a diverse array of polymer nanocomposite materials. Numerous areas of materials science can benefit from utilizing this facile and effective method to extend the reliability of polymer materials.
Kβ X-ray Emission Spectra Analysis Using Bayesian Optimization
The Kβ X-ray emission spectrum of 3 d transition metals is rich with electronic and structural information due to strong exchange interactions with the valence shell of the metal, and has become crucial for understanding their spin and oxidation states. The spectrum is commonly treated using crystal-field multiplet theory, a semi-empirical theory that uses tunable parameters to control the strength of the effects present in X-ray emission spectroscopy (XES). However, determining the experimental values of these parameters remains a challenge. We present a methodology that applies Bayesian optimization to crystal-field multiplet theory to determine parameter values. The algorithm is tested on the X-ray emission spectra of a collection of Mn, Co, and Ni oxides. We are able to find optimal values for the four most impactful parameters: Slater−Condon reduction factors F dd , F pd , and G pd , and crystal field splitting 10 Dq . The algorithm produces significantly improved accuracy compared to current analysis methods, and probes interparameter dependencies by modeling the error landscape. This advancement enhances XES analysis by offering an approach of obtaining quantitative electronic structural information on 3 d transition metal valence shells, facilitating applications across various scientific fields.
Experimental Quantification of Spin–Phonon Coupling in Molecular Qubits Using Inelastic Neutron Scattering
Electronic spin superposition states enable nanoscale sensing through their sensitivity to the local environment, yet their sensitivity to vibrational motion also limits their coherence times. In molecular spin systems, chemical tunability and atomicscale resolution are accompanied by a dense, thermally accessible phonon spectrum that introduces efficient spin relaxation pathways. Despite extensive theoretical work, there is little experimental consensus on which vibrational energies dominate spin relaxation or how molecular structure controls spin−phonon coupling (SPC). We present a fully experimental method to quantify SPC coefficients by combining temperature-dependent vibrational spectra from inelastic neutron scattering with spin relaxation rates measured by electron paramagnetic resonance. We apply this framework to two model S = 1/2 systems, copper(II) phthalocyanine (CuPc) and copper(II) octaethylporphyrin (CuOEP). Two distinct relaxation regimes emerge: below 40 K, weakly coupled lattice modes below 50 cm −1 dominate, whereas above 40 K, optical phonons above ∼185 cm −1 become thermally populated and drive relaxation with SPC coefficients nearly 3 orders of magnitude larger. Structural distortions in CuOEP that break planar symmetry soften the crystal lattice and enhance anharmonic scattering but also raise the energy of stretching modes at the molecular core where the spins reside. This redistributes vibrational energy toward the molecular periphery and out of plane, ultimately reducing SPC relative to CuPc and enabling room-temperature spin coherence in CuOEP. Although our method does not provide mode-specific SPC coefficients, it quantifies contributions from distinct spectral regions and establishes a broadly applicable, fully experimental link between crystal structure, lattice dynamics, and spin relaxation.
Bridging the time scale in exascale computing of chemical systems (Final Technical Report)
This report summarizes the work carried out with support of the United States Department of Energy under Award DE-SC0019441. The theme of this project was to develop and apply methods that allowed for the acceleration of atomistic calculations, particularly in challenging areas such as multiphase systems, electrified interfaces, uncertainty estimation, and applications requiring chemical accuracy, which tend to be applications where simulation time is severely bottlenecked by the computational time requirements. Much of the focus was on the application of emerging machine-learning methodologies, although a wide range of methodologies were employed. This report has two major sections. The first focuses on the methodological advances themselves. Within this part, we report a number of major advances, a few examples of which are described here. We report the first machine-learning scheme for the acceleration of electronically grand-canonical calculations (that is, those applicable to electrochemistry). We report new methods of performing transfer learning, in which physics-based priors can be used to provide predictions, often with uncertainty estimates, of images well outside of training sets; we also offer ways to fine-tune these transfer-learning models. We provide a new systematic means to generate and apply minimal training data sets to very large (10,000’s of atoms) systems, with only small training sets appropriate for electronic structure. We developed new methodologies to integrate surface vibrations into surface adsorption calculations. We made advances to the applicability of diffusion Monte Carlo methods to allow (learned) force prediction, finite-size error correction, and force-free means of searching for transition states. We integrated machine-learned atomistic predictions into mechanism generation codes. Additionally, we released new software including AmpTorch, a modernized version of our original atomistic machine-learning code Amp. The second part of this report focuses on the scientific applications that accompanied, and were often enabled by, the methodological advances described earlier. A few examples follow, but full details are in the individual chapters of the report. For example, we developed a general theory of phonon-induced friction on molecular adsorbates. We showed fundamentally how solvent influences the adsorption and desorption process and how it differs from the processes typically involved at the solid–gas interface, making aqueous-phase and electrocatalysis different from traditional thermocatalysis. We examined how metal–insulator and magnetic transitions can be probed, and accelerated exciton dynamics via Frenkel Hamiltonian parameters. We showed that the nearsighted force-training approach, developed within this project, can predict both the stability and reactivity of large nanoparticles, and can also lead to insights on catalyst coverage on binding energies and entropies. These applied studies, which generally integrated with our method development, allowed us to push forward the theoretical understanding of several reaction classes.
Discovering invariant spatial features in electron energy loss spectroscopy images on the mesoscopic and atomic levels
Over the last two decades, Electron Energy Loss Spectroscopy (EELS) imaging with a scanning transmission electron microscope has emerged as a technique of choice for visualizing complex chemical, electronic, plasmonic, and phononic phenomena in complex materials and structures. The availability of the EELS data necessitates the development of methods to analyze multidimensional data sets with complex spatial and energy structures. Traditionally, the analysis of these data sets has been based on analysis of individual spectra, one at a time, whereas the spatial structure and correlations between individual spatial pixels containing the relevant information of the physics of underpinning processes have generally been ignored and analyzed only via the visualization as 2D maps. Here, we develop a machine learning-based approach and workflows for the analysis of spatial structures in 3D EELS data sets using a combination of dimensionality reduction and multichannel rotationally invariant variational autoencoders. This approach is illustrated for the analysis of both the plasmonic phenomena in a system of nanowires and in the core excitations in functional oxides using low loss and core-loss EELS, respectively. The code developed in this manuscript is open sourced and freely available and provided as a Jupyter notebook for the interested reader.
Multi–length scale characterization of point defects in thermally oxidized, proton irradiated iron oxides
A key for the success of safe nuclear power generation system is to consider structural materials that are economical, meet mechanical property needs, possess good corrosion resistance, and are radiation tolerant. Nevertheless, fundamental insights that elucidate the details of radiation damage on materials corrosion performance are lacking. This includes the behavior of surface oxides which often regulate corrosion. For example, it is unclear how non-equilibrium point defects, oxide structure, mass transport in oxides, and subsequent oxidation behavior are altered by the radiation. Here, in this work, some of the effects of proton irradiation on the attributes of point defects, iron oxide microstructures, and the physical nature of the oxidation product were correlated with corrosion behavior. Iron oxides, fabricated by thermal oxidation in air at 400°C and 800°C for 1 h, were subjected to 200 keV, 0.03 dpa (displacements per atom) of proton irradiation, and subjected to corrosion reactivity assessment using AC and DC electrochemical methods. Experimental methods that target materials information at different length scales, such as positron annihilation spectroscopy (atomistic), transmission electron microscopy (mesoscopic), and electrochemical methods (macroscopic), were coupled to shed light on the impact of radiation-induced defect modifications and structural changes in oxides on corrosion reactivity which ultimately affects durability in harsh environments.