Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional calculation”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Fundamental Insights into Cathode Stability: Linking Compositional Tuning and Local Coordination in Complex Metal Oxides under Aqueous Transformations

Compositional tuning of complex metal oxides in Li-ion battery materials influences their performance as well as their end-of-life behavior, in particular, the tendency to release toxic metal cations in aqueous solution. We modeled ternary variants of a parent LiCoO 2 delafossite structure by varying the metal identity and relative amounts. This yielded ten model formulations of Li(A 4/6 B 1/6 C 1/6 )O 2 , where the material is enriched with the A metal and doped with B and C, with Ni, Mn, Co, Fe, Al, V, and Ti as constituent metals. To assess their stability in aqueous conditions, metal release energetics were calculated using a combination of Density Functional Theory calculations and thermodynamics. Metal release in ternary oxides is dictated by subtle variations in the coordination environment of the leaving group. To identify governing chemical features across diverse compositions with varying local coordination environments, we leverage random forest regression and descriptor importance analysis. A key result is that metal–oxygen orbital hybridization, quantified using a projected density-of-states-derived descriptor, H d/p , provides a physically grounded measure of interaction strength that governs metal release energetics. This refined perspective goes beyond conventional oxidation state considerations and offers more robust insights for materials science. Finally, we model defect surface-bound O 2 dimer formation as a proxy for reactive oxygen species (ROS) generation. The results show that Ni-rich compositions more readily stabilize spin-polarized O 2 dimers, corroborating experimental reports of an increased ROS-driven biological response. In conclusion, our results establish a compositional and electronic basis for metal release and surface oxygen reactivity that form a rationale for complex metal oxide design principles.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Photoelectron spectra of early 3 d -transition metal dioxide molecular anions from GW calculations

Photoelectron spectra of early 3 d -transition metal dioxide anions, Sc O 2 - , Ti O 2 - , V O 2 - , Cr O 2 - , and Mn O 2 - , are calculated using semilocal and hybrid density functional theory (DFT) and many-body perturbation theory within the GW approximation using one-shot perturbative and eigenvalue self-consistent formalisms. Different levels of theory are compared with each other and with available photoelectron spectra. We show that one-shot GW with a PBE0 starting point ( G 0 W 0 @PBE0) consistently provides very good agreement for all experimentally measured binding energies (within 0.1 eV–0.2 eV or less). We attribute this to the success of PBE0 in mitigating self-interaction error and providing good quasiparticle wave functions, which renders a first-order perturbative GW correction effective. One-shot GW calculations with a Perdew–Burke–Ernzerhof (PBE) starting point do poorly in predicting electron removal energies by underbinding orbitals with typical errors near 1.5 eV. A higher exact exchange amount of 50% in the DFT starting point of one-shot GW does not provide very good agreement with experiment by overbinding orbitals with typical errors near 0.5 eV. While not as accurate as G 0 W 0 @PBE0, the G -only eigenvalue self-consistent GW scheme with W fixed to the PBE level provides a reasonably predictive level of theory (typical errors near 0.3 eV) to describe photoelectron spectra of these 3 d -transition metal dioxide anions. Adding eigenvalue self-consistency also in W , on the other hand, worsens the agreement with experiment overall. Overall, our findings on the performance of various GW methods are discussed in the context of our previous studies on other transition metal oxide molecular systems.

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

Comparison of first principles and semi-empirical models of the structural and electronic properties of Ge{sub 1−x}Sn{sub x} alloys

