Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “electronic structure method”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Electronic Structure and Mechanical Properties of Solvated Montmorillonite Clay Using Large-Scale DFT Method

Montmorillonite clay (MMT) has been widely used in engineering and environmental applications as a landfill barrier and toxic waste repository due to its unique property as an expandable clay mineral that can absorb water easily. This absorption process rendered MMT to be highly exothermic due to electrostatic interactions among molecules and hydrogen bonds between surface atoms. A detailed study of a large supercell model of structural clay enables us to predict long-term nuclear waste storage. Herein, a large solvent MMT model with 4071 atoms is studied using ab initio density functional theory. The DFT calculation and analysis clarify the important issues, such as bond strength, solvation effect, elasticity, and seismic wave velocities. These results are compared to our previous study on crystalline MMT (dry). The solvated MMT has reduced shear modulus (G), bulk modulus (K), and Young’s modulus (E). We observe that the conduction band (CB) in the density of states (DOS) of solvated MMT model has a single, conspicuous peak at –8.5 eV. Moreover, the atom-resolved partial density of states (PDOS) summarizes the roles played by each atom in the DOS. These findings illuminate numerous potential sophisticated applications of MMT clay.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Unveiling the Redox Noninnocence of Metallocorroles: Exploring K-Edge X-ray Absorption Near-Edge Spectroscopy with a Multiconfigurational Wave Function Approach

X-ray absorption near-edge spectroscopy (XANES) is an advanced technique for probing the local electronic structure of catalysts, effectively identifying the noninnocent nature of ligands in transition-metal complexes. Metallocorroles with noninnocent corrole rings exhibit unusual electronic structures that challenge traditional density functional theory (DFT) methods, necessitating more rigorous approaches to describe electron correlation accurately. We explored K-edge XANES spectra of Fe, Mn, and Co metallocorroles using TDDFT and wave function-based methods. This is the first investigation employing multireference methods, specifically RASSCF, RASPT2, and MC-PDFT, to analyze the redox noninnocent nature of metallocorroles reflected in their XANES spectra. We quantified the noninnocent character of the corrole and the oxidation states of the metals, capturing more than singly excited excitations responsible for the pre-edge peak. Our findings demonstrate the importance of these advanced computational techniques for accurately predicting XANES spectra, providing a reliable understanding of the electronic properties of such complexes. In conclusion, this study offers a new strategy for investigating ligand redox noninnocence via integrated experimental and computational XANES.

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

Force-Free Identification of Minimum-Energy Pathways and Transition States for Stochastic Electronic Structure Theories

Here, the accurate mapping of potential energy surfaces (PESs) is crucial to our understanding of the numerous physical and chemical processes mediated by atomic rearrangements, such as conformational changes and chemical reactions, and the thermodynamic and kinetic feasibility of these processes. Stochastic electronic structure theories, e.g., Quantum Monte Carlo (QMC) methods, enable highly accurate total energy calculations that in principle can be used to construct the PES. However, their stochastic nature poses a challenge to the computation and use of forces and Hessians, which are typically required in algorithms for minimum-energy pathway (MEP) and transition state (TS) identification, such as the nudged elastic band (NEB) algorithm and its climbing image formulation. Here, we present strategies that utilize the surrogate Hessian line-search method, previously developed for QMC structural optimization, to efficiently identify MEP and TS structures without requiring force calculations at the level of the stochastic electronic structure theory. By modifying the surrogate Hessian algorithm to operate in path-orthogonal subspaces and at saddle points, we show that it is possible to identify MEPs and TSs by using a force-free QMC approach. We demonstrate these strategies via two examples, the inversion of the ammonia (NH 3 ) molecule and the nucleophilic substitution (S N 2) reaction F – + CH 3 F → FCH 3 + F – . We validate our results using Density Functional Theory (DFT)- and Coupled Cluster (CCSD, CCSD(T))-based NEB calculations. We then introduce a hybrid DFT-QMC approach to compute thermodynamic and kinetic quantities, free energy differences, rate constants, and equilibrium constants that incorporates stochastically optimized structures and their energies, and show that this scheme improves upon DFT accuracy. Our methods generalize straightforwardly to other systems and other high-accuracy theories that similarly face challenges computing energy gradients, paving the way for highly accurate PES mapping, transition state determination, and thermodynamic and kinetic calculations at significantly reduced computational expense.

