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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Understanding the influence of solvents on the Pt-catalyzed hydrodeoxygenation of guaiacol

Bio-oil derived from fast pyrolysis of lignocellulosic biomass needs to be deoxygenated to become a substitute for petroleum fuels. Here, we study the hydrodeoxygenation mechanism of guaiacol, a bio-oil model compound derived from the lignin fraction of biomass, on Pt(1 1 1) terrace sites in the presence of water, diethyl ether, 1-butanol, and n-hexane as solvent. Using first-principles periodic density functional theory (DFT) calculations and mean-field microkinetic reactor modeling, a detailed reaction mechanism is investigated targeting various products such as catechol, phenol, anisole, benzene, cyclohexanone, and cyclohexanol. Solvent phase DFT outcomes are mostly similar to that of the vapor phase; however, microkinetic modeling results suggest that rate controlling species and transition states differ somewhat in the various reaction environments. Catechol was found to be the major aromatic product across all reaction environments. Over Pt(1 1 1), unsaturated monooxygenate production from catechol is unlikely.

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

Trends in Siting of Metals in Heterometallic Nd–Yb Metal–Organic Frameworks and Molecular Crystals

Several studies suggest that metal ordering within metal–organic frameworks (MOFs) is important for understanding how MOFs behave in relevant applications; however, these siting trends can be difficult to determine experimentally. To garner insight into the energetic driving forces that may lead to nonrandom ordering within heterometallic MOFs, we employ density functional theory (DFT) calculations on several bimetallic metal–organic crystals composed of Nd and Yb metal atoms. We also investigate the metal siting trends for a newly synthesized MOF. Our DFT-based energy of mixing results suggest that Nd will likely occupy sites with greater access to electronegative atoms and that local homometallic domains within a mixed-metal Nd–Yb system are favored. We also explore the use of less computationally extensive methods such as classical force fields and cluster expansion models to understand their feasibility for large system sizes. Here, this study highlights the impact of metal ordering on the energetic stability of heterometallic MOFs and crystal structures.

energy of mixing↗
DOE OSTI2022

Sulfated Zirconium Metal–Organic Frameworks as Well-Defined Supports for Enhancing Organometallic Catalysis

Understanding heterogeneous catalysts is a challenging pursuit due to surface site nonuniformity and aperiodicity in traditionally-used materials. One example is sulfated metal oxides, which function as highly-active catalysts and as supports for organometallic complexes. These applications are due to traits such as acidity, ability to act as a weakly coordinating ligand, and aptitude for promoting transformations via radical cation intermediates. Research is ongoing about the structural features of sulfated metal oxides that imbue the aforementioned properties, such as sulfate geometry and coordination. To better understand these materials, metal–organic frameworks (MOFs) have been targeted as structurally-defined analogs. Composed of inorganic nodes and organic linkers, MOFs possess features such as high porosity and crystallinity, which make them attractive for mechanistic studies of heterogeneous catalysts. In this work, a Zr 6 -based MOF NU-1000 is sulfated and characterized using atomically-precise techniques such as single crystal X-ray diffraction (SCXRD) in addition to diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS). Here, the dynamic nature of the sulfate binding motif is found to transition from monodentate, to bidentate, to tridentate depending on the degree of hydration, as supported by density functional theory (DFT) calculations. Heightened Brønsted acidity compared to the parent MOF was observed upon sulfation, and probed through trimethylphosphine oxide (TMPO) physisorption, ammonia sorption, in-situ ammonia DRIFTS, and DFT studies. With the support structure benchmarked, an organoiridium complex was chemisorbed onto the sulfated MOF node and the efficacy of this supported catalyst was demonstrated for stoichiometric and catalytic activation of benzene-d6 and toluene with structure-activity relationships derived.

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

Optimizing Binding Energies of Key Intermediates for CO 2 Hydrogenation to Methanol over Oxide-Supported Copper

Rational optimization of catalytic performance has been one of the major challenges in catalysis. In this paper, we report a bottom-up study on the ability of TiO 2 and ZrO 2 to optimize the CO 2 conversion to methanol on Cu, using combined density functional theory (DFT) calculations, kinetic Monte Carlo (KMC) simulations, in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) measurements, and steady-state flow reactor tests. Furthermore, the theoretical results from DFT and KMC agree with in situ DRIFTS measurements, showing that both TiO 2 and ZrO 2 help to promote methanol synthesis on Cu via carboxyl intermediates and the reverse water–gas-shift (RWGS) pathway; the formate intermediates, on the other hand, likely act as a spectator eventually. The origin of the superior promoting effect of ZrO 2 is associated with the fine-tuning capability of reduced Zr 3+ at the interface, being able to bind the key reaction intermediates, e.g. *CO 2 , *CO, *HCO, and *H 2 CO, moderately to facilitate methanol formation. Our study demonstrates the importance of synergy between theory and experiments to elucidate the complex reaction mechanisms of CO 2 hydrogenation for the realization of a better catalyst by design.

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

