Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density Functional Theory”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Machine learning-accelerated path integral molecular dynamics simulations of reactive organic electrolytes

Hydrogen bonded electrolytes that exhibit accelerated proton transport via sequential reactive hops have drawn interest for their promise in clean energy applications. Molecular dynamics simulations of these electrolytes offer the opportunity to uncover microscopic mechanistic details that could be used to design and tune the properties of candidate electrolyte technologies. However, accurately modeling the proton transfer reactions and transport properties that give rise to high charge conductivites in these electrolytes proves computationally challenging because of the need to perform lengthy condensed phase simulations, treating both the electronic and nuclear degrees of freedom quantum mechanically. In this paper, we demonstrate that such a modeling task can be efficiently achieved with the use of density functional theory (DFT)-trained machine learning potentials (MLP) to accelerate path integral molecular dynamics (PIMD) simulations. We highlight the practical utility of this approach by using it to benchmark how closely PIMD simulations employing different DFT exchange–correlation functionals reproduce the composition-dependent densities, diffusion coefficients, and electrical conductivities of mixtures consisting of imidazole and levulinic acid. Even with the speedup afforded by our MLPs, PIMD simulations remain quite expensive. Furthermore, in order to render PIMD more computationally tractable, we introduce and benchmark the accuracy of a ring polymer contraction approach that leverages a computationally efficient short-range MLP to accelerate our PIMD simulations by an additional factor of four.

Chemical bonding↗
DOE OSTI2024

Solvation Environment and Interface Dynamics of Li 2 B 12 H 12 and Li 2 B 12 F 12 Electrolytes Uncovered by Theory and Operando Optical and FTIR Spectroelectrochemistry

In this work, we evaluated two closo-borate salts (Li 2 B 12 H 12 and Li 2 B 12 F 12 ) in propylene carbonate from theoretical and experimental perspectives to understand how the coordination environment influences their spectroscopic and electrochemical properties. The coordination environments of the closo-borate salts were modeled via density functional theory (DFT) and molecular dynamics (MD). Vibrational spectra calculated from the predicted coordination environments are in agreement with experimentally measured steady-state FTIR data. This theoretical investigation also suggested that Li 2 B 12 F 12 would possess a higher ionic conductivity than Li 2 B 12 H 12 , which was corroborated experimentally. Additionally, an electrochemical cell was designed and fabricated that enabled operando optical and FTIR spectroelectrochemical (OP-IR-SEC) measurements. This allowed for the simultaneous measurement of the relative changes of species at a lithium electrode-liquid electrolyte interface and the visualization of lithium plating at the electrode surface. This technique could provide new chemical insights and potentially link optical changes at the electrode–electrolyte interface to specific chemical species in similar electrochemical systems. The Li 2 B 12 F 12 electrolyte was found to have a higher thermal stability, which may find utility in applications for batteries that are subject to high-temperature conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Low-temperature etching of silicon oxide and silicon nitride with hydrogen fluoride

Etching of high aspect ratio features into alternating SiO 2 and SiN layers is an enabling technology for the manufacturing of 3D NAND flash memories. In this paper, we study a low-temperature or cryo plasma etch process, which utilizes HF gas together with other gas additives. Compared with a low-temperature process that uses separate fluorine and hydrogen gases, the etching rate of the SiO 2 /SiN stack doubles. Both materials etch faster with this so-called second generation cryo etch process. Pure HF plasma enhances the SiN etching rate, while SiO 2 requires an additional fluorine source such as PF 3 to etch meaningfully. Further, the insertion of H 2 O plasma steps into the second generation cryo etch process boosts the SiN etching rate by a factor of 2.4, while SiO 2 etches only 1.3 times faster. We observe a rate enhancing effect of H 2 O coadsorption in thermal etching experiments of SiN with HF. Ammonium fluorosilicate (AFS) plays a salient role in etching of SiN with HF with and without plasma. AFS appears weakened in the presence of H 2 O. Density functional theory calculations confirm the reduction of the bonding energy when NH 4 F in AFS is replaced by H 2 O.

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

Weak decays in superheavy nuclei

