Search NASA⌕ Search

SEARCH · Search NASA

Results for “DFT”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Thermodynamic properties of the Nd-Bi system via emf measurements, $\mathrm{DFT}$ calculations, machine learning, and $\mathrm{CALPHAD}$ modeling

Thermodynamic properties of the Nd-Bi system were investigated using a combination of experimental measurements, first-principles calculations based on density functional theory (DFT), data mining and machine learning (DM + ML) predictions, and calculation of phase diagrams (CALPHAD) modeling. The electromotive force (emf) of Nd-Bi alloys in molten LiCl-KCl-NdCl 3 at 773–973 K was measured via coulometric titration of Nd into Bi for the determination of thermochemical properties such as activity coefficients and solubilities of Nd in Bi. A new peritectic reaction of [liquid + NdBi 2 = Nd 3 Bi 7 ] at 774 K was confirmed using differential scanning calorimetry, structural (X-ray diffraction), and microstructural (scanning electron microscopy) analyses. The unknown crystal structure of NdBi2 was suggested to be a mixture of the anti-La 2 Sb configuration and the La 2 Te-type configuration based on ML predictions for over 26,000 data-mined AB 2 -type configurations together with DFT-based verifications. Using the newly acquired experimental data and DFT-based calculations, the thermodynamic description of the Nd-Bi system was remodeled, and a more complete Nd-Bi phase diagram was calculated, including the Nd 3 Bi 7 compound, invariant transition reactions, and liquidus temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

CO adsorption on MnO(100): Experimental benchmarks compared to DFT

In this study, CO adsorption on the MnO(100) surface was studied using temperature programmed desorption (TPD) and density functional theory (DFT). TPD results show that CO is weakly-bound on MnO(100), with an experimental adsorption energy of -35.6 ± 2.1 kJ/mol at terrace sites in the low coverage limit. PBE simulations suggest that CO adsorption causes an implausible (2 × 2) surface reconstruction. PBE+U simulations show no signs of surface reconstruction, and provide an accurate estimate of the adsorption energy (-36.4 kJ/mol) when combined with the DFT-D3 method with Becke-Jonson damping to correct for van der Waals interactions. This simulation also shows that CO adsorbs C-down onto the Mn 2+ terrace site in a tilted geometry, which is also observed experimentally and computationally on the similarly-structured NiO(100) transition metal oxide surface. TPD results for large doses show a plateauing of the coverage at about 0.4 monolayers of CO at 85 K, with a defect coverage equivalent to 0.08 monolayers. Adsorption associated with defect sites is indicated by a high-temperature desorption tail which is not satisfactorily explained by DFT simulations of simple step or oxygen vacancy defects.

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

DFT-based mesoscopic interactions of N adatoms on Cu(100)

The first-principles density functional theory (DFT) method is used to directly calculate the interactions of N atoms of various configurations on Cu(100) surfaces. This allows us to derive a mesoscopic substratemediated interaction model. We find that, for these Cu(100)-supported two-dimensional N clusters, oscillating substrate-mediated interactions obtained from an elastic eigenvector approach including a mesoscopic interacting range (called the mesoscopic model) have reasonable agreement with the DFT results. The DFT periodic boundary conditions are properly mapped into the mesoscopic model. The model is used to analyze the energies of c(2×2) square islands independent of their sizes. Previous experimental studies have shown that, in the low coverage region, N atoms agglomerate to isolated square islands of about 5 nm × 5 nm in size with a c(2×2) structure. With increasing coverage, those islands form regular patterns but do not touch below about 0.4 monolayers. From previous experiments and calculations, the square island pattern shows significant elastic effects which are now confirmed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

Interaction of Methane with Calcium Ferrite in the Chemical Looping Partial Oxidation Application: Experimental and DFT Study