We present and compare three distinct atomistic models—based on first principles and semi-empirical approaches—of the structural and electronic properties of Ge{sub 1−x}Sn{sub x} alloys. Density functional theory calculations incorporating Heyd–Scuseria–Ernzerhof (HSE), local density approximation (LDA) and modified Becke–Johnson (mBJ) exchange-correlation functionals are used to perform structural relaxation and electronic structure calculations for a series of Ge{sub 1−x}Sn{sub x} alloy supercells. Based on HSE calculations, a semi-empirical valence force field (VFF) potential and sp{sup 3}s{sup ∗} tight-binding (TB) Hamiltonian are parametrised. Comparing the HSE, LDA+mBJ and VFF+TB models, and using the HSE results as a benchmark, we demonstrate that: (1) LDA+mBJ calculations provide an accurate first principles description of the electronic structure at reduced computational cost, (2) the VFF potential is sufficiently accurate to circumvent the requirement to perform first principles structural relaxation, and (3) VFF+TB calculations provide a good quantitative description of the alloy electronic structure in the vicinity of the band edges. Our results also emphasise the importance of Sn-induced band mixing in determining the nature of the conduction band structure of Ge{sub 1−x}Sn{sub x} alloys. The theoretical models and benchmark calculations we present inform and enable predictive, computationally efficient and scalable atomistic calculations for disordered alloys and nanostructures. This provides a suitable platform to underpin further theoretical investigations of the properties of this emerging semiconductor alloy.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Even Faster Exact Exchange for Solids via Tensor Hypercontraction

Hybrid density functional theory (DFT) remains intractable for large periodic systems due to the demanding computational cost of exact exchange. We apply the tensor hypercontraction (THC) (or interpolative separable density fitting) approximation to periodic hybrid DFT calculations with Gaussian-type orbitals using the Gaussian plane wave approach. This is done to lower the computational scaling with respect to the number of basis functions (N) and k-points (N k ) at a fixed system size. Additionally, we propose an algorithm to fit only occupied orbital products via THC (i.e., a set of points, N ISDF ) to further reduce computation time and memory usage. This algorithm has linear scaling cost with k-points, no explicit dependence of N ISDF on basis set size, and overall cubic scaling with unit cell size. Significant speedups and reduced memory usage may be obtained for moderately sized k-point meshes, with additional gains for large k-point meshes. Adequate accuracy can be obtained using THC-oo-K for self-consistent calculations. In conclusion, we perform illustrative hybrid density function theory calculations on the benzene crystal in the basis set and thermodynamic limits to highlight the utility of this algorithm.

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

Pt-Assisted Carbon Remediation of Mo2C Materials for CO Disproportionation

Using the CO disproportionation (Boudouard) reaction as a probe reaction, an in-depth analysis of temperature-programmed pulse response data shows that the addition of Pt to Mo2C mitigates deactivation of Mo active sites by acting as a carbon collector. CO2 production on Mo2C and Pt/Mo2C materials is dependent on both the activation energy and the CO surface concentration. Detailed plane-wave density functional theory calculations of the CO adsorption and disproportion reactions on Mo2C-supported Pt nanoparticles (NPs) are reported. The Mo2C was modeled by the ß-Mo2C (100) surface, and the Pt/Mo2C interface was modeled by the addition of 12 Pt atoms to the Mo2C (100) surface (12Pt@Mo2C). The potential energy surfaces of the Boudouard reaction were calculated on pure Mo2C, 12Pt@Mo2C, and Pt (111) surfaces. CO dissociation readily occurs on the Mo2C (100) surface, but not on the Pt (111) surface, with the former being exothermic and the latter being endothermic. At the Pt/Mo2C interface, CO dissociation is still exothermic, but with a larger energy barrier. The Boudouard reaction takes place on the Mo2C region, where CO2 is formed from a surface O atom dissociated from one CO molecule in reaction with another CO molecule, leaving one C atom on the surface. C adsorption is preferential on the Pt site in comparison to the Mo site. The supported Pt domains can collect the remaining C atoms, facilitating further CO2 formation on the active Mo sites. A Bader charge analysis shows that the surface metal-carbon bond is a mixture of covalent and ionic bonds, whereas the surface metal-oxygen bond is ionic. Electron localization function (ELF) and partial charge density calculations agree well with the Bader charge analysis. These computational results are consistent with experimental observations of the interaction of CO with Mo2C nanotube supported Pt domains in the transient regime under far from equilibrium conditions. The Boudouard reaction is an important side reaction, and the unexpected role found for Pt as a carbon collector, with Mo serving as a disproportionation site, provides a unique vantage point for understanding carbon and coke formation on catalytic materials.