Iyer, Gopal R.↗
DOE OSTI2025

How the choice of exchange–correlation functional affects DFT-based simulations of the hydrated electron

Hydrated electrons are anionic species that are formed when an excess electron is introduced into liquid water. Building an understanding of how hydrated electrons behave in solution has been a long-standing effort of simulation methods, of which density functional theory (DFT) has come to the fore in recent years. The ability of DFT to model the reactive chemistry of hydrated electrons is an attractive advantage over semi-classical methodologies; however, relatively few density functional approximations (DFAs) have been used for the hydrated electron simulations presented in the literature. Here, we simulate hydrated electron systems using a series of exchange–correlation (XC) functionals spanning Jacob’s ladder. We calculate a variety of experimental and other observables of the hydrated electron and compare the XC functional dependence for each quantity. We find that the formation of a stable localized hydrated electron is not necessarily limited to hybrid XC functionals and that some hybrid functionals produce delocalized hydrated electrons or electrons that react with the surrounding water at an unphysically fast rate. Here, we further characterize how different DFAs impact the solvent structure and predicted spectroscopy of the hydrated electron, considering several methods for calculating the hydrated electron’s absorption spectrum for the best comparison between structures generated using different density functionals. None of the dozen or so DFAs that we investigated are able to correctly predict the hydrated electron’s spectroscopy, vertical detachment energy, or molar solvation volume.

Ab-initio molecular dynamics↗
DOE OSTI2025

Relationship between atomic and electronic structure in Ln-bearing oxides

Metastable states of matter are of great interest as they offer the promise of novel functionality. They are often a natural consequence of exposure to nonequilibrium environments. In oxides, metastability can take the form of new polymorphs, chemical disorder, and even amorphization. While significant attention has been given to the impact those changes have on the atomic properties of the material, the corresponding changes in the electronic structure have received less attention. Here, using density functional theory, we consider how the electronic structure varies with potential metastable structures in two classes of lanthanide-bearing oxides—pyrochlores and interlanthanide sesquioxides. We find that the changes depend strongly on both the crystal structure and crystal chemistry of the compound with, for example, disordering and amorphization either increasing or decreasing the bandgap depending on the chemistry. For the 𝐴 2 ⁢𝐵 2 ⁢O 7 pyrochlores, we find different dependencies of the bandgap on the 𝐴 = 𝐿⁢𝑛 cations as the 𝐵 cation is changed, which we relate to the nature of the density of states at the conduction band minimum for different 𝐵 chemistries. Our calculations are validated by electron energy loss spectroscopy measurements for two pyrochlore compounds in which amorphization does reduce the bandgap, consistent with our calculations on these two compounds. In conclusion, our results highlight the relationship between atomic and electronic structure and how radiation can be used to modify and potentially control the electronic properties of oxides.

36 MATERIALS SCIENCE↗
NASA NTRS1988

Spatially resolved ballistic electron spectroscopy of subsurface interfaces

Buried interfaces have traditionally been inaccessible to direct investigation by surface-sensitive techniques. A unique and novel electronic probe, which is sensitive to subsurface electronic structure, has been utilized to probe Schottky barrier interfaces. The method, ballistic electron emission microscopy (BEEM), is based on scanning tunnelling microscopy (STM) techniques. A theoretical treatment has been developed to describe the ballistic electron spectra obtained by BEEM and this treatment is applied to data obtained on the Au-Si Schottky barrier interface system. Excellent agreement between experiment and theory is obtained. In addition, the treatment predicts nanometer spatial resolution for interface imaging by BEEM.

Bell, L. D.↗
NASA NTRS1984

Sensitivity analysis of discrete structural systems: A survey

Methods for calculating sensitivity derivatives for discrete structural systems are surveyed, primarily covering literature published during the past two decades. Methods are described for calculating derivatives of static displacements and stresses, eigenvalues and eigenvectors, transient structural response, and derivatives of optimum structural designs with respect to problem parameters. The survey is focused on publications addressed to structural analysis, but also includes a number of methods developed in nonstructural fields such as electronics, controls, and physical chemistry which are directly applicable to structural problems. Most notable among the nonstructural-based methods are the adjoint variable technique from control theory, and the Green's function and FAST methods from physical chemistry.

Adelman, H. M.↗
NASA NTRS1988

