Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DFT Calculations”

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 1,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Exploring two-dimensional van der Waals heavy-fermion material: Data mining theoretical approach

Abstract The discovery of two-dimensional (2D) van der Waals (vdW) materials often provides interesting playgrounds to explore novel phenomena. One of the missing components in 2D vdW materials is the intrinsic heavy-fermion systems, which can provide an additional degree of freedom to study quantum critical point (QCP), unconventional superconductivity, and emergent phenomena in vdW heterostructures. Here, we investigate 2D vdW heavy-fermion candidates through the database of experimentally known compounds based on dynamical mean-field theory calculation combined with density functional theory (DFT+DMFT). We have found that the Kondo resonance state of CeSiI does not change upon exfoliation and can be easily controlled by strain and surface doping. Our result indicates that CeSiI is an ideal 2D vdW heavy-fermion material and the quantum critical point can be identified by external perturbations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Mechanistic insight of KBiQ 2 (Q = S, Se) using panoramic synthesis towards synthesis-by-design

Solid-state synthesis has historically focused on reactants and end products; however, knowledge of reaction pathways, intermediate phases and their formation may provide mechanistic insight of solid-state reactions. With an increased understanding of reaction progressions, design principles can be deduced, affording more predictive power in materials synthesis. In pursuit of this goal, in situ powder X-ray diffraction is employed to observe crystalline phase evolution over the course of the reaction, thereby constructing a “panoramic” view of the reaction from beginning to end. We conducted in situ diffraction studies in the K–Bi–Q (Q = S, Se) system to understand the formation of phases occurring in this system in the course of their reactions. Powder mixtures of K 2 Q to Bi 2 Q 3 in 1:1 and 1.5:1 ratios were heated to 800 °C or 650 °C, while simultaneously collecting diffraction data. Three new phases, K 3 BiS 3 , β-KBiS 2 , and β-KBiSe 2 , were discovered. Panoramic synthesis showed that K 3 BiQ 3 serves an important mechanistic role as a structural intermediate in both chalcogen systems (Q = S, Se) in the path to form the KBiQ 2 structure. Thermal analysis and calculations at the density functional theory (DFT) level show that the cation-ordered β-KBiQ 2 polymorphs are the thermodynamically stable phase in this compositional space, while Pair Distribution Function (PDF) analysis shows that all α-KBiQ 2 structures have local disorder due to stereochemically active lone pair expression of the bismuth atoms. The formation of the β-KBiQ 2 structures, both of which crystallize in the α-NaFeO 2 structure type, show a boundary where the structure can be disordered or ordered with regards to the alkali metal and bismuth. A cation radius tolerance for six-coordinate cation site sharing of ($\frac{r^+}{r^{3+}}$) ~ 1.3 is proposed. The mechanistic insight the panoramic synthesis technique provides in the K–Bi–Q system is progress towards the overarching goal of synthesis-by-design.

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

Diffusion Monte Carlo evaluation of disiloxane linearisation barrier

The disiloxane molecule is a prime example of silicate compounds containing the Si-O-Si bridge. The molecule is of significant interest within the field of quantum chemistry, owing to the difficulty in theoretically predicting its properties. Herein, the linearisation barrier of disiloxane is investigated using a fixed-node diffusion Monte Carlo (FNDMC) approach, which is one of the most reliable ab initio methods in accounting for the electronic correlation. Calculations utilizing the density functional theory (DFT) and the coupled cluster method with single and double substitutions, including noniterative triples (CCSD(T)) are carried out alongside FNDMC for comparison. It is concluded that FNDMC successfully predicts the disiloxane linearisation barrier and does not depend on the completeness of the basis-set as much as DFT or CCSD(T), thus establishing its suitability.

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

Alloying multiple halide perovskites on the same sublattice in search of stability and target band gaps