First-principles calculated decomposition pathways for LiBH 4 nanoclusters

Here, we analyze thermodynamic stability and decomposition pathways of LiBH 4 nanoclusters using grand-canonical free-energy minimization based on total energies and vibrational frequencies obtained from density-functional theory (DFT) calculations. We consider (LiBH 4 ) n nanoclusters with n = 2 to 12 as reactants, while the possible products include (Li) n , (B) n , (LiB) n , (LiH) n , and Li 2 B n H n ; off-stoichiometric LinBnHm (m ≤ 4n) clusters were considered for n = 2, 3, and 6. Cluster ground-state configurations have been predicted using prototype electrostatic ground-state (PEGS) and genetic algorithm (GA) based structural optimizations. Free-energy calculations show hydrogen release pathways markedly differ from those in bulk LiBH 4 . While experiments have found that the bulk material decomposes into LiH and B, with Li 2 B 12 H 12 as a kinetically inhibited intermediate phase, (LiBH 4 ) n nanoclusters with n ≤ 12 are predicted to decompose into mixed Li n B n clusters via a series of intermediate clusters of Li n B n H m (m ≤ 4n). The calculated pressure-composition isotherms and temperature-pressure isobars exhibit sloping plateaus due to finite size effects on reaction thermodynamics. Generally, decomposition temperatures of free-standing clusters are found to increase with decreasing cluster size due to thermodynamic destabilization of reaction products.

36 MATERIALS SCIENCE↗
DOE OSTI2017

First-principles equation-of-state table of beryllium based on density-functional theory calculations

Beryllium has been considered a superior ablator material for inertial confinement fusion (ICF) target designs. An accurate equation-of-state (EOS) of beryllium under extreme conditions is essential for reliable ICF designs. Based on density-functional theory (DFT) calculations, we have established a wide-range beryllium EOS table of density ρ = 0.001 to 500 g/cm 3 and temperature T = 2000 to 10 8 K. Our first-principle equation-of-state (FPEOS) table is in better agreement with the widely used SESAME EOS table (SESAME 2023) than the average-atom INFERNO and Purgatorio models. For the principal Hugoniot, our FPEOS prediction shows ~10% stiffer than the last two models in the maximum compression. Although the existing experimental data (only up to 17 Mbar) cannot distinguish these EOS models, we anticipate that high-pressure experiments at the maximum compression region should differentiate our FPEOS from INFERNO and Purgatorio models. Comparisons between FPEOS and SESAME EOS for off-Hugoniot conditions show that the differences in the pressure and internal energy are within ~20%. By implementing the FPEOS table into the 1-D radiation–hydrodynamic code LILAC, we studied in this paper the EOS effects on beryllium-shell–target implosions. Finally, the FPEOS simulation predicts higher neutron yield (~15%) compared to the simulation using the SESAME 2023 EOS table.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Comparing Classical and Machine Learning Force Fields for Modeling Deformation of Metal–Organic Frameworks Relevant for Direct Air Capture

Deformation of metal–organic frameworks (MOFs) induced by adsorbate molecules can affect adsorption properties such as capacity and selectivity, but most computational studies of MOFs assume framework rigidity to simplify calculations. Although flexible force fields (FFs) for MOFs have been parametrized for specific materials, the generality of FFs for reliably modeling adsorbate-induced deformation to accuracy nearing that of density functional theory (DFT) has not been established. This work confirms using DFT calculations that adsorbate-induced deformation can affect CO 2 and H 2 O adsorption energies in a considerable fraction of MOFs promising for direct air capture (DAC). We then benchmark the efficacy of several general-purpose FFs in describing adsorbate-induced deformation for DAC against DFT. Our results show that current classical FFs are insufficient for describing MOF deformation, especially in cases of interest for DAC where strong interactions exist between adsorbed molecules and MOF frameworks. Some emerging machine learning force fields (MLFFs) we tested, particularly CHGNet, MACE-MP-0, and Equiformer V2, appear to be more promising than the classical FF for emulating the deformation behavior described by DFT. The best performing FF (CHGNet), however, fails to achieve the accuracy required for practical predictions with a mean absolute adsorption energy error of 0.124 eV.

