Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy density functional”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Extracting correlation effects from momentum-resolved electron energy loss spectroscopy: Synergistic origin of the dispersion kink in Bi 2.1 Sr 1.9 CaCu 2 O 8 + x

We employ momentum-resolved electron energy loss spectroscopy (M-EELS) on Bi 2.1 Sr 1.9 CaCu 2 O 8 + x to resolve the issue of the kink feature in the electron dispersion widely observed in the cuprates. To this end, we utilize the GW approximation to relate the density response function measured in M-EELS to the self-energy, isolating contributions from phonons, electrons, and the momentum dependence of the effective interaction to the decay rates. The phononic contributions, present in the M-EELS spectra due to electron-phonon coupling, lead to kink features in the corresponding single-particle spectra at energies between 40 and 80 meV, independent of the doping level. We find that a repulsive interaction constant in momentum space is able to yield the kink attributed to phonons in ARPES. Hence, our analysis of the M-EELS spectra points to local repulsive interactions as a factor that enhances the spectroscopic signatures of electron-phonon coupling in cuprates. We conclude that the strength of the kink feature in cuprates is determined by the combined action of electron-phonon coupling and electron-electron interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1984

Electron escape depth variation in thin SiO2 films measured with variable photon energy

A double crystal monochromator at the Stanford Synchrotron Radiation Laboratory is used to study the Si/SiO2 interface, using photon energies of hv = 1950-3700 eV. This photon energy range allows interfaces to be observed through oxide layers 50 A thick or more. Variations in electron escape depth and/or oxide density as a function of distance from the interface are observed over the entire kinetic energy range (100-3600 eV). These differences are attributed to a strained oxide layer near the interface.

Hecht, M. H.↗
DOE OSTI2022

A Universal Machine Learning Model for Elemental Grain Boundary Energies

The grain boundary (GB) energy has a profound influence on the grain growth and properties of polycrystalline metals. Here, we show that the energy of a GB, normalized by the bulk cohesive energy, can be described purely by four geometric features. By machine learning on a large computed database of 361 small Σ (Σ<10) GBs of more than 50 metals, we develop a model that can predict the grain boundary energies to within a mean absolute error of 0.13 J m –2 . More importantly, this universal GB energy model can be extrapolated to the energies of high Σ GBs without loss in accuracy. These results highlight the importance of capturing fundamental scaling physics and domain knowledge in the design of interpretable, extrapolatable machine learning models for materials science.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Effect of Zn Addition on Phase Evolution in AlCrFeCoNiZn High–Entropy Alloy

The addition of Zn to AlCrFeCoNi high-entropy alloy (HEA) poses intriguing questions as to how it would affect phase evolution. Herein, the phase evolution in AlCrFeCoNiZn is studied using a combination of experimental techniques (X-ray diffraction, scanning electron microscopy, energy-dispersive spectroscopy, and differential scanning calorimetry) and computational (density-functional theory [DFT], calculation of phase diagrams, and machine-learning) methods. Mechanically alloyed and spark-plasma-sintered AlCrFeCoNiZn assumes a metastable single-phase, body-centered-cubic (BCC) structure that undergoes diffusion-controlled phase separation upon subsequent heat treatment to form separate (Al, Cr)-rich, (Fe, Co)-rich, and (Zn, Ni)-rich phases. The formation of (Al, Cr)-rich phase, not reported previously in AlCrFeCoNi-based HEAs, is attributed to strong clustering tendency of Cr–Zn and Cr–Ni pairs, combined with the strong ordering of Zn–Ni pair, driving out Cr that in turn combines with Al to form a (Al, Cr)-rich phase. In the DFT results, the formation of thermodynamically stable L1 2 phase is shown wherein Cr–Fe–Zn [Al–Ni-Co] preferably occupy1a (000) [3c (0 ½ ½)] positions. Furthermore, the sluggish diffusional transformation to L1 2 phase from BCC precursors is attributed to the small stacking-fault energy of AlCrFeCoNiZn. The equilibrated HEA exhibits a high microhardness of 8.24 GPa with an elastic modulus of 184 GPa.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Design of refractory multi-principal-element alloys for high-temperature applications

Abstract Refractory multi-principal-element alloys (RMPEAs) exhibit high specific strength at elevated temperatures ( T ). However, current RMPEAs lack a balance of room-temperature (RT) ductility, high- T strength, and high- T creep resistance. Using density-functional theory methods, we scanned composition space using four criteria: (1) formation energies for operational stability: $$-150\le {E}_{{\rm {f}}}$$ − 150 ≤ E f ≤ +70 meV per atom; (2) higher strength found via interstitial electron density with Young’s moduli E > 250 GPa; (3) inverse Pugh ratio for ductility: G / B < 0.57; and (4) high melting points: T m > 2500 °C. Using rapid bulk alloy synthesis and characterization, we validated theory and down-selected promising alloy compositions and discovered Mo 72.3 W 12.8 Ta 10.0 Ti 2.5 Zr 2.5 having well-balanced RT and high- T mechanical properties. This alloy has comparable high- T compressive strength to well-known MoNbTaW but is more ductile and more creep resistant. It is also superior to a commercial Mo-based refractory alloy and a nickel-based superalloy (Haynes-282) with improved high- T tensile strength and creep resistance.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ordering in GaAs co-doped with Bi and N

