Search NASA⌕ Search

SEARCH · Search NASA

Results for “density functional theory (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 OSTI2023

Understanding the Effect of Single Atom Cationic Defect Sites in an Al 2 O 3 (012) Surface on Altering Selenate and Sulfate Adsorption: An Ab Initio Study

Adsorption is a promising under-the-sink selenate remediation technique for distributed water systems. Recently it was shown that adsorption induced water network rearrangement control adsorption energetics on the α-Al 2 O 3 (012) surface. Here, we aim to elucidate the relative importance of the water network effects and surface cation identity on controlling selenate and sulfate adsorption energy using density functional theory calculations. Density functional theory (DFT) calculations predicted the adsorption energies of selenate and sulfate on nine transition metal cations (Sc–Cu) and two alkali metal cations (Ga and In) in the α-Al 2 O 3 (012) surface under simulated acidic and neutral pH conditions. We find that the water network effects had a larger impact on the adsorption energy than the cationic identity. However, cation identity secondarily controlled adsorption. Most cations decreased the adsorption energy, weakening the overall performance, the larger Sc and In cations enabled inner-sphere adsorption in acidic conditions because they relaxed outward from the surface, providing more space for adsorption. Additionally, only Ti induced Se selectivity over S by reducing the adsorbing selenate to selenite but not reducing the sulfate. Altogether, this study indicates that tuning water network structure will likely have a larger impact than tuning cation–selenate interactions for increasing adsorbate effectiveness.

Adsorption↗
NASA NTRS1995

Using Density Functional Theory (DFT) for the Calculation of Atomization Energies

The calculation of atomization energies using density functional theory (DFT), using the B3LYP hybrid functional, is reported. The sensitivity of the atomization energy to basis set is studied and compared with the coupled cluster singles and doubles approach with a perturbational estimate of the triples (CCSD(T)). Merging the B3LYP results with the G2(MP2) approach is also considered. It is found that replacing the geometry optimization and calculation of the zero-point energy by the analogous quantities computed using the B3LYP approach reduces the maximum error in the G2(MP2) approach. In addition to the 55 G2 atomization energies, some results for transition metal containing systems will also be presented.

Bauschlicher, Charles W., Jr.↗
DOE OSTI2025

Towards a Deeper Fundamental Understanding of (Al,Sc)N Ferroelectric Nitrides

Density functional theory (DFT) calculations, within the virtual crystal alloy approximation, are performed, along with the development of a Landau-type model employing a symmetry-allowed analytical expression of the internal energy and having parameters determined from first principles, to investigate properties and energetics of Al1-xScxN ferroelectric nitrides in their hexagonal forms. These DFT computations and this model predict the existence of two different types of minima, namely, the fourfold-coordinated wurtzite (WZ) polar structure and a five-fold coordinated paraelectric hexagonal phase (denoted as H5), for any Sc composition up to 40%. The H5 minimum progressively becomes the lowest-energy state within hexagonal symmetry as the Sc concentration increases from 0 to 0.4. Furthermore, the model points to several key findings. Examples include the crucial role of the coupling between polarization and strains to create the WZ minimum, in addition to polar and elastic energies, and that the origin of the H5 state overcoming the WZ phase as the global minimum within hexagonal symmetry when increasing the Sc composition mostly lies in the compositional dependency of only two parameters-one linked to the polarization and another one being purely elastic in nature. Other examples are that forcing Al1-xScxN systems to have no or a weak change in lattice parameters when heating them allows us to reproduce their finite-temperature polar properties well and that a value of the axial ratio close to that of the ideal WZ structure implies a large polarization at low temperatures but not necessarily at high temperatures because of the ordered-disordered character of the temperature-induced formation of the WZ state. Such findings should allow for a better fundamental understanding of (Al,Sc)N ferroelectric nitrides, which may be used to design efficient devices having, e.g., low operating voltages.

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

In Situ Studies of Methanol Decomposition Over Cu(111) and Cu2O/Cu(111): Effects of Reactant Pressure, Surface Morphology, and Hot Spots of Active Sites

The dissociative adsorption of methanol was investigated on Cu(111) and ultrathin Cu 2 O films. We employed synchrotron-based Ambient Pressure X-ray Photoelectron Spectroscopy (AP-XPS) and Scanning Tunneling Microscopy (STM) to study the dynamics of gas–solid interactions, and calculations based on Density Functional Theory (DFT) were used to examine the reaction path. C 1s XPS spectra revealed that methanol underwent dissociative adsorption on plain Cu(111) to form methoxy (CH 3 O), formaldehyde (H 2 CO), and formate (HCOO) at a pressure range of 0.5–10 mTorr, with these species remaining on the surface after evacuation. This was accompanied by the appearance of a low coverage (~0.05 ML) of O ads in the O 1s which can be considered a highly active site for methanol activation. The high activity is apparent by a coverage of 0.8 ML of methoxy at room temperature. STM was unable to image these species at room temperature as they were highly mobile on metallic copper. In contrast, for CH 3 OH on Cu 2 O/Cu(111), STM showed clear hot spots for reaction and a complex array of adsorption structures. On the oxide substrate, there was decomposition of methanol to H 2 CO, CH 3 O, HCOO, and hydrocarbon species (CH x ) due to the subsequent interactions of methanol with lattice oxygen. Cu(111) remained entirely saturated with decomposition products under 10 mTorr of methanol (θ ≈ 0.97 ML), whereas the Cu 2 O overlayer was saturated at a much lower coverage (θ ≈ 0.30 ML). STM revealed rows and step edges of Cu 2 O decorated with decomposition products and metallic Cu islands ~5 nm in size. The difference in activity between Cu(111) and Cu 2 O/Cu(111) is attributed to the significant amount of O present on the oxide surface. Additionally, Density Functional theory (DFT) calculations described the XPS measurements well, showing a likely methanol dissociation to *CH 3 O and therefore a surface reduction. More importantly, the DFT results revealed that it was the chemisorbed oxygen on Cu 2 O/Cu(111) which oxidized the dissociated *CH 3 O to *HCOO and eventually CO 2 , while the reaction only involving upper oxygen on the Cu 2 O hexagonal ring led to the formation of H 2 CO.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Band Energy Dependence of Defect Formation in the Topological Semimetal Cd3As2

Cadmium Arsenide (Cd3As2) is a prototypical Dirac semimetal that manifests topological properties in a 3D bulk material. In defect-free Cd3As2, the Fermi level lies at a minimum in the density of states at the Dirac point, but experimentally it forms with excess electron carriers and an elevated EF, thereby masking the topological features. To computationally study the self-doping of Cd3As2, we combine density functional theory (DFT) calculations for defect formation energies with quasi-particle self-consistent GW (QSGW) electronic structure calculations. We demonstrate an innate dependence of the point defect formation energies on carrier concentrations and use the QSGW calculated density of states to extrapolate formation energies to arbitrary electron concentrations. This approach allows the quantitative modeling of thermodynamic defect equilibria in topological semimetals and is used to predict how Cd3As2 growth conditions affect the position of EF relative to the Dirac point.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗
DOE OSTI2026

TiV_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys titanium-vanadium (Ti-V). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Ti and V. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2026

NbTi_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys niobium-titanium (Nb-Ti). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP), https://github.com/ORNL/ScalableWorkflow_VASP_Calculations. Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Nb and Ti. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2026

