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

Nonadiabatic couplings from a variational excited state method based on constrained DFT

Excited Costrained Density Functional Theory (XCDFT) is a variational excited state method that extends ground state DFT to the computation of low-lying excited states. It borrows much of the machinery of Constrained DFT (CDFT) with a crucial difference: the constraint imposes a population of one electron in the Hilbert space spanned by the virtuals of a reference ground state. Here, we present theory and implementation for evaluating nonadiabatic coupling vectors (NACVs) between the first excited state computed with XCDFT and the ground state. Our NACVs are computed analytically using density functional perturbation theory with a formalism that is general enough that could be applied to CDFT diabatic states. We showcase the new method with pilot NACV calculations for the conical intersection in H 3 , the avoided crossing in selenoacrolein, and the NACV magnitudes in azobenzene. Despite complications from the nonorthogonality of the wavefunctions, XCDFT’s energy surfaces and NACVs reproduce benchmark values and respect known sum rules within a reasonable degree. This shows that XCDFT is a viable method for nonadiabatic dynamics simulations.

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

Subspace recursive Fermi-operator expansion strategies for large-scale DFT eigenvalue problems on HPC architectures

Quantum mechanical calculations for material modeling using Kohn–Sham density functional theory (DFT) involve the solution of a nonlinear eigenvalue problem for N smallest eigenvector-eigenvalue pairs, with N proportional to the number of electrons in the material system. Here, these calculations are computationally demanding and have asymptotic cubic scaling complexity with the number of electrons. Large-scale matrix eigenvalue problems arising from the discretization of the Kohn–Sham DFT equations employing a systematically convergent basis traditionally rely on iterative orthogonal projection methods, which are shown to be computationally efficient and scalable on massively parallel computing architectures. However, as the size of the material system increases, these methods are known to incur dominant computational costs through the Rayleigh–Ritz projection step of the discretized Kohn–Sham Hamiltonian matrix and the subsequent subspace diagonalization of the projected matrix. This work explores the potential of polynomial expansion approaches based on recursive Fermi-operator expansion as an alternative to the subspace diagonalization of the projected Hamiltonian matrix to reduce the computational cost. Subsequently, we perform a detailed comparison of various recursive polynomial expansion approaches to the traditional approach of explicit diagonalization on both multi-node central processing unit and graphics processing unit architectures and assess their relative performance in terms of accuracy, computational efficiency, scaling behavior, and energy efficiency.

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

A DFT-based kinetic Monte Carlo simulation of multiphase oxide-metal thin film growth

Functional thin films of nanoscale metal pillars in oxide or nitride matrices known as vertically aligned nanocomposite (VAN) have gained much interest owing to their unique strain-coupled and highly anisotropic properties. So far, the deposition of these films has been explored mostly experimentally. In this work, a density functional theory (DFT)-based kinetic Monte Carlo simulation model using Bortz–Kalos–Lebowitz algorithm was developed to understand the growth of VAN films deposited by pulsed laser technique on mismatching substrates. The model has been parameterized and applied to understand the kinetics of growth thin films consisting of Au pillars in CeO2 matrix deposited on SrTiO3 substrates. The effects of pulsed laser deposition (PLD) conditions including the pulse frequency, deposition flux, and substrate temperature were explored. The simulations indicate that the Au pillar size and shape exhibit significant dependence on the PLD conditions. Namely, increasing the temperature increases the average pillar size and lowers the pillar density, and vice versa. In addition, the simulations revealed that increasing the deposition rate results in lowering the average pillar size and increasing the density. Particularly, the DFT results suggest that Au pillar size can be tuned during the initial growth of the first monolayer due to the significantly low activation barrier. Our analysis showed that the relationship between the average pillar size and pillar density is influenced by the kinetics. Furthermore, autocorrelation analysis showed that pillars self-organize in quasi-ordered patterns at certain windows of the deposition conditions, which is attributed to the complex nature of the chemical interactions in the system, the kinetics, and the deposition parameters.

Physics↗
DOE OSTI2025

Xenon–metal pair formation in UO 2 investigated using DFT + U

