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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Interaction of molecular nitrogen with vanadium oxide in the absence and presence of water vapor at room temperature: Near-ambient pressure XPS

Interactions of N2 at oxide surfaces are important for understanding electrocatalytic nitrogen reduction reaction (NRR) mechanisms. Interactions of N2 at the polycrystalline vanadium oxide/vapor interface were monitored at room temperature and total pressures up to 10−1 Torr using Near-Ambient Pressure X-ray Photoelectron Spectroscopy (NAP-XPS). The oxide film was predominantly V(IV), with V(III) and V(V) components. XPS spectra were acquired in environments of both pure N2 and equal pressures of N2 and H2O vapor. In pure N2, broad, partially resolved N1s features were observed at binding energies of 401.0 and 398.7 eV, with a relative intensity of ∼3:1, respectively. These features remained upon subsequent pumpdown to 10−9 Torr. The observed maximum N surface coverage was ∼1.5 × 1013 cm−2—a fraction of a monolayer. In the presence of equal pressures of H2O, the adsorbed N intensity at 10−1 Torr is ∼25% of that observed in the absence of H2O. The formation of molecularly adsorbed H2O was also observed. Density functional theory-based calculations suggest favorable absorption energies for N2 bonding to both V(IV) and V(III) cation sites but less so for V(V) sites. Hartree–Fock-based cluster calculations for N2–V end-on adsorption show that experimental XPS doublet features are consistent with the calculated shake-up and normal, final ionic configurations for N2 end-on bonding to V(III) sites but not V(IV) sites. The XPS spectra of vanadium oxide transferred in situ between electrochemical and UHV environments indicate that the oxide surfaces studied here are stable upon exposure to the electrolyte under NRR-relevant conditions.

Balogun, K.↗
DOE OSTI2024

Density Functional Tight-Binding Models for Band Structures of Transition-Metal Alloys and Surfaces across the d -Block

First-principles electronic structure simulations are an invaluable tool for understanding chemical bonding and reactions. While machine-learning models such as interatomic potentials significantly accelerate the exploration of potential energy surfaces, electronic structure information is generally lost. Particularly in the field of heterogeneous catalysis, simulated electron band structures provide fundamental insights into catalytic reactivity. This ab initio knowledge is preserved in semiempirical methods such as density functional tight binding (DFTB), which extend the accessible computational length and time scales beyond first-principles approaches. In this paper here we present Shell-Optimized Atomic Confinement (SOAC) DFTB electronic-part-only parametrizations for bulk and surface band structures of all d-block transition metals that enable efficient predictions of electronic descriptors for large structures or high-throughput studies on complex systems outside the computational reach of density functional theory.

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

Photoionization of the I 4d and valence orbitals of methyl iodide

Here, the photoabsorption and photoionization dynamics of the I 4d and valence orbitals in methyl iodide have been studied both experimentally and theoretically. Synchrotron radiation has been employed to measure the total ion yield in the vicinity of the I 4d ionization thresholds. The observed structure, due to excitations into Rydberg or valence states, has been assigned using transition energies and relative intensities computed with time-dependent density functional theory within the Tamm–Dancoff approximation. Photoelectron spectra, recorded with plane polarized radiation in two polarization geometries, have allowed the effect of autoionization on the valence electron angular distributions to be investigated. The spectra obtained at photon energies of 50.62 and 52.34 eV, coinciding respectively with the I 4d 5/2 → σ* and 4d 3/2 → σ* transitions, reveal, in addition to valence shell photoelectron bands, features not associated with simple photoionization of the parent molecule. High resolution photoelectron spectra of the I 4d main-lines display structure resulting from spin–orbit coupling and molecular field splitting. The binding energies of the five states contributing to the (I 4d) -1 ionization have been determined. The iodine (in CH 3 I) N 45 VV Auger spectrum has been measured and the observed structure has been assigned using the core hole binding energies derived in the present work together with established ionization energies of the doubly charged ion. The experimentally determined Auger electron angular distributions have been discussed in relation to the theoretical angular distribution parameter characterizing the spatial alignment of molecular axes in the (I 4d) -1 state.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Defect Energetics in Pseudo-Cubic Mixed Halide Lead Perovskites from First-Principles