adsorption↗
DOE OSTI2022

Electron Counting and High-Pressure Phase Transformations in Metal Hexaborides

Pressure-induced structural transitions of the alkaline earth hexaborides, CaB 6 , SrB 6 , and BaB 6 , are studied theoretically using electron counting rules and density functional theory calculations. Here, we demonstrate the applicability of gas-phase borane electron counting methods to solid-state metal borides under pressure and validate the assumptions of the rules by density functional theory (DFT) calculations. All three compounds share ambient-pressure and high-pressure structures, but BaB 6 differs from CaB 6 and SrB 6 at intermediate pressures. The unique BaB 6 phase is shown to break electron counting rules, while all other phases obey them. This anomaly is resolved by DFT, which reveals B–Ba covalency and unusual B–B π bonding under pressure. The relationships between structure and bonding can help us to understand the exotic behavior of lanthanide hexaborides and design new borides with desirable properties. Developing electron counting procedures for solids will enhance materials discovery efforts with chemical intuition.

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

Bulk Fermi surfaces of the Dirac type-II semimetallic candidate NiTe 2

Here, we present a study on the Fermi surface of the Dirac type-II semimetallic candidate NiTe 2 via the temperature and angular dependence of the de Haas–van Alphen (dHvA) effect measured in single crystals grown through Te flux. In contrast to its isostructural compounds like PtSe 2 , band-structure calculations predict NiTe 2 to display a tilted Dirac node very close to its Fermi level that is located along the Γ to A high-symmetry direction within its first Brillouin zone. The angular dependence of the dHvA frequencies is found to be in agreement with the first-principles calculations when the electronic bands are slightly shifted with respect to the Fermi level (ε F ), and therefore provide support for the existence of a Dirac type-II node in NiTe 2 . Nevertheless, we observed mild disagreements between experimental observations and density functional theory (DFT) calculations as, for example, nearly isotropic and light experimental effective masses. This indicates that the dispersion of the bands is not well captured by DFT. Despite the coexistence of Dirac-like fermions with topologically trivial carriers, samples of the highest quality display an anomalous and large either linear or sublinear magnetoresistivity. Finally, this suggests that Lorentz invariance breaking Dirac-like quasiparticles dominate the carrier transport in this compound.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Hydride Accessibility and Reactivity in the Configurational and Stoichiometric Space of β-Ga 2 O 3 for CO 2 Hydrogenation

Understanding how surface species evolve under reaction conditions is essential for improving catalyst design for efficient CO 2 hydrogenation. This work combines systematic DFT calculations with grand canonical sampling to investigate the stability and reactivity of Ga–H species on β-Ga 2 O 3 across a range of reaction conditions. Initial DFT studies reveal that when Ga–H species are present, they facilitate formate formation via a low-barrier pathway, largely independent of the surface termination or hydrogen site. However, grand canonical sampling shows that under a broad range of reaction conditions─especially at high oxygen chemical potentials associated with high water content─Ga–H species are thermodynamically inaccessible. Furthermore, adsorbed water molecules can block reactive sites, inhibiting CO 2 activation even when hydrides are present. These findings suggest that the lack of accessible hydride species, rather than their intrinsic reactivity, could contribute to reduced catalytic performance of β-Ga 2 O 3 under more oxidizing, high-conversion conditions.

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

Kohn–Sham accuracy from orbital-free density functional theory via Δ-machine learning

Here, we present a Δ-machine learning model for obtaining Kohn–Sham accuracy from orbital-free density functional theory (DFT) calculations. In particular, we employ a machine-learned force field (MLFF) scheme based on the kernel method to capture the difference between Kohn–Sham and orbital-free DFT energies/forces. We implement this model in the context of on-the-fly molecular dynamics simulations and study its accuracy, performance, and sensitivity to parameters for representative systems. We find that the formalism not only improves the accuracy of Thomas–Fermi–von Weizsäcker orbital-free energies and forces by more than two orders of magnitude but is also more accurate than MLFFs based solely on Kohn–Sham DFT while being more efficient and less sensitive to model parameters. We apply the framework to study the structure of molten Al 0.88 Si 0.12 , the results suggesting no aggregation of Si atoms, in agreement with a previous Kohn–Sham study performed at an order of magnitude smaller length and time scales.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Ground State Properties of the Wide Band Gap Semiconductor Beryllium Sulfide (BeS)