A recent experimental study on a spent uranium dioxide (UO 2 ) fuel sample from Belgium Reactor 3 identified a unique pair structure formed by the noble metal phase (NMP) and fission gas [xenon (Xe)] precipitate. However, the fundamental mechanism behind this structure remains unclear. The present study aims to provide an understanding of the interaction between five different metal precipitates [molybdenum (Mo), ruthenium (Ru), palladium (Pd), technetium (Tc), and rhodium (Rh)] and the Xe fission gas atoms in UO 2 , by using density functional theory (DFT) in combination with the Hubbard U correction to compute the formation energies involved. All DFT + U calculations were performed with occupation matrix control to ensure antiferromagnetic ordering of UO 2 . The calculated formation and binding energies of the Xe and solid fission products in the NMP reveal that these metal precipitates form stable pair structures with Xe. Notably, the formation energy of Xe–metal pairs is lower than that of the isolated single defects in all instances, with Pd and Mo showing the most favorable binding energy, likely accounting for the observed pair structure formation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

How the choice of exchange–correlation functional affects DFT-based simulations of the hydrated electron

Hydrated electrons are anionic species that are formed when an excess electron is introduced into liquid water. Building an understanding of how hydrated electrons behave in solution has been a long-standing effort of simulation methods, of which density functional theory (DFT) has come to the fore in recent years. The ability of DFT to model the reactive chemistry of hydrated electrons is an attractive advantage over semi-classical methodologies; however, relatively few density functional approximations (DFAs) have been used for the hydrated electron simulations presented in the literature. Here, we simulate hydrated electron systems using a series of exchange–correlation (XC) functionals spanning Jacob’s ladder. We calculate a variety of experimental and other observables of the hydrated electron and compare the XC functional dependence for each quantity. We find that the formation of a stable localized hydrated electron is not necessarily limited to hybrid XC functionals and that some hybrid functionals produce delocalized hydrated electrons or electrons that react with the surrounding water at an unphysically fast rate. Here, we further characterize how different DFAs impact the solvent structure and predicted spectroscopy of the hydrated electron, considering several methods for calculating the hydrated electron’s absorption spectrum for the best comparison between structures generated using different density functionals. None of the dozen or so DFAs that we investigated are able to correctly predict the hydrated electron’s spectroscopy, vertical detachment energy, or molar solvation volume.

Ab-initio molecular dynamics↗
DOE OSTI2024

Point defects in CdTe and CdSeTe alloy: a first principles investigation with DFT+U

CdTe and its alloy CdSeTe are widely used in optoelectronic devices, such as radiation detectors and solar cells, due to their superior electrical properties. However, the formation of defects and defect complexes in these materials can significantly affect their performance. As a result, understanding the defect formation and recombination processes in CdTe and CdSeTe alloy is of great importance. In recent years, density functional theory (DFT) calculations have emerged as a powerful tool for investigating the properties of defects in semiconductors. In this paper, we use DFT+U calculations to comprehensively study the properties of intrinsic defects as well as extrinsic defects induced by commonly used dopants, such as Cu and group V elements, in CdTe and CdSeTe alloy. This work provides insights into the effects of these defects on the electrical and optical properties of the material.

14 SOLAR ENERGY↗
DOE OSTI2025

Collection: TD-DFT and EOM-CCSD Calculations for the GDB-9-Ex Dataset

We present two datasets that contain quantum chemical electronic structure calculations for organic molecules from the GDB-9-Ex dataset. The “GDB-9-Ex_TD-DFT-PBE0” dataset contains calculations performed using the time-dependent density functional theory (TD-DFT) first principles method, and the “GDB-9-Ex_EOMCCSD” dataset contains calculations performed using the equation-of-motion coupled cluster (EOM-CCSD) method. Both types of calculations were performed using the ORCA software and provided ultraviolet-visible spectra with a high level of accuracy.

