Search NASA⌕ Search

SEARCH · Search NASA

Results for “DFT 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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Mechanistic Study of the Direct Hydrodeoxygenation of m -Cresol over WO x -Decorated Pt/C Catalysts

In this report hydrodeoxygenation (HDO) of m-cresol to produce toluene over carbon-supported Pt and Pt-WOx catalysts was studied. In stark contrast to Pt/C that exhibits only modest selectivity and low stability for this reaction, Pt-WOx/C was found to be unusually active and selective to toluene with greater than 94% selectivity to this product while exhibiting little to no deactivation under a wide range of reaction conditions. Reactivity studies in combination with density functional theory (DFT) calculations for the adsorption and reaction of m-cresol on structurally optimized WOx-decorated Pt(111) structures indicate that the HDO reaction on Pt-WOx/C proceeds via a direct hydrogenolysis of the C–O bond in m-cresol adsorbed on oxygen vacancy (or redox) sites on WOx species. The DFT results also indicate that Pt helps stabilize the WOx film while facilitating oxygen vacancy formation.

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

Oxygen Vacancies Alter Methanol Oxidation Pathways on NiOOH

A thorough comprehension of the mechanism underlying the methanol oxidation reaction (MOR) on Ni-based catalysts is critical for future electrocatalytic design and development. However, the mechanism of MOR on these materials remains a matter of controversy. Herein, we combine in situ surface-enhanced infrared absorption spectroscopy (SEIRAS) and density functional theory (DFT) calculations to identify the active sites and determine the mechanism of MOR on monometallic Ni-based catalysts in alkaline media. The SEIRAS results show that formate and (bi)carbonate are formed after the commencement of the MOR with potential-dependent relative distributions. These spectroscopic results are in good agreement with the DFT-computed reaction profiles over an oxygen vacancy, suggesting that the MOR mainly proceeds through the formate-involving pathway, in which the early consumption of methanol yields formate as the major product, while increasing potential drives further oxidation of formate to (bi)carbonate. We also find a parallel pathway for the generation of (bi)carbonate at high potentials that bypasses the formation of formate. The two main pathways are thermodynamically more feasible than the one predominantly reported in the literature for MOR on NiOOH that involves CHO and/or CO as key intermediates. These DFT results are supported by spectroscopic evidence showing that no band associated with CHO or CO can be detected by SEIRAS, which is attributed to the nature of the oxygen vacancies as the active sites, suppressing deep dehydrogenation of CH 2 O to CHO. Furthermore, this work thus shows the promising role of defect engineering in promoting the electrocatalytic MOR activity and selectivity.

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

Theoretical Prediction and Interpretation of 237 Np Mössbauer Isomer Shifts

Different theoretical approaches for the calculation of 237 Np Mössbauer isomer shifts are investigated. The traditional contact density route is compared to a previously proposed alternative approach that uses energy derivatives with respect to the nuclear radius. Both approaches yield similar results as long as suitable basis sets augmented with large exponents and relativistic methods are used. Density functional theory (DFT) calculations do not show a strong dependency of the 237 Np isomer shift on the chosen functional. Wavefunction calculations show that dynamic electron correlation can be important when covalent bonding influences the isomer shift. Effects from spin–orbit coupling are small. The isomer shifts of ionic solids and Np(III) organometallic complexes are largely governed by the oxidation state of Np. Isomer shifts of organometallic Np(IV) complexes are strongly affected by donation bonding. Furthermore, detailed analysis of the wavefunction results with different active spaces demonstrates that correlation among the outer core Np and occupied ligand frontier orbitals contributes significantly to isomer shifts of Np(IV) compounds.

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

Quasiperiodic potassium adlayer on decagonal Al–Ni–Co quasicrystal

