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

Guiding Principles for the Rational Design of Hybrid Materials: Use of DFT Methodology for Evaluating Non‐Covalent Interactions in a Uranyl Tetrahalide Model System

Abstract Together with the synthesis and experimental characterization of 14 hybrid materials containing [UO 2 X 4 ] 2− (X=Cl − and Br − ) and organic cations, we report on novel methods for determining correlation trends in their formation enthalpy (ΔH f ) and observed vibrational signatures. Δ H f values were analyzed through isothermal acid calorimetry and a Density Functional Theory+Thermodynamics (DFT+T) approach with results showing good agreement between theory and experiment. Three factors (packing efficiency, cation protonation enthalpy, and hydrogen bonding energy [ ]) were assessed as descriptors for trends in Δ H f . Results demonstrated a strong correlation between and Δ H f , highlighting the importance of hydrogen bonding networks in determining the relative stability of solid‐state hybrid materials. Lastly, we investigate how hydrogen bonding networks affect the vibrational characteristics of uranyl solid‐state materials using experimental Raman and IR spectroscopy and theoretical bond orders and find that hydrogen bonding can red‐shift U≡O stretching modes. Overall, the tightly integrated experimental and theoretical studies presented here bridge the trends in macroscopic thermodynamic energies and spectroscopic features with molecular‐level details of the geometry and electronic structure. This modeling framework forms a basis for exploring 3D hydrogen bonding as a tunable design feature in the pursuit of supramolecular materials by rational design.

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

Microkinetic Modeling of the Transient CO 2 Methanation with DFT–Based Uncertainties in a Berty Reactor

The transient operation of methanation reactors can become desirable when coupled with fluctuating renewable energies in Power-to-Gas scenarios, which requires suitable kinetic approaches that can describe the transient catalytic phenomena. A combined experimental and theoretical investigation of the transient CO 2 methanation is conducted using concentration forcing to derive a suitable microkinetic model. Methanation experiments are performed with a Ni/SiO 2 catalyst in a Berty-type reactor at industrially relevant conditions. The microkinetics are based on previous work and were automatically constructed for the Ni(111) facet using the Reaction Mechanism Generator. A feasible set of energetic parameters of the microkinetic models was identified in a theory-constrained optimization procedure within the DFT uncertainty space that can accurately reproduce the experimental results on a first-principles basis. The microkinetic model unravels that the formation of H 2 O* and CH 3 * control the activity and selectivity of Ni(111) under the investigated conditions.

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

Resolving the Chemical Formula of Nesquehonite via NMR Crystallography, DFT Computation, and Complementary Neutron Diffraction

Nesquehonite is a magnesium carbonate mineral relevant to carbon sequestration envisioned for carbon capture and storage of CO 2 . Its chemical formula remains controversial today, assigned as either a hydrated magnesium carbonate [MgCO 3 · 3H 2 O], or a hydroxy- hydrated- magnesium bicarbonate [Mg(HCO 3 )OH · 2H 2 O]. Here, the resolution of this controversy is central to understanding this material‘s thermodynamic, phase, and chemical behavior. In an NMR crystallography study, using rotational-echo double-resonance 13 C{ 1 H} (REDOR), 13 C- 1 H distances are determined with precision, and the combination of 13 C static NMR lineshapes and density functional theory (DFT) calculations are used to model different H atomic coordinates. [MgCO 3 · 3H 2 O] is found to be accurate, and evidence from neutron powder diffraction bolsters these assignments. Refined H positions can help understand how H-bonding stabilizes this structure against dehydration to MgCO 3 . More broadly, these results illustrate the power of NMR crystallography as a technique for resolving questions where X-ray diffraction is inconclusive.

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

Single vs double atom catalyst for N 2 activation in nitrogen reduction reaction: A DFT perspective

Abstract Ammonia synthesis through electrochemical reduction of nitrogen molecules is a promising strategy to significantly reduce the energy consumption in traditional industrial process. Detailed mechanism study of multistep complex nitrogen reduction reaction is prerequisite for the design of highly efficient catalyst. Stable atomically dispersed catalyst with unique geometric and electronic structure is suitable for the mechanism clarification of such a complex reaction. In this study, d ‐block transition‐metal (TM) anchored C 2 N single layer catalyst is investigated by the density functional theory (DFT) calculation. Both single TM‐anchored single atom catalyst (SAC) and double TM‐anchored double atom catalyst (DAC) exhibit good thermodynamic stability in atomically dispersed catalyst. In the case of SACs, IVB metals (Ti, Zr, Hf) exhibit the highest reactivity and lowest overpotential. While in the case of DACs, Cr─Cr system leads to the NH 3 formation, but V─V system leads to the N 2 H 4 formation. The SACs show much lower overpotential and stronger activation of N 2 molecule than the DACs due to the different activation mechanisms: traditional σ ‐donation/ π ‐backdonation N 2 activation mechanism is found in SACs, while a new π ‐donation/π‐backdonation N 2 activation mechanism is found in the DACs. The present work demonstrates that the different catalytic effect for NRR between SAC and DAC and their corresponding electronic structure origin, which gives more insight into the single atom catalyst. image