Single-component halide perovskites (HPs) rarely satisfy all the necessary criteria for optoelectronic applications, such as achieving an optimal band gap while maintaining high chemical and structural stability. Alloying halide perovskites has emerged as a promising strategy, not only to enhance stability but also to fine-tune their electronic and optical properties. In this work, we explore multiple degrees of freedom in alloy design, considering different substitution sublattices sites (A, B, or X in ABX3 perovskites), various chemical species (isovalent and hetero-valent elements), and multi-component compositions on a given sublattice. Using first-principles calculations based on density functional theory (DFT), we investigate how compositional variations influence the electronic (band gap) and structural properties (mixing enthalpy) of HP alloys. Our approach employs the polymorphous cell model, allowing full local relaxation which breaks local symmetry while preserving global cubic symmetry—an essential framework for accurately modeling HPs. Our results reveal that X-site mixing (halogen substitution) primarily affects the valence band maximum, allowing target band gap engineering. Additionally, variations in halogen radii introduce internal strain through octahedral distortions, influencing the mixing enthalpy. A-site substitution, while not directly contributing to the band edge states, modifies structural stability via volume effects, indirectly impacting the band gap. B-site alloying plays a dominant role in band gap modulation, leading to either positive or negative band gap bowing. Specifically, isovalent B-site mixing (Sn–Pb) induces strong positive bowing, where the alloy band gap is smaller than the average gap of parent compounds, whereas hetero-valent mixing (Cd–Pb) results in pronounced negative bowing. As an aside, we investigate the competition between the excess energy of disordered alloys vs. that of long-range ordered double perovskites of the same compositions, seeking examples of ordered phases emerging from disordered alloys. Furthermore, our findings provide fundamental insights into the electronic and structural behavior of HP alloys, offering valuable design principles for the development of stable and efficient materials for next-generation photovoltaic and optoelectronic devices.

14 SOLAR ENERGY↗
DOE OSTI2023

Metal–insulator transition in composition-tuned nickel oxide films

Abstract Thin films of the solid solution Nd 1 − x La x NiO 3 are grown in order to study the expected 0 K phase transitions at a specific composition. We experimentally map out the structural, electronic and magnetic properties as a function of x and a discontinuous, possibly first order, insulator–metal transition is observed at low temperature when x = 0.2. Raman spectroscopy and scanning transmission electron microscopy show that this is not associated with a correspondingly discontinuous global structural change. On the other hand, results from density functional theory (DFT) and combined DFT and dynamical mean field theory calculations produce a 0 K first order transition at around this composition. We further estimate the temperature-dependence of the transition from thermodynamic considerations and find that a discontinuous insulator–metal transition can be reproduced theoretically and implies a narrow insulator–metal phase coexistence with x . Finally, muon spin rotation ( µ SR) measurements suggest that there are non-static magnetic moments in the system that may be understood in the context of the first order nature of the 0 K transition and its associated phase coexistence regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Revealing symmetry-broken superconducting configurations by density functional theory

A coherent theory for the superconductivity of both conventional and unconventional superconductors is currently lacking. Here we show that superconductivity arises from the formation of a symmetry-broken superconducting configuration (SCC) due to atomic perturbation of the normal conducting configuration (NCC). This electron–phonon interaction creates straight one-dimensional tunnels (SODTs) for charge density of electrons and/or holes as revealed by the calculations based on density functional theory (DFT). The SODTs act as resistance-free superhighways and are correlated to the Cooper pairs in the Bardeen–Cooper–Schrieffer (BCS) theory. The formation of SODTs implies that the electron–phonon interaction in the BCS theory can be represented by the difference in charge densities between SCC and NCC predicted by DFT. The present work highlights that in conventional superconductors, SODTs are embedded within the bulk materials and are easily destroyed by phonon vibrations, resulting in a low critical superconducting temperature (T C ). Conversely, in unconventional superconductors such as YBa 2 Cu 3 O 7 (YBCO 7 ), SODTs are protected by a layered pontoon structure with very weak bonding to the bulk materials, maintaining SODTs’ stability at higher temperatures and leading to a much higher T C . The present approach is validated for 14 conventional superconductors of 18 pure elements and MgB 2 examined in this work, including the presently predicted superconductivity in Cu, Ag, Au, Sb, Bi, and MgB 2 at 0 K and 0 GPa, and one unconventional superconductor of YBCO 7 . Our discovery indicates that DFT can be a practical tool for predicting superconductors, enabling a systematic search for new superconducting materials in the future.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Charged point defect benchmark of Hybrid and GGA-PBE