Quasiperiodicity in free-electron-like metals is a subject of significant interest within the scientific community. In this work, utilizing scanning tunneling microscopy (STM), low energy electron diffraction (LEED), and density functional theory (DFT), we demonstrate the formation of a quasiperiodic potassium monolayer on the tenfold surface of decagonal Al–Ni–Co quasicrystal. A dispersed growth comprising of isolated K adatoms is observed at sub-monolayer coverage, which coalesce with increasing coverage and forms pentagonal and decagonal quasiperiodic motifs. LEED demonstrates distinct rings of spots displaying decagonal symmetry. Furthermore, our DFT calculation using the W-approximant surface that closely resembles d-Al–Ni–Co shows that sizable adsorbate-substrate interaction makes the potassium adatoms bind to the quasiperiodically dispersed favorable adsorption sites resulting in the formation of quasiperiodic potassium monolayer. Notably, the experimental motifs obtained from STM measurements align remarkably well with the DFT predictions, underscoring the intricate relationship between the electronic structure of the substrate that drives the quasiperiodic growth of the potassium adlayer.

alkali metal↗
DOE OSTI2026

Machine learning approach for vibronically renormalized electronic band structures

Here, we present a machine learning (ML) method for efficient computation of vibrational thermal expectation values of physical properties from first principles. Our approach is based on the nonperturbative frozen phonon formulation in which stochastic Monte Carlo algorithm is employed to sample configurations of nuclei in a supercell at finite temperatures based on a first-principles phonon model. A deep-learning neural network is trained to accurately predict physical properties associated with sampled phonon configurations, thus bypassing the time-consuming ab initio calculations. To incorporate the point-group symmetry of the electronic system into the ML model, group-theoretical methods are used to develop a symmetry-invariant descriptor for phonon configurations in the supercell. We apply our ML approach to compute the temperature dependent electronic energy gap of silicon based on density functional theory (DFT). We show that, with less than a hundred DFT calculations for training the neural network model, an order of magnitude larger number of sampling can be achieved for the computation of the vibrational thermal expectation values. Our work highlights the promising potential of ML techniques for finite temperature first-principles electronic structure methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Capturing multireference excited states by constrained-density-functional theory

The computation of excited electronic states with commonly employed (approximate) methods is challenging, typically yielding states of lower quality than the corresponding ground state for a higher computational cost. In this work, we present a mean-field method that extends the previously proposed excited constrained-density-functional theory (XCDFT) from single Slater determinants to ensemble one-body reduced density matrices for computing low-lying excited states. The method still retains an associated computational complexity comparable to a semilocal density functional theory (DFT) calculation while at the same time it is capable of approaching states with multireference character. Further, we benchmark the quality of this method on well-established test sets, finding good descriptions of the electronic structure of multireference states and maintaining an overall accuracy for the predicted excitation energies comparable to semilocal time-dependent DFT.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Density functional theory calculations of stability and diffusion mechanisms of impurity atoms in Ge crystals

Ge-based substrates are being developed for applications in advanced nano-electronic devices because of their higher intrinsic carrier mobility than Si. The stability and diffusion mechanism of impurity atoms in Ge are not well known in contrast to those of Si. Systematic studies of the stable sites of 2nd to 6th row element impurity atoms in Ge crystal were undertaken with density functional theory (DFT) and compared with those in Si crystal. It was found that most of the impurity atoms in Ge were stable at substitutional sites, while transition metals in Si were stable at interstitial sites and the other impurity atoms in Si were stable at substitutional sites. Furthermore, DFT calculations were carried out to clarify the mechanism responsible for the diffusion of impurity atoms in Ge crystals. The diffusion mechanism for 3d transition metals in Ge was found to be an interstitial-substitutional diffusion mechanism, while in Si this was an interstitial diffusion mechanism. The diffusion barriers in the proposed diffusion mechanisms in Ge and Si were quantitatively verified by comparing them to the experimental values in the literature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Synthesis, Spectra, and Theoretical Investigations of 1,3,5-Triazines Compounds as Ultraviolet Rays Absorber Based on Time-Dependent Density Functional Calculations and three-Dimensional Quantitative Structure-Property Relationship