Scanning Tunneling Microscopy methods for spectroscopic imaging of subsurface interfaces

A new method for spatially-resolved, spectroscopic investigation of subsurface interface structure has been developed. The method, Ballistic Electron Emission Microscopy (BEEM), is based on Scanning Tunneling Microscopy (STM) techniques. BEEM combines STM vacuum tunneling with unique ballistic electron spectroscopy capabilities. BEEM enables, for the first time, direct imaging of subsurface interface electronic properties with nanometer spatial resolution. STM topographic images of surface structure and BEEM images of subsurface properties are obtained simultaneously. BEEM capabilities are demonstrated by investigation of important metal-semiconductor interfaces.

Bell, L. D.↗
DOE OSTI2020

Determination of excited state molecular structures from time-resolved gas-phase X-ray scattering

We present a comprehensive investigation of a recently introduced method to determine transient structures of molecules in excited electronic states with sub-ångstrom resolution from time-resolved gas-phase scattering signals. The method, which is examined using time-resolved X-ray scattering data measured on the molecule N-methylmorpholine (NMM) at the Linac Coherent Light Source (LCLS), compares the experimentally measured scattering patterns against the simulated patterns corresponding to a large pool of molecular structures to determine the full set of structural parameters. In addition, we examine the influence of vibrational state distributions and find the effect negligible within the current experimental detection limits, despite that the molecules have a comparatively high internal vibrational energy. Here, the excited state structures determined using three structure pools generated using three different computational methods are in good agreement, demonstrating that the procedure is largely independent of the computational chemistry method employed as long as the pool is sufficiently expansive in the vicinity of the sought structure and dense enough to yield good matches to the experimental patterns.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2023

Ceramic Capacitor Grain Size Analysis Using Electron Backscatter Diffraction (EBSD)

The overall scope of this NESC study was to identify (or develop) a tool, methodology, or process that can inspect a ceramic capacitor and consistently return repeatable grain size distribution results. It is important to note that the development of a machine learning technique was considered as an option for this study, but only if an existing automated tool was not available. This report explains why the Electron Backscatter Diffraction analysis method for grain structure and size distribution was deemed a suitable candidate for this work.

Electron Backscatter Diffraction↗
DOE OSTI2021

Strategies for solving the excited-state self-consistent-field problem for highly excited and multiply ionized states

The dynamics of molecules exposed to intense x-ray radiation involve a large number of multiply ionized and highly excited electronic configurations. To model these dynamics a reliable and efficient electronic structure model is imperative. Employing the Hartree-Fock-Slater electronic structure model in combination with the maximum overlap method, we quantify the associated convergence failures when calculating electronic states of carbon monoxide with multiple vacancies in the core and valence levels. Here, we characterize these cases and describe strategies to overcome the convergence problems. The described techniques not only eliminate all convergence issues for CO but also result in a significant reduction of convergence failures for simulations of the x-ray-induced multiple ionization dynamics of the phenol molecule.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

MISPR : an open-source package for high-throughput multiscale molecular simulations

Computational tools provide a unique opportunity to study and design optimal materials by enhancing our ability to comprehend the connections between their atomistic structure and functional properties. However, designing materials with tailored functionalities is complicated due to the necessity to integrate various computational-chemistry software (not necessarily compatible with one another), the heterogeneous nature of the generated data, and the need to explore vast chemical and parameter spaces. The latter is especially important to avoid bias in scattered data points-based models and derive statistical trends only accessible by systematic datasets. Here, we introduce a robust high-throughput multi-scale computational infrastructure coined MISPR (Materials Informatics for Structure–Property Relationships) that seamlessly integrates classical molecular dynamics (MD) simulations with density functional theory (DFT). By enabling high-performance data analytics and coupling between different methods and scales, MISPR addresses critical challenges arising from the needs of automated workflow management and data provenance recording. The major features of MISPR include automated DFT and MD simulations, error handling, derivation of molecular and ensemble properties, and creation of output databases that organize results from individual calculations to enable reproducibility and transparency. In this work, we describe fully automated DFT workflows implemented in MISPR to compute various properties such as nuclear magnetic resonance chemical shift, binding energy, bond dissociation energy, and redox potential with support for multiple methods such as electron transfer and proton-coupled electron transfer reactions. The infrastructure also enables the characterization of large-scale ensemble properties by providing MD workflows that calculate a wide range of structural and dynamical properties in liquid solutions. MISPR employs the methodologies of materials informatics to facilitate understanding and prediction of phenomenological structure–property relationships, which are crucial to designing novel optimal materials for numerous scientific applications and engineering technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Excited-State and Nonadiabatic Molecular Dynamics Methods with Broad Applicability