Pt-assisted carbon remediation, Pt/Mo2C interface,↗
DOE OSTI2015

Orbital-free extension to Kohn-Sham density functional theory equation of state calculations: Application to silicon dioxide

The liquid regime equation of state of silicon dioxide SiO 2 is calculated via quantum molecular dynamics in the density range of 5 to 15 g/cc and with temperatures from 0.5 to 100 eV, including the α-quartz and stishovite phase Hugoniot curves. Below 8 eV calculations are based on Kohn-Sham density functional theory (DFT), and above 8 eV a new orbital-free DFT formulation, presented here, based on matching Kohn-Sham DFT calculations is employed. Recent experimental shock data are found to be in very good agreement with the current results. Finally both experimental and simulation data are used in constructing a new liquid regime equation of state table for SiO 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2018

Intra-molecular Charge Transfer and Electron Delocalization in Non-fullerene Organic Solar Cells

Two types of electron acceptors were synthesized by coupling two kinds of electron-rich cores with four equivalent perylene diimides (PDIs) at the a position. With fully aromatic cores, TPB and TPSe have pi-orbitals spread continuously over the whole aromatic conjugated backbone, unlike TPC and TPSi, which contain isolated PDI units due to the use of a tetrahedron carbon or silicon linker. Density functional theory calculations of the projected density of states showed that the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) for TPB are localized in separate regions of space. Further, the LUMO of TPB shows a greater contribution from the orbitals belonging to the connective core of the molecules than that of TPC. Overall, the properties of the HOMO and LUMO point at increased intra-molecular delocalization of negative charge carriers for TPB and TPSe than for TPC and TPSi and hence at a more facile intra-molecular charge transfer for the former. The film absorption and emission spectra showed evidences for the inter -molecular electron delocalization in TPB and TPSe, which is consistent with the network structure revealed by X-ray diffraction studies on single crystals of TPB. These features benefit the formation of charge transfer states and/or facilitate charge transport. Thus, higher electron mobility and higher charge dissociation probabilities under J sc condition were observed in blend films of TPB:PTB7-Th and TPSe:PTB7-Th than those in TPC:PTB7Th and TPSi:PTB7-Th blend films. As a result, the J sc and fill factor values of 15.02 mA/cm 2 , 0.58 and 14.36 mA/cm 2 , 0.55 for TPB- and TPSe-based solar cell are observed, whereas those for TPC and TPSi are 11.55 mA/cm 2 , 0.47 and 10.35 mA/cm 2 , 0.42, respectively.

14 SOLAR ENERGY↗
DOE OSTI2020

Computational and Experimental Mechanistic Insights into the Ethanol-to-Butanol Upgrading Reaction over MgO

The mechanism of ethanol upgrading to higher products is still under debate, especially regarding intermediate species and hydrogenation and dehydrogenation steps. In this work, we conducted a combined theoretical and experimental approach to contribute to this discussion. For such, detailed electronic structure density functional theory calculations (aiming at probing density of states, infrared spectra, geometric parameters, charge densities, and reaction energetics) and diffuse reflectance infrared Fourier transform spectroscopy experiments were carried out revealing the relevance of an appropriate combination of reactive surface sites to support the formation of several intermediates that are formed in the C-C coupling over MgO. The roles of Mg and O sites were also studied under an electronic perspective and different geometrical arrangements. We found that a kink configuration was the most adequate for ethanol to 1-butanol upgrading. Our calculations also gave us arguments to propose distinct reaction routes, whose mutual predominance would depend upon reaction temperature. At temperatures up to 573 K, the so-called β-route, which goes through scission of a Cβ-H bond and formation of an oxametallacycle-like intermediate, would dominate the coupling, whereas at higher temperatures, up to 673 K, a more usual Guerbet mechanism, via an aldol coupling step and then consecutive hydrogenations, would be expected. Here, the theoretical conclusions were followed by a careful experimental strategy using sequential experimental planning techniques in order to estimate accurate parameters with the lowest possible experimental load. Information from these different sources were coupled to develop a mathematical model for the rate of the ethanol upgrading reaction, using a Langmuir-Hinshelwood-Hougen-Watson approach. The developed and statistically validated model adequately described the experimental data at 673 K and 1.1 bar total pressure for ethanol partial pressures in the range from 0 to 20 kPa.