Introducing only a few atomic percent of Bi or N in GaAs has a large effect on the band gap of the material. Specifically Bi doped GaAs shows potential for local band gap engineering in optoelectronic applications. The incorporation of Bi and N into GaAs is difficult due to strain effects. In this work we study the ordering of these dopants at the atomic scale in order to get a better understanding of the behavior of these dopants in the host lattice. Cross-sectional scanning tunneling microscopy is used to find the exact position of Bi and N dopants in the GaAs matrix, allowing us to study both their nearest neighbor pair occurrences and pair correlation functions. An attractive interaction between Bi dopants at short ranges (1-2 nm) is found and a similar effect is observed between N dopants. Here, we find a repulsive interaction with a similar length scale between Bi and N dopants. A similar repulsion is found in the Bi-N nearest neighbor pairs. Density functional theory is used to calculate the different nearest neighbor pair energies and test these results to the experimental pair occurrences. It is concluded from the experimental and theoretical results that the growth conditions and N inclusion greatly affects the Bi distribution in GaAs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Assessing the accuracy of compound formation energies with quantum Monte Carlo

Accurately predicting the formation energy of a compound, which describes its thermodynamic stability, is a key challenge in materials physics. Here, we employ a many-body quantum Monte Carlo (QMC) with single-reference trial functions to compute the formation energy of two electronically disparate compounds, the intermetallic VPt2 and the semiconductor CuI, for which standard density functional theory (DFT) predictions using both the Perdew-Burke-Ernzerhof (PBE) and the strongly constrained and appropriately normed (SCAN) density functional approximations deviate markedly from available experimental values. For VPt2, we find an agreement between QMC, SCAN, and PBE0 estimates, which therefore remain in disagreement with the much less exothermic experimental value. For CuI, the QMC result agrees with neither SCAN nor PBE pointing towards DFT exchange-correlation biases, likely related to the localized Cu 3d electrons. Compared to the behavior of some density functional approximations within DFT, spin-averaged QMC exhibits a smaller but still appreciable deviation when compared to experiment. The QMC result is slightly improved by incorporating spin-orbit corrections for CuI and solid I 2 , so that experiment and theory are brought into imperfect but reasonable agreement within about 120 meV/atom.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1987

Electron beam experiments at high altitudes

Experiments with the electron gun on the SCATHA satellite produced evidence of beam-plasma interactions, and heating of the low energy electrons around the satellite. These experiments were conducted near geosynchronous orbit, in the dusk bulge, and plasma sheet, with one short operation in the lobe regions, providing a range of ambient plasma densities. The electron gun was operated at 50 eV, with beam currents of 1, 10, and 100 micro-A. Data from electrostatic analyzers and the dc electric field experiment show that the satellite charged to near the beam energy in sunlight, if the beam current was sufficient. Higher ambient densities required higher beam currents. The electrostatic analyzers showed distribution functions which had peaks, or plateaus, at energies greater than the satellite potential. These measurements indicate heating of the ambient plasma at several Debye lengths from the satellite, with the heated plasma then accelerated into the satellite. It is likely that the ambient plasma is in fact the photoelectron sheath generated by the satellite.

Olsen, R. C.↗
NASA NTRS1988

Electron beam experiments at high altitudes

Experiments with the electron gun on the SCATHA satellite produced evidence of beam-plasma interactions, and heating of the low energy electrons around the satellite. These experiments were conducted near geosynchronous orbit, in the dusk bulge, and plasma sheet, with one short operation in the lobe regions, providing a range of ambient plasma densities. The electron gun was operated at 50 eV, with beam currents of 1, 10, and 100 micro-A. Data from electrostatic analyzers and the DC electric field experiment show that the satellite charged to near the beam energy in sunlight, if the beam current was sufficient. Higher ambient densities required higher beam currents. The electrostatic analyzers showed distribution functions which had peaks, or plateaus, at energies greater than the satellite potential. These measurements indicate heating of the ambient plasma at several Debye lengths from the satellite, with the heated plasma then accelerated into the satellite. It is likely that the ambient plasma is in fact the photoelectron sheath generated by the satellite.

Olsen, R. C.↗
DOE OSTI2021

Nickel B-site substitution in bulk Sr 1-x Ca x FeO 3 perovskite oxygen carriers: Benefits and limitations