Superheavy nuclei represent the extreme atoms and nuclides known at the limit of mass and charge. The observed superheavy nuclei are all proton-rich; they decay primarily by emitting 𝛼 particles and by fission with a possible small electron capture (EC) branch. Here, due to the huge atomic numbers and associated relativistic effects, EC decays of superheavy systems are expected to differ from what is known in lighter nuclei. In this letter, using the quantified relativistic nuclear density functional theory and the quasiparticle random-phase approximation with the interaction optimized to experimental EC/𝛽 ± -decay half-lives, and Gamow-Teller resonance energies, we study the EC/𝛽 ± -decays in 𝑍=101–118 nuclei. Both allowed (1 + ) and first-forbidden (0 − ,1 − and 2 − ) transitions are considered. We show that the first-forbidden 1 − transitions dominate the decay rates in almost all studied nuclei. For proton-rich nuclei, EC dominates over 𝛽 + decay. Based on calculations with two relativistic energy density functionals, we identify 45 candidate nuclei in which a competition between weak decays and 𝛼 decay and spontaneous fission is expected.

A ≥ 220↗
DOE OSTI2024

Solar enhanced oxygen evolution reaction with transition metal telluride

The photo-enhanced electrocatalytic method of oxygen evolution reaction (OER) shows promise for enhancing the effectiveness of clear energy generation through water splitting by using renewable and sustainable source of energy. However, despite benefits of photoelectrocatalytic (PEC) water splitting, its uses are constrained by its low efficiency as a result of charge carrier recombination, a large overpotential, and sluggish reaction kinetics. Here, we illustrate that Nickel telluride (NiTe) synthesized by hydrothermal methods can function as an extremely effective photo-coupled electrochemical oxygen evolution reaction (POER) catalyst. In this study, NiTe was synthesized by hydrothermal method at 145°C within just an hour of reaction time. In dark conditions, the NiTe deposited on carbon cloth substrate shows a small oxygen evolution reaction overpotential (261 mV) at a current density of 10 mA cm –2 , a reduced Tafel slope (65.4 mV dec -1 ), and negligible activity decay after 12 h of chronoamperometry. By virtue of its enhanced photo response, excellent light harvesting ability, and increased interfacial kinetics of charge separation, the NiTe electrode under simulated solar illumination displays exceptional photoelectrochemical performance exhibiting overpotential of 165 mV at current density of 10 mA cm -2 , which is about 96 mV less than on dark conditions. In addition, Density Functional Theory investigations have been carried out on the NiTe surface, the results of which demonstrated a greater adsorption energy for intermediate -OH on the catalyst site. Since the -OH adsorption on the catalyst site correlates to catalyst activation, it indicates the facile electrocatalytic activity of NiTe owing to favorable catalyst activation. DFT calculations also revealed the facile charge density redistribution following intermediate -OH adsorption on the NiTe surface. This work demonstrates that arrays of NiTe elongated nanostructure are a promising option for both electrochemical and photoelectrocatalytic water oxidation and offers broad suggestions for developing effective PEC devices.

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

Spin dynamics of the centrosymmetric skyrmion material GdRu 2 Si 2

Magnetic skyrmion crystals are traditionally associated with non-centrosymmetric crystal structures; however, it has been demonstrated that skyrmion crystals can be stabilized by competing interactions in centrosymmetric crystals. To understand and optimize the physical responses associated with topologically nontrivial skyrmion textures, it is important to quantify their magnetic interactions by comparing theoretical predictions with spectroscopy data. Here, we present neutron diffraction and spectroscopy data on the centrosymmetric skyrmion material GdRu 2 Si 2 and show that the key spectroscopic features can be explained by magnetic interactions calculated using density functional theory. We further show that the recently proposed 2-q “topological spin stripe” structure yields better agreement with our data than a 1-q helical structure and identify how the magnetic structure evolves with temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Increased Defect Resistance and Ordering in MnBi 2 (Se 1– x Te x ) 4 via Accurate Diffusion Monte Carlo

Stabilizing materials and controlling defect formation remain key challenges in materials science, particularly for theory, where small energy differences must be resolved for accurate predictions. Here, we applied state-of-the-art theoretical methods to topological materials, focusing on MnBi 2 Te 4 (MBT), which is a promising intrinsic magnetic topological insulator. Antisite defects in MBT alter its electronic structure and magnetism, degrading topological properties and causing experimental inconsistencies. Using diffusion Monte Carlo and density functional theory, we investigated the thermodynamic stability and defect formation in MBT, MnBi 2 Se 4 (MBS), and MnBi 2 (Se 1–x Te x ) 4 . We found that MnBi 2 Se 2 Te 2 can be stable at finite temperatures, with higher defect formation energies due to stronger Mn–Se bonding and reduced internal strain. Se preferentially substitutes Te near Mn instead of Te in the outer layer, encouraging long-range ordering when incorporated. For MnBi 2 (Se 1–x Te x ) 4 , cluster expansion phase diagrams revealed solid solution behavior when x <0.5 and phase separation for larger x. MBT and MBS are topological insulators; therefore, the MnBi 2 (Se 1–x Te x ) 4 family could offer tunable topological behavior and improved stability.