09 BIOMASS FUELS↗
DOE OSTI2024

Regioselective On-Surface Synthesis of [3]Triangulene Graphene Nanoribbons

The integration of low-energy states into bottom-up engineered graphene nanoribbons (GNRs) is a robust strategy for realizing materials with tailored electronic band structure for nanoelectronics. Low-energy zero-modes (ZMs) can be introduced into nanographenes (NGs) by creating an imbalance between the two sublattices of graphene. This phenomenon is exemplified by the family of [n]triangulenes (n ϵ $\mathbb{N}$). Here, we demonstrate the synthesis of [3]triangulene-GNRs, a regioregular one-dimensional (1D) chain of [3]triangulenes linked by five-membered rings. Hybridization between ZMs on adjacent [3]triangulenes leads to the emergence of a narrow band gap, E g,exp ~ 0.7 eV, and topological end states that are experimentally verified using scanning tunneling spectroscopy. Tight-binding and first-principles density functional theory calculations within the local density approximation corroborate our experimental observations. Our synthetic design takes advantage of a selective on-surface head-to-tail coupling of monomer building blocks enabling the regioselective synthesis of [3]triangulene-GNRs. Detailed ab initio theory provides insights into the mechanism of on-surface radical polymerization, revealing the pivotal role of Au-C bond formation/breakage in driving selectivity.

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

Experimental determination of the valence band offsets of ZnGeN 2 and (ZnGe) 0.94 Ga 0.12 N 2 with GaN

A predicted type-II staggered band alignment with an approximately 1.4 eV valence band offset at the ZnGeN 2 /GaN heterointerface has inspired novel band-engineered III-N/ZnGeN 2 heterostructure-based device designs for applications in high performance optoelectronics. We report on the determination of the valence band offset between metalorganic chemical vapor deposition grown (ZnGe) 1- x Ga 2 x N 2 , for x = 0 and 0.06, and GaN using x-ray photoemission spectroscopy. The valence band of ZnGeN 2 was found to lie 1.45–1.65 eV above that of GaN. This result agrees well with the value predicted by first-principles density functional theory calculations using the local density approximation for the potential profile and quasiparticle self-consistent GW calculations of the band edge states relative to the potential. For (ZnGe) 0.94 Ga 0.12 N 2 the value was determined to be 1.29 eV, ~10%–20% lower than that of ZnGeN 2 . The experimental determination of the large band offset between ZnGeN 2 and GaN provides promising alternative solutions to address challenges faced with pure III-nitride-based structures and devices.

Physics↗
DOE OSTI2019

Superconductivity in Y 7 R u 4 InG e 12

We report a type-II intermetallic superconductor Y 7 Ru 4 InGe 12 with a transition temperature ( T c ) of ~5.8 K, which is confirmed by zero resistivity, diamagnetic magnetic susceptibility, and specific-heat jump. Single crystals of Y 7 Ru 4 InGe 12 were grown from a reactive indium flux. Y 7 Ru 4 InGe 12 crystallizes in the tetragonal space group P 4/m and features a [Ru 4 InGe 12 ] polyanionic network with Y atoms located in three different channels. The upper critical fields of Y 7 Ru 4 InGe 12 at 0 K are determined to be ~5.3 and 2.4 T along the c axis and the ab plane, respectively. The estimated coherence length along the c axis (~174 Å) is much larger than the estimated mean free path along the c axis (~29 Å) in Y 7 Ru 4 InGe 12 , placing its superconductivity in the so-called dirty regime. The compound exhibits a large superconducting specific-heat jump ΔC/γe T c ≈ 2.4, significantly well above the weak-coupling Barden-Cooper-Schrieffer theoretical value of 1.43 and pointing to a strong-coupling scenario in Y 7 Ru 4 InGe 12 . As a result, density-functional-theory calculations show that the density of states in Y 7 Ru 4 InGe 12 exhibits a broad peak near the Fermi level which mainly derives from Y–4 d ,Ru–4 d , and Ge–4 p states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Equations of state and stability of MgSiO 3 perovskite and post-perovskite phases from quantum Monte Carlo simulations