In this study, experimental studies using mass spectroscopy, thermogravimetric analysis, in situ X-ray diffraction (XRD), and theoretical density functional theory calculations were carried out to investigate the mechanism of methane partial oxidation over a CaFe 2 O 4 oxygen carrier. The XRD analysis revealed the formation of the calcium ferrite phases during methane reduction occurred in the sequence CaFe 2 O 4 →CaFe 3 O 5 + Ca 2 Fe 2 O 5 →FeO + Fe + CaO→Fe + CaO. Each of these phases was modeled using periodic density functional theory to study the activity of these phases for the reactivity with CH 4 , CO, and H 2 along with the adsorbed surface species to help elucidate the effect that phase changes during reduction have on the chemical looping products CO and H 2 . Mass spectroscopic analysis revealed the formation of H 2 O and CO 2 occurred prior to the formation of CO and H 2 , and these products may be due to the removal of lattice oxygen bound in the discrete primary phase, which contributes to a small portion (~16%) of the total transferable oxygen in the solid. The density functional theory (DFT) results suggest that the gaseous CO and H 2 are not thermodynamically favorable to react with the Ca 2 Fe 2 O 5 where the release of CO and H 2 led to an increase in selectivity for these products. The DFT results suggest that some of the adsorbed atomic hydrogen combined to form adsorbed H 2 and then gaseous H 2 on the CaFe 2 O 4 surface, while the carbon deposition and some adsorbed atomic hydrogen react to form the CaFe 3 O 5 and Ca 2 Fe 2 O 5 phases simultaneously. DFT simulations suggest that both the gaseous H 2 and adsorbed atomic H are less thermodynamically favorable to drive the deep reduction of Ca 2 Fe 2 O 5 to CaO and Fe0. The deep reduction of the calcium ferrite therefore favors the generation of syngas products (CO and H 2 ) through mechanisms of methane dehydrogenation and selective oxidation of carbon to CO. The significant difference in mechanisms for partial oxidation of methane to syngas with the CaFe 2 O 4 carrier can be correlated to the differences in the oxide phases present during the reduction process, concentration of oxygen species, and bonding in the lattice structure.

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

Cu(II) Stability and UV-Induced Electron Transfer in a Metal–Organic Hybrid: An EPR, DFT, and Crystallographic Characterization of Copper-Doped Zinc Creatininium Sulfate

Single-crystal X-ray diffraction and electron paramagnetic resonance (EPR) spectroscopic experiments, complemented by quantum chemical DFT calculations, were carried out on the copper-doped metal–organic hybrid and Tutton salt analogue zinc creatininium sulfate to determine its crystal structure, to characterize the electronic structure of the doped Cu(II) binding site, and to propose a pathway for an excited-state, proton-coupled electron transfer (PCET) process in UV-exposed crystals. The crystal structure is isomorphous to that of cadmium creatininium sulfate, which has the transition ion, not in direct coordination with the creatinine, but forming a hexahydrate complex, which is bridged to a creatininium through an intervening sulfate ion. The EPR g (2.446, 2.112, 2.082) and copper hyperfine (A Cu : -327, -59.6, 10.8 MHz) tensor parameters are consistent with doped copper replacing host zinc in the metal–hexahydrate complex. These parameters are similar to those observed for copper hexahydrate in doped Tutton salt systems at low temperature, where the unpaired electron occupies mainly the copper 3d x 2 –y 2 orbital. At room temperature in the Tutton systems, vibration couplings stemming from a dynamic Jahn–Teller effect cause tensor averaging which results in a reduction in their maximum g-tensor and hyperfine tensor values. However, like for the doped isomorphous Cd creatinine crystal, the Cu(II) EPR exhibits little, or no room temperature averaging compared to its low temperature pattern. Samples exposed to 254 nm UV light generate a carbon-centered free radical species, characterized by an isotropic g-tensor (g = 2.0029) and an alpha-proton hyperfine coupling (-24 -14 +4 G). These parameters identify it as a creatinine radical cation formed by the oxidative release of one of its C2 methylene hydrogens. DFT calculations confirm the unpaired electronic structures of both the Cu(II) site and free radical. The growth in radical concentration with an increase in the UV exposure time coincides with a decrease in the copper EPR signal, indicating a coupled light-induced oxidation reduction process. A comparison of the crystal structure with the EPR parameters and DFT results provides evidence for a UV-induced PCET.

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

