Search NASA⌕ Search

SEARCH · Search NASA

Results for “DFT calculations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Probing Active-Site Relocation in Cu/SSZ-13 SCR Catalysts during Hydrothermal Aging by in situ EPR Spectroscopy, Kinetic Studies, and DFT Calculations

Cu/SSZ-13 selective catalytic reduction catalysts display activity loss after mild hydrothermal aging without support degradation or loss of isolated Cu ion content. By applying electron paramagnetic resonance spectroscopy to compare hydrated, dehydrated and NH3-saturated catalysts before and after hydrothermal aging, as well as density functional theory calculations, it is found that hydrothermal aging induces Cu relocation, leading to stronger Cu-support interactions. This is particularly evidenced by ZCuOH conversion to Z2Cu, which leads to the generation of a unique structure with two Cu(II) ions positioned in a double 6-membered ring prism. The Cu(II) ions that do not undergo such dramatic changes also appear to stay closer to framework windows as evidenced by the overall increased anisotropy during EPR measurements. In situ EPR measurements under SCR conditions demonstrate that, in comparison to the fresh catalyst, higher percentages of Cu species stay as Cu(II) in the aged catalyst. This suggests that the reduction half-cycle of the SCR redox, i.e., Cu(II) ? Cu(I), becomes slower for the hydrothermally aged catalyst.

42 ENGINEERING↗
DOE OSTI2022

TTDFT: A GPU accelerated Tucker tensor DFT code for large-scale Kohn-Sham DFT calculations

We present the Tucker tensor DFT (TTDFT) code which uses a tensor-structured algorithm with graphic processing unit (GPU) acceleration for conducting ground-state DFT calculations on large-scale systems. The Tucker tensor DFT algorithm uses a localized Tucker tensor basis computed from an additive separable approximation to the Kohn-Sham Hamiltonian. The discrete Kohn-Sham problem is solved using Chebyshev filtered subspace iteration method that relies on matrix-matrix multiplications of a sparse symmetric Hamiltonian matrix and a dense wavefunction matrix, expressed in the localized Tucker tensor basis. These matrix-matrix multiplication operations, which constitute the most computationally intensive step of the solution procedure, are GPU accelerated providing ~8-fold GPU-CPU speedup for these operations on the largest systems studied. In conclusion, the computational performance of the TTDFT code is presented using benchmark studies on aluminum nano-particles and silicon quantum dots with system sizes ranging up to ~7,000 atoms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Locating anionic hydrogen in Ba 3 (Yb,Lu) 2 O 5 H 2 : A combined approach of X-ray diffraction, crystal chemistry, and DFT calculations