In this study, we have performed quantum Monte Carlo (QMC) simulations and density functional theory calculations to study the equations of state of MgSiO 3 perovskite (Pv, bridgmanite) and post-perovskite (PPv) up to the pressure and temperature conditions of the base of Earth's lower mantle. The ground-state energies were derived using QMC simulations and the temperature-dependent Helmholtz free energies were calculated within the quasiharmonic approximation and density functional perturbation theory. The equations of state for both phases of MgSiO 3 agree well with experiments, and better than those from generalized gradient approximation calculations. The Pv-PPv phase boundary calculated from our QMC equations of state is also consistent with experiments, and better than previous local density approximation calculations. Lastly, we discuss the implications for double crossing of the Pv-PPv boundary in the Earth.

58 GEOSCIENCES↗
DOE OSTI2016

A multi-scale study of the adsorption of lanthanum on the (110) surface of tungsten

In this study, we utilize a multi-scale approach to studying lanthanum adsorption on the (110) plane of tungsten. The energy of the system is described from density functional theory calculations within the framework of the cluster expansion method. It is found that including two-body figures up to the sixth nearest neighbor yielded a reasonable agreement with density functional theory calculations as evidenced by the reported cross validation score. The results indicate that the interaction between the adsorbate atoms in the adlayer is important and cannot be ignored. The parameterized cluster expansion expression is used in a lattice gas Monte Carlo simulation in the grand canonical ensemble at 773 K and the adsorption isotherm is recorded. Implications of the obtained results for the pyroprocessing application are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Optimized Auxiliary Functions for Robust Mitigation of Finite-Size Errors in Periodic Hybrid Density Functional Theory

When calculating properties of periodic systems at the thermodynamic limit (TDL), the dominant source of finite size error (FSE) arises from the long-range Coulomb interaction, and can manifest as a slowly converging quadrature error when approximating an integral in the reciprocal space by a finite sum. The singularity subtraction (SS) method offers a systematic approach for reducing this quadrature error and thus the FSE. Here, in this work, we first investigate the performance of the SS method in the simplest setting, aiming at reducing the FSE in exact exchange calculations by subtracting the Coulomb contribution with a single, adjustable Gaussian auxiliary function. We demonstrate that a simple fitting method can robustly estimate the optimal Gaussian width and leads to rapid convergence toward the TDL. Furthermore, we suggest new forms of the auxiliary function, whose optimal parameters could also be determined through least-squares fitting. For a range of semiconductors and insulators, the proposed auxiliary functions achieve robust, millihartree-level accuracy in hybrid density functional theory calculations, including cases with sparse k-meshes and large basis sets.