We report oxygen (O 2 ) storage materials often rely on the presence of cobalt (Co) to reduce the thermodynamic penalty and increase the kinetics necessary for efficient O 2 storage and release. In this work, we investigate nickel (Ni) as an alternative B-site dopant in Sr 1-x Ca x FeO 3 to identify Co-free carriers that still show improved kinetics at low temperatures. In fact, we show a substantial increase in the reversible O2 release rate through mild Ni B-site substitution (y = 0.06) in select Sr 1-x Ca x Fe 1-y Ni y O 3 systems at 400 to 500 °C, reaching 2.00 wt.% O 2 release up to approximately 75% faster than Ni-free systems. To explain the role of Ni in these systems, we use density functional theory to calculate the O 2 vacancy (V O ) formation energy from separate metal-oxygen (M-O) bonding and relaxation components. We computationally show elongated Ni-O bonds are directly responsible for the decrease in V O upon Ni substitution.

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

Computational catalyst discovery: Active classification through myopic multiscale sampling

We report the recent boom in computational chemistry has enabled several projects aimed at discovering useful materials or catalysts. We acknowledge and address two recurring issues in the field of computational catalyst discovery. First, calculating macro-scale catalyst properties is not straightforward when using ensembles of atomic-scale calculations [e.g., density functional theory (DFT)]. We attempt to address this issue by creating a multi-scale model that estimates bulk catalyst activity using adsorption energy predictions from both DFT and machine learning models. The second issue is that many catalyst discovery efforts seek to optimize catalyst properties, but optimization is an inherently exploitative objective that is in tension with the explorative nature of early-stage discovery projects. In other words, why invest so much time finding a “best” catalyst when it is likely to fail for some other, unforeseen problem? We address this issue by relaxing the catalyst discovery goal into a classification problem: “What is the set of catalysts that is worth testing experimentally?” Here, we present a catalyst discovery method called myopic multiscale sampling, which combines multiscale modeling with automated selection of DFT calculations. It is an active classification strategy that seeks to classify catalysts as “worth investigating” or “not worth investigating” experimentally. Our results show an ~7–16 times speedup in catalyst classification relative to random sampling. These results were based on offline simulations of our algorithm on two different datasets: a larger, synthesized dataset and a smaller, real dataset.

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

An Accurate, Affordable Density Functional Tight-Binding Model for Excited State Hydrocarbon Polymer Molecular Dynamics

We have developed a density functional tight-binding model for hydrocarbon excited state dynamics by referencing to high-level electronic structure theory and incorporating a many-body repulsive energy. We then validate our model against n-octane geometry optimizations, bond dissociation scans, and vibrational frequencies. Our model is approximately 10 3 times more efficient than hybrid time-dependent density functional theory calculations with comparable accuracy. In conclusion, our efforts enable longer timescale excited state simulations of photochemistry and scattering of incident radiation.

Chemical calculations↗
DOE OSTI2021

Accurate Thermochemistry of Complex Lignin Structures via Density Functional Theory, Group Additivity, and Machine Learning

A molecular-level understanding of lignin structures and bond dissociation energies could facilitate depolymerization technologies. Still, this information is currently limited due to the lack of databases and the simplification of surrogate models. Here, substitution effects on seven common linkages in lignin polymers are systematically investigated. An automated reaction network generator is employed to create a database of structures. A new group additivity (GA) model based on principal component analysis (PCA) descriptors is introduced and trained on gas-phase density functional theory data of 4100 species at the M06-2X/6-311++G(d,p) level. Hydrogen bonds, local steric, and nonaromatic ring contributions are also incorporated. Lastly, we improve the accuracy of the group additivity model to reach the G4 theory by computing a data set of 770 species at this level and using a data fusion approach.

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

A non-perturbative pairwise-additive analysis of charge transfer contributions to intermolecular interaction energies

Energy decomposition analysis (EDA) based on absolutely localized molecular orbitals (ALMOs) decomposes the interaction energy between molecules into physically interpretable components like geometry distortion, frozen interactions, polarization, and charge transfer (CT, also sometimes called charge delocalization) interactions. In this work, a numerically exact scheme to decompose the CT interaction energy into pairwise additive terms is introduced for the ALMO-EDA using density functional theory. Unlike perturbative pairwise charge-decomposition analysis, the new approach does not break down for strongly interacting systems, or show significant exchange–correlation functional dependence in the decomposed energy components. Both the energy lowering and the charge flow associated with CT can be decomposed. Complementary occupied–virtual orbital pairs (COVPs) that capture the dominant donor and acceptor CT orbitals are obtained for the new decomposition. It is applied to systems with different types of interactions including DNA base-pairs, borane-ammonia adducts, and transition metal hexacarbonyls. While consistent with most existing understanding of the nature of CT in these systems, we report the results also reveal some new insights into the origin of trends in donor–acceptor interactions.

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