MnBi2Te4↗
DOE OSTI2024

Ghost states and surface structures of the charge density wave kagome metal ScV 6 Sn 6

In this work, we investigate the high-temperature phase of the kagome metal ScV 6 Sn 6 using scanning tunneling microscopy/spectroscopy (STM/S) and density functional theory calculations. STM topographic images of the cleaved sample reveal two distinct surface terminations: flat islands with Sn termination and trenches terminated by kagome layers with Sn as the outermost atomic layer. STS measurements on the Sn-terminated and kagome-terminated surfaces show significant differences, in particular the presence of large density of states near the Fermi level in the former case. Our first-principles calculations reveal that the charge density on the kagome-terminated surface gives rise to “ghost states” which show intensity away from surface atoms, arising due to hybridization of orbitals above the surface. These states can obscure the intrinsic properties of the surface, potentially leading to misattribution of the surface termination. This underscores the need for careful interpretation in STM studies, especially when discerning surface states of localized states. Understanding the surface structure of this versatile quantum material provides essential information for interpreting surface-sensitive experiments, tailoring material properties, engineering interfaces, and controlling stability and reactivity. This knowledge paves the way for further exploration and potential applications of kagome lattice materials in various fields, including quantum computing, topological physics, and advanced electronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Evaluating the gradients of localized diabatic state energies and couplings at minimum cost

We calculate the gradients of Boys diabatic state energies and couplings when the electronic vector space is generated by configuration interaction singles. Our approach follows the Lagrangian approach of Paz and Glover (rather than direct differentiation of the adiabatic-to-diabatic approaches that have been published previously). The result is that we achieve a dramatic increase in savings over previous approaches, and the present approach should be immediately useful to scientists focused on electronic relaxation, especially chemists studying electron transfer who wish to go beyond the Condon approximation. Here, a future extension to time-dependent density functional theory in the Tamm–Dancoff approximation is clear.

Chemical physics↗
DOE OSTI2025

Pressure-induced structural and dielectric changes in liquid water at room temperature

Understanding the pressure-dependent dielectric properties of water is crucial for a wide range of scientific and practical applications. In this study, we employ a deep neural network trained on density functional theory data to investigate the dielectric properties of liquid water at room temperature across a pressure range of 0.1–1000 MPa. We observe a nonlinear increase in the static dielectric constant ɛ 0 with increasing pressure, a trend that is qualitatively consistent with experimental observations. This increase in ɛ 0 is primarily attributed to the increase in water density under compression, which enhances collective dipole fluctuations within the hydrogen-bonding network as well as the dielectric response. Despite the increase in ɛ 0 , our results reveal a decrease in the Kirkwood correlation factor G K with increasing pressure. Furthermore, this decrease in G K is attributed to pressure-induced structural distortions in the hydrogen-bonding network, which weaken dipolar correlations by disrupting the ideal tetrahedral arrangement of water molecules.

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

Assessing the Performance of Exchange‐Correlation Density Functionals in Describing the Iron‐Catalyzed Ammonia Synthesis System

Density functional theory (DFT) has been widely employed for elucidating mechanistic aspects of heterogeneous catalysis. However, the accuracy of DFT calculations relies heavily on selecting an exchange-correlation (XC) functional that correctly describes the electronic structure of materials involved in the reactions. This study assesses the accuracy of several XC density functionals for modeling the iron-catalyzed ammonia synthesis system. In the assessment of functional accuracy, experimental references are compared to DFT-calculated values for the formation energy of gas-phase ammonia and nitrogen, bulk Fe/Fe-nitride (γ′-Fe 4 N) lattice constants and cohesive/formation energies, and nitrogen and ammonia binding energies on Fe(100), Fe(111), Fe(110), and γ′-Fe 4 N(111). It is observed that the experimental value for each of these descriptors is accurately modeled by at least one functional. RPBE alone provides reliable estimates for both the lattice constant and cohesive energy of Fe and γ′-Fe 4 N. Temperature-programmed desorption experiments led to estimates for N and NH 3 adsorption across several Fe-based facets that are best captured by RPBE. These results highlight the importance of choosing an appropriate XC functional that accurately describes experimental systems and offer insights into effectively modeling the interactions between nitrogen and ammonia on Fe-based surfaces.