Benchmarking DFT Accuracy in Predicting O 1s Binding Energies on Metals

X-ray photoelectron spectroscopy (XPS) is a powerful tool for probing the electronic structure and composition of materials, particularly metals and metal oxides of relevance to solar cells and catalysis. Density functional theory (DFT) is often used to support XPS peak assignments, but its reliability for predicting oxygen species is not well established. Here, we compile a large data set of experimental oxygen binding energies and evaluate corresponding DFT predictions. We find that as the binding energies of metal-bound atomic oxygen species increase, especially above ≈530 eV, there is a general decrease in the accuracy of DFTpredicted values. Thus, high-binding-energy atomic oxygen species, such as those proposed as active for selective Ag-catalyzed epoxidation, are less well represented. The chemical nature of the oxygen species also influences accuracy, with molecularly bound species more reliably captured across the entire range of energies. These findings illustrate the limitations of DFT for interpreting XPS spectra and provide a benchmark for improving computational methods.

Adsorption↗
DOE OSTI2020

DFT simulation of the X-ray diffraction pattern of aluminum-ion-intercalated graphite used as the cathode material of the aluminum-ion battery

Many breakthroughs have been achieved in rechargeable aluminum-ion battery technologies in recent years. Most recently, operando X-ray diffraction (XRD) combined with density functional theory (DFT) calculations was reported to study the chloroaluminate anion (AlCl 4 -)-intercalated graphite cathode of the battery. However, there are quite a few discrepancies between the measured and simulated XRD patterns. This work is focused on the simulation of XRD patterns of graphite intercalation compounds (GICs) with DFT calculations. Our results reveal that both the geometry of AlCl 4 - in graphite and the gallery height of GICs are dependent on the intercalant density. At low intercalant density, the gallery height keeps constant, but at high intercalant densities, the gallery height is linearly related to the intercalant density. Our simulated XRD patterns are highly consistent with the measured operando XRD patterns. Not only do the angles of the peaks match very well, but also the relative intensities and the corresponding electrode capacities show reasonable agreement with the experimental results. The DFT simulation of the XRD pattern provides significant information on the stage index and the charge capacity of the GIC electrode.

25 ENERGY STORAGE↗
DOE OSTI2020

Evolution of lithium ordering with (de)-lithiation in β-LiVOPO 4 : insights through solid-state NMR and first principles DFT calculations

The lithium-ion battery cathode material β-VOPO 4 is capable of intercalating more than one Li ion per transition metal ion due to the accessibility of both the V 5+ /V 4+ and V 4+ /V 3+ redox couples at ~4.5 V and ~2.3 V vs. Li, respectively, giving a theoretical capacity greater than ~300 mA h g -1 . The ability to perform full and reversible two Li-ion intercalation in this material, however, has been a matter of debate and the poor crystallinity of the fully lithiated phase has thus far precluded its complete structural characterisation by conventional diffraction-based methods. In this work, 7 Li and 31 P NMR spectroscopy, in combination with first principles DFT calculations, indicate that chemical lithiation results in a single phase β-Li 2 VOPO 4 exhibiting a complex Li ordering scheme with lithium ions occupying multiple disordered environments. 2D NMR 7 Li correlation experiments were used to deduce the most likely Li ordering for the β-Li 2 VOPO 4 phase from amongst several DFT optimised structures. In contrast, electrochemically lithiated β-Li 2- x VOPO 4 discharged to 1.6 V exhibits, in addition to β-Li 2 VOPO 4 , a β-Li 1.5 VOPO 4 phase. Additionally, the existence of β-Li 1.5 VOPO 4 is not reflected in the flat galvanostatic charge and discharge curves nor is evident from diffraction-based methods due to the very close structural similarity between the β-Li 1.5 VOPO 4 phase and β-Li 2 VOPO 4 phases. We demonstrate that solid state NMR spectroscopy, in combination with DFT results, provides a powerful tool for identifying intermediate states formed during charge/discharge of these complex phosphates as these phases can be distinguished from the end member phases primarily by the nature of the lithium ordering.

