Search NASA⌕ Search

SEARCH · Search NASA

Results for “density functional calculation”

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

Structural transition and possible pressure-induced superconductivity in a suboxide La 5 Pb 3 O

Here, in this study, we report a structural phase transition and its interplay with superconductivity in the suboxide La 5⁢ Pb 3 ⁢O. Upon cooling through 𝑇 𝑡 = 225K, La 5⁢ Pb 3⁢ O transforms from a high-temperature 𝐼⁢4/mcm to a low-temperature 𝑃⁢4/ncc structure in which La-Pb dimerization along the 𝑐 axis occurs. This transition is accompanied by anomalies in the temperature dependence of electrical resistivity and specific heat. High-pressure electrical transport measurements reveal that hydrostatic pressure suppresses the structural transition and possibly induces superconductivity with a maximum superconducting temperature of 10 K. Density functional theory calculations show minimal changes in the electronic density of states and no gap opening at 𝐸 𝐹 across 𝑇 𝑡 , suggesting that the transition is driven by bonding effects rather than Fermi surface instability. These findings establish La 5⁢ Pb 3 ⁢O as a promising platform for exploring the interplay between weak structural transitions and superconductivity.

Yan, Jiaqiang [Oak Ridge National Laboratory (ORNL↗
DOE OSTI2025

Direct Deoxygenation of Phenol over Fe-Based Bimetallic Surfaces Using On-the-Fly Surrogate Models

We present an accelerated nudged elastic band (NEB) study of phenol direct deoxygenation (DDO) on Fe-based bimetallic surfaces using a recently developed Gaussian process regression (GPR) calculator. Our test calculations demonstrate that the GPR calculator achieves up to 3 times speedup compared to conventional density functional theory calculations while maintaining high accuracy, with energy barrier errors below 0.015 eV. Using GPR-NEB, we systematically examine the DDO mechanism on pure Fe(110) and surfaces modified with Co and Ni in both top and subsurface layers. Our results show that subsurface Co and Ni substitutions preserve favorable thermodynamics and kinetics for both C–O bond cleavage and C–H bond formation, comparable to those on the pure Fe(110) surface. In contrast, top-layer substitutions generally increase the C–O bond cleavage barrier, render the step endothermic, and result in significantly higher reverse reaction rates, making DDO unfavorable on these surfaces. This work demonstrates the effectiveness of GRR-accelerated transition state searches for complex surface reactions and provides insights into rational design of bimetallic catalysts for selective deoxygenation.

Aromatic compounds↗
DOE OSTI2025

Regioisomeric Engineering for Multicharge and Spin Stabilization in Two-Electron Organic Catholytes

Developing multicharge and spin stabilization strategies is fundamental to enhancing the lifetime of functional organic materials, particularly for long-term energy storage in multiredox organic redox flow batteries. Current approaches are limited to the incorporation of electronic substituents to increase or decrease the overall electron density or bulky substituents to sterically shield reactive sites. With the aim to further expand the molecular toolbox for charge and spin stabilization, we introduce regioisomerism as a scaffold-diversifying design element that considers the collective and cumulative electronic and steric contributions from all of the substituents based on their relative regioisomeric arrangements. Through a systematic study of regioisomers of near-planar aromatic cyclic triindoles and nonplanar nonaromatic cyclic tetraindoles, we demonstrate that this regioisomeric engineering strategy significantly enhances the H-cell cycling stability in the above two new classes of 2e – catholytes, even when current strategies failed to stabilize the multicharged species. Density functional theory calculations reveal that the strategy operates by redistributing the charge and spin densities while highlighting the role of aromaticity in charge stabilization. The most stable 2e – catholyte candidate was paired with a viologen derivative anolyte to achieve a proof-of-concept all-organic flow battery with 1.26–1.49 V, 98% capacity retention, and only 0.0117% fade/h and 0.00563% fade/cycle over 400 cycles (192 h), which is the highest capacity retention ever reported over 400 cycles in a multielectron all-organic flow battery setup. We anticipate regioisomeric engineering to be a promising strategy complementary to conventional electronic and steric approaches for multicharge and spin stabilization in other functional organic materials.

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

Ferromagnetic Interlayer Exchange Coupling in a Few Layers of CrSBr on a Gold Thin Film

Abstract The two-dimensional character of van der Waals magnets allows for efficient control of their properties via proximity effects and electrical stimuli, making them promising candidates for application in spin electronics. We use spin-polarized low-energy electron microscopy to directly image the magnetic texture of thin CrSBr on top of a Au film, discovering a ferromagnetic ground state. We argue that the stabilization of the ferromagnetic ordering─as compared to the conventional antiferromagnetic one─is obtained via electron transfer from the Au film to the CrSBr flakes, in agreement with ab initio density functional theory calculations. Reflected-electron spectroscopy shows clear differences in the unoccupied density of states between a few layers of CrSBr on Au and bulk CrSBr, pointing toward electronic band structure modification in thin CrSBr. This work sheds light on the possibility of tuning the magnetic properties of two-dimensional magnets via substrate engineering.

Bosma, Rixt [University of Groningen , , ,]↗
DOE OSTI2025

An Accurate, Affordable Density Functional Tight-Binding Model for Excited State Hydrocarbon Polymer Molecular Dynamics

We have developed a density functional tight-binding model for hydrocarbon excited state dynamics by referencing to high-level electronic structure theory and incorporating a many-body repulsive energy. We then validate our model against n-octane geometry optimizations, bond dissociation scans, and vibrational frequencies. Our model is approximately 10 3 times more efficient than hybrid time-dependent density functional theory calculations with comparable accuracy. In conclusion, our efforts enable longer timescale excited state simulations of photochemistry and scattering of incident radiation.

Chemical calculations↗
DOE OSTI2025

Planar Defect Layers Template a High-Pressure InBi Polymorph

The short- and long-range order of III–V materials under high pressure has long been the subject of debate, with advancements in structural characterization leading to significant revisions to the accepted structural models. Despite these revisions, previous high-pressure structural assignments in the In–Bi system include the site-disordered β-Sn structure type, a structure type demonstrated to be nonexistent in analogous III–V systems. While X-ray diffraction is consistent with site disordering in InBi at high pressure, cluster expansion calculations indicate that disordering requires temperatures above 3000 K. Here, we propose InBi as a model material for studying unique high-pressure planar defects due to its highly anisotropic stress-dependent properties and structure. Specifically, we identify two sets of planar defects that mimic the diffraction pattern of a site disordered β-Sn structure type and are compatible with the calculated disorder barrier. We derive these defects by symmetry relations over crystallographic transitions. Density functional theory calculations of the proposed defects suggest that these defects are stabilized by diminishing interlayer separations with pressure. Further, we find that one of the proposed defects closely resembles a bulk high-pressure phase of InBi, InBi-ϵ, and we assert that the proposed defects order upon heating, acting as a template for InBi-ϵ growth. The proposed defects and their electronic structure provide a basis for the trend of superconducting critical temperature with increasing pressure. These methods for identifying defects are generalizable to other materials with reports of site disorder at high pressure, prompting a broader search for related high-pressure defects.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Point defect energetics in gallium arsenide, a comprehensive density functional theory study

In materials, point defects often control or modify functional properties. To predict the performance of materials intended for application in optoelectronic devices, it is imperative to understand the properties of those point defects. For the first time, all six intrinsic defects of GaAs, a key optoelectronics material, and their charge transition levels are calculated using density functional theory with the HSE06 functional. For comparison, both PBE and r 2 SCAN calculations are also carried out. The HSE06 results are found to be in better agreement with experimental data than previous calculations. In conclusion, the importance of using the exact electron exchange present in hybrid functionals and larger supercells to accurately determine defect levels and ground state defect configurations is demonstrated.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Surface dependence of electronic growth of Cu(111) on MoS 2

Scanning tunneling microscopy shows that copper deposited at room temperature onto a freshly exfoliated MoS 2 surface forms Cu(111) clusters with periodic preferred heights of 5, 8, and 11 atomic layers. These height intervals correlate with Fermi nesting regions along the necks of the bulk Cu Fermi surface, indicating a connection between physical and electronic structures. Density functional theory calculations of freestanding Cu(111) films support this as well, predicting a lower density of states at the Fermi level for these preferred heights. This is consistent with other noble metals deposited on MoS 2 that exhibit electronic growth, in which the metal films self-assemble as nanostructures minimizing quantum electronic energies. Here, we have discovered that it is critical for the metal deposition to begin on a clean MoS 2 surface. If copper is deposited onto an already Cu coated surface, even if the original film displays electronic growth, the resulting Cu film lacks quantization. Instead, the preferred heights of the Cu clusters simply increase linearly with the amount of Cu deposited upon the surface. We believe this is due to different bonding conditions during the initial stages of growth. Newly deposited copper would bond strongly to the already present copper clusters, rather than the weak bonding, which exists to the van der Waals terminated surface of MoS 2 . The stronger bonding with previously deposited clusters hinders additional Cu atoms from reaching their lowest quantum energy state. The interface characteristics of the van der Waals surface enable surface engineering of self-assembled structures to achieve different applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Metal Identity and Coordination Environment Modulate Single-Atom Catalyst Stability During Electrocatalysis

A major hurdle to the implementation of single-atom catalysts (SACs) in real-world systems is a poor understanding of their stability under operating conditions, which is particularly relevant due to the high surface free energy of SACs. Here, we evaluated the aggregation behavior of a suite of SACs varied by metal identity (Fe, Co, Ni, and Cu) during electrocatalytic nitrate reduction using in situ X-ray absorption spectroscopy. The metal center had significant influence on reconstruction, where under identical applied reductive potentials, SACs underwent varying levels of reconstruction, ranging from no discernible change to complete reduction into metallic nanoparticles. Such in situ experiments revealed Cu SACs to be the most susceptible to aggregation, prompting a deeper investigation into how coordination environment (O-, B-, and N-graphene) affected Cu SAC aggregation. We further conducted density functional theory calculations to elucidate the relationship between Cu SAC structure and stability. This work deconvolutes the relationship between SAC architecture and stability, which is essential to evaluate and explain for the realization of SACs for electrocatalysis.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Targeted Chemical Looping Materials Discovery by an Inverse Design

Chemical looping with oxygen uncoupling (CLOU) materials is actively sought for combustion of carbonaceous materials to achieve complete conversion and capture of carbon dioxide. These materials may play a vital role in reducing atmospheric carbon via negative carbon output. However, there is no one‐size‐fits‐all approach as different operating conditions and feedstocks may require different CLOU materials. As a result, the exploration and discovery of high‐performance CLOU materials can be a slow process. To address this challenge, a high‐throughput inverse machine learning workflow that identifies optimum materials from perovskite oxides for a given set of targets is developed—temperature and Gibbs free energy of oxygen formation. The model is trained on high‐throughput density functional theory calculations of CLOU materials and inverts the materials design process using a genetic algorithm to produce realistic substituted SrFeO 3‐δ compositions as output. Using the inverse model, it is able to identify several interesting new families of CLOU materials: Sr 1‐ x A x Fe 1‐ y B y O 3‐δ (e.g., A = Ca or K; B = Mg, Bi, Mn, Ni, Co, Cu, or Zn). These materials have shown promising properties, and some of them even outperform the benchmark material in terms of oxygen release kinetics under relevant CLOU operating conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Heavy boron doping effects on biaxially tensile strained germanium (>1.5%) investigated via structural characterization, effective lifetime assessment and atomistic modeling

Highly tensile strained germanium (ε-Ge) represents an essential material system for emerging electronic and photonics applications. Moreover, adjusting the doping levels to moderate or high concentrations can effectively tailor the properties of ε-Ge for specific applications. This article combines experimental characterization with a theoretical framework to examine the effects of heavy elemental boron (B) doping on pseudomorphic sub-50 nm ε-Ge. High resolution X-ray diffractometry is used to validate tensile strain levels of 1.53% and 1.68% in Ge epilayers, surpassing the indirect-to-direct band gap crossover point at ∼1.5% biaxial tensile strain. Cross-sectional transmission electron microscopy revealed visual evidence of stacking faults and surface roughening in 1.68% ε-Ge, although a coherent and abrupt Ge/III–V heterointerface is observed, devoid of interfacial misfit dislocations. Effective lifetime measurements demonstrated approximately twofold enhancement in 1.53% B-doped ε-Ge (N B ∼7 × 10 19 cm −3 ) compared to its unstrained B-doped counterpart, while no such improvement was observed in 1.68% B-doped ε-Ge. This lack of enhancement is attributed to the presence of stacking faults and surface roughness within the ε-Ge epilayer. Through density functional theory calculations, we independently demonstrate that substitutional B atoms induce local deformation of Ge–Ge bonds in both unstrained Ge and ε-Ge epilayers, resulting in an additive tensile strain. This phenomenon could potentially lead to dynamic reduction and overcoming of the critical layer thickness for the system, facilitating the nucleation and subsequent glide of 90° leading Shockley partial dislocations, thereby generating stacking faults. In essence, these findings establish an upper limit on the B-doping concentration that can be achieved in highly ε-Ge epilayers, and collectively, offer valuable insights into the significance of heavy doping in Ge-based heterostructures. As such, this study delineates a fundamental constraint for integrating heavily doped ε-Ge in high-performance optoelectronic systems, necessitating precise strain-doping co-optimization to avoid performance degradation.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Ultrafast Optically Induced Perturbation of Oxygen Octahedral Rotations in Multiferroic BiFeO 3 Thin Films

The functional properties of complex oxides, including magnetism and ferroelectricity, are closely linked to subtle structural distortions. Ultrafast optical excitations provide the means to manipulate structural features and ultimately to affect the functional properties of complex oxides with picosecond-scale precision. Here, we report that the lattice expansion of multiferroic BiFeO 3 following above-bandgap optical excitation leads to distortion of the oxygen octahedral rotation (OOR) pattern. The continuous coupling between OOR and strain was probed using time-resolved X-ray free-electron laser diffraction with femtosecond time resolution. Density functional theory calculations predict a relationship between the OOR and the elastic strain consistent with the experiment, demonstrating a route to employing this approach in a wider range of systems. Ultrafast control of the functional properties of BiFeO 3 thin films is enabled by this approach because the OOR phenomena are related to ferroelectricity, and via the Fe-O-Fe bond angles, the superexchange interaction between Fe atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Benchmarking Density Functional Theory Methods for Efficient Calculations of a Strongly Correlated Li 1– x Ni 1– y O 2−δ System

Transition metal oxides (TMOs), such as LiNiO 2 , are promising candidates for energy storage and electronic devices due to their unique electronic properties, exceptional physical and chemical characteristics, and ability to adopt multiple oxidation states. However, accurately predicting their properties using mean-field density functional theory (DFT) is challenging due to the presence of strongly correlated d-electrons and the complex interplay between their structural, electronic, and magnetic responses. These challenges are further exacerbated by the need to model defects, surfaces, and interfaces, which require computationally efficient, large-scale simulations. To address these issues, we carry out a benchmark study on the Li 1–x NiO 2 system, evaluating the performance of several popular functionals. Our findings demonstrate that combining SCAN functional relaxation with single-step HSE calculations provides a practical and scalable computational strategy. This approach balances accuracy and efficiency, enabling high-throughput simulations of strongly correlated TMOs and improved predictive modeling capability of TMOs for practical applications.

25 ENERGY STORAGE↗
DOE OSTI2025

Defluorination Mechanisms and Real-Time Dynamics of Per- and Polyfluoroalkyl Substances on Electrified Surfaces

Per- and polyfluoroalkyl substances (PFAS) are persistent environmental contaminants found in groundwater sources and a wide variety of consumer products. In recent years, electrochemical approaches for the degradation of these harmful contaminants have garnered a significant amount of attention due to their efficiency and chemical-free modular nature. However, these electrochemical processes occur in open, highly non-equilibrium systems, and a detailed understanding of PFAS degradation mechanisms in these promising technologies is still in its infancy. To shed mechanistic insight into these complex processes, we present the first constant-electrode potential (CEP) quantum calculations of PFAS degradation on electrified surfaces. These advanced CEP calculations provide new mechanistic details about the intricate electronic processes that occur during PFAS degradation in the presence of an electrochemical bias, which cannot be gleaned from conventional density functional theory calculations. We complement our CEP calculations with large-scale ab initio molecular dynamics simulations in the presence of an electrochemical bias to provide time scales for PFAS degradation on electrified surfaces. Taken together, our CEP-based quantum calculations provide critical reaction mechanisms for PFAS degradation in open electrochemical systems, which can be used to prescreen candidate material surfaces and optimal electrochemical conditions for remediating PFAS and other environmental contaminants.

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

Defect levels and self-compensation in iodine-doped CdTe single crystals

This study systematically documents defect levels in n-type iodine-doped cadmium telluride (CdTe:I) crystals as measured by thermoelectric effect spectroscopy, Hall-effect measurements, and photoluminescence and calculated by density functional theory. The primary donor, I Te , was identified with an activation energy of ∼0.05 eV (measured)/∼0.13 eV (calculated). Deep acceptor states, V Cd and I Te –V Cd complex, exhibited activation energies of ∼0.1 eV (measured)/∼0.19 eV (calculated) and ∼0.12 eV (measured)/∼0.24 eV (calculated), respectively. Self-compensation and Fermi level pinning were observed at ∼0.7 eV (measured)/∼1.08 eV (calculated), contributing to high resistivity in as-grown samples. Post-growth Cd annealing effectively removes compensation centers, releases iodine donors, and significantly reduces resistivity, resulting in an ideal n-type solar cell material.

Cd annealing↗
DOE OSTI2024

Antimony stable isotope fractionation during adsorption onto birnessite: A molecular perspective from X-ray absorption spectroscopy and density functional theory

Sorption of antimony (Sb) onto birnessite significantly influences the fate of Sb in oceanic and terrestrial environments and fractionates Sb isotopes. Nevertheless, little is known about Sb isotopic fractionation during its adsorption on birnessite. Here, in this study, we show the value of Δ 123 Sb adsorbed-aqueous increases from −0.398 to −0.332 ‰ in 1 h and then decreases and stabilizes at −0.384 ‰ in 72 h. The enrichment of the light Sb isotope is predominantly due to the distortion of the octahedral symmetry. X-ray absorption spectroscopy results indicate Sb first forms a double-corner-sharing complex on birnessite and then transforms to a double-edge-sharing complex during adsorption. The optimized bond distances for double-corner-sharing (3.37 Å) and double-edge-sharing (2.90 Å) complexes calculated using density functional theory (DFT) fits well with the structure (3.41 and 3.00 Å) revealed by X-ray absorption spectroscopy, respectively. The fractionation derived from reduced partition function ratios calculated using DFT aligns well with the experimental results. Therefore, the variation in Sb isotopic fractionation during adsorption is attributed to the evolving structure of Sb complexes on birnessite. Our results demonstrate the isotopic fractionation of Sb during adsorption on birnessite and provide a molecular-scale understanding of Sb behavior, contributing to the correct reconstruction of the Sb isotope composition of ancient seawater using ferromanganese crusts and nodules, and efforts to trace Sb migration in epigenetic mining environments.

Adsorption↗
DOE OSTI2025

Selective Hydrogenation of Furfural Acetone over a Cu Catalyst: Combined Theoretical and Experimental Study

The selective hydrogenation of biomass-derived hydroxymethyl furfural (HMF)-acetone-HMF (HAH) presents an alternative route to producing highly functionalized polyesters and polyurethanes. While HAH undergoes furan ring hydrogenation over Pd, Ru, and Ni catalysts, furan ring hydrogenation is not observed over Cu. Herein, we combine reaction kinetics experiments and density functional theory calculations to elucidate the selective hydrogenation behavior of HAH over Cu catalysts. We identified furfural acetone (FAc) as a suitable surrogate for modeling HAH hydrogenation over Cu surfaces and performed reaction kinetics experiments between temperatures of 313–393 K and a H 2 partial pressure of 55 bar. Similar to the behavior of HAH, hydrogenation of FAc follows a consecutive two-step hydrogenation pathway over Cu and does not undergo furan ring hydrogenation. The apparent activation energy barriers for hydrogenation of the aliphatic double bond (0.58 eV) and carbonyl (0.43 eV) of FAc measured in a continuous flow reactor setup are consistent with those reported in prior studies for batch HAH hydrogenation. Reaction orders with respect to each reactant, including H 2 and FAc, were determined to be nearly one. Density functional theory (DFT; GGA-PBE-D3) calculations on Cu(111) showed that the hydrogenation of the aliphatic double bond of FAc is more facile than the hydrogenation of the furan ring, which displays weak interactions with the Cu surface. We determined an apparent activation energy barrier for FAc hydrogenation (0.58 eV) that agreed with our DFT predictions (highest barrier for FAc hydrogenation of 0.57 eV). Finally, our DFT calculations further show that weak interactions between the furan ring and Cu surface are responsible for the selective hydrogenation behavior.

Biomass↗