Data used by the publication "High-throughput calculations of charged point defect properties with semi-local density functional theory - performance benchmarks for materials screening applications." This work presented an in-depth benchmark analysis of automated, semi-local point defect calculations with a-posteriori corrections, compared to 245 “gold standard” hybrid calculations previously published. We considered three different a-posteriori correction sets for semi-local calculations, implemented in a fully automated workflow, and consider the qualitative and quantitative differences for four different categories of defect information: thermodynamic transition levels, formation energies, fermi levels, and dopability limits. We highlighted the type of qualitative information about point defect properties that can be extracted from high-throughput calculations based on semi-local DFT methods, while also demonstrating the limits of quantitative accuracy that can be achieved by these approaches.

Broberg, Danny↗
DOE OSTI2024

Sulfur-Doped Carbon Support Boosts CO2RR Activity of Ag Electrocatalysts

For presentation at the 70th AVS International Symposium and Exhibition. In this work, we show that the activity of Ag electrocatalysts for electrochemical CO2 to CO conversion is improved when supported on sulfur-doped (S-doped) carbon materials. S-doped carbon support was created by treating the heavily sputtered, highly oriented pyrolytic graphite (HOPG) in H2S at elevated temperatures, as confirmed by the S 2p X-ray photoelectron spectroscopy (XPS) peak. Scanning tunneling microscopy (STM) images indicated that Ag nanoparticles supported on S-doped HOPG had similar size distributions as those supported on sulfur-free (S-free) HOPG. While both catalysts reached > 90% CO Faradaic efficiency (FECO) at E = -1.3 V vs. the reversible hydrogen electrode (RHE) in the CO2 reduction reaction (CO2RR), Ag catalysts supported on S-doped HOPG demonstrated 70% higher CO turnover frequency (TOFCO = 3.4 CO/atomAg/s) than those supported on S-free HOPG (TOFCO = 2.0 CO/atomAg/s). Preliminary calculations based on density functional theory (DFT) indicated a more favorable energetic pathway of CO2-to-CO at the C-S-Ag interface, tentatively consistent with experiments. These results hint at a new approach to design active and selective electrocatalysts for CO2 conversion.

Deng, Xingyi↗
NASA NTRS2018

Aladyn - Adaptive Neural Network Molecular Dynamics Simulation Code: Computational Materials Mini-Application

This report provides an overview and commands description of the Computational Materials mini-application, Aladyn. Aladyn is a simple molecular dynamics code written in FORTRAN 2008, which is designed to demonstrate the use of adaptive neural networks (ANNs) in atomistic simulations. The role of ANNs is to reproduce the very complex energy landscape resulting from the atomic interactions in materials with the accuracy of quantum mechanics-based energy calculations. The ANN is trained on a large set of atomic structures calculated using the density functional theory (DFT) method. The Aladyn code is being released to serve as a training testbed for students and professors in academia to explore possible optimization algorithms for parallel computing on multicore central processing unit (CPU) computers or computers utilizing many core architectures based on graphic processing units (GPUs). The effort is related to the High Performance Computing Incubator (HPCI) project at NASA Langley Research Center.

Yamakov, Vesselin I.↗
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 OSTI2022

Organic Photovoltaic Efficiency Predictor: Data-Driven Models for Non-Fullerene Acceptor Organic Solar Cells

In the design of organic solar cells, there has been a need for materials with high power conversion efficiencies. Scharber’s model is commonly used to predict efficiency; however, it exhibits poor performance with new non-fullerene acceptor (NFA) devices, since it was designed for fullerene-based devices. In this work, an empirical model is proposed that can be a more accurate alternative for NFA organic solar cells. Additionally, many screening studies use computationally expensive methods. A model based on using semiempirical simplified time-dependent density functional theory (sTD-DFT) as an alternative method can accelerate the calculations and yield a similar accuracy. The models presented in this paper, termed organic photovoltaic efficiency predictor (OPEP) models, have shown significantly lower errors than previous models, with OPEP/B3LYP yielding errors of 1.53% and OPEP/sTD- DFT of 1.55%. As a result, the proposed computational models can be used for the fast and accurate screening of new high-efficiency NFAs/donor pairs.

14 SOLAR ENERGY↗
DOE OSTI2022

Predicting band gaps and band-edge positions of oxide perovskites using density functional theory and machine learning