25 ENERGY STORAGE↗
DOE OSTI2022

DFT-based QM/MM with Particle-Mesh Ewald for Direct, Long-Range Electrostatic Embedding

In this work, we present a DFT-based, QM/MM implementation with long-range electrostatic embedding achieved by direct real-space integration of the particle mesh Ewald (PME) computed electrostatic potential. The key transformation is the interpolation of the electrostatic potential from the PME grid to the DFT quadrature grid, from which integrals are easily evaluated utilizing standard DFT machinery. We provide benchmarks of the numerical accuracy with choice of grid size and real-space corrections, and demonstrate that good convergence is achieved while introducing nominal computational overhead. Furthermore, the approach requires only small modification to existing software packages, as is demonstrated with our implementation in the OpenMM and Psi 4 software. After presenting convergence benchmarks, we evaluate the importance of long-range electrostatic embedding in three solute/solvent systems modeled with QM/MM. Water and BMIM/BF 4 ionic liquid were considered as "simple" and "complex" solvents respectively, with water and p-phenylenediamine (PPD) solute molecules treated at QM level of theory. While electrostatic embedding with standard real-space truncation may introduce negligible error for simple systems such as water solute in water solvent, errors become more significant when QM/MM is applied to complex solvents such as ionic liquids. An extreme example is the electrostatic embedding energy for oxidized PPD in BMIM/BF 4 for which real-space truncation produces severe error even at 2-3 nm cutoff distances. This latter example illustrates that utilization of QM/MM to compute redox potentials within concentrated electrolytes/ionic media requires carefully chosen long-range electrostatic embedding algorithms, with our presented algorithm providing a general and robust approach.

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

Assessment of DFT functionals for a minimal nitrogenase [Fe(SH) 4 H]- model employing state-of-the-art ab initio methods