Mehta, Kshitij [Oak Ridge National Laboratory (ORN↗
DOE OSTI2020

DFT-Based Calculation of Dissolution Activation Energy and Kinetics of Ni–Cr Alloys

A density functional theory investigation of the process of dissolution of Ni and Cr adatoms from model Ni–Cr(111) alloy surfaces is presented, both in vacuum and with explicit water molecules. The goal is to understand how the electronic structures solved using DFT can provide insights as to changes in valence, energy and coordination of the adatoms in the process of dissolving from the alloy surface. It is found that nearby Cr solute atoms increase the dissolution activation energy of Ni. Cr adatoms have a similar dissolution activation energy as Ni adatoms, except for one particular surface configuration that has a much smaller dissolution activation energy, which might promote selective dissolution of Cr in Ni–Cr alloys. To interpret the first-principles modeling results for the potential energy trajectory along the dissolution reaction coordinate, we provide a thermodynamic breakdown of the various terms that contribute to the dissolution activation and total reaction energy. The role of addition of Cr is explored on the dissolution kinetics of a Ni–Cr alloys of varying composition by constructing a model for the corrosion current density with DFT-based activation energies. Lastly, dissolution resistance index is proposed as a quantifiable descriptor of the dissolution resistance of corrosion resistant alloys.

36 MATERIALS SCIENCE↗
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.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A Filon-like integration strategy for calculating exact exchange in periodic boundary conditions: a plane-wave DFT implementation

An efficient and accurate approach for calculating exact exchange and other two-electron integrals has been developed for periodic electronic structure methods. Traditional approaches used for integrating over the Brillouin zone in band structure calculations, e.g. trapezoidal or Monkhorst-Pack, are not accurate enough for two-electron integrals. This is because their integrands contain multiple singularities over the double integration of the Brillouin zone, which with simple integration methods lead to very inaccurate results. A common approach to this problem has been to replace the Coulomb interaction with a screened Coulomb interaction that removes singularities from the integrands in the two-electron integrals, albeit at the inelegance of having to introduce a screening factor which must precomputed or guessed. Instead of introducing screened Coulomb interactions in an ad hoc way, the method developed in this work derives an effective screened potential using a Filon-like integration approach that is based only on the lattice parameters. This approach overcomes the limitations of traditionally defined screened Coulomb interactions for calculating two-electron integrals, and makes chemistry many-body calculations tractable in periodic boundary conditions. This method has been applied to several systems for which conventional DFT methods do not work well, including the reaction pathways for the addition of H 2 to phenol and Au$^{–}_{20}$ nanoparticle, and the electron transfer of a charge trapped state in the Fe(II) containing mica, annite.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Physi-sorption of H 2 on pure and boron-doped graphene monolayers: A dispersion-corrected DFT study

Carbon based materials are of interest as potential candidates for H 2 storage. Earlier work has been inconclusive on the effect of boron doping on the energy of H 2 binding. However, earlier work has been inconclusive on the definitive effect of boron doping on the energy of H 2 binding, i.e. isosteric heats of adsorption (Qst). In this work, we completed a systematic DFT study to evaluate this effect, and found that doping graphene with boron provides only minor enhancement in H 2 binding. More importantly, the presence of the electron deficient boron into a graphene ring introduces a defect, such as terminal hydrogen or distortion from planarity, which creates hydrogen adsorption sites with slightly increased Qst. The increase is from ~ 5 kJ/mol H 2 for the pure carbon matrix to ~6 - 7 kJ/mol for the boron doped system. The more strongly bond H 2 is located near the defect and shows little direct interaction with the boron. Most significant enhancement is found in systems where H 2 is confined between layers at a distance of about 7Å. In this case, the H 2 binding nearly doubles, to ~10 kJ/mol for both pure, undistorted graphene, and the 2% boron doped system, containing tetrahedral C atom distortion. Interestingly, at higher doping levels of boron, we found that the Qst decreases compared to the all carbon materials. These finding suggest that interplanar nanoconfinement may be more effective way to enhance H 2 binding than boron doping.

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

The potential application of VS 2 as an electrode material for Mg ion battery: A DFT study

Herein, by means of first-principle calculations based on density functional theory (DFT), the electrochemical properties of monolayer VS 2 (M-VS 2 ), double-layer VS 2 (D-VS 2 ) and bulk VS 2 (B-VS 2 ) as electrode materials for Mg-ion batteries (MIBs) were comprehensively explored. The computation results reveal that Mg atom can strongly bind with the three different forms of VS 2 . All of the Mg adsorbed VS 2 systems demonstrate metallic characteristics, which indicates a good electronic conductivity. In addition, crystal orbital hamiltonian population shows that the stability of V-S bond is weakened after adsorption of Mg atom. The low diffusion barriers of Mg give rise to the high rate performance of VS 2 in MIBs. More interestingly, the three types of VS 2 display same storage ability for Mg cations, which can adsorb 0.5 Mg atoms for VS 2 , producing maximum theoretical capacity 233 mA h g -1 for MIBs. Furthermore, the average working voltages results suggest that M-VS 2 can be employed as anode materials, while D-VS 2 , B-VS 2 can be used as cathode materials for MIBs.

25 ENERGY STORAGE↗
DOE OSTI2022

Acetate formation on metals via CH 4 carboxylation by CO 2 : A DFT study

Here, density functional theory calculations have been performed to investigate CH 4 activation and coupling to CO 2 forming C 2 carboxylates such as acetate on the close-packed (111) or (0001) surfaces of ten mid-to-late transition and coinage metals (Co, Ni, Cu, Ru, Rh, Pd, Ag, Ir, Pt, and Au). Consistent with the literature, the activation energy (E a ) for the initial C-H bond scission in CH 4 is mild, being ca. 1 eV or less on all but the coinage metals, of which Ag exhibits the highest E a at 2.13 eV, followed by Au and Cu. E a for the CH 3 -CO 2 coupling step is 0.8 ~ 1.1 eV on Co, Ru, Rh, and Ag, 1.2 ~ 1.5 eV on Ni, Cu, Pd, and Ir, and 1.8 ~ 2.1 eV on Pt and Au. While the two E a are comparable for several metals in terms of DFT total energies, free energy analysis indicates CH 3 -CO 2 coupling to be much more rate-limiting than CH 4 activation. Overcoming it would require over 800 K even on the most active of the metals considered, Ru, which makes the formation of acetate not feasible on the monometallic metal surfaces. Instead, we propose that single atom alloys based on early transition metals doped into a host metal such as Ni(111) could be viable catalysts. The dopant sites serve to stabilize the transition state of C-C coupling while Ni sites continue to activate CH 4 , thereby significantly lowering the required temperature for acetate formation.

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

Screening of transition metals for PFAS adsorption: A comparative DFT investigation

Transition metals are promising candidates for catalytic destruction of per- and polyfluoroalkyl substances (PFAS). This study employed density functional theory (DFT) calculations to systematically investigate the adsorption of PFAS, which is a critical step in the catalytic process, onto a group of transition metals, namely Cu, Pd, Pt, and Rh. The investigation considered the influence of different PFAS characteristics, including chain length, functional group (i.e., sulfonic vs. carboxylic), and protonation state (i.e., protonated, anionic, and deprotonated radical) on the adsorption mechanism. Overall, the findings indicated that the adsorption of PFAS on transition metals exhibited thermodynamically favorable energetics, with Rh showing the most favorable adsorption for PFAS. However, the adsorption strength was significantly dependent on the protonation state of PFAS. Cu was the highest electron donor, while Pt was the highest electron acceptor. Results were compared to experimental studies in the literature for qualitative validation.

Adsorption↗
DOE OSTI2022

Degradation and fate of 2,4-dinitroanisole (DNAN) and its intermediates treated with Mg/Cu bimetal: Surface examination with XAS, DFT, and LDI-MS

In this article, a novel Mg-based bimetal reagent (Mg/Cu) was used as an enhanced reductive system to degrade insensitive munition 2,4-dinitroanisole (DNAN), a contaminant found in energetic-laden waste. Degradation of DNAN was significantly impacted by dissolved oxygen and studied in anoxic and oxic bimetal systems (i.e., purging with N 2 , air, or O 2 gas). Degradation occurred through sequential nitroreduction: first one nitro group was reduced (ortho or para) to form short-lived intermediates 2-amino-4-nitroanisole or 4-amino-2-nitroanisole (2-ANAN or 4-ANAN), and then subsequent reduction of the other nitro group formed 2,4-diaminoanisole (DAAN). The nitro-amino intermediates demonstrated regioselective reduction in the ortho position to 2-ANAN; Regioselectivity was also impacted by the anoxic/oxic environment. Under O 2 -purging DNAN degradation rate was slightly enhanced, but most notably O 2 significantly accelerated DAAN generation. DAAN also further degraded only in the oxygenated Mg/Cu system. Adsorption of DNAN byproducts to the reagent occurred regardless of anoxic/oxic condition, resulting in a partition of carbon mass between the adsorbed phase (27%-35%) and dissolved phase (59%-72%). Additional surface techniques were applied to investigate contaminant interaction with Cu. Density functional theory (DFT) calculations identified preferential adsorption structures for DNAN on Cu with binding through two O atoms of one or both nitro groups. X-ray absorption spectroscopy (XAS) measurements determined the oxidation state of catalytic metal Cu and formation of a Cu-O-N bond during treatment. Laser desorption ionization mass spectrometry (LDI-MS) measurements also identified intermediate 2-ANAN adsorbed to the bimetal surface.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

DFT-based investigation of solvation of Nd(III)/Yb(III) cations in room-temperature ionic-liquid

We report a DFT study on the solvation of bare rare-earth (RE) cations [Nd(III) and Yb(III)] in the room-temperature ionic-liquid (RTIL) based on 1,3-dimethyl-imidazolium + , PF 6 - , [MMI][PF 6 ]. In the present article, we aim to investigate whether the chosen RTIL can preferentially solvate any of the RE cations under study. First, the RTIL binding to the RE cation is studied in the gas phase, and then explored with continuum solvation model employing RTIL parameters. First, we considered the formation of the first solvation shell, where the PF 6 - anions directly bond with the Ln(III) cation. We then further surrounded the first solvation shell with [MMI]+ cations to form the second solvation shell. The stability of PF 6 - coordinated complexes is studied by calculating the binding energy values. The number of PF 6 - anions for both Ln(III) cations starts with three ligands and goes until the maximum binding energy in the RTIL medium is reached. The coordination number and binding modes of PF 6 - depend on the size of the RE cation. Comparing the binding energy values between the RE cations for the first solvation shell complexes, it is noticed that Yb-complexes have relatively greater binding energy than Nd-complexes. We have further performed a Boltzmann population analysis to find out the most dominant conformer among each category of PF 6 - coordinated complexes. A study using Natural Population Analysis shows that after complexation with PF 6 - the charge on the Yb-centers lower to a greater extent than Nd-centers, implying stronger bonding between Yb and PF 6 - , which is in harmony with the binding energy analysis. Overall, the computational investigation provides an in-depth understanding of the complexation phenomena of RE cations in RTIL, which can be useful for the practical separation of RE cations using RTIL.

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

The effects of vanadium substitution on one-dimensional tunnel structures of cryptomelane: Combined TEM and DFT study

Rational design of energy storage materials requires mechanistic understanding which relies on atomic scale characterization of structure and defects. In this paper we utilize vanadium (V)-substituted one-dimensional cryptomelane as a model material to study structural features critical to electrochemical reversibility of the cathode in lithium-ion batteries. We found the substitution of manganese (Mn) atoms in the tunnel wall with V enhances structural uniformity in the parent KMn 8 O 16 material by inhibiting formation of irregular (2 × 3 and 2 × 1) manganese oxide tunnel structures. Upon lithiation, both transmission electron microscopy (TEM) and density functional theory (DFT) reveal significant structural evolution, eventually resulting in phase separation and formation of highly reduced manganese oxide. During this process, Mn 4+ accepts electrons donated by intercalated Li, while the V 4+ substituent is not electrochemically reduced but acts to stabilize the tunnel framework structure. As a result, the V substituted material shows higher operating voltage and improved electrochemical reversibility. Our results highlight the benefits of V-substitution in cryptomelane, and the significance of integrating characterization with simulation for understanding and developing design principles for lithium ion battery materials.

25 ENERGY STORAGE↗
DOE OSTI2022

Elucidating facet dependent electronic and electrochemical properties of Cu 2 O nanocrystals using AFM/SCEM and DFT

Cuprous oxide (Cu 2 O) has extensively been studied owing to its excellent optical, magnetic, and catalytic properties. Many of these properties are facet-dependent and have not been well elucidated. This work synthesized cubic, cuboctahedral, octahedral, and rhombic dodecahedral shaped Cu 2 O nanocrystals of ~300 nm in size to evaluate the facet-dependent electrochemical activities. Here, cyclic voltammetry (CV) and electrochemical impedance spectroscopy (EIS) were firstly used to reveal the average electrochemical activities at the ensemble level. Atomic force microscopy-scanning electrochemical microscopy (AFM-SECM) was further used to assess the electrochemical activities of different Cu 2 O nanocrystals at the facet level. Hexaammineruthenium (III) chloride ({Ru(NH 3 ) 6 }Cl 3 ) was employed as the probe molecules that reacted with four different Cu 2 O nanocrystals under —400mV and yielded ~300 pA current between the probing tip and the nanocrystal surface. The tip-current mapping results indicate that rhombic dodecahedral Cu 2 O exhibits higher electrocatalytic activity than other shaped Cu 2 O, due to the presence of dominant exposed facet of {110} as indicated by the relatively high tip current. Density-functional theory (DFT) calculations confirmed the facet dependence of local surface energy and electronic structure of Cu 2 O nanocrystals. Besides electrochemical activity, the surface work function and adsorptive properties were both observed to vary with the shape and dominant exposed facets of Cu 2 O. This study presented a unique experimental and computational chemistry approach to analyze surface electrochemical properties of Cu 2 O crystals at a crystalline facet level.

36 MATERIALS SCIENCE↗