By a combination of x-ray diffraction, structural chemistry, and DFT calculations, the presence and location of anionic hydrogen in the two new, layered lanthanide oxyhydrides, Ba 3 Ln 2 O 5 H 2 (Ln ​= ​Yb, Lu) is inferred. Single crystals of the compounds have been synthesized from a molten barium flux with the addition of small amounts of BaH 2 . These phases crystallize in space group I4/mmm (#139, Z ​= ​2) with lattice parameters a ​= ​4.3336(2) Å and c ​= ​22.7197(6) Å, and a ​= ​4.3291(1) Å and c ​= ​22.597(1) Å, respectively. The Ba 3 Ln 2 O 5 H 2 phases comprise two different structural moieties: a perovskite double layer of stoichiometry Ba 2 Ln 2 O 5 H – formed by corner-connected LnO 5 tetragonal bi-pyramids with a terminating hydrogen anion, and a puckered rocksalt-type (BaH) + layer that is stretched along the c-axis. DFT calculations were used to arrive at hydrogen positions that minimize energy and are consistent with structural chemistry principles. Furthermore, the calculations show that the valence band edge is dominated by oxygen 2p orbitals with hydrogen 1s states admixed. The conduction band is formed by barium 5d-orbitals and Lu (Yb) 5d-orbitals. These are characteristics of materials with anionic H – . These new phases are isostructural with the Ba 3 Ln 2 O 5 Cl 2 (Ln ​= ​Gd–Lu) family of compounds with the chlorine atom in the same apical position as the hydrogen atom. Finally, steric effects limit the size of the lanthanide ion for Ba 3 Ln 2 O 5 H 2 .

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

Fast Correction of Errors in the DFT‐Calculated Energies of Gaseous Nitrogen‐Containing Species

Abstract Modeling adsorption phenomena on surfaces by DFT calculations often involves substantial errors, resulting in inaccurate predictions of catalytic activities. Such errors partly stem from the inaccurate description of the energetics of free molecules. Herein, we use a semiempirical group‐additivity method to correct the DFT‐calculated heats of formation of 106 carbon‐ and nitrogen‐containing gaseous compounds belonging to 15 different chemical families. PBE, PW91, RPBE and BEEF‐vdW initially yield mean absolute errors (MAEs) with respect to experiments in the range of 0.32–0.75 eV. After correcting the systematic errors, the overall MAEs decrease to ∼0.05 eV. Additionally, upon applying the corrections to three types of reaction enthalpies, the resulting MAEs are below 0.10 eV. These functional‐group corrections can be used in (electro)catalysis to correct the gas‐phase references necessary to evaluate equilibrium potentials and adsorption energies, predict error cancellation, and assess conflicting experimental data.

Urrego‐Ortiz, Ricardo↗
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 OSTI2020

Identifying the Molecular Edge Termination of Exfoliated Hexagonal Boron Nitride Nanosheets with Solid-State NMR Spectroscopy and Plane-Wave DFT Calculations

Hexagonal boron nitride nanosheets (h-BNNS), the isoelectronic analog to graphene, have received interest over the past decade due to their high thermal oxidative resistance, high bandgap, catalytic activity, and low cost. The functional groups that terminate boron and nitrogen zigzag and/or armchair edges directly affect their chemical, physical, and electronic properties. However, an understanding of the molecular edge termination present in h-BNNS is lacking. Here, high-resolution magic-angle spinning (MAS) solid-state NMR (SSNMR) spectroscopy, and plane-wave density-functional theory (DFT) calculations are used to determine the molecular edge termination in exfoliated h-BNNS. 1 H → 11 B cross-polarization MAS (CPMAS) SSNMR spectra of h-BNNS revealed multiple hydroxyl/oxygen coordinated boron edge sites that were not detectable in direct excitation experiments. A dynamic nuclear polarization (DNP)-enhanced 1 H → 15 N CPMAS spectrum of h-BNNS displayed four distinct 15 N resonances while a 2D 1 H{ 14 N} dipolar-HMQC spectrum acquired with fast MAS revealed three distinct 14 N environments. Plane-wave DFT calculations were used to construct model edge structures and predict the corresponding 11 B, 14 N and 15 N SSNMR spectra. Comparison of the experimental and predicted SSNMR spectra confirms that zigzag and armchair edges with both amine and boron hydroxide/oxide termination are present. The detailed characterization of h-BNNS molecular edge termination will prove useful for many material science applications. The techniques outlined here should also be applicable to understand the molecular edge terminations in other 2D materials

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

Mechanistic Study of 1,2-Dichloroethane Hydrodechlorination on Cu-Rich Pt–Cu Alloys: Combining Reaction Kinetics Experiments with DFT Calculations and Microkinetic Modeling

Cu-rich Pt–Cu bimetallic catalysts are among the most promising candidates for actively catalyzing the hydrodechlorination of 1,2-dichloroethane (1,2-DCA) toward ethylene production. Combining reaction kinetics experiments with density functional theory (DFT) calculations and mean-field microkinetic modeling, we present a systematic mechanistic study for 1,2-DCA hydrodechlorination on Cu-rich Pt–Cu alloy catalysts. Our DFT (PBE+(TS+SCS)) results suggest that increasing Cu content in the Pt–Cu alloy destabilizes C 2 -species adsorption while stabilizing the binding of atomic chlorine. The reaction energetics of all the elementary steps in the 1,2-DCA reaction network were calculated on a Pt 1 Cu 3 (111) model surface. The DFT results were then used to construct a microkinetic model, and the model-predicted reaction rates were compared with our reaction kinetics experimental results on a Cu-rich SiO 2 -supported Pt–Cu alloy catalyst through a parameter estimation procedure. Both the reaction kinetics experiments and the microkinetic model after parameter adjustments yielded 100% selectivity to ethylene. The microkinetic model pointed to a reaction pathway involving two sequential chlorine-removal steps on the Pt–Cu alloy catalyst, a mechanism distinct from the one previously identified on pure Pt/SiO 2 catalysts, which involved an initial hydrogen-removal step. Adjustments to the DFT-derived parameters indicate the possible formation of chlorine-induced Cu-enriched surface sites during 1,2-DCA hydrodechlorination conditions, sites that are more active than those encountered in the bulk Pt 1 Cu 3 (111) alloy surface. Furthermore, our study offers valuable initial insights on the 1,2-DCA hydrodechlorination reaction mechanism and the nature of the active sites on PtCu bimetallic catalysts.

1,2-Dichloroethane↗
DOE OSTI2026

pyRMG: A framework for high-throughput, large-cell DFT calculations on supercomputers

Exascale computing delivers the raw power to simulate ever larger and more chemically realistic systems, but realizing this potential requires codes that can efficiently use thousands of processors. Our real-space multigrid (RMG) density functional theory (DFT) code’s grid-decomposition approach scales nearly linearly with the number of graphics processing units (GPUs), even for simulations exceeding thousands of atoms. This scalability makes RMG a compelling tool for high-throughput DFT studies of materials that would otherwise be bottlenecked in other codes (for example, by global fast Fourier transforms in plane-wave DFT). However, the limited workflow infrastructure for RMG has thus far constrained its adoption to a small user community. In this work, we present pyRMG, a Python package designed to streamline the setup and execution of RMG DFT calculations. Built on the pymatgen and ASE (Atomic Simulation Environment) computational materials science Python packages, pyRMG automates input generation and convergence checking, and it integrates with modern job schedulers (e.g., Flux) on leadership-class platforms such as Frontier and Perlmutter. Here, we demonstrate pyRMG for a high-throughput study of strain effects in 2D 2L-Bi 2 Se 3 /2L-NbSe 2 heterostructures, which offers chemical insights into this system and shows that RMG-based workflows can converge with limited user intervention.

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

Investigation of the High-Pressure Behaviors of Amblygonite by Single-Crystal X-ray Diffraction, Raman Spectroscopy, and DFT Calculations

In the present study, we extensively explored the high-pressure behaviors and vibrational properties of amblygonite LiAlPO 4 F with elevated pressures up to 34.3 GPa based on single-crystal X-ray diffraction measurements, Raman spectroscopy, and DFT calculations. The compressibility and elastic properties of amblygonite were determined first. Specifically, the obtained isothermal bulk modulus of LiAlPO 4 F is 128(4) GPa and the triclinic phase exhibited anisotropic compression with axial compressibility β c > β a > β b with a ratio of 1.11:1.00:1.20. The Raman spectra showed no indication of phase transformation and were used to obtained mode Grüneisen parameters. The average Grüneisen parameter for PO 4 tetrahedral sites was smaller than for the LiO 4 F sites. Our results provide new insights into the phase stability and elastic properties of lithium-fluorite granites at extreme conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

DFT calculations and theory do not support enantiospecificity in NMR J-coupling constants

Distinguishing the enantiomers of small organic molecules is an industrially relevant problem with important implications for the health of the population. In a recent publication, Bouchard and co-workers have suggested that large differences in indirect spin-spin (J) coupling constants between enantiomers are possible. A close inspection of their work revealed significant flaws in their density functional theory (DFT) calculations and that the reported effects disappear with appropriate care. We thus conclude that enantiospecificity in spin-spin coupling constants has not be demonstrated either experimentally or theoretically.

Perras, Frédéric A. [Ames Laboratory (AMES), Ames,↗
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 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↗
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 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 OSTI2022

In-situ synthesized N-doped ZnO for enhanced CO 2 sensing: Experiments and DFT calculations

Chemiresistive CO 2 sensing is attractive due to low cost and ease of chip-level integration. Our previous studies (Yong Xia, 2021) showed the well-developed ZnO material fabricated by in-situ annealing exhibited good CO 2 sensing performance. Here, we have expanded on those studies, including CO 2 cyclic tests under both dry air and N 2 background whereby a much higher response to CO 2 in N 2 background was observed. In this study, detailed density functional theory calculations were conducted to understand the behavior. The results indicated nitrogen doping is mainly responsible for the observed response. In the presence of pre-adsorbed O 2 , N-doped ZnO can no longer interact with CO 2 , which agrees well with the observation of higher response in N 2 background. Furthermore, density of states analysis showed N sp 2 hybridized orbital and N 2p orbital of the N dopant mixed with sp 2 hybridized orbital of C atom and 2p orbitals of C/O atoms in CO 2 to form σ and π bonds, respectively. However, they mixed with O 2s/2p orbitals of O atom in O 2 when pre-adsorbed O 2 was present, hindering CO 2 interaction with N-doped ZnO, and resulting in limited response in air. The illustrated mechanism does not only further the understanding of metal oxide-based CO 2 sensing, but also guide the design of new functional materials for CO 2 sensing or capture.

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

DFT Calculation of Thermal Expansion within Debye-Grüneisen Framework Made Simple and Accurate

The density functional theory (DFT) based Debye-Grüneisen approach is a computational method used to calculate the thermodynamic properties of materials by relating the volume dependence of the Debye frequency to the Grüneisen parameter. In this study, we propose a new model to describe the volume dependence of the Grüneisen parameter. This model offers a simple yet accurate approach to predict the thermal expansion of materials within the framework of the Debye-Grüneisen theory. Specifically, for pure elements or species used as a reference system, we calibrated the parameters to reproduce the experimental coefficient of thermal expansion, heat capacity, and entropy. For compounds with arbitrary compositions, the parameters are determined by taking the arithmetic or geometric average of those obtained from the pure elements or species based on their composition. We then compared the model prediction with experiments. This allowed for a comprehensive characterization of the thermal expansion behavior of various materials.

Wang, Yi↗
DOE OSTI2020

Experimental and DFT Calculated IR Spectra of Guests in Zeolites: Acyclic Olefins and Host-Guest Interactions

We performed experimental and periodic density functional theory (DFT) IR spectroscopy to investigate the adsorption of acyclic olefins over both acidic and non-acidic zeolites. Two conjugated polyenes, 2,4-dimethyl-1,3-pentadiene (I) and 2,6-dimethyl-2,4,6-octatriene (II) were studied to probe organic intermediates that can form during methanol conversion, and that can lead to deactivating species known collectively as “coke.” We computed vibrational spectra using zeolite-adsorbed and gas-phase models for both neutral and protonated forms of I and II, and compared these DFT results to diffuse reflectance IR Fourier transform (DRIFT) spectra of zeolite-guest systems. Our experimental and computational results are precise enough to pinpoint that the gauche s-cis conformation of species I predominates during adsorption over de-aluminated zeolite beta. Computed zeolite-adsorbed spectra of the protonated species I and II best represent the DRIFT spectra obtained after the adsorption of the olefins on HMOR at 20 ?, with computed bands at 1543 and 1562 cm–1 for molecules I+ and II+, respectively, attributed to the allylic stretching mode, ?(C=C–C+). These computed band frequencies are within 6 cm–1 of experimental data and confirm that the interaction between neutral acyclic olefins and acidic zeolites leads to protonation of the olefin. Comparison of computed spectra of the protonated species in the gas phase to those in the zeolite indicates that the electrostatic interaction between alkenyl and alkadienyl cations and negative zeolite framework does not significantly impact the position of the allylic stretching bands. These results highlight that computed spectroscopy and thermodynamics coupled with experimental spectra can be used for elucidating complex mixtures in zeolites, and certain spectral features of adsorbed olefins can be accurately modeled by gas-phase calculations. ACKNOWLEDGMENTS BM and SMA were supported by the National Science Foundation under Award # CBET-1512442. Acknowledgment is made to the Donors of the American Chemical Society Petroleum Research Fund for support (or partial support) of this research (#589978-ND5, EDH and FCJ). MDB was supported by the US Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences (BES), Division of Materials Science and Engineering. CJM was supported by the DOE, Office of Science, BES, Division of Chemical Sciences, Geosciences, and Biosciences. Pacific Northwest National Laboratory (PNNL) is a multiprogram national laboratory operated by Battelle for the U.S. DOE. Computing resources were generously provided through PNNL’s institutional computing (PIC).

Manookian, Babgen↗