We report the results from self-consistent calculations of electronic, transport, and bulk properties of beryllium sulfide (BeS) in the zinc-blende phase, and employed an ab-initio local density approximation (LDA) potential and the linear combination of atomic orbitals (LCAO). We obtained the ground state properties of zb-BeS with the Bagayoko, Zhao, and Williams (BZW) computational method, as enhanced by Ekuma and Franklin (BZW-EF). Our findings include the electronic energy bands, the total (DOS) and partial (pDOS) densities of states, electron and hole effective masses, the equilibrium lattice constant, and the bulk modulus. The calculated band structure clearly shows that zb-BeS has an indirect energy band gap of 5.436 eV, from Γ to a point between Γ and X, for an experimental lattice constant of 4.863 Å. This is in excellent agreement with the experiment, unlike the findings of more than 15 previous density functional theory (DFT) calculations that did not perform the generalized minimization of the energy functional, required by the second DFT theorem, which is inherent to the implementation of our BZW-EF method.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Highly tunable band inversion in AB 2 X 4 (A=Ge, Sn, Pb; B=As, Sb, Bi; X=Se, Te) compounds

Topological materials have been discovered so far largely by searching for existing compounds in crystallographic databases, but there are potentially new topological materials with desirable features that have not been synthesized. One of the desirable features is high tunability resulting from the band inversion with a very small direct band gap, which can be tuned by changes in pressure or strain to induce a topological phase transition. Here, using density-functional theory (DFT) calculations, we have studied the septuple layered AB 2 X 4 series compounds, where A=(Ge, Sn and Pb), B=(As, Sb and Bi), and X=(Se and Te). With the DFT thermodynamic stability validated by the already-reported compounds in these series, we predict stable Se compounds, which are not found in crystallographic database. Among them, we find that GeBi 2 Se 4 and GeSb 2 Se 4 having a small direct band gap at the Z point are very close to a strong topological insulator, which can be tuned by a moderate pressure to induce the band inversion. Importantly, the topological features with the small direct band gap are well isolated in both momentum and energy windows, which offers high tunability for studying the topological phase transition.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Effect of applied potential on metal surfaces: Surface energy, Wulff shape and charge distribution

Here we use grand canonical density functional theory to predict the surface energies, Wulff shapes, charge distributions and catalytically active sites of different metal surfaces under electrochemical conditions. We propose a method for computing surface energies from grand canonical density functional theory (GC-DFT) calculations of periodic slab models and use it to compute the surface energies of the facets of Pt, Cu, and Ag crystals to predict their Wulff shapes under electrochemical conditions. GC-DFT predicts that, for the pure metals studied, solvation only slightly affects the Wulff shape while applied potentials considerably affect the surface energies and corresponding Wulff shapes. We used Bader charge analysis of GC-DFT computed electron densities to investigate the effect of applied potential on the distribution of electron density over the atoms of the surfaces of Pt, Cu, Ag, and the 75–25 Ag-Pt and Au-Ni alloys. This analysis shows that, under an applied potential, the electron density is unevenly distributed over the surface atoms and that the charges of atoms more exposed to solvent are more sensitive to bias. Our results show that the most sensitive atom to bias can be used to identify the most favorable adsorption site and thus, the active sites of electrochemical reactions, which is computationally less demanding than calculating the adsorption energies on all possible adsorption sites.

36 MATERIALS SCIENCE↗
DOE OSTI2013

125Te NMR chemical-shift trends in PbTe–GeTe and PbTe–SnTe alloys

Complex tellurides, such as doped PbTe, GeTe, and their alloys, are among the best thermoelectric materials. Knowledge of the change in 125Te NMR chemical shift due to bonding to dopant or “solute” atoms is useful for determination of phase composition, peak assignment, and analysis of local bonding. We have measured the 125Te NMR chemical shifts in PbTe-based alloys, Pb 1-x Ge x Te and Pb 1-x Sn x Te, which have a rocksalt-like structure, and analyzed their trends. For low x, several peaks are resolved in the 22-kHz MAS 125Te NMR spectra. A simple linear trend in chemical shifts with the number of Pb neighbors is observed. No evidence of a proposed ferroelectric displacement of Ge atoms in a cubic PbTe matrix is detected at low Ge concentrations. The observed chemical shift trends are compared with the results of DFT calculations, which confirm the linear dependence on the composition of the first-neighbor shell. The data enable determination of the composition of various phases in multiphase telluride materials. They also provide estimates of the 125Te chemical shifts of GeTe and SnTe (+970 and +400±150 ppm, respectively, from PbTe), which are otherwise difficult to access due to Knight shifts of many hundreds of ppm in neat GeTe and SnTe.

125Te chemical shift calculation↗
DOE OSTI2023

Understanding of Active Sites and Interconversion of Pd and PdO during CH 4 Oxidation