A series of 1,3,5-triazines were synthesized and their UV absorption properties were tested. The computational chemistry methods were used to construct quantitative structure-property relationship (QSPR), which was used to computer aided design of new 1,3,5-triazines ultraviolet rays absorber compounds. The experimental UV absorption data are in good agreement with those predicted data using the Time-dependent density functional theory (TD-DFT) [B3LYP/6–311 + G(d,p)]. A suitable forecasting model (R > 0.8, P < 0.0001) was revealed. Predictive three-dimensional quantitative structure-property relationship (3D-QSPR) model was established using multifit molecular alignment rule of Sybyl program, which conclusion is consistent with the TD-DFT calculation. The exceptional photostability mechanism of such ultraviolet rays absorber compounds was studied and confirmed as principally banked upon their ability to undergo excited-state deactivation via an ultrafast excited-state proton transfer (ESIPT). The intramolecular hydrogen bond (IMHB) of 1,3,5-triazines compounds is the basis for the excited state proton transfer, which was explored by IR spectroscopy, UV spectra, structural and energetic aspects of different conformers and frontier molecular orbitals analysis.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Assessing the stability of Pd-exchanged sites in zeolites with the aid of a high throughput quantum chemistry workflow

Abstract Cation exchanged-zeolites are functional materials with a wide range of applications from catalysis to sorbents. They present a challenge for computational studies using density functional theory due to the numerous possible active sites. From Al configuration, to placement of extra framework cation(s), to potentially different oxidation states of the cation, accounting for all these possibilities is not trivial. To make the number of calculations more tractable, most studies focus on a few active sites. We attempt to go beyond these limitations by implementing a workflow for a high throughput screening, designed to systematize the problem and exhaustively search for feasible active sites. We use Pd-exchanged CHA and BEA to illustrate the approach. After conducting thousands of explicit DFT calculations, we identify the sites most favorable for the Pd cation and discuss the results in detail. The high throughput screening identifies many energetically favorable sites that are non-trivial. Lastly, we employ these results to examine NO adsorption in Pd-exchanged CHA, which is a promising passive NO x adsorbent (PNA) during the cold start of automobiles. The results shed light on critical active sites for NO x capture that were not previously studied.

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

Parameterizing empirical interatomic potentials for predicting thermophysical properties via an irreducible derivative approach: the case of ThO 2 and UO 2

The accuracy of classical physical property predictions using molecular dynamics simulations is determined by the quality of the interatomic potentials. Here we introduce a training approach for empirical interatomic potentials (EIPs) which is well suited for capturing phonons and phonon-related properties. Our approach is based on direct comparisons of the second- and third-order irreducible derivatives (IDs) between an EIP and the Born–Oppenheimer potential within density functional theory (DFT) calculations. IDs fully exploit space group symmetry and allow for training without redundant information. We demonstrate the fidelity of our approach in the context of ThO 2 and UO 2 , where we optimize parameters of an embedded-atom method potential in addition to core–shell interactions. Our EIPs provide thermophysical properties in good agreement with DFT and outperform widely utilized EIPs for phonon dispersion and thermal conductivity predictions. Reasonable estimates of thermal expansion and formation energies of Frenkel pairs are also obtained.

empirical interatomic potential↗
DOE OSTI2025

Toward machine learning interatomic potentials for modeling uranium mononitride

Uranium mononitride (UN) is a promising accident-tolerant fuel because of its high fissile density and high thermal conductivity. In this study, we developed the first machine learning interatomic potentials for reliable atomic-scale modeling of UN at finite temperatures. We constructed a training set using density functional theory (DFT) calculations that was enriched through an active learning procedure, and two neural network potentials were generated. Both potentials successfully reproduce key thermophysical properties of interest, such as temperature-dependent lattice parameter, specific heat capacity, and bulk modulus. We also evaluated the energy of stoichiometric defect reactions and defect migration barriers and found close agreement with DFT predictions, demonstrating that our potentials can be used for modeling defects in UN. Additional tests provide evidence that our potentials are reliable for simulating diffusion, noble gas impurities, and radiation damage.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Conditional probability density functional theory