Density functional theory (DFT) within the local or semilocal density approximations, i.e., the local density approximation (LDA) or generalized gradient approximation (GGA), has become a workhorse in the electronic structure theory of solids, being extremely fast and reliable for energetics and structural properties, yet remaining highly inaccurate for predicting band gaps of semiconductors and insulators. The accurate prediction of band gaps using first-principles methods is time consuming, requiring hybrid functionals, quasiparticle GW, or quantum Monte Carlo methods. Efficiently correcting DFT-LDA/GGA band gaps and unveiling the main chemical and structural factors involved in this correction is desirable for discovering novel materials in high-throughput calculations. In this direction, we, in this study, use DFT and machine learning techniques to correct band gaps and band-edge positions of a representative subset of ABO 3 perovskite oxides. Relying on the results of HSE06 hybrid functional calculations as target values of band gaps, we find a systematic band-gap correction of ~1.5 eV for this class of materials, where ~1 eV comes from downward shifting the valence band and ~0.5 eV from uplifting the conduction band. The main chemical and structural factors determining the band-gap correction are determined through a feature selection procedure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Wannier–Koopmans method calculations for transition metal oxide band gaps

The widely used density functional theory (DFT) has a major drawback of underestimating the band gaps of materials. Wannier–Koopmans method (WKM) was recently developed for band gap calculations with accuracy on a par with more complicated methods. WKM has been tested for main group covalent semiconductors, alkali halides, 2D materials, and organic crystals. Here we apply the WKM to another interesting type of material system: the transition metal (TM) oxides. TM oxides can be classified as either with d 0 or d 10 closed shell occupancy or partially occupied open shell configuration, and the latter is known to be strongly correlated Mott insulators. We found that, while WKM provides adequate band gaps for the d0 and d10 TM oxides, it fails to provide correct band gaps for the group with partially occupied d states. This issue is also found in other mean-field approaches like the GW calculations. We believe that the problem comes from a strong interaction between the occupied and unoccupied d-state Wannier functions in a partially occupied d-state system. We also found that, for pseudopotential calculations including deep core levels, it is necessary to remove the electron densities of these deep core levels in the Hartree and exchange–correlation energy functional when calculating the WKM correction parameters for the d-state Wannier functions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Towards Elucidating Structure–Spectra Relationships in Rhamnogalacturonan II: Computational Protocols for Accurate 13C and 1H Shifts for Apiose and Its Borate Esters

Apiose is a naturally occurring, uncommon branched-chain pentose found in plant cell walls as part of the complex polysaccharide Rhamnogalacturonan II (RG-II). The structural elucidation of the three-dimensional structure of RG-II by nuclear magnetic resonance (NMR) spectroscopy is significantly complicated by the ability of apiose to cross-link via borate ester linkages to form RG-II dimers. Here, we developed a computational approach to gain insight into the structure–spectra relationships of apio–borate complexes in an effort to complement experimental assignments of NMR signals in RG-II. Our protocol involved structure optimizations using density functional theory (DFT) followed by isotropic magnetic shielding constant calculations using the gauge-invariant atomic orbital (GIAO) approach to predict chemical shifts. We evaluated the accuracy of 23 different functional–basis set (FBS) combinations with and without implicit solvation for predicting the experimental 1 H and 13 C shifts of a methyl apioside and its three borate derivatives. The computed NMR predictions were evaluated on the basis of the overall shift accuracy, relative shift ordering, and the ability to distinguish between dimers and monomers. We demonstrate that the consideration of implicit solvation during geometry optimizations in addition to the magnetic shielding constant calculations greatly increases the accuracy of NMR chemical shift predictions and can correctly reproduce the ordering of the 13 C shifts and yield predictions that are, on average, within 1.50 ppm for 13 C and 0.12 ppm for 1 H shifts for apio–borate compounds.

BASIC BIOLOGICAL SCIENCES,BIOMASS FUELS↗
DOE OSTI2020

Robust, Accurate, and Efficient: Quantum Embedding Using the Huzinaga Level-Shift Projection Operator for Complex Systems

Using wave function (WF) in density functional theory (DFT) embedding methods provides a framework for performing localized, high-accuracy WF calculations on a system, while not incurring the full computational cost of the WF calculation on the full system. To effectively partition a system into localized WF and DFT subsystems, we utilize the Huzinaga level-shift projection operator within an absolutely localized basis. In this work, we study the ability of the absolutely localized Huzinaga level-shift projection operator method to study complex WF and DFT partitions, including partitions between multiple covalent bonds, a double bond, and transition-metal–ligand bonds. We find that our methodology can accurately describe all of these complex partitions. Additionally, we study the robustness of this method with respect to the WF method, specifically where the embedded systems were described using a multiconfigurational WF method. We found that the method is systematically improvable with respect to both the number of atoms in the WF region and the size of the basis set used, with energy errors less than 1 kcal/mol. Additionally, we calculated the adsorption energy of H2 to a model of an iron metal–organic framework (Fe-MOF-74) to within 1 kcal/mol compared to CASPT2 calculations performed on the full model while incurring only a small fraction of the full computational cost. In conclusion, this work demonstrates that the absolutely localized Huzinaga level-shift projection operator method is applicable to very complex systems with difficult electronic structures.

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