Owing to the increasing popularity of lead-based hybrid perovskites for photovoltaic (PV) applications, it is crucial to understand their defect energetics and its influence on their optoelectronic properties. In this work, we simulate various point defects in pseudocubic structures of mixed iodide-bromide and bromide-chloride methylammonium lead perovskites with the general formula MAPbI 3-y Br y or MAPbBr 3-y Cl y (where y is between 0 and 3), and use first-principles based density functional theory computations to study their relative formation energies and charge transition levels. We identify vacancy defects and Pb on MA antisite defect as the lowest energy native defects in each perovskite. Here, we observe that while the low energy defects in all MAPbI 3-y Br y systems only create shallow transition levels, the Br or Cl vacancy defects in the Cl-containing pervoskites have low energy and form deep levels which become deeper for higher Cl content. We examine the structures and density of states of pure and defect-containing perovskite systems to obtain an understanding of the nature of defect levels. Further, we study extrinsic substitution by different elements at the Pb site in MAPbBr 3 , MAPbCl 3 , and the 50-50 mixed halide perovskite, MAPbBr 1.5 Cl 1.5 , and identify some transition metals that create lower energy defects than the dominant intrinsic defects and also create midgap charge transition levels.

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

Solvent and A-Site Cation Control Preferred Crystallographic Orientation in Bromine-Based Perovskite Thin Films

Preferred crystallographic orientation in polycrystalline films is desirable for efficient charge carrier transport in metal halide perovskites and semiconductors. However, the mechanisms that determine the preferred orientation of halide perovskites are still not well understood. In this work, we investigate crystallographic orientation in lead bromide perovskites. We show that the solvent of the precursor solution and organic A-site cation strongly affect the preferred orientation of the deposited perovskite thin films. Specifically, we show that the solvent, dimethylsulfoxide, influences the early stages of crystallization and induces preferred orientation in the deposited films by preventing colloidal particle interactions. Additionally, the methylammonium A-site cation induces a higher degree of preferred orientation than the formamidinium counterpart. We use density functional theory to show that the lower surface energy of the (100) plane facets in methylammonium-based perovskites, compared to the (110) planes, is the reason for the higher degree of preferred orientation. In contrast, the surface energy of the (100) and (110) facets is similar for formamidinium-based perovskites, leading to lower degree of preferred orientation. Furthermore, we show that different A-site cations do not significantly affect ion diffusion in bromine-based perovskite solar cells but impact ion density and accumulation, leading to increased hysteresis. Our work highlights the interplay between the solvent and organic A-site cation which determine crystallographic orientation and plays a critical role in the electronic properties and ionic migration of solar cells.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Mode-Selective Vibrational Energy Transfer Dynamics in 1,3,5-Trinitroperhydro-1,3,5-triazine (RDX) Thin Films

The coupling of inter- and intramolecular vibrations plays a critical role in initiating chemistry during the shock-to-detonation transition in energetic materials. In this work, we report on the subpicosecond to subnanosecond vibrational energy transfer (VET) dynamics of the solid energetic material 1,3,5-trinitroperhydro-1,3,5-triazine (RDX) by using broadband, ultrafast infrared transient absorption spectroscopy. Experiments reveal VET occurring on three distinct time scales: subpicosecond, 5 ps, and 200 ps. The ultrafast appearance of signal at all probed modes in the mid-infrared suggests strong anharmonic coupling of all vibrations in the solid, whereas the long-lived evolution demonstrates that VET is incomplete, and thus thermal equilibrium is not attained, even on the 100 ps time scale. Density functional theory and classical molecular dynamics simulations provide valuable insights into the experimental observations, revealing compression-insensitive time scales for the initial VET dynamics of high-frequency vibrations and drastically extended relaxation times for low-frequency phonon modes under lattice compression. Mode selectivity of the longest dynamics suggests coupling of the N–N and axial NO2 stretching modes with the long-lived, excited phonon bath.

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