We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ΔE estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree–Fock included (0%–54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ΔE is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn–Sham orbitals gives results within 4–11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ~5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.

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

Accuracy of DFT computed oxygen-vacancy formation energies and high-throughput search of solar thermochemical water-splitting compounds

The enthalpy change involved in metal oxide reduction is a key quantity in various processes related to energy conversion and storage, and is of particular interest for computational prediction. Often this prediction involves the simulation of a high temperature reduction process with a 0K methodology like density functional theory (DFT), and it is not infrequent for the high temperature and 0K stable crystal structures to differ. This introduces a conundrum with regards to the choice of crystal structure to utilize in the computation, with approaches in the literature varying and experimental validation remaining scarce. In this work we address both the crystal structure conundrum and the experimental validation, and then apply the insights we gain to guide a high-throughput search for new materials for solar thermochemical water-splitting applications. By computing the DFT+U oxygen vacancy formation energy (ΔE vf ) of a selection of ABO 3 compounds and comparing different crystal structures for each composition, we highlight the issues that arise when the structure utilized in the computation is dynamically unstable at 0K, namely the presence of an artificial lowering of ΔE vf , and the lack of convergence of ΔE vf with cell size. We solve these limitations by identifying and employing a suitable surrogate dynamically stable structure. We then validate the predictive power of our calculations against appositely generated experimental measurements of reduction enthalpy for a series of Hubbard U values, finding an accuracy ranging between 0.2-0.6 eV/O. In light of such conclusions, we revise and expand a previous a high-throughput DFT study on ABO 3 perovskite oxides. As a result, we provide a list of candidate STCH materials, highlight trends with redox-active cation and structural distortion, and identify Mn 4+ , Mn 3+ and Co 3+ as the most promising redox-active cations.

08 HYDROGEN↗
DOE OSTI2025

Strong Correlation DMRG and DFT

This project developed new ways to improve computer simulations of materials where electrons interact strongly with each other, a challenge for today’s most widely used method, density functional theory (DFT). We used an exact numerical method, the density matrix renormalization group (DMRG), to create highly accurate reference results for simple model systems, and used these to test DFT, prove when it will converge, and even train machine-learned functionals. We also invented new kinds of localized basis functions (“gausslets” and “multi-sliced gausslets”) and a “sliced-basis” approach that make high-accuracy simulations faster and more practical. These methods were applied to extended hydrogen systems, enabling the direct derivation of accurate low-energy models from first-principles calculations. We also introduced a new formalism, Conditional-Probability DFT, which could bypass traditional approximations. The tools and results from this work, including open-source software releases, will help scientists design and understand complex quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1986

VLSI Architectures for Computing DFT's

Simplifications result from use of residue Fermat number systems. System of finite arithmetic over residue Fermat number systems enables calculation of discrete Fourier transform (DFT) of series of complex numbers with reduced number of multiplications. Computer architectures based on approach suitable for design of very-large-scale integrated (VLSI) circuits for computing DFT's. General approach not limited to DFT's; Applicable to decoding of error-correcting codes and other transform calculations. System readily implemented in VLSI.

Truong, T. K.↗
NASA NTRS1986

A VLSI pipeline design of a fast prime factor DFT on a finite field

A conventional prime factor discrete Fourier transform (DFT) algorithm is used to realize a discrete Fourier-like transform on the finite field, GF(q sub n). A pipeline structure is used to implement this prime factor DFT over GF(q sub n). This algorithm is developed to compute cyclic convolutions of complex numbers and to decode Reed-Solomon codes. Such a pipeline fast prime factor DFT algorithm over GF(q sub n) is regular, simple, expandable, and naturally suitable for VLSI implementation. An example illustrating the pipeline aspect of a 30-point transform over GF(q sub n) is presented.

Truong, T. K.↗
NASA NTRS1986

On the number of multiplications necessary to compute a length-2 exp n DFT

The number of multiplications necessary and sufficient to compute a length-2 exp n DFT is determined. The method of derivation is shown to apply to the multiplicative complexity results of Winograd (1980, 1981) for a length-p exp n DFT, for p an odd prime number. The multiplicative complexity of the one-dimensional DFT is summarized for many possible lengths.

Heideman, M. T.↗
DOE OSTI2025

PFOA catalytic reaction mechanisms on zerovalent iron (Fe 0 ): A DFT investigation

In this study, density functional theory (DFT) calculations with dispersion corrections are employed to investigate catalytic reactions of perfluorooctanoic acid (PFOA) onto zerovalent iron (Fe 0 ). The main goal of this investigation is to explain the reaction mechanisms, including bond dissociation energies, activation barriers, and rate-determining steps (RDSs), of PFOA degradation. Additionally, we aim to reveal the catalytic effects of Fe 0 by comparing the reaction energy profile for the degradation of PFOA both in solution (i.e., isolated) and after adsorption on Fe 0 . Along with investigating previously reported pathways, this study proposes a new PFOA degradation pathway on Fe0. Overall, the results revealed that Fe 0 plays a key role in the catalytic degradation of PFOA, significantly affecting the RDSs and increasing defluorination rates. For instance, the activation barrier for the decarboxylation step significantly decreased from 3.08 eV to 1.11 eV over the Fe0 surface, while the activation barrier for F dissociation decreased from ∼5 eV to 0.54 eV. Additionally, the proposed pathway in this study (i.e., successive defluorination) was found to be the most favorable pathway for PFOA degradation on Fe 0 . Taken together, the DFT results are consistent with reported experimental findings and provide critical insights that will help advance the design of catalysts for PFAS degradation.

Catalysis↗
DOE OSTI2021

DFT Mechanism Studies: Biomimetic 1,4-NADH Chemoselective, Co-factor Regeneration with [Cp*Rh(bpy)H] + , in Tandem with the Biocatalysis Pathways of a Core Model of the (HLADH)-Zn(II) Mediated Enzyme, in the Enantioselective Reduction of Achiral Ketones to Chiral S-Alcohols

In this study, Quantum Chemical (QC) calculations, utilizing Density Functional Theory (DFT), were performed to investigate the mechanistic aspects of the chemoselective catalyzed reaction of [Cp*Rh(bpy)H] + with the biomimetic NAD + analogues, N-benzylnicotinamide triflate, 1, and β-nicotinamide ribose-5'-methyl phosphate, 2, in the conversion to their 1,4-NADH analogues, 1,4-dihydro-N-benzylnicotinamide, 4, and β-1,4-dihydronicotinamide-5'-ribose methyl phosphate, 5. This reaction was in tandem with the 1,4-NADH dependent HLADH-Zn(II)- catalyzed reduction of achiral ketones to chiral S-alcohols. The [Cp*Rh(bpy)H] + complex, and not its equilibrium tautomer, [η 4 -Cp*HRh(bpy)] + , was found to control the hydride transfer during the biomimetic NAD + /1,4-NADH conversion, through the non-covalent interactions of the biomimetic co-factors with [Cp*Rh(bpy)H] + . The thermodynamics and kinetics for the chiral reduction of the Zn(II) bound ketones, 2-pentanone and 4-phenyl-2-butanone, with co-factor, 4, catalyzed by Zn(SCH 3 ) 2 (Imidazole), a core model of the Zn(II)-based catalytic center of HLADH, was also investigated by the evaluation of two possible reaction pathways: (1) formation of a ZnH from the C4-H hydride transfer of co-factor, 4, followed by reaction of the postulated ZnH with the bound 2-pentanone or 4-phenyl-2-butanone substrate, and (2), the direct C4-H transfer to the bound achiral ketone substrates, to provide the dominant chiral alcohols, S-2-pentanol or S-4-phenyl-2-butanol. The latter pathway was found most viable, and DFT calculations also revealed an essential η 2 -coordination of the 5,6 double bond of co-factor, 4, to the HLADH-Zn(II) metal ion center, upon imidazole decomplexation, providing an asymmetric differentiation of S-η 2 -5,6-1,4-NADH-Zn(II) binding. A proposed new paradigm for the Zn(II)'s non-innocent role in the HLADH-Zn(II) biocatalysis reduction mechanism, for enantioselective hydride transfer to a Zn(II) bound ketone, providing S-alcohols.

1,4 NADH co-factors↗
DOE OSTI2021

Bayesian automated weighting of aggregated DFT, MD, and experimental data for candidate thermodynamic models of aluminum with uncertainty quantification

Atomic-scale modeling methods such as density functional theory (DFT) and molecular dynamics (MD) can predict the thermodynamic properties of materials at a lower cost than experimental measurements. However, their regular usage in thermodynamic model construction is hampered by the lack of quantitative agreement with experimental measurements and the lack of uncertainty estimates on the data. To make regular usage of this atomistic simulation data, it is important to assess whether the atomistic simulation datasets, by themselves or in combination with experimental measurements, result in the same physics-informed models best supported by experimental measurements alone. Here, models of aluminum thermodynamic properties are discussed using three data sources: atomistic calculations (DFT and MD), experiments, and a combination of atomistic calculations and experiments. The study shows that, after ensuring self-consistency in predicting key invariant points, both experimental measurements and atomistic calculations can significantly contribute to an optimal model.

36 MATERIALS SCIENCE↗