Quiton, Stephen Jon [University of California, Ber↗
DOE OSTI2020

Rethinking the magnetic properties of lepidocrocite: A density functional theory and cluster expansion study

The iron oxyhydroxide lepidocrocite (g-FeOOH) is an abundant mineral critical to a number of chemical and technological applications. Of particular interest is the ground state and nite temperature magnetic order, and the subsequent impact this has upon crystal properties. The magnetic properties, investigated in this work are governed primarily through superexchange interactions, and have been calculated using density functional theory and cluster expansion methods. Quantication of these exchange terms has facilitated the determination of the ground state magneto-crystalline structure and subsequent calculation of its lattice constants, elastic moduli, cohesive enthalpy, and electronic density of states. Further, using a collinear magnetic conguration model, the magnetic heat capacity versus temperature has been studied and the Neel temperature obtained.

Pope, Daniel J.↗
DOE OSTI2015

Manganese-Based Molecular Electrocatalysts for Oxidation of Hydrogen

Oxidation of H2 (1 atm) is catalyzed by the manganese electrocatalysts [(P2N2)MnI(CO)(bppm)]+ and [(PNP)MnI(CO)(bppm)]+ (P2N2= 1,5-dibenzyl-3,7-diphenyl-1,5-diaza-3,7-diphosphacyclooctane; PNP = (Ph2PCH2)2NMe); bppm = (PArF2)2CH2, and ArF = 3,5-(CF3)2C6H3). In fluorobenzene solvent using 2,6-lutidine as the exogeneous base, the turnover frequency for [(P2N2)MnI(CO)(bppm)]+ is 3.5 s-1 with an estimated overpotential of 590 mV. For [(PNP)MnI(CO)(bppm)], in fluorobenzene solvent using N-methylpyrrolidine as the exogeneous base, the turnover frequency is 1.4 s-1 with an estimated overpotential of 700 mV. Density functional theory calculations suggest that the slow step in the catalytic cycle is proton transfer from the oxidized 17-electron manganese hydride, e.g., [(P2N2)MnIIH(CO)(bppm)]+ to the pendant amine. The computed activation barrier for intramolecular proton transfer from the metal to the pendant amine is 20.4 kcal/mol in [(P2N2)MnIIH(CO)(bppm)]+ and 21.3 kcal/mol in [(PNP)MnI(CO)(bppm)]. The high barrier appears to result from both the unfavorability of metal-to-nitrogen proton transfer (thermodynamically uphill by 6.6 pKa units, 9 kcal/mol), as well as the relatively long manganese-nitrogen separation in the MnIIH complexes.

Hydrogen, oxidation, electrocatalysis, proton tran↗
DOE OSTI2019

Tuning Optical Properties of Conjugated Molecules by Lewis Acids: Insights from Electronic Structure Modeling

Understanding and controlling the optoelectronic properties of organic semiconductors at the molecular level remains a challenge due to the complexity of chemical structures and intermolecular interactions. A common strategy to address this challenge is to utilize both experimental and computational approaches. In this contribution, we show that density functional theory (DFT) calculation is a useful tool to provide insights into the bonding, electron population distribution, and optical transitions of adducts between conjugated molecules and Lewis acids (CM-LA). Adduct formation leads to relevant modifications of key properties, including a red shift in optical transitions and an increase in charge carrier density and charge mobility, compared to the parent conjugated molecules. We show that electron density transfer from the CM to the LA, which was hypothesized to cause the experimental red shift in absorption spectra upon LA binding, can be quantified and interpreted by population analysis. Experimental red shifts in optical transitions for all molecular families can also be predicted by time dependent DFT calculations with different density functionals. Finally, these detailed insights help to optimize a priori design guidelines for future applications.

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

Electronic and magnetic properties of epitaxial perovskite SrCrO3(001)

We have investigated the intrinsic properties of SrCrO3 epitaxial thin films synthesized by molecular beam epitaxy. We find compelling evidence that SrCrO3 is a correlated metal. X-ray photoemission valence band and O K-edge x-ray absorption spectra indicate a strongly hybridized Cr3d-O2p state crossing the Fermi level, leading to metallic behavior. Comparison between valence band spectra near the Fermi level and the densities of states calculated using density functional theory (DFT) also suggests the presence of coherent and incoherent states and points to a strong electron-electron correlation effects. The magnetic susceptibility can be described by Pauli paramagnetism at temperatures above 100 K, but reveals antiferromagnetic behavior at lower temperatures resulting from orbital ordering as suggested by Ortega-San-Martin et al. [Phys. Rev. Lett. 99, 255701 (2007)].

Environmental Molecular Sciences Laboratory↗