Density-Functional Tight-Binding Parameters for Bulk Zirconium: A Case Study for Repulsive Potentials

Density-functional tight-binding (DFTB) parameters are presented for the simulation of the bulk phases of zirconium. Electronic parameters were obtained using a band structure fitting strategy, while two-center repulsive potentials were created by particle swarm optimization. As objective functions for the repulsive potential fitting, we employed the Birch–Murnaghan equations of state for hexagonal close-packed (HCP), body-centered cubic (BCC) and ω phases of Zr from density-functional theory (DFT). When fractional atomic coordinates are not allowed to change in the generation of the equation-of-state curves, long-range repulsive DFTB potentials are able to almost perfectly reproduce equilibrium structures, relative DFT energies of the bulk phases, and bulk moduli. However, the same potentials lead to artifacts in the DFTB potential energy surfaces when atom positions in the unit cell are allowed to fully relax during the change of unit cell parameters. Finally, conventional short-range repulsive DFTB potentials, while inferior in their ability to reproduce DFT bulk energetics, are able to correctly reproduce the qualitative shape of the DFT potential energy surfaces, including the location of global minima, and can therefore be considered more transferable.

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

Targeted Chemical Looping Materials Discovery by an Inverse Design

Chemical looping with oxygen uncoupling (CLOU) materials is actively sought for combustion of carbonaceous materials to achieve complete conversion and capture of carbon dioxide. These materials may play a vital role in reducing atmospheric carbon via negative carbon output. However, there is no one‐size‐fits‐all approach as different operating conditions and feedstocks may require different CLOU materials. As a result, the exploration and discovery of high‐performance CLOU materials can be a slow process. To address this challenge, a high‐throughput inverse machine learning workflow that identifies optimum materials from perovskite oxides for a given set of targets is developed—temperature and Gibbs free energy of oxygen formation. The model is trained on high‐throughput density functional theory calculations of CLOU materials and inverts the materials design process using a genetic algorithm to produce realistic substituted SrFeO 3‐δ compositions as output. Using the inverse model, it is able to identify several interesting new families of CLOU materials: Sr 1‐ x A x Fe 1‐ y B y O 3‐δ (e.g., A = Ca or K; B = Mg, Bi, Mn, Ni, Co, Cu, or Zn). These materials have shown promising properties, and some of them even outperform the benchmark material in terms of oxygen release kinetics under relevant CLOU operating conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Deformation Behaviors in Single BCC‐Phase Refractory Multi‐Principal Element Alloys under Dynamic Conditions

The mechanical behavior and microstructural evolution of a BCC-phase NbTaTiV refractory multi-principal element alloy (RMPEA) is studied over a wide range of strain rates (10 −3 to 10 3 s −1 ) and temperatures (room temperature to 850 °C). The mechanical property of present RMPEA shows less strain-rate dependence and strong resistance to softening at high temperatures. Under high strain-rate loading, the formation of thin type-I twins is observed, which could lead to an increase in strain-hardening rates. However, this hardening mechanism competes with adiabatic heating effects, resulting in the deterrence of strain-hardening behaviors. In contrast, substantial strain-hardening occurs at cryogenic temperatures due to the formation of twins, which act as stronger barriers to dislocation motion and interact with each other. To further understand the different strain-hardening behaviors, density functional theory (DFT) calculations predict relatively low stacking fault energies and high twinning stress for the NbTaTiV RMPEA.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dissecting Critical Factors for Electrochemical CO 2 Reduction on Atomically Precise Au Nanoclusters

Here, this work investigates the critical factors impacting electrochemical CO 2 reduction reaction (CO 2 RR) using atomically precise Au nanoclusters (NCs) as electrocatalysts. First, the influence of size on CO 2 RR is studied by precisely controlling NC size in the 1–2.5 nm regime. We find that the electrocatalytic CO partial current density increases for smaller NCs, but the CO Faradaic efficiency (FE) is not directly associated with the NC size. This indicates that the surface-to-volume ratio, i.e. the population of active sites, is the dominant factor for determining the catalytic activity, but the selectivity is not directly impacted by size. Second, we compare the CO 2 RR performance of Au 38 isomers (Au 38 Q and Au 38 T) to reveal that structural rearrangement of identical size NCs can lead to significant changes in both CO 2 RR activity and selectivity. Au 38 Q shows higher activity and selectivity towards CO than Au 38 T, and density functional theory (DFT) calculations reveal that the average formation energy of the key *COOH intermediate on the proposed active sites is significantly lower on Au 38 Q than Au 38 T. These results demonstrate how the structural isomerism can impact stabilization of reaction intermediates as well as the overall CO 2 RR performance of identical size Au NCs. Overall, this work provides important structure–property relationships for tailoring the NCs for CO 2 RR.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