Strongly Correlated States of Transition Metal Spin Defects: The Case of an Iron Impurity in Aluminum Nitride

We investigate the electronic properties of an exemplar transition metal impurity in an insulator, with the goal of accurately describing strongly correlated defect states. Here, we consider iron in aluminum nitride, a material of interest for hybrid quantum technologies, and we carry out calculations with quantum embedding methods, density matrix embedding theory (DMET) and quantum defect embedding theory (QDET), and with spin-flip time-dependent density functional theory (TDDFT). We show that both DMET and QDET accurately describe the ground state and low-lying excited states of the defect and that TDDFT yields photoluminescence spectra in agreement with experiments. In addition, we provide a detailed discussion of the convergence of our results as a function of the active space used in the embedding methods, thus defining a protocol to obtain converged data directly comparable with experiments.

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

Low-temperature etching of silicon oxide and silicon nitride with hydrogen fluoride

Etching of high aspect ratio features into alternating SiO 2 and SiN layers is an enabling technology for the manufacturing of 3D NAND flash memories. In this paper, we study a low-temperature or cryo plasma etch process, which utilizes HF gas together with other gas additives. Compared with a low-temperature process that uses separate fluorine and hydrogen gases, the etching rate of the SiO 2 /SiN stack doubles. Both materials etch faster with this so-called second generation cryo etch process. Pure HF plasma enhances the SiN etching rate, while SiO 2 requires an additional fluorine source such as PF 3 to etch meaningfully. Further, the insertion of H 2 O plasma steps into the second generation cryo etch process boosts the SiN etching rate by a factor of 2.4, while SiO 2 etches only 1.3 times faster. We observe a rate enhancing effect of H 2 O coadsorption in thermal etching experiments of SiN with HF. Ammonium fluorosilicate (AFS) plays a salient role in etching of SiN with HF with and without plasma. AFS appears weakened in the presence of H 2 O. Density functional theory calculations confirm the reduction of the bonding energy when NH 4 F in AFS is replaced by H 2 O.

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

Magnetoelastic interactions in SrCu 2 (BO 3 ) 2 studied by Raman scattering experiments and first principles calculations

Dynamic and static crystal lattice properties of SrCu 2 (BO 3 ) 2 are studied by means of Raman scattering, magnetostriction, and thermal expansion measurements in magnetic fields to 45 T. Raman experiments versus temperature reveal that some phonon modes show an unusual behavior: their frequencies soften (modes at 200 and 450 cm –1 ) while others harden (modes at 385 and 478 cm –1 ) when decreasing the temperature below 15 K. Magneto-Raman experiments show that their field dependence correlates with their respective temperature dependencies; e.g., modes that are hardened with increasing temperature also harden with applied magnetic fields and modes that become softer with temperature also soften with applied fields. We use density functional theory to successfully model and compute the energies of these modes, classifying them into two types: pantograph (modes that soften when decreasing the temperature) and nonpantograph. We conclude that the former involves the modification of the intradimer exchange interaction J and the latter the interdimer J'. Lastly, dilatometry is used to correlate field-dependent Raman modes to the closing of the spin gap as well as fractional-magnetization stripe states M = 1/4 M s and M = 1/3 M s , where M s is the saturation magnetization.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Structure sensitivity in adsorbate-induced adatom formation on FCC transition-metal surfaces