Electronically excited states play a pivotal role in the capture, transformation, storage, and emission of radiative energy. Despite much recent progress for small model systems, excited state properties and dynamics in many larger molecular devices and materials have remained elusive for experiment and theory. Nonadiabatic transitions between electronic states turn radiative energy into heat and are particularly difficult to control. This project aims to develop new electronic structure and nonadiabatic molecular dynamics (NAMD) methods to simulate, predict, and ultimately control the flow and transformation of excitation energy in a wide-range of light-driven molecular devices such as dyes, photocatalysts, and fluorescence emitters.

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

Computational Spectroscopy of the Cr–Cr Bond in Coordination Complexes

In this work, we report the accurate computational vibrational analysis of the Cr–Cr bond in dichromium complexes using second-order multireference complete active space methods (CASPT2), allowing direct comparison with experimental spectroscopic data both to facilitate interpreting the low-energy region of the spectra and to provide insights into the nature of the bonds themselves. Recent technological development by the authors has realized such computation for the first time. Accurate simulation of the vibrational structure of these compounds has been hampered by their notorious multiconfigurational electronic structure that yields bond distances that do not correlate with bond order. Some measured Cr–Cr vibrational stretching modes, ν(Cr 2 ), have suggested weaker bonding, even for so-called ultrashort Cr–Cr bonds, while others are in line with the bond distance. Here, we optimize geometries and compute ν(Cr 2 ) with CASPT2 for three well-characterized complexes, Cr 2 (O 2 CCH 3 ) 4 (H 2 O) 2 , Cr 2 (mhp) 4 , and Cr 2 (dmp) 4 . We obtain CASPT2 harmonic ν(Cr 2 ) modes in good agreement with experiment at 282 cm –1 for Cr 2 (mhp) 4 and 353 cm –1 for Cr 2 (dmp) 4 , compute 50 Cr and 54 Cr isotope shifts, and demonstrate that the use of the so-called IPEA shift leads to improved Cr–Cr distances. Additionally, normal mode sampling was used to estimate anharmonicity along ν(Cr 2 ), leading to an anharmonic mode of 272 cm –1 for Cr 2 (mhp) 4 and 333 cm –1 for Cr 2 (dmp) 4 .

36 MATERIALS SCIENCE↗
DOE OSTI2022

Coarse-Grained Density Functional Theory Predictions via Deep Kernel Learning

Scalable electronic predictions are critical for soft materials design. Recently, the Electronic Coarse-Graining (ECG) method was introduced to renormalize all-atom quantum chemical (QC) predictions to coarse-grained (CG) resolutions using deep neural networks (DNNs). While DNNs can learn complex representations that prove challenging for kernel-based methods, they are susceptible to overfitting and the overconfidence of uncertainty estimations. Here, we develop ECG within a GPU-accelerated Deep Kernel Learning (DKL) framework to enable CG QC predictions using range-separated hybrid density functional theory (DFT), obtaining a 107 speedup relative to naive all-atom QC. By treating the predicted electronic properties as random Gaussian Processes, DKL incorporates CG mapping degeneracy by learning the distribution of electronic energies as a function of CG configuration. DKL-ECG accurately reproduces molecular orbital energies from range-separated DFT while facilitating efficient training via active learning using the uncertainties provided by DKL. Further, we show that while active learning algorithms enable efficient sampling of a more diverse configurational space relative to random sampling, all explored query methods exhibit comparable performance for the examined system. We attribute this result to the significant overlap of the feature space and output property distributions across multiple temperatures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Fe Kα XANES, Fe Kβ HERFD XANES and EPMA flank method determinations of the oxidation state of Fe in garnet