Using ground state and excited state density functional theory to decipher 3d dopant defects in GaN

Abstract Using ground state density functional theory (DFT) and implementing an occupation-constrained DFT (occ-DFT) for self-consistent excited state calculations, we decipher the electronic structure of the Mn dopant and other 3 d defects in GaN across the band gap. Our analysis, validated with broad agreement with defect levels (ground-state calculations) and photoluminescence data (excited-state calculations), mandates reinterpretation and reassignment of 3 d defect data in GaN. The Mn Ga defect is determined to span stable charge states from (1−) in n -type GaN through (2+) in p -type GaN. The Mn(2+) is predicted to be a d 2 ground state spin triplet defect with a singlet excited state, isoelectronic with the defect associated with the 1.19 eV photoluminescence in n -type GaN. The combined analysis of defect levels and excited states invites reassessment of all d 2 -capable dopants in GaN. We demonstrate that the 1.19 eV defect, a candidate defect for optically controlled quantum applications, cannot be the Cr(1+) assumed in literature and instead must be the V(0). The combined ground-state/excited-state DFT analysis is shown to be able to chemically fingerprint defects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Adsorption and activation of CO 2 on Pt/CeOx/TiO 2 (110): Role of the Pt-CeO x interface

The adsorption and dissociation of CO 2 on TiO 2 (110), CeO x /TiO 2 (110) and Pt/CeO x /TiO 2 (110) surfaces has been examined using Ambient Pressure X-ray Photoelectron Spectroscopy (AP-XPS). The substrates under study exhibited different degrees of complexity which were tested for the binding of the adsorbate and the cleavage of C-O bonds. The surfaces were prepared by depositing CeO x (0.1 ML) onto TiO 2 (110) to form a mixed oxide support, onto which Pt nanoparticles (0.2 ML) were deposited. This configuration yields a complex set of interfaces between metal and oxides and we have systematically titrated the active role of each component (Pt 4f, Ce 3d and Ti 2p regions) and the arising surface intermediates (C 1s and O 1s regions). CO 2 barely bonds to stoichiometric TiO 2 (110). It heals oxygen vacancies of this oxide surface (CO 2,gas → CO gas + O a ) and does not form stable carbonates. A stable carbonate was seen upon adsorption of CO 2 on CeO x /TiO 2 (110) and on this type of substrate the adsorbate also removed O vacancies leading to the oxidation of Ti 3+ and Ce 3+ sites. Pt nanoparticles dispersed on CeO x /TiO 2 (110) were highly effective for the binding and dissociation of CO 2 , with the formation of CO 3 , CO, C and CH x species on the Pt/CeO x /TiO 2 (110) system. The results of theoretical calculations based on density-functional theory (DFT) show that Pt/CeOx/TiO2(110) binds CO2 much stronger than surfaces of bulk platinum {(111), (100), (110)} or other late transition metals. On a Pt-CeO x interface, the molecule adsorbs with a bent configuration (~ 130° O-C-O bond angle) and with a substantial elongation (~ 0.1 Å) of the C-O bonds, facilitating its transformation into high value chemicals.

36 MATERIALS SCIENCE↗
NASA NTRS1998

C2 Fragmentation Energy of C60 Revisited: Theory Disagrees with Most Experiments

Following our earlier work on the subject, we have carried out density functional theory (DFT) and second-order Moller-Plesset perturbation theory (MP2) calculations of the dissociation energy of the reaction C60 yields C58 + C2 using polarized basis sets and geometries optimized with DFT methods. The present theoretical results support an electronic fragmentation energy D(sub e) around 10-11 eV in disagreement with most experimental results that place the dissociation energy D(sub o) (including zero point energy) around 7-8 eV. The plausible errors remaining in the theoretical calculations are unlikely to account for this big difference (2-4 eV).

Boese, A. Daniel↗