Recent surface science discoveries reveal that adsorbates may induce in situ sub-nanometer cluster formation on transition-metal surfaces. To elucidate the structure sensitivity behind this phenomenon, we performed density functional theory calculations to construct an adatom formation energy database for eight fcc metals and 26 adsorbates commonly involved in catalytic reactions. We show that the adatom formation on (100) surfaces is generally easier than that on (111) surfaces. Many adsorbate/metal pairs exist, mostly on Pd, Ni, Rh, Pt, and Ir, for which adsorbates might induce adatom formation under near-ambient conditions on the (100) facet, but not on (111), highlighting the role of more open facets in adatom formation on metals intrinsically harder than Ag, Au, and Cu. Furthermore, our study offers a new perspective towards understanding structure sensitivity in heterogeneous thermal- and electro-catalytic systems such as methane steam reforming, Fischer-Tropsch synthesis, and ammonia decomposition.

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

Platinum Graphene Catalytic Condenser for Millisecond Programmable Metal Surfaces

Accelerating catalytic chemistry and tuning surface reactions require precise control of the electron density of metal atoms. In this work, nanoclusters of platinum were supported on a graphene sheet within a catalytic condenser device that facilitated electron or hole accumulation in the platinum active sites with negative or positive applied potential, respectively. The catalytic condenser was fabricated by depositing on top of a p-type Si wafer an amorphous HfO 2 dielectric (70 nm), on which was placed the active layer of 2–4 nm platinum nanoclusters on graphene. A potential of ±6 V applied to the Pt/graphene layer relative to the silicon electrode moved electrons into or out of the active sites of Pt, attaining charge densities more than 1% of an electron or hole per surface Pt atom. At a level of charge condensation of ±10% of an electron per surface atom, the binding energy of carbon monoxide to a Pt(111) surface was computed via density functional theory to change 24 kJ mol –1 (0.25 eV), which was consistent with the range of carbon monoxide binding energies determined from temperature-programmed desorption (ΔBE CO of 20 ± 1 kJ mol –1 or 0.19 eV) and equilibrium surface coverage measurements (ΔBE CO of 14 ± 1 kJ mol –1 or 0.14 eV). Impedance spectroscopy indicated that Pt/graphene condensers with potentials oscillating at 3000 Hz exhibited negligible loss in capacitance and charge accumulation, enabling programmable surface conditions at amplitudes and frequencies necessary to achieve catalytic resonance.

Inorganic carbon compounds↗
DOE OSTI2020

Exploring Ca–Ce–M–O (M = 3d Transition Metal) Oxide Perovskites for Solar Thermochemical Applications

Solar thermochemical (STC) processes hold promise as efficient ways to generate renewable fuels, fuel precursors, or chemical feedstocks using concentrated sunlight. Specifically, one actively researched approach is the two-step STC cycle, which uses a redox-active, off-stoichiometric, transition-metal oxide material to split water and/or CO 2 , generating H 2 and/or CO, respectively, or syngas (a combination of H 2 and CO). Identifying novel metal oxides that yield larger reduction extents (practically achievable off-stoichiometries) than the state-of-the-art CeO 2 is critical. Here, we explore the chemical space of Ca–Ce–M–O (M = 3d transition metal, except Cu and Zn) metal oxide perovskites, with Ca and/or Ce occupying the A site and M occupying the B site within an ABO 3 framework, as potential STC candidates. We use density functional theory (DFT)-based calculations and systematically evaluate the oxygen vacancy (VaO) formation energy (≈ enthalpy of reduction in an STC cycle), electronic properties, thermodynamic stability of CaMO 3 , CeMO 3 , and Ca 0.5 Ce 0.5 MO 3 perovskites, and the VaO formation energy within Ca 0.5 Ce 0.5 Ti 0.5 Mg 0.5 O 3 perovskite. We consider only Ca and/or Ce on the A site because of their similar size and the potential redox activity of Ce 4+ . If both Ce and M exhibit simultaneous reduction with Va O formation, the resulting perovskite could exhibit a larger entropy of reduction than a single cation reduction. The increased entropy produces increased reduction for fixed temperature, partial pressure of oxygen, and reduction enthalpy, and therefore increased STC efficiency. Importantly, we identify Ca 0.5 Ce 0.5 MnO 3 , Ca 0.5 Ce 0.5 FeO 3 , and Ca 0.5 Ce 0.5 VO 3 to be promising candidates based on their Va O formation energy and thermodynamic (meta)stability. Moreover, based on our calculated on-site magnetic moments, electron density of states, and electron density differences between pristine and defective structures, we find Ca 0.5 Ce 0.5 MnO 3 to exhibit simultaneous reduction of both Ce 4+ (A-site) and Mn 3+ (B-site), highlighting a particularly promising candidate for STC applications with a predicted higher entropy of reduction than CeO 2 . Lastly, we extract metrics that govern the trends in Va O formation energies, such as standard reduction potentials, and provide pointers for further experimental and theoretical studies, which will enable the design of improved materials for the STC cycle.