The ferric to total iron ratios (Fe 3+ /ΣFe) of garnets can be paired with thermodynamic mineral activity models to quantify the oxygen fugacity of garnet-bearing rocks. However, techniques with a high analytical and spatial resolution are necessary to distinguish differences in garnet Fe 3+ /ΣFe ratios at the percent level and to accurately measure garnets that are zoned or contain inclusions. We acquired conventional Fe Kα and high-resolution energy fluorescence detection (HERFD) Fe Kβ X-ray absorption near edge structure (XANES) spectra and electron microprobe flank method analyses on a suite of 27 peridotitic and eclogitic garnets with Fe 3+ /ΣFe ratios previously determined by Mössbauer spectroscopy to evaluate the precision of each technique. We examined variations in the energy and intensity of three XANES spectral features as a function of Fe 3+ /ΣFe ratios: 1) the intensity ratio of two-post edge features (I-ratio; Fe Kα only); 2) the energy of the Fe edge at 90% normalized intensity (E 0.9 ; Fe Kα only) and 3) the pre-edge centroid energy (Fe Kα and HERFD Fe Kβ). In accordance with previous work, we find the energies of garnet pre-edge centroids are relatively insensitive to Fe 3+ /ΣFe ratios. The I-ratios of peridotitic and eclogitic garnets are offset from each other at low Fe 3+ /ΣFe ratios (≤0.13); I-ratio garnet XANES calibrations are composition-specific. The E 0.9 feature is independent of garnet major element composition in spectra that have been corrected for the effects of self-absorption. We produce two Fe Kα garnet XANES calibrations based on variations in the E0.9 feature; one calibration with all garnet reference materials included (Fe 3+ /ΣFe up to 1.0; “all garnet calibration”) and another calibration specific to garnets with low Fe 3+ /ΣFe ratios (“low ferric calibration”). Fe 3+ /ΣFe ratios calculated from the mean of up to 25 flank method measurements on eight garnet reference materials fall within 4% absolute of a one-to-one correlation with Fe 3+ /ΣFe ratios measured by Mössbauer. The standard error of the mean Fe 3+ /ΣFe ratio calculated from flank method approaches the Mössbauer-determined Fe 3+ /ΣFe ratio within estimated error (3%) after three analyses. Flank method precision is enhanced at higher beam current; however, the precision of the flank method does not approach the precision of XANES under any microprobe analytical condition tested here. Garnet reference materials detailed here are available by request to the Smithsonian Institution.

58 GEOSCIENCES↗
DOE OSTI2023

Femtosecond Electronic and Hydrogen Structural Dynamics in Ammonia Imaged with Ultrafast Electron Diffraction

Directly imaging structural dynamics involving hydrogen atoms by ultrafast diffraction methods is complicated by their low scattering cross-sections. Here we demonstrate that megaelectronvolt ultrafast electron diffraction is sufficiently sensitive to follow hydrogen dynamics in isolated molecules. In a study of the photodissociation of gas phase ammonia, we simultaneously observe signatures of the nuclear and corresponding electronic structure changes resulting from the dissociation dynamics in the time-dependent diffraction. Both assignments are confirmed by ab initio simulations of the photochemical dynamics and the resulting diffraction observable. While the temporal resolution of the experiment is insufficient to resolve the dissociation in time, our results here represent an important step towards the observation of proton dynamics in real space and time.

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

Mechanical and Thermal Properties for Uranium and U–6Nb Alloy from First-Principles Theory

Elasticity, lattice dynamics, and thermal expansion for uranium and U–6Nb alloy (elastic moduli) are calculated from density functional theory that is extended to include orbital polarization (DFT+OP). Introducing 12.5 at.% of niobium, substitutionally, in uranium softens all the cii elastic moduli, resulting in a significantly softer shear modulus (G). Combined with a nearly invariant bulk modulus (B), the quotient B/G increases dramatically for U–6Nb, suggesting a more ductile material. Lattice dynamics from a harmonic model coupled with a DFT+OP electronic structure is applied for α uranium, and the obtained phonon density of states compares well with inelastic neutron-scattering measurements. The Debye temperature associated with the lattice dynamics falls within the range of experimentally observed Debye temperatures and it also validates our quasi-harmonic (QH) phonon model. The QH Debye–Grüneisen phonon method is combined with a DFT+OP electronic structure and used to explore the anisotropic thermal expansion in α uranium. The anomalous negative thermal expansion (contraction) of the b lattice parameter of the α-phase orthorhombic cell is relatively well reproduced from a free-energy model consisting of QH-phonon and DFT+OP electronic structure contributions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