Search NASA⌕ Search

SEARCH · Search NASA

Results for “density functional theory 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Optimized Auxiliary Functions for Robust Mitigation of Finite-Size Errors in Periodic Hybrid Density Functional Theory

When calculating properties of periodic systems at the thermodynamic limit (TDL), the dominant source of finite size error (FSE) arises from the long-range Coulomb interaction, and can manifest as a slowly converging quadrature error when approximating an integral in the reciprocal space by a finite sum. The singularity subtraction (SS) method offers a systematic approach for reducing this quadrature error and thus the FSE. Here, in this work, we first investigate the performance of the SS method in the simplest setting, aiming at reducing the FSE in exact exchange calculations by subtracting the Coulomb contribution with a single, adjustable Gaussian auxiliary function. We demonstrate that a simple fitting method can robustly estimate the optimal Gaussian width and leads to rapid convergence toward the TDL. Furthermore, we suggest new forms of the auxiliary function, whose optimal parameters could also be determined through least-squares fitting. For a range of semiconductors and insulators, the proposed auxiliary functions achieve robust, millihartree-level accuracy in hybrid density functional theory calculations, including cases with sparse k-meshes and large basis sets.

Quiton, Stephen Jon [University of California, Ber↗
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 OSTI2021

Computation of the Thermal Expansion Coefficient of Graphene with Gaussian Approximation Potentials

Direct experimental measurement of thermal expansion coefficient without substrate effects is a challenging task for two-dimensional (2D) materials, and its accurate estimation with large-scale ab initio molecular dynamics is computationally very expensive. Machine learning-based interatomic potentials trained with ab initio data have been successfully used in molecular dynamics simulations to decrease the computational cost without compromising the accuracy. In this work, we investigated using Gaussian approximation potentials to reproduce the density functional theory-level accuracy for graphene within both lattice dynamical and molecular dynamical methods, and to extend their applicability to larger length and time scales. Two such potentials are considered, GAP17 and GAP20. GAP17, which was trained with pristine graphene structures, is found to give closer results to density functional theory calculations at different scales. Further vibrational and structural analyses verify that the same conclusions can be deduced with density functional theory level in terms of the reasoning of the thermal expansion behavior, and the negative thermal expansion behavior is associated with long-range out-of-plane phonon vibrations. Thus, it is argued that the enabled larger system sizes by machine learning potentials may even enhance the accuracy compared to small-size-limited ab initio molecular dynamics.

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

Selective-Area Remote Epitaxy of ZnO Microrods Using Multilayer–Monolayer-Patterned Graphene for Transferable and Flexible Device Fabrications

Here, selective-area remote epitaxy (SA-REpi) is demonstrated for fabricating mechanically releasable position-controlled ZnO microrod (MR) arrays from donor wafers in an arrayed form. Intaglio-patterned graphene, consisting of basal single-layer graphene (SLG) overlayered with multilayer graphene (MLG) patterned with perforated holes, is transferred onto a GaN/Al 2 O 3 wafer on which the hydrothermal synthesis is performed for growing ZnO MRs. The basal SLG area exposed through the MLG pattern yields ZnO MRs, whereas the MLG plateau inhibits the growth. The noncovalent remote epitaxial heterointerface enables the release of the MR overlayer in the arrayed form, and the original source wafer is refurbished for reproducibly repeating the SA-REpi. Density-functional theory calculations suggest that localized surface charge density is induced on the surface of SLG by the underlying GaN across ultrathin SLG, which possibly provides a driving force for precursor adatoms and the following remote epitaxy of ZnO. In contrast, the induction of the charge density redistribution does not clearly occur through MLG; so, that keeps the surface of MLG nearly charge-neutral. The diameter and spacing of ZnO MRs are controlled in a designed way by changing the pattern geometries. High-resolution scanning transmission electron microscopy reveals the remote heteroepitaxial relationship at an atomic level. The remote epitaxy is expected to provide an ideal platform to transfer the addressable spatial arrays of nano- or micro-architecture semiconductor components to arbitrary target surfaces directly after the growth without the assembly procedures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

One-bond 13 C– 13 C spin-coupling constants in saccharides: a comparison of experimental and calculated values by density functional theory using solid-state 13 C NMR and X-ray crystallography

Methyl aldohexopyranosides were 13 C-labeled at contiguous carbons, crystallized, and studied by single-crystal X-ray crystallography and solid-state 13 C nuclear magnetic resonance (NMR) spectroscopy to examine the degree to which density functional theory (DFT) can calculate one-bond 13C–13C spin-coupling constants ( 1 J CC ) in saccharides with sufficient accuracy to permit their use in MA'AT analysis, a newly-reported hybrid DFT/NMR method that provides probability distributions of molecular torsion angles in solution (Zhang et al., J. Phys. Chem. B, 2017, 121, 3042–3058; Meredith et al., J. Chem. Inf. Model., 2022, 62, 3135–3141). Experimental 1 J CC values in crystalline samples of the doubly 13 C-labeled compounds were measured by solid-state 13 C NMR and compared to those calculated from five different DFT models: (1) 1 J CC values calculated from single structures identical to those observed in crystalline samples by X-ray crystallography (all atom refinement); (2) 1 J CC values calculated from the single structures in (1) but after Hirshfeld atom refinement (HAR); (3) 1 J CC values calculated from the single structures in (1) after DFT-optimization of hydrogen atoms only; and (4 and 5) 1 J CC values calculated in rotamers of torsion angle θ 2 (C1–C2–O2–O2H) or ω(C4–C5–C6–O6) from which either specific or generalized parameterized equations were obtained and used to calculate 1 J CC values in the specific θ 2 or ω rotamers observed in crystalline samples. Good qualitative agreement was observed between calculated 1 J CC values and those measured by solid-state 13 C NMR regardless of the DFT model, but in no cases were calculated 1 J CC values quantitative, differing (over-estimated) on average by 4–5% from experimental values. These findings, and those reported recently from solution NMR studies (Tetrault et al., J. Phys. Chem. B 2022, 126, 9506–9515), indicate that improvements in DFT calculations are needed before calculated 1 J CC values can be used directly as reliable constraints in MA'AT analyses of saccharides in solution.

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

Exceptional magnetic and magnetoelastic behavior of rare-earth non-centrosymmetric Sm 7 Pd 3

Magnetic compounds possessing an intrinsic combination of near-zero magnetization with high magnetic anisotropy are highly desirable for spinronic applications and memory recording. A comprehensive study of Sm 7 Pd 3 binary compound uncovered a unique combination of strong magnetoelastic behavior, very low net magnetization, and exceptionally high magnetic coercivity. The temperature-dependent X-ray synchrotron powder diffraction study indicates the abrupt changes in the compound's lattice parameters at the magnetic ordering temperature of T C =169 K, although the crystal structure remains non-centrosymmetric hexagonal Th 7 Fe 3 -type down to 6 K. Density functional theory calculations confirm high intrinsic magnetocrystalline anisotropy of Sm 7 Pd 3 , which explains the extremely large coercivity of the polycrystalline sample, up to H cr = 130 kOe at 2 K. In conclusion, this discovery brings to life a novel class of highly anisotropic materials that are distinctly different from known spintronic materials, making them interesting future systems for magnetic memory research.

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

Phase Boundaries, Isotope Effect and Superconductivity of Lithium Under Hydrostatic Conditions

We present theoretical and experimental studies of superconductivity and low temperature structural phase boundaries in lithium. We mapped the structural phase diagram of 6Li and 7Li under hydrostatic conditions between 5 – 55 GPa and within the temperature range of 15 – 75 K, observing the FCC → hR1 → cI16 phase transitions. 6Li and 7Li show some differences at the structural boundaries, with a potential shift of the phase boundaries of 6Li to lower pressures. Density functional theory calculations and topological analysis of the electron density elucidates the superconducting properties and interatomic interactions within these phases of lithium.

Deemyad, Shanti↗
DOE OSTI2022

A new perspective on the initial hydrogenation of $\mathrm{TiFe_{0.9}M_{0.1}}$ ($\mathrm{M = V, Cr, Fe, Co, Ni}$) gained from surface oxide analysis and nucleation thermodynamics

Despite the promise of TiFe-based alloys as low-cost solid-state hydrogen storage materials with mild operating conditions and reasonable hydrogen capacity, their initial hydrogenation process is difficult, hindering broad utilization. The effect of alloying element on the initial hydrogenation kinetics of TiFe alloys, TiFe 0.9 M 0.1 (M = V, Cr, Fe, Co and Ni), was evaluated by analyzing changes to the passivating surface oxide layer that inhibits hydrogen permeation, as well as the ease of initial-stage hydrogen absorption into the underlying matrix. X-ray photoelectron spectroscopy and atom probe tomography revealed key variations in surface oxide compositions and thinning of the passivating oxide layer compared to pure TiFe, which suggests suppressed oxide growth by alloying-induced elemental redistribution. At the same time, density functional theory calculations predicted exothermic formation of hydride nuclei when alloying with V or Cr, as well as a reduced nucleation barrier when alloying with Co or Ni. Overall, these results are consistent with the observed experimental trend of the activation kinetics. We propose that improvements in activation kinetics of TiFe with alloying arises from the combined effect of reduced passivating oxide thickness and easier hydride nucleation, offering a starting point for alloy design strategies towards further improvement.

25 ENERGY STORAGE↗
DOE OSTI2021

Mechanistic Investigation of Dimethylmercury Formation Mediated by a Sulfide Mineral Surface

Mercury (Hg) pollution is a global environmental problem, but many key transformations among 29 various Hg species are not well understood. Although anaerobic microorganisms are the primary 30 producers of monomethylmercury (MMeHg), the abiotic formation of dimethylmercury (DMeHg) 31 from MMeHg have been suggested to account for a large portion of DMeHg in the ocean. It has 32 been shown that abiotic formation of DMeHg from MMeHg can be facilitated by reduced sulfur 33 groups on mineral surfaces. However, a detailed mechanistic understanding of this reaction is 34 lacking. Here we perform density functional theory (DFT) calculations to explore mechanisms of 35 DMeHg formation on the 110 surface of a CdS(s) (hawleyite) nanoparticle. We show that 36 coordination of MMeHg substituents to adjacent reduced sulfur groups protruding from the surface 37 facilitate DMeHg formation and that the reaction proceeds through direct transmethylation from 38 one MMeHg substituent to another. Solvation effects were found to play an important role in the 39 surface reconstruction of the nanoparticle and in decreasing the energetic barrier for DMeHg 40 formation. We also investigate DMeHg formation from MMeHg mediated by organic dithiols and 41 compare the results to those mediated by the mineral surface. In contrast to the surface reaction, 42 the calculations suggest that dithiol-mediated formation of DMeHg instead proceeds through a 43 mechanism in which Hg–Hg interactions activate a methyl group from one MMeHg substituent 44 for intramolecular transfer to the adjacent MMeHg substituent. These findings refine previously 45 proposed mechanisms of abiotic DMeHg formation and fill gaps in our understanding of global 46 Hg cycling in the environment.

Lian, Peng↗
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 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

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 OSTI2025

Large bandgap observed on the surfaces of EuZn 2 As 2 single crystals

EuM 2 As 2 (M = Zn, Cd, In, Sn etc.) is an excellent material system for studying magnetism-tuned topological properties. However, discrepancies exist between experimental data and theoretical calculations regarding the bulk and surface bandgaps. In this work, cleaved EuZn 2 As 2 crystals are studied using scanning tunneling microscopy/spectroscopy and density functional theory calculations. Triangular-shaped defect-induced modifications in the local density of states help distinguish between Eu-terminated and AsZn-terminated surfaces. While large bandgaps (~1.5 eV at 77 K) are observed on both pristine surfaces, the bandgap size is found to be highly sensitive to local heterogeneity, tending to decrease. By combining experimental observations with theoretical simulations, we conclude that the reduced bandgap in heterogeneous regions arises from Zn vacancies and/or substitution by As atoms, both impacting greater in the Eu surface electronic properties than those in the AsZn surface. This demonstrates the intimate relationship between the electronic structure and magnetism in EuZn 2 As 2 .

interfaces↗
DOE OSTI2020

Localization in the SCAN meta-generalized gradient approximation functional leading to broken symmetry ground states for graphene and benzene

Density functional theory calculations play a central role in understanding chemical and solid-state systems. Progress depends on density functionals that accurately reproduce both energies, for thermochemistry, and properly describe ground states and other properties that are of interest. The Cr dimer, benzene and graphene are especially important benchmark systems for quantum chemistry and condensed matter physics. The Strongly Constrained and Appropriately Normed (SCAN) functional, which is an advanced meta-generalized gradient approximation functional that significantly improves molecular energies is shown to perform poorly for the Cr dimer. This is connected with its poor performance for itinerant solid-state magnets and is a consequence of over localization of electrons, thus illustrating an analogy between the Cr dimer and itinerant magnets. The Cr dimer is a notoriously difficult system for density functionals. However, we additionally find that SCAN predicts an incorrect symmetry broken ground state for 2D graphene and for the benzene molecule, which is surprising considering that ground states of these are known to be well described even by the simplest local density approximation. We show that SCAN overly favors localized spin polarized states, which is a serious deficiency of this approach. Thus, the challenge of finding density functionals that accurately treat both localized and delocalized electronic systems remains.

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

Ultrafast x-ray pump x-ray probe transient absorption spectroscopy: A computational study and proposed experiment probing core-valence electronic correlations in solvated complexes

Femtosecond X-ray pump X-ray probe experiments are currently possible at free electron lasers such as the linac coherent light source (LCLS), which opens new opportunities for studying solvated transition metal-ligand complexes. In order to make the most effective use of these kinds of experiments, it is necessary to determine which chemical properties an X-ray pump pulse will measure. In this paper, we have combined electron cascade calculations and excited-state time-dependent density functional theory calculations to predict the initial state prepared by an X-ray pump and the subsequent X-ray probe spectra at the Fe K-edge in the solvated model transition metal complex, K 4 Fe II (CN) 6 . We find several key spectral features that report on the ligand-field splitting, solvent-solute interactions, and the 3p and 3d electron interactions. We then show how these features could be measured in an experiment.

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

A Noncentrosymmetric Polymorph of LuRuGe

In this work, we report a new polymorph of LuRuGe, obtained in indium flux. This phase exhibits the noncentrosymmetric ZrNiAl-type structure with the space group $Ρ\bar 62m$ as determined by single-crystal X-ray diffraction. This polymorph can convert into another centrosymmetric polymorph (TiNiSi-type structure, space group Pnma ) at high temperatures. In this work, we performed electrical transport, magnetization, and specific heat measurements on this new phase. It shows metallic behavior with a Hall sign change from negative at 2 K to positive at 125 K. LuRuGe exhibits Pauli paramagnetism as the ground state with no local magnetic moments from either the Ru or Lu site. The Debye temperature Θ = 348 K and electronic coefficient γ e = 3.6 mJ K –2 mol –1 are extracted from the low-temperature specific heat data in LuRuGe. We also carried out first-principles density functional theory calculations to map out the electronic band structure and density of states. There are several electronic bands crossing the Fermi level, supporting a multiband scenario consistent with the Hall sign change. The density of states around the Fermi level is mainly from Ru 4d and Ge 4p electrons, indicating a strong hybridization between those atomic orbitals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