Qian, Yumin↗
DOE OSTI2025

Examining the Impact of Local Constraint Violations on Energy Computations in DFT

ABSTRACT This work examines the impact of locally imposed constraints in Density Functional Theory (DFT). Using a metric referred to as the extent of violation index (EVI), we quantify how well exchange‐correlation functionals adhere to local constraints. Applying EVIs to a diverse set of molecules for GGA functionals reveals constraint violations, particularly for semi‐empirical functionals. We leverage EVIs to explore potential connections between these violations and errors in chemical properties. While no correlation is observed for atomization energies, a significant statistical correlation emerges between EVIs and total energies. Similarly, the analysis of reaction energies suggests weak positive correlations for specific constraints. However, definitive conclusions about error cancellation mechanisms cannot be made at this time. These observations revealed by EVIs may be useful for consideration when designing future generations of semilocal functionals.

Khanna, Vaibhav [Department of Chemistry Universit↗
DOE OSTI2024

Effect of XC functionals and dispersion corrections on the DFT‐computed structural and vibrational properties of SrCl 2 –NaCl and ZrF 4 –LiF

Density functional theory (DFT) calculations were performed to examine the impact of exchange–correlation (XC) functionals and van der Waals corrections (specifically the D3 method) on the structural and vibrational properties of the SrCl 2 –NaCl and ZrF 4 –LiF salt systems. Multiple XC functionals, including the local density approximation (LDA), the generalized gradient approximation using the Perdew–Burke–Ernzerhof (PBE) model, and its modified form suitable for solids (PBEsol), the dispersion-corrected PBE-D3 and PBEsol-D3, were considered. Of these functionals, LDA was found to exhibit the highest degree of error, while PBEsol and PBE-D3 displayed the least error. Underestimated lattice parameters compared with experimental values were observed to result in higher force constants, leading to an overprediction of vibrational frequencies. Conversely, an overestimation of lattice parameters was associated with lower vibrational frequencies. The methodology presented in this study yielded results that are in good agreement with experiment, irrespective of the method (finite differences vs. density functional perturbation theory) employed for calculating infrared and Raman spectra. It was further demonstrated that for alkali halides with weak Raman scattering, utilizing a supercell constructed from primitive cells better predicts Raman features than does the use of conventional cells.

Raman↗
DOE OSTI2025

Hydride Migration within RhH 2 Ag 19 Superatom: A Combined Neutron Diffraction and DFT Analysis

Abstract An investigation combining neutron diffraction and DFT allows determining the most likely hydride migration pathway within the icosahedral metal framework of [RhH 2 Ag 19 {S 2 P(O n Pr) 2 } 12 ] (RhH 2 Ag 19 ). Starting from the experimentally derived solid‐state structures, a computational analysis is able to reveal an energetically favorable migration pathway with a maximum energy barrier of 4.2 kcal mol −1 . The two hydrides migrate simultaneously within the Rh@Ag 12 icosahedral core, traversing several positional isomers. This study expands the understanding of hydride dynamics in nanoclusters and provides critical insights into the structural flexibility of the superatom framework. These findings have significant implications for hydrogen storage, catalysis, and the design of advanced hydride‐containing materials.

Chemistry↗
DOE OSTI2019

Structure-antioxidant activity relationships, QSAR, DFT calculation, and mechanisms of flavones and flavonols

DFT calculations were performed on nine flavones and flavonols to explain their high antioxidant activity and variations in their activity. Conformational analysis showed that only flavonols with 3-OH directed toward the B-ring are nonplanar; however, flavones and resulted radicals are planar. Hydroxyl group eligible for dissociation is the one with ortho OH directed toward it; otherwise, one of B-ring (in 2′ or 4′-position) hydroxyl group. There are two main factors responsible for stabilizing the resulted radicals and lowering the bond dissociation energy and hence there were found well correlated with the experimental activity. First, driving force resulted from the conversion of nonplanar flavonols to planar radicals accompanied by resonance toward the carbonyl group and H-bond formation with 3-OH and second, radical stabilization by H-bond with ortho hydroxyl group with resonance toward carbonyl or pyrone oxygen. All resonance and H-bonds were confirmed by spin density, bond length, and molecular orbital calculations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Monte-Carlo modeling of phonon thermal transport using DFT-based anisotropic dispersion relations over the full Brillouin zone

In this work, we present a Monte Carlo (MC) approach to solve the phonon Boltzmann transport equation (BTE) in which the anisotropic phonon dispersion relations over the full Brillouin zone (BZ) are used. In this approach, the discretization of the BZ used to compute the phonon relaxation time places constraints on the direction of scattered phonons in the real-space simulation domain. The phonon dispersion and phonon relaxation times are calculated using the density functional theory (DFT) approach. The modified MC approach is validated by a close examination of its ability to simulate phonon transport in both the ballistic and diffusive regimes for multiple materials including GaAs, InAs, ThO 2 , and α-U. In doing so, the phonon thermal conductivities from 100 K to 1000 K are calculated and compared with traditional non-transport solution of the phonon BTE. It is found that the phonon thermal conductivities of α-U and ThO 2 obtained from MC simulations using isotropic dispersion are larger than the values obtained using anisotropic phonon dispersion relations over the full BZ. The effect of phonon-defect scattering on the thermal conductivity of ThO 2 is also studied as an application of the current MC approach and found to agree with previously computed values in the literature. The MC solver developed here has been parallelized as a step to demonstrate its potential to solving computationally intensive phonon thermal transport problems at the mesoscale.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Identifying factors controlling the selective ethane dehydrogenation on Pt-based catalysts from DFT based micro-kinetic modeling

Herein, a detailed EDH reaction mechanism on Pt(111), Pt(211), Pt 3 Sn(111) and Pt 3 Sn(211) surfaces has been systematically computed with DFT. A mean-field micro-kinetic model has been established to simulate the product distribution and selectivity on each surface at experimentally relevant conditions. We have uncovered a number of different factors controlling the catalytic performance of Pt-based EDH catalysts, thus leading to a better understanding of experimentally detected phenomenon. In this study, the (111) and (211) model surfaces of Pt and Pt 3 Sn catalysts were chosen to represent the terrace and step-edge sites when simulating the EDH reaction mechanisms.

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

Understanding of three different polyvinylpyrrolidone (PVP) based battery binders blends on graphene surfaces from first principles via DFT simulations

Binders play a crucial role in binding electrodes, current collectors, and conductive agents during battery manufacturing. The selection of the binders and their blends considerably affect various properties such as mechanical durability, adhesion performance, ionic/electronic conductivities, and solid electrolyte interface stability for lithium-ion batteries. In this work, polyvinylidene di-fluoride (PVDF), polyacrylic acid (PAA), polyvinylpyrrolidone (PVP), and lithiated polyacrylic acid (Li-PAA) were selected as the model binders, and their blends were studied on graphene surfaces. Further, among them, PVDF is hardly soluble in water unlike the other binders, so in practice, the blend of PVP:PVDF is not thought to be a proper combination. However, the PVP:PVDF binder has also been included because the purpose of this work is to calculate several electronic properties such as binding energies, intermolecular interaction energies, bond critical points, electron density, and the Laplacian of electron density at critical bonding regions of all these blends at the molecular level via density functional theory (DFT) simulations in order to evaluate and compare how is the interaction strength and bonding type of three different lithium-ion batteries (LIBs) binder blends on graphene surfaces. Most stable binder pairs and their binding mechanisms on graphene surfaces were studied. The delocalization of lithium in Li-PAA was studied with the presence of water.

25 ENERGY STORAGE↗
DOE OSTI2024

Insights into the interaction of nitrobenzene and the Ag(111) surface: A DFT study

Here, this study explores the potential of nitrobenzene as an anolyte material for nonaqueous redox flow batteries (RFBs) by theoretically examining its low-coverage adsorption behavior on neutral and charged Ag(111) model electrode surfaces. At the low coverage limit, DFT calculations show a preference for nitrobenzene to adsorb parallel to the surface, with the benzene ring and nitro group centered over HCP sites. Interactions between nitrobenzene and the surface were analyzed using induced charge density analysis, Bader charge analysis, and projected density of states (PDOS). It was found that nitrobenzene adsorbs primarily through van der Waals interactions with the surface. As nitrobenzene accumulates negative charge, the strength of adsorption diminishes. Understanding the electrode-electrolyte interface is crucial for enhancing RFB electrochemical performance, and this study sheds light on nitrobenzene's interaction with a model Ag electrode.

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

Insights into Copper Sulfide Formation from Cu and S K edge XAS and DFT studies

An understanding of the fundamentals of the reaction between CuO with trace amounts of H 2 S to form CuS products is critical for the optimal utilization of this process in sulfur removal applications. Unfortunately, CuS is a complex material, featuring various Cu 2-x S compounds (with 0 ≤ x ≤ 1), distorted crystal phases, and varying electronic structures and coordination environments of Cu and S ions. In this work, we combine ex situ and in situ X-ray absorption spectroscopy (XAS) at S and Cu K edges, fixed bed sorption experiments, DFT simulations, and other characterization techniques to speciate the CuS products formed at different temperatures (298–383 K) and from CuO sorbents with different crystallite sizes (2.8–40 nm). The results of our analysis identify the formation of a distorted CuS layer at the surface of CuO crystals with disulfide groups with shorter Cu–S bonds and higher delocalization of the positive charge of the Cu center into (S 1– )2. This distorted CuS layer dominates the XAS signal at lower temperatures (298–323 K) and at the initial stages of sulfidation at higher temperatures (353 and 383 K) where conversion is low (<40%). First-principles atomistic simulations confirm the thermodynamic favorability of the formation of surface (S 1– )2 on both CuO (111) and ($\bar{1}11$) surfaces, providing further support for our experimental observations. Furthermore, these simulations reveal that the presence of disulfide bonds stabilized surface hydroxyl groups, leading to lower Gibbs Free Energies of their surface migration.

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

Improving Predictions of Spin-Crossover Complex Properties through DFT Calculations with a Local Hybrid Functional

We conducted a study on the performance of the local hybrid exchange-correlation functional PBE0r for a set of 95 experimentally-characterized iron spin crossover (SCO) complexes. The PBE0r functional is a variant of PBE0 where the exchange correction is restricted to on-site terms formulated within the basis of local orbitals. We determine the free parameters of the PBE0r functional against experimental data and other hybrid functionals. With a Hartree-Fock (HF) exchange factor of 4%, the PBE0r functional accurately reproduces the electronic and free energy trends predicted in prior DFT studies for these 95 complexes using the B3LYP functional. Larger values of HF exchange stabilize high-spin states. The PBE0r-predicted bond lengths tend to exceed the experimental bond lengths, and bond lengths are less sensitive to HF exchange. The predicted SCO transition temperatures T 1/2 from PBE0r correlate moderately with the experimental transition temperatures, showing a slight improvement compared to the previous modB3LYP-predicted T 1/2 . Furthermore, this study suggests the PBE0r functional as computationally cost-effective and offers the possibility of simulating larger complexes with accuracy comparable to other global hybrid functionals, provided the HF exchange parameter is carefully optimized.

25 ENERGY STORAGE↗
DOE OSTI2025

Exploring the Coexistence of Spin States in [Fe(tpy-Ph) 2 ] 2+ Complexes on Au(111) Using DFT Calculations

In this work, we systematically study the electronic structure and stability of spin states of the [Fe-(tpy-ph) 2 ] 2+ molecule in both the gas phase and on a Au(111) substrate using density functional theory + U (DFT+ U ) calculations. We find that the stability of the Fe 2+ ion’s spin states predicted by the computations is significantly influenced by the Hubbard U parameter. In the gas phase, the low-spin (LS, S = 0) state is found to be energetically favorable for U (Fe) ≤ 3 eV, whereas the high-spin (HS, S = 2) state is stabilized for U (Fe) > 3 eV. Interaction with the Au(111) substrate is found to elevate the critical U for the spin-state transition to 3.5 eV. Additionally, we perform L-edge X-ray absorption spectroscopy (XAS) calculations for both HS and LS states. The calculated XAS suggests that the HS state more closely aligns with the experimental observations, indicating the potential coexistence of the HS state as the initial state during the X-ray excitation process. These findings enrich our understanding of spin-state dynamics in [Fe(tpy-Ph) 2 ] 2+ .

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

Concentration-Dependent Solvation Structure and Dynamics of Aqueous Sulfuric Acid Using Multinuclear NMR and DFT

Sulfuric acid is a ubiquitous compound for industrial processes and aqueous sulfate solutions also play a critical role in electrolytes for many prominent battery chemistries. While the thermodynamic literature for it is quite well-developed, comprehensive studies of the solvation structure, and particularly molecular-scale dynamical and transport properties, are less available. This study applies a multinuclear NMR approach to the elucidation of the solvation structure and dynamics over a wide temperature (-10 – 50?C) and concentration (0 – 18 M) range, combining 17O shift, linewidth, and T1 relaxation measurements, 33S shift and linewidth measurements, and 1H PFG-NMR measurements of the proton self-diffusivity. In conjunction, these results indicate a crossover between two regimes of solvation structure and dynamics, occurring above the concentration associated with the deep eutectic point (~ 4.5 M), with the high-concentration regime dominated by strong water-sulfate correlation. This description was borne out in detail by the activation energy trends with increasing concentration derived from the relaxation of both the H2O/H3O+ and H2SO4/HSO4-/SO42- 17O resonances and the 1H self-diffusivity. However, the 17O chemical shift difference between the H2O/H3O+ and H2SO4/HSO4-/SO42- resonances across the entire temperature range is nevertheless strikingly linear. A computational approach coupling MD simulations and DFT NMR shift calculations to reproduce this trend is presented, which will be the subject of further development. This combination of multinuclear, dynamical NMR and computational methods, and the results furnished by this study, will provide a platform for future studies on battery electrolytes where aqueous sulfate chemistry plays a central role in the solution structure.

Bazak, Jonathan D.↗
DOE OSTI2023

DFT-Based Calculation of Molecular Hyperpolarizability and SFG Intensity of Symmetric and Asymmetric Stretch Modes of Alkyl Groups

Vibrational sum frequency generation (SFG) spectroscopy has been extensively used for obtaining structural information of molecular functional groups at two-dimensional (2D) interfaces buried in the gas or liquid medium. Although the SFG experiment can be done elegantly, interpreting the measured intensity in terms of molecular orientation with respect to the lab coordinate is quite complicated. One of the main reasons is the difficulty of determining the hyperpolarizability tensors of even simple molecules that govern their SFG responses. The single-bond polarizability derivative (SBPD) model has been proposed to estimate the relative magnitude of SFG-active hyperpolarizability by assuming that the perturbation associated to each vibration is negligible. In this study, density functional theory (DFT) was used to calculate the polarizability and dipole derivative tensors of the CH 3 stretch mode of CH 3 I, CH 3 CH 2 I, CH 3 OH, CH 3 CH 2 OH. Then, the hyperpolarizability tensors of symmetric and asymmetric vibration modes were calculated considering the Boltzmann distribution of representative conformers, which allowed theoretical calculation of their SFG intensities at all polarization combinations as a function of the tilt angle of the CH 3 group with respect to the surface normal direction. Then the ratios of the calculated SFG intensities for the CH 3 peaks used in experimental studies for the CH 3 tilt angle determination were compared. This comparison clearly showed the effect of vibrational coupling among neighboring functional groups. As a result, this study presents new parameters that can be used in determining the average tilt angle of the CH 3 group at the 2D interface with SFG measurements as well as limitations of the method.

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

Nature of Terminating Hydroxyl Groups and Intercalating Water in Ti 3 C 2 T x MXenes: A Study by 1 H Solid-State NMR and DFT Calculations

Since the discovery of two-dimensional transition metal carbides, referred to as MXenes, research efforts targeted their applications in energy storage, as lithium-ion batteries and supercapacitors. This interest is attributable to MXenes’ large volumetric capacitance, high rate handling capability and stable cycling performance, which largely rely on surface chemistry provided by the terminating groups, such as -OH, -O, and -F. However, the atomic-scale characterization of these surface terminations challenges the diffraction methods. Solid-state (SS)NMR spectroscopy, especially 1 H SSNMR, is a promising approach for scrutinizing the surface terminations and the intercalated water in an atomistic-scale, yet only a few SSNMR studies in MXenes have been reported to date offering conflicting results and limited understanding of -OH terminations. Here, we used 1 H SSNMR experiments in concert with the DFT calculations of NMR parameters to identify multiple types of -OH groups, residing on the external and internal surfaces, in a commonly studied MXene, Ti 3 C 2 T x . The study also identifies bulk-like water trapped between the MXene flakes and interfacial water stranded on the surface. Lastly, two-dimensional 1 H -1 H correlation spectra elucidated the water-surface interactions and the mechanism of de-intercalation of water upon annealing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