adsorption↗
DOE OSTI2025

Self-consistent microscopic calculations for electron captures on nuclei in core-collapse supernovae

Calculations for electron capture rates on nuclei with atomic numbers between 𝑍 = 20 and 𝑍 = 52 are performed in a self-consistent finite-temperature covariant energy density functional theory within the relativistic quasiparticle random-phase approximation. Electron captures on these nuclei contribute most to reducing the electron fraction during the collapse phase of core-collapse supernovae. The rates include contributions from allowed (Gamow-Teller) and first-forbidden (FF) transitions, and it is shown that the latter become dominant at high stellar densities and temperatures. Temperature-dependent effects such as Pauli unblocking and transitions from thermally excited states are also included. The new rates are implemented in a spherically symmetric one-dimensional simulation of the core-collapse phase. The results indicate that the increase in electron capture rates, due to inclusion of FF transitions, leads to reductions of the electron fraction at nuclear saturation density, the peak neutrino luminosity, and enclosed mass at core bounce. The new rates reaffirm that the most relevant nuclei for the deleptonization situate around the 𝑁 = 50 and 82 shell closures, but, compared to previous simulations, nuclei are less proton rich. Here, the new rates developed in this work are available, and will be of benefit to improve the accuracy of multidimensional supernova simulations.

Electron & muon capture↗
DOE OSTI2025

Thermodynamics of Liquid Uranium from Atomistic and Ab Initio Modeling

We present thermodynamic properties for liquid uranium obtained from classical molecular dynamics (MD) simulations and the first-principles theory. The coexisting phases method incorporated within MD modeling defines the melting temperature of uranium in good agreement with the experiment. The calculated melting enthalpy is in agreement with the experimental range. Classical MD simulations show that ionic contribution to the total specific heat of uranium does not depend on temperature. The density of states at the Fermi level, which is a crucial parameter in the determination of the electronic contribution to the total specific heat of liquid uranium, is calculated by ab initio all electron density functional theory (DFT) formalism applied to the atomic configurations generated by classical MD. The calculated specific heat of liquid uranium is compared with the previously calculated specific heat of solid γ-uranium at high temperatures. The liquid uranium cannot be supercooled below T sc ≈ 800 K or approximately about 645 K below the calculated melting point, although, the self-diffusion coefficient approaches zero at T D ≈ 700 K. Uranium metal can be supercooled about 1.5 times more than it can be overheated. The features of the temperature hysteresis are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Enhanced stability and ordering of 2D MnBi2Se2Te2