HfTi_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-titanium (Hf-Ti). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Ti. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2025

HfZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-zirconium (Hf-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

HfTa_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-tantalum (Hf-Ta). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Ta. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect BCC lattice sites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,040 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Origin of Unusual Acidity and Li+ Diffusivity in a Series of Water-in-Salt Electrolytes.

Superconcentrated aqueous electrolytes ("water-in-salt" electrolytes, or WiSEs) enable various aqueous battery chemistries beyond the voltage limits imposed by the Pourbaix diagram of water. However, their detailed structural and transport properties remain unexplored and could be better understood through added studies. Here, we report on our observations of strong acidity (pH 2.4) induced by lithium bis(trifluoromethane sulfonyl)imide (LiTFSI) at superconcentration (at 20 mol/kg). Multiple nuclear magnetic resonance (NMR) and pulsed-field gradient (PFG) diffusion NMR experiments, density functional theory (DFT) calculations, and molecular dynamics (MD) simulations reveal that such acidity originates from the formation of nanometric ion-rich structures. The experimental and simulation results indicate the separation of water-rich and ion-rich domains at salt concentrations >= 5 m and the acidity arising therefrom is due to deprotonation of water molecules in the ion-rich domains. As such, the ion-rich domain is composed of hydrophobic -CF3 (of TFSI-) and hydrophilic hydroxyl (OH-) groups. At 20 m concentration, the tortuosity and radius of water diffusion channels are estimated to be similar to 10 and similar to 1 nm, respectively, which are close to values obtained from hydrated Nafion membranes that also have hydrophobic polytetrafluoroethylene (PTFE) backbones and hydrophilic channels consisting of SO3- ion cluster networks providing for the transport of ions and water. Thus, we have discovered the structural similarity between WiSE and hydrated Nafion membranes on the nanometer scale.

Han, Kee Sung↗
DOE OSTI2020

Electronic Structures and Magnetism of Zr-, Th-, and U-based Metal-Organic Frameworks (MOFs) by Density Functional Theory

Metal-organic frameworks (MOFs) have recently gained wide interest as candidate materials for nuclear waste immobilization. While the fundamental thermodynamic properties, such as the substitution energies determine the favorability of radionuclide sequestration by utilization of a MOF matrix, the studies of MOF electronic structure reveal the role of d-, and/or f-electrons on changes in physical properties of actinide-containing materials. We use density functional theory (DFT) calculations to investigate the electronic structures of Zr-, Th-, and U-MOFs, including their electronic band structures and, where appropriate, their magnetic properties. We employ various DFT methods including DFT+U, collinear spin-polarization, spin-orbit coupling, and different flavors of exchange-correlation functionals to assess the robustness to the specific exchange-correlation functional. Unlike the Zr-, and Th-MOFs, the U-MOF is found to be sensitive to electron localization and spin; hence we explore the magnetic structure of the U-MOF in further detail.

Metal-Organic Framework, MOF, Actinides, DFT↗
DOE OSTI2021

Catalyst Design in Nitrate Removal

Based on the volcano plot developed by Dr. Goldsmith group (Report linked in submission), we utilized DFT (density functional theory) calculations to search for bimetallic materials in the application of catalysts in aqueous nitrate removal. The calculations are conducted via the high-throughput automated workflow package developed by our group (Github linked in submission) using VASP commercial first-principles calculation software.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Validation of the Cossee–Arlman mechanism for propylene oligomerization on Ni/UiO-66

Steady state rate expressions can be derived to distinguish the Cossee–Arlman and metallacycle mechanisms postulated for propylene oligomerization on nickel-based catalysts based on product selectivities, where product selectivities for the former are a function of olefin pressure because sequential coordination and insertion steps lead to independent mechanistic pathways for different hexene isomers. In contrast, the metallacycle mechanism presents pressure-independent product selectivities due to successive coordination prior to the kinetically relevant steps in each mechanism. In this work, steady state propylene oligomerization rates and selectivities were measured in the absence of an activator on nickel functionalized UiO-66 metal organic framework (MOF), Ni/UiO-66, to validate the Cossee–Arlman mechanism for light olefin oligomerization. In situ NO titrations reveal that ~5% of nickel sites were active during the reaction, and thus, not all nickel sites are relevant for catalysis. Propylene dimerization was first order in propylene pressure from 5 to 500 kPa with an apparent activation energy of ~20 kJ mol -1 from 453 to 493 K. Calculated apparent activation energies with density functional theory (DFT) calculations on cluster models of Ni/UiO-66 are in agreement with experiment to corroborate the Cossee–Arlman mechanism. Selectivities of hexene products and the ratio of hexene product selectivities on Ni/UiO-66 are in accordance with selectivity expressions derived from the Cossee–Arlman mechanism. In conclusion, analysis of product selectivities can be used more extensively to demarcate the Cossee–Arlman and metallacycle mechanisms for olefin oligomerization on metal-based catalysts.

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

Correct symmetry in density functional theory calculations of δ-Pu [Slides]

Density functional theory (DFT) calculations in the delta-phase of Pu have long relied on the antiferromagnetic (AFM) spin structure to approximate the material’s unique electronic structure. The resulting tetragonal distortion away from the correct cubic symmetry was believed to be a necessary consequence of allowing an ordered spin structure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Chromium enhances the mechanical performance of 3d transition metal high entropy alloys

A comprehensive examination of the compositional effects on the deformation behavior of high entropy alloys (HEAs) was conducted through nanoindentation, indentation creep, and stress relaxation experiments at ambient temperature. The alloys investigated included NiCoFe, NiCoCr, NiCoCrFe, and NiCoCrFeMn, all of which underwent identical thermomechanical processing. From the experimental results, NiCoCr exhibited the highest maximum shear stress for dislocation nucleation (13 GPa) and nanoindentation hardness (3.8 GPa) among the alloys tested. It also showed the lowest stress sensitivity (0.02) and activation volume (3–4 b 3 ) under creep and stress relaxation conditions, indicative of its superior plastic flow properties. Density functional theory (DFT) calculations further revealed an uneven charge density and larger bonding-length variation in NiCoCr due to the presence of Cr, which increased lattice distortion, impeding dislocation movement. Furthermore, in this study, post-deformation microscopy using transmission electron microscopy (TEM) revealed a high density of stacking faults and twins in NiCoCr that effectively enhanced its creep strength. The results highlight the significance of specific elemental effects, particularly from Cr in this study, over configurational entropy effect (i.e., compositional complexity) in governing the deformation microstructure and mechanical properties of HEAs. Finally, these insights will be instrumental for the design and optimization of advanced alloys for load-bearing applications.

36 - MATERIALS SCIENCE↗
DOE OSTI2021

High-pressure study of the low- Z rich superconductor Be 22 Re

With T c ~ 9.6K, Be 22 Re exhibits one of the highest critical temperatures among Be-rich compounds.We have carried out a series of high-pressure electrical resistivity measurements on this compound to 30 GPa. The data show that the critical temperature T c is suppressed gradually at a rate of dT c /dP = –0.05K/GPa. Using density functional theory (DFT) calculations of the electronic and phonon density of states (DOS) and the measured critical temperature, we estimate that the rapid increase in lattice stiffening in Be 22 Re overwhelms a moderate increase in the electron-ion interaction with pressure, resulting in the decrease in T c . Furthermore, high-pressure x-ray diffraction measurements show that the ambient pressure crystal structure of Be 22 Re persists to at least 154 GPa.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