14 SOLAR ENERGY↗
DOE OSTI2024

Zentropy theory for accurate prediction of free energy, volume, and thermal expansion without fitting parameters

Based on statistical mechanics, a macroscopically homogeneous system, i.e., a single phase in the present context, is composed of many independent configurations that the system embraces. The macroscopical properties of the system are determined by the properties and statistical probabilities of those configurations with respect to external conditions. The volume of a single phase is thus the weighted sum of the volumes of all configurations. Consequently, the derivative of the volume to temperature of a single phase depends on both the derivatives of the volumes of every configuration to temperature and the derivatives of their statistical probabilities to temperature, with the latter introducing nonlinear emergent behaviors. It is shown that the derivative of the volume to the temperature of the single phase can be negative, i.e., negative thermal expansion, due to the symmetry-breaking non-ground-state configurations with smaller volumes than that of the ground-state configuration and the rapid increase of the statistical probabilities of the former, and negative thermal expansion can be predicted without fitting parameters from the zentropy theory that combines quantum mechanics and statistical mechanics with the free energy of each configuration predicted from quantum mechanics and the partition function of each configuration calculated from its free energy.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ethane Dehydrogenation on Single and Dual Centers of Ga-modified γ-Al 2 O 3

Density-functional theory calculations and microkinetic analysis are used to investigate the efficacy of Ga-modified γ-Al 2 O 3 (110) surfaces for the catalytic dehydrogenation of ethane and elucidate the synergy between Ga and Al sites. The model surfaces are modified by either Ga grafting or doping. We consider and analyze numerous active sites and rank them using microkinetic analysis. The kinetic parameters obtained from microkinetic modelling are compared with experimental values for ethane dehydrogenation over Ga 2 O 3 -Al 2 O 3 mixed oxides prepared by co-precipitation. The dominant reaction pathway proceeds via heterolytic C-H bond dissociation to a surface proton and a metal-carbanion intermediate that undergoes β-hydride elimination. We find that grafted Ga sites are catalytically inactive. In contrast, Ga-doped sites exhibit fivefold enhancement in catalytic activity when compared to the sites on pristine Al 2 O 3 , owed to the synergy between neighboring AlIII and GaIV sites. Furthermore, we model and investigate the effect of surface hydroxylation, demonstrate how surface water interferes with the aforementioned synergy between AlIII and GaIV sites and discuss the implications for the catalytic activity of the modified surfaces. Here, an increase in the partial pressure of H 2 O significantly increases the apparent activation energies of dehydrogenation and interestingly changes the most active site.

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

Accurate Free Energies for Complex Condensed-Phase Reactions Using an Artificial Neural Network Corrected DFTB/MM Methodology

Semiempirical methods like density functional tight-binding (DFTB) allow extensive phase space sampling, making it possible to generate free energy surfaces of complex reactions in condensed-phase environments. Such a high efficiency often comes at the cost of reduced accuracy, which may be improved by developing a specific reaction parametrization (SRP) for the particular molecular system. Thiol–disulfide exchange is a nucleophilic substitution reaction that occurs in a large class of proteins. Its proper description requires a high-level ab initio method, while DFT-GAA and hybrid functionals were shown to be inadequate, and so is DFTB due to its DFT-GGA descent. We develop an SRP for thiol–disulfide exchange based on an artificial neural network (ANN) implementation in the DFTB+ software and compare its performance to that of a standard SRP approach applied to DFTB. Furthermore, as an application, we use both new DFTB-SRP as components of a QM/MM scheme to investigate thiol–disulfide exchange in two molecular complexes: a solvated model system and a blood protein. Demonstrating the strengths of the methodology, highly accurate free energy surfaces are generated at a low cost, as the augmentation of DFTB with an ANN only adds a small computational overhead.

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