This dataset includes density functional theory and quantum Monte Carlo calculations on MnBi(Se(1-x)Te(x))4 system for its phase stability and defect formation energies. The data supports publication "Increased Defect Resistance and Ordering in MnBi2(Se(1–x)Te(x))4 via Accurate Diffusion Monte Carlo" (https://doi.org/10.1021/acs.jpcc.5c04299).

36 MATERIALS SCIENCE↗
DOE OSTI2024

Calibration and validation of the foundation for a multiphase strength model for tin

In this work, the Common Model of Multi-phase Strength and Equation of State (CMMP) model was applied to tin. Specifically, calibrations of the strength-specific elements of the CMMP foundation were developed with a combination of experiments and theory, and then the model was validated experimentally. The first element of the foundation is a multi-phase analytic treatment of the melt temperature and the shear modulus for the solid phases. These models were parameterized for each phase based on ab initio calculations using the software VASP (Vienna Ab initio Simulations Package) based on density functional theory. The shear modulus model for the ambient phase was validated with ultrasonic sound speed measurements as a function of pressure and temperature. The second element of the foundation is a viscoplastic strength model for the phase, upon which strength for inaccessible higher-pressure phases can be scaled as necessary. The stress–strain response of tin was measured at strain rates of 10 -3 to 3 x 10 3 s -1 and temperatures ranging from 87 to 373 K. The Preston–Tonks–Wallace (PTW) strength model was fit to that data using Bayesian model calibration. For validation, six forward and two reverse Taylor impact experiments were performed at different velocities to measure large plastic deformation of tin at strain rates up to ⁠10 5 s -1 . The PTW model accurately predicted the deformed shapes of the cylinders, with modest discrepancies attributed to the inability of PTW to capture the effects of twinning and dynamic recrystallization. Some material in the simulations of higher velocity Taylor cylinders reached the melting temperature, thus testing the multiphase model because of the presence of a second phase, the liquid. In simulations using a traditional modeling approach, the abrupt reduction of strength upon melt resulted in poor predictions of the deformed shape and non-physical temperatures. With CMMP, the most deformed material points evolved gradually to a mixed solid–liquid but never a fully liquid state, never fully lost strength, stayed at the melt temperature as the latent heat of fusion was absorbed, and predicted the deformed shape well.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ordering in GaAs co-doped with Bi and N

Introducing only a few atomic percent of Bi or N in GaAs has a large effect on the band gap of the material. Specifically Bi doped GaAs shows potential for local band gap engineering in optoelectronic applications. The incorporation of Bi and N into GaAs is difficult due to strain effects. In this work we study the ordering of these dopants at the atomic scale in order to get a better understanding of the behavior of these dopants in the host lattice. Cross-sectional scanning tunneling microscopy is used to find the exact position of Bi and N dopants in the GaAs matrix, allowing us to study both their nearest neighbor pair occurrences and pair correlation functions. An attractive interaction between Bi dopants at short ranges (1-2 nm) is found and a similar effect is observed between N dopants. Here, we find a repulsive interaction with a similar length scale between Bi and N dopants. A similar repulsion is found in the Bi-N nearest neighbor pairs. Density functional theory is used to calculate the different nearest neighbor pair energies and test these results to the experimental pair occurrences. It is concluded from the experimental and theoretical results that the growth conditions and N inclusion greatly affects the Bi distribution in GaAs.

36 MATERIALS SCIENCE↗
DOE OSTI2024

An Impedance Humidity Sensor Based on CVD Grown WSe 2 2D Films

2D materials-based planar devices have been sought after for gas and humidity sensor applications owing to their high sensitivity due to their ultrathin nature and increased surface area. Herein, 2D WSe 2 films-based planner devices were fabricated to evaluate their performance for potential application in humidity sensors. The structure and morphology of the WSe 2 sensing films were studied by Raman spectroscopy, photoluminescence, optical microscopy, and scanning electron microscopy (SEM) techniques. The relative humidity-dependent sensing performances of the 2D WSe 2 films were evaluated using current–voltage (I–V), relative impedance, and complex total impedance (electrochemical impedance spectroscopy, EIS) techniques. This sensor showed a very stable and repeatable performance investigated over the period of 18 months. The response and recovery times of the WSe 2 -based impedance sensors were found to be 40 and 30 s, respectively. The WSe 2 sensor exhibits a 2 ± 1% minimum hysteresis at a lower humidity level (50% RH) and around 8 ± 2% maximum hysteresis at relative humidity level (70% RH). This value of hysteresis is in the acceptable range during the initial investigation of any material to explore its potential for application in humidity sensors. The density functional theory (DFT) analysis was performed with the interaction from water molecules of the WSe 2 sample for a better understanding of the sensing mechanism of devices.

2D- WSe2↗
DOE OSTI2025

Supporting Data for "Constraints on magnetism and correlations in RuO2 from lattice dynamics and Mössbauer spectroscopy"

Contents of this DOI are data for the research paper "Constraints on magnetism and correlations in RuO2 from lattice dynamics and Mossbauer spectroscopy" by the authors: George Yumnam, et. al. The dataset contains data from inelastic x-ray scattering measurements of RuO2 single crystals performed at the Advanced Photon Source in Argonne National Lab. This dataset also contains data from inelastic neutron scattering experiments of RuO2 powder performed at the ARCS spectrometer at Spallation Neutron Source at ORNL. Supporting theoretical calculations based on density functional theory with r2SCAN and DFT+U methods is also included in this dataset. Please see the README.txt files for more information of the dataset and file contents. For comments and questions, please contact: George Yumnam (yumnamg@ornl.gov) and/or Raphael P Hermann (hermannrp@ornl.gov)

density functional theory↗