Pd-based catalysts are widely used in the oxidation of CH 4 and have a significant impact on global warming. However, understanding their active sites remains controversial, because interconversion between Pd and PdO occurs consecutively during the reaction. Understanding the intrinsic active sites under reaction conditions is critical for developing highly active and selective catalysts. In this study, we demonstrated that partially oxidized palladium (PdO x ) on the surface plays an important role for CH 4 oxidation. Regardless of whether the initial state of Pd corresponds to oxides or metallic clusters, the topmost surface is PdO x , which is formed during CH4 oxidation. A quantitative analysis using CO titration, diffuse reflectance infrared Fourier-transform spectroscopy, X-ray diffraction, and scanning transmission electron microscopy demonstrated that a surface PdO layer was formed on top of the metallic Pd clusters during the CH 4 oxidation reaction. Furthermore, the time-on-stream test of CH 4 oxidation revealed that the presence of the PdO layer on top of the metallic Pd clusters improves the catalytic activity. Our periodic density functional theory (DFT) calculations with a PdO x slab and nanoparticle models aided the elucidation of the structure of the experimental PdO particles, as well as the experimental C-O bands. The DFT results also revealed the formation of a PdO layer on the metallic Pd clusters. This study helps achieve a fundamental understanding of the active sites of Pd and PdO for CH 4 oxidation and provides insights into the development of active and durable Pd-based catalysts through molecular-level design.

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

Structural stability, vibrational, and bonding properties of potassium 1, 1′-dinitroamino-5, 5′-bistetrazolate: An emerging green primary explosive

Potassium 1,1′-dinitroamino-5,5′-bistetrazolate (K{sub 2}DNABT) is a nitrogen rich (50.3% by weight, K{sub 2}C{sub 2}N{sub 12}O{sub 4}) green primary explosive with high performance characteristics, namely, velocity of detonation (D = 8.33 km/s), detonation pressure (P = 31.7 GPa), and fast initiating power to replace existing toxic primaries. In the present work, we report density functional theory (DFT) calculations on structural, equation of state, vibrational spectra, electronic structure, and absorption spectra of K{sub 2}DNABT. We have discussed the influence of weak dispersive interactions on structural and vibrational properties through the DFT-D2 method. We find anisotropic compressibility behavior (b<a<c) from pressure dependent structural properties. The predicted equilibrium bulk modulus reveals that K{sub 2}DNABT is softer than toxic lead azide and harder than the most sensitive cyanuric triazide. A complete assignment of all the vibrational modes has been made and compared with the available experimental results. The calculated zone center IR and Raman frequencies show a blue-shift which leads to a hardening of the lattice upon compression. In addition, we have also calculated the electronic structure and absorption spectra using recently developed Tran Blaha-modified Becke Johnson potential. It is found that K{sub 2}DNABT is a direct band gap insulator with a band gap of 3.87 eV and the top of the valence band is mainly dominated by 2p-states of oxygen and nitrogen atoms. K{sub 2}DNABT exhibits mixed ionic (between potassium and tetrazolate ions) and covalent character within tetrazolate molecule. The presence of ionic bonding suggests that the investigated compound is relatively stable and insensitive than covalent primaries. From the calculated absorption spectra, the material is found to decompose under ultra-violet light irradiation.

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

Density Matrix Implementation of the Fermi–Löwdin Orbital Self-Interaction Correction Method

The Fermi–Löwdin orbital self-interaction correction (FLOSIC) method effectively provides a transformation from canonical orbitals to localized Fermi–Löwdin orbitals which are used to remove the self-interaction error in the Perdew–Zunger (PZ) framework. This transformation is solely determined by a set of points in space, called Fermi–Löwdin descriptors (FODs), and the occupied canonical orbitals or the density matrix. In this work, we provide a detailed workflow for the implementation of the FLOSIC method for removal of self-interaction error in DFT calculations in an orbital-by-orbital basis that takes advantage of the unitary invariant nature of the FLOSIC method. In this way, it is possible to cast the self-consistent energy minimization at fixed FODs in the same manner than standard Kohn–Sham with one additional term in the Kohn–Sham Hamiltonian that introduces the PZ self-interaction correction. Each energy minimization iteration is divided in two substeps, one for the density matrix and one for the FODs. Expressions for the effective Kohn–Sham matrix and FOD gradients are provided such that its implementation is suitable for most electronic structure codes. Here, we analyze the convergence characteristics of the algorithm and present applications for the evaluation of NMR shielding constants and real-time time-dependent DFT simulations based on the Liouville–von Neumann equation to calculate excitation energies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