Here we present conditional probability (CP) density functional theory (DFT) as a formally exact theory. In essence, CP-DFT determines the ground-state energy of a system by finding the CP density from a series of independent Kohn-Sham (KS) DFT calculations. By directly calculating CP densities, we bypass the need for an approximate XC energy functional. In this paper we discuss and derive several key properties of the CP density and corresponding CP-KS potential. Illustrative examples are used throughout to help guide the reader through the various concepts and theory presented. We explore a suitable CP-DFT approximation and discuss exact conditions, limitations, and results for selected examples.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Potassium-induced effect on structure and chemical activity of Cu x O/Cu(111) (x≤2) surface: A combined scanning tunneling microscopy and density functional theory study

Potassium (K) plays an essential role in promoting catalytic reaction in many established industrial catalytic processes. Here, we report a combined study using scanning tunneling microscopy (STM) and density functional theory (DFT) in understanding the effect of depositing K on the atomic and electronic structures as well as chemical activities of Cu x O/Cu(111) (x≤2). The DFT calculations observe a pseudomorphic growth of K on Cu x O/Cu(111) up to 0.19 monolayer (ML) of coverage, where K binds the surface via strong ionic interaction with chemisorbed oxygen and the relatively weak electrostatic interactions with copper ions, lower and upper oxygen on the Cu x O rings. The simulated STM pattern based on the DFT results agrees well with the experimental observations. The deposited K displays great impact on the surface electronic structure of Cu x O/Cu(111), which induces significant reduction in work function and leads to a strong electron polarization on the surface. The promotion of K on the surface binding properties is selective. It varies depending on the nature of adsorbates. According to our results, K has little effect on surface acidity, while it enhances the surface basicity significantly. As a consequence, the presence of K does not help for CO adsorption on Cu x O/Cu(111), but being able to accelerate the activation of CO 2 . Thus, such promotion strongly depends on the combinations from both geometric and electronic effects. Our results highlight the origin of promoting effect of alkalis in the design of catalysts for the complex reactions.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2023

Phonons in Copper Diphosphide (CuP 2 ): Raman Spectroscopy and Lattice Dynamics Calculations

Copper diphosphide (CuP 2 ) is an emerging binary semiconductor with promising properties for energy conversion and storage applications. While functionality and possible applications of CuP 2 have been studied, there is a curious gap in the investigation of its vibrational properties. In this work, we provide a reference Raman spectrum of CuP 2 , with a complete analysis of all Raman active modes from both experimental and theoretical perspectives. Raman measurements have been performed on polycrystalline CuP 2 thin films with close to stoichiometric composition. Detailed deconvolution of the Raman spectrum with Lorentzian curves has allowed identification of all theoretically predicted Raman active modes (9A g and 9B g ), including their positions and symmetry assignment. Furthermore, calculations of the phonon density of states (PDOS), as well as the phonon dispersions, provide a microscopic understanding of the experimentally observed phonon lines, in addition to the assignment to the specific lattice eigenmodes. We further provide the theoretically predicted positions of the infrared (IR) active modes, along with the simulated IR spectrum from density functional theory (DFT). Overall good agreement is found between the experimental and DFT-calculated Raman spectra of CuP 2 , providing a reference platform for future investigations on this material.

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

Identification of a catalytic iron-hydride at the H-cluster of [FeFe]-hydrogenase