Extreme sensitivity of higher-order interatomic force constants and thermal conductivity to the energy surface roughness of exchange-correlation functionals

In this Letter, we report that the fourth-order interatomic force constants (4th-IFCs) are significantly sensitive to the energy surface roughness of exchange-correlation (XC) functionals in density functional theory calculations. This sensitivity, which is insignificant for the second- (2nd-) and third-order (3rd-) IFCs, varies for different functionals in different materials and can cause misprediction of thermal conductivity by several times of magnitude. As a result, when calculating the 4th-IFCs using the finite difference method, the atomic displacement needs to be taken large enough to overcome the energy surface roughness, in order to accurately predict phonon lifetime and thermal conductivity. We demonstrate this phenomenon on a benchmark material (Si), a high-thermal conductivity material (BAs), and a low thermal conductivity material (NaCl). For Si, we find that the LDA, PBE, and PBEsol XC functionals are all smooth to the 2nd- and 3rd-IFCs but all rough to the 4th-IFCs. This roughness can lead to a prediction of nearly one order of magnitude lower thermal conductivity. For BAs, all three functionals are smooth to the 2nd- and 3rd-IFCs, and only the PBEsol XC functional is rough for the 4th-IFCs, which leads to a 40% underestimation of thermal conductivity. For NaCl, all functionals are smooth to the 2nd- and 3rd-IFCs but rough to the 4th-IFCs, leading to a 70% underprediction of thermal conductivity at room temperature. In conclusion, with these observations, we provide general guidance on the calculation of 4th-IFCs for an accurate thermal conductivity prediction.

42 ENGINEERING↗
DOE OSTI2024

Influence of excess silicon on polytype selection during metal-mediated epitaxy of GaN nanowires

For this work, we have examined the origins of polytype selection during metal-mediated molecular-beam epitaxy of GaN nanowires (NWs). High-angle annular dark-field scanning transmission electron microscopy reveals [111]-oriented zinc blende (ZB) NWs and [0001]-oriented wurtzite (WZ) NWs, with SixNy at the interface between individual NWs and the Si (001) substrate. Quantitative energy dispersive x-ray spectroscopy reveals a notably higher Si concentration of 7.0% ± 2.3% in zinc blende (ZB) NWs than 2.3% ± 1.2% in wurtzite (WZ) NWs. Meanwhile, density functional theory calculations show that incorporation of 8 at. % Si on the Ga sublattice inverts the difference in formation energies between WZ and ZB GaN, such that the ZB polytype of GaN is stabilized. This identification of Si and other ZB polytype stabilizers will enable the development of polytype heterostructures in a wide variety of WZ-preferring compounds.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Machine learned features from density of states for accurate adsorption energy prediction

Materials databases generated by high-throughput computational screening, typically using density functional theory (DFT), have become valuable resources for discovering new heterogeneous catalysts, though the computational cost associated with generating them presents a crucial roadblock. Hence there is a significant demand for developing descriptors or features, in lieu of DFT, to accurately predict catalytic properties, such as adsorption energies. Here, we demonstrate an approach to predict energies using a convolutional neural network-based machine learning model to automatically obtain key features from the electronic density of states (DOS). The model, DOSnet, is evaluated for a diverse set of adsorbates and surfaces, yielding a mean absolute error on the order of 0.1 eV. In addition, DOSnet can provide physically meaningful predictions and insights by predicting responses to external perturbations to the electronic structure without additional DFT calculations, paving the way for the accelerated discovery of materials and catalysts by exploration of the electronic space.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