Hydrogenases couple electrochemical potential to the reversible chemical transformation of H 2 and protons, yet the reaction mechanism and composition of intermediates are not fully understood. In this Communication we describe the biophysical properties of a hydride-bound state (H hyd ) of the [FeFe]-hydrogenase from Chlamydomonas reinhardtii . The catalytic H-cluster of [FeFe]-hydrogenase consists of a [4Fe-4S] subcluster ([4Fe-4S] H ) linked by a cysteine thiol to an azadithiolate-bridged 2Fe subcluster ([2Fe] H ) with CO and CN- ligands. Mossbauer analysis and density functional theory (DFT) calculations show that H hyd consists of a reduced [4Fe-4S] H + coupled to a diferrous [2Fe] H with a terminally bound Fe-hydride. The existence of the Fe-hydride in Hhyd was demonstrated by an unusually low Mossbauer isomer shift of the distal Fe of the [2Fe] H subcluster. As a result, a DFT model of H hyd shows that the Fe-hydride is part of a H-bonding network with the nearby bridging azadithiolate to facilitate fast proton exchange and catalytic turnover.

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

Effect of interatomic potential on the energetics of hydrogen and helium-vacancy complexes in bulk, or near surfaces of tungsten

Hydrogen (H) trapping by helium-vacancy (He$-$V) complexes in bulk and the near surface region of tungsten (W) have been investigated by molecular statics calculations that evaluate two different W$-$H interatomic potentials, which use the same W$-$He, He$-$He and He$-$H potentials. One of the W$-$H potentials is a bond-order potential (BOP) developed by Juslin et al., while the other is an embedding atom method (EAM) potential developed by Wang et al. Both potentials overestimate the H binding energies to He clusters in bulk W, as compared to DFT calculations, but properly predict the functional form of the H binding energies to He clusters with increasing number of He and H. The BOP simulations reveal that H binding energies to He x V complexes generally increase with increasing number of He. However, the EAM results indicate that the H binding energy as a function of number of He depends on the number of H, and the H binding energies change slightly at high He content. Compared with available DFT data, both BOP and EAM underestimate the H binding energies to He x V 2 H m complexes. Here, the BOP reproduces the He formation energy below a W surface, while the EAM potential better reproduces the H formation energy and the interactions between H and He$-$V complexes. Based on these comparisons, we determine that the EAM potential is more accurate than BOP for large-scale molecular dynamics simulations of W$-$He$-$H interactions. The EAM potential predicts that the difference in the average binding energies of H to stable He$-$V complexes near the W surface is less than 0.2 eV and the difference decreases with increasing He content. Thus, the EAM potential indicates that the effect of surfaces on H binding energies to large He$-$V complexes below the W surfaces can be ignored.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Predicting Device Parameters for Dye-Sensitized Solar Cells from Electronic Structure Calculations to Reproduce Experiment.

Given that improvements to the power-conversion efficiency (PCE) of dye-sensitized solar cells (DSSCs) have slowed in recent years, a means to accurately predict the device parameters yielded by trial dyes in silico, without having to synthesize them, would be extremely valuable to speed up the design process. Currently, the best-performing methods of calculating device parameters rely on a set of experimentally determined kinetic coefficients. In practice, it is very difficult to measure these kinetic parameters accurately, limiting the overall accuracy of such predictive methods. This work proposes a model to obtain key parameters such as J(SC), V-OC, and PCE using only the results from density functional theory (DFT) and time-dependent DFT calculations, noting that rates of electron-transfer steps are ultimately linked to the electronic structure of the dye center dot center dot center dot TiO2 working electrode. Six organic DSSC dyes from dissimilar chemical classes (L0, L1, L2, WS-2, WS-92, and C281) were chosen to demonstrate the power of this approach. Their a priori known experimentally determined device performance metrics served to validate our predictions. The greatest absolute error in our predicted PCE values was 0.36% relative to the experiment, while the greatest fractional error was 0.042. This indicates that the proposed model offers a dramatic improvement on previous predictive methods for DSSC device parameters, both in accuracy and in consistency. Moreover, such a predictive model has great potential to be applied to other photovoltaic applications, further enabling the design of novel, highly efficient photoactive materials.

density functional theory↗