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 37 records · Page 2

Exploring Anomalous Photoelectron Angular Distributions in the Photoelectron Spectra of Gd 3 O 3 – : Study of Gd 3 O 2 – and Gd 3 O 3 – Using Photoelectron Spectroscopy and Density Functional Theory Calculations

Anion photoelectron (PE) spectra of lanthanide oxide clusters obtained previously have exhibited anomalous photoelectron angular distributions which were attributed to strong PE–valence electron (PEVE) interactions. Here, to further explore this effect, we have obtained the PE spectra of Gd 3 O 2 – and Gd 3 O 3 – , two clusters that have similarly complex electronic structures but contrasting symmetries. The spectra exhibit manifolds of detachment transitions at similar binding energies in a 0.5 eV window of energy. The electron affinity of Gd 3 O 2 is measured to be 1.29 ± 0.05 eV, and that of Gd 3 O 3 is 1.31 ± 0.05 eV. As seen in previous studies on lanthanide oxide cluster anions in lower than conventional oxidation states, transitions in spectra obtained lower photon energies are more congested than those obtained with higher photon energy, a signature of strong PEVE interactions. While the detachment transitions have predominantly parallel photoelectron angular distributions (PAD), the PAD varies across the manifold of transitions in the PE spectrum of Gd 3 O 3 – in a way that suggests four different subgroups of transitions. Results of calculations on Gd 3 O 2 – suggest kite or V-shape structures with antiferromagnetic coupling between one of the 4f 7 subshells with the two others. Calculations on Gd 3 O 3 – more definitively point to ring structures with a nearly isoenergetic ferromagnetically coupled high spin (24-tet) state and a dectet state in which one of the 4f 7 subshells is antiferromagnetically coupled with the other two. Taking these results as qualitative, we propose that strong mixing between the unperturbed states predicted computationally leads to overlapping transitions with different PADs.

anions↗

Quantum Monte Carlo Calculations of Chemical Binding and Reactions

The auxiliary field quantum Monte Carlo method developed by the PIs has been shown to provide the most accurate description of strongly correlated electronic systems, from molecules to solids. Unlike other explicitly many‐body approaches, the quantum Monte Carlo method scales as a low order polynomial of systems size, similar to mean‐field methods such as density functional theory. However, the auxiliary field quantum Monte Carlo algorithm is significantly more expensive than traditional density functional calculations. This creates a bottleneck for applications to extended systems, such as large molecules and solids. One principal objective of this proposal was to develop new auxiliary field quantum Monte Carlo computational strategies to achieve improved scaling with system size, using downfolding and localization schemes, without sacrificing the predictive power of the calculations. A second goal is to extend the reach of auxiliary field quantum Monte Carlo to calculate excited states. This final report summarizes what has been achieved during the course the project toward these goals.

97 MATHEMATICS AND COMPUTING↗

Electronic structure Fermi liquid theory of high T(sub c) superconductors: Comparison with experiments

For years, there has been controversy on whether the normal state of the Cu-oxide superconductors is a Fermi liquid or some other exotic ground state. However, some experimentalists are clarifying the nature of the normal state of the high T(sub c) superconductors by surmounting the experimental difficulties in producing clean, well characterized surfaces so as to obtain meaningful high resolved photoemission data, which agrees with earlier positron-annihilation experiments. The experimental work on high resolution angle resolved photoemission by Campuzano et al. and positron-annihilation studies by Smedskjaer et al. has verified the calculated Fermi surfaces in YBa2Cu3O7 superconductors and has provided evidence for the validity of the energy band approach. Similar good agreement was found for Bi2Sr2CaCu2O8 by Olson et al. As a Fermi liquid (metallic) nature of the normal state of the high T(sub c) superconductors becomes evident, these experimental observations have served to confirm the predictions of the local density functional calculations and hence the energy band approach as a valid natural starting point for further studies of their superconductivity.

Freeman, A. J.↗

Unoccupied electronic resonances of Sc adsorbed on W(001) by k-resolved inverse photoemission

Scandium adsorbed on the (001) face of tungsten has been studied using Auger-electron spectroscopy, low-energy electron diffraction, k-resolved inverse photoelectron spectroscopy, work-function measurements, and relativistic-electronic-structure calculations. We find that the work function of W(001) does not decrease monotonically as a function of scandium adsorption but reaches a minimum value of Phi = 3.25 eV at 50% coverage. For a complete monolayer, an order 1 X 1 Sc overlayer is formed. Inverse photoemission in the isochromat mode was used to map the unoccupied energy levels of this overlayer along the bar-Gamma(bar-Mu) symmetry direction of the surface Brillouin zone. Local-density-functional calculations using the muffin-tin orbital method were performed for a 1 X 1 Sc overlayer on W and are compared with the experimental two-dimensional band structure.

Lamouri, A.↗

Magnetic interactions and excitations in SrMnSb 2

The magnetic interactions in the antiferromagnetic (AFM) Dirac semimetal candidate SrMnSb 2 are investigated using ab initio linear response theory and inelastic neutron scattering (INS). Our calculations reveal that the first two nearest in-plane couplings (J 1 and J 2 ) are both AFM in nature, indicating a significant degree of spin frustration, which aligns with experimental observations. The orbital resolution of exchange interactions shows that J 1 and J 2 are dominated by direct and superexchange, respectively. In a broader context, a rigid-band model suggests that electron doping fills the minority spin channel and results in a decrease in the AFM coupling strength for both J 1 and J 2 . To better compare with INS measurements, we calculate the spin-wave spectra within a linear spin-wave theory, utilizing the computed exchange parameters. Although the calculated spin-wave spectra somewhat overestimate the magnon bandwidth, they exhibit overall good agreement with measurements from INS experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Coexisting paramagnetic spins and long-range magnetic order in Ba4⁢(Ru0.92⁢Ir0.08)3⁢O10

We investigate the effect of dilute Ir substitution on the magnetism of the trimer-based ruthenate Ba4⁢Ru3⁢O10 using neutron diffraction, magnetic susceptibility measurements, first-principles calculations, and atomistic simulations. Neutron diffraction shows that Ir doping preserves the zigzag antiferromagnetic structure and the ordered-moment magnitude of the parent compound, in which the moments reside exclusively on the two outer Ru(2) sites of each Ru3⁢O12 trimer, while the central Ru(1) site remains nonmagnetic. The Néel temperature is reduced from ∼105 to 84.0(1) K upon 8% Ir substitution, while magnetic susceptibility reveals a pronounced low-temperature Curie-like upturn, indicating the coexistence of paramagnetic spins with long-range antiferromagnetic order. Density-functional calculations show that Ir preferentially occupies the central Ru(1) site, which disrupts the exchange pathways of the affected Ru3⁢O12 trimer state. Atomistic simulations incorporating this paramagnetic dilution reproduce the suppressed ordering temperature and the coexistence of ordered and paramagnetic components.

Islam, Farhan [Ames Laboratory]↗

Fine-tuning catalytic selectivity by modulating catalyst-environment interactions: CO 2 hydrogenation over Pd-based catalysts

Capturing catalytic behaviors under operational conditions is pivotal to gaining a mechanistic understanding and promoting the design of robust catalysts. The challenge lies in the difficulty of monitoring real-time surface dynamics driven by catalyst-environment interactions. Here, in this work, we introduce a framework based on density functional calculations and kinetic modeling. This framework significantly improves the accuracy of theoretical models’ descriptions of experimental observations by quantifying environmental impacts on surface phases and active sites. CO 2 hydrogenation over Pd-based catalysts is taken as a showcase. The observed selectivity variations of Pd and Pd-M bimetallic catalysts strongly correlate with hydrogen coverage maintained under typical CO 2 hydrogenation conditions. By reducing the amount of surface hydrogen, the selectivity tuned effectively from formic acid toward CO and methanol. This study not only deepens the comprehension of dynamics of active sites under active chemical conditions but also introduces an alternative opportunity for catalytic tuning by modulating catalyst-environment interactions.

03 NATURAL GAS↗

Self-compensation of group-V acceptors in CdTe

Cadmium Telluride is at the core of an important thin-film technology for photovoltaics that is already commercially available, yet the CdTe-based solar cell efficiency remains limited at 22%, well below the theoretical limit of ~30%. Increasing the hole concentration is crucial for higher efficiency, and group-V elements such as As, P, and Sb are potential dopants as they are shallow acceptors. Nevertheless, group-V doped p-type CdTe often exhibits low doping activation, and the compensation source remains debated. Here, we performed hybrid density functional calculations with spin-orbit coupling to investigate possible sources of hole compensation in group-V doped CdTe. First, regarding possible self-compensation of the group-V dopants, we find that the formation of AX centers is unlikely since they are found to be unstable relative to the shallow acceptor forms. However, if the group-V dopants come in during growth (such as dimer molecules As2, P2, and Sb2), we find that the impurity atoms, which would occupy nearest neighbor sites, maintain the V-V bonds, limiting the hole density. For the native defects, our study reveals that Cd interstitial is the lowest energy donor defect in p-type CdTe. Still, it has a small migration barrier of 0.5 eV, making it highly mobile and unstable at room temperature. The Te vacancy is the next lowest formation energy donor. The migration barrier of 1.4 eV indicates that the Te vacancy is stable at room temperature. The antisite CdTe is also a donor, with low formation energy and stable at room temperature, potentially limiting the hole concentration. Our results, therefore, shed light on possible compensation centers and some guidance on how to avoid them.

14 SOLAR ENERGY↗

Modeling an S N 2 Reaction Mechanism for Hydrolysis of Siloxane Linkages with Density Functional Theory under Basic Conditions and Implications for Dissolution of Quartz

An extended silicate molecular cluster was used to perform density functional theory calculations to determine if an S N 2 mechanism, where OH – attacks a Q 1 Si coupled with Na + charge balancing an adjacent siloxane bridging O atom, could explain the observed energy of activation of quartz dissolution under basic conditions.

chemical calculations↗

Computational studies of metal-metal and metal-ligand interactions

Accurate calculations on the bond length, dipole moment, and harmonic frequency of CO are presented, using large basis sets and high levels of electron correlation. The geometric structure, force constants and binding energies of Cr(C0)6 and Cr(CO)5 are computed using large basis sets and high levels of electron correlation. The molecule 04(+) is studied, using large basis sets and high levels of electron correlation, including the CASSCF, CASSI and CASPT2 methods. Binding energies, geometries and frequencies are computed. Symmetry breaking is a particular problem for the antisymmetric stretch, which is addressed using the CASSI method. The symmetry breaking problem in 04(+) has also been studied using the Brueckner coupled cluster method. This gives results in good agreement with CASSI. A multi-region numerical integration scheme is investigated for use in Density Functional Calculations. This scheme is found to give comparable results to a widely used scheme based on the Euler-Maclaurin technique.

Barnes, Leslie A.↗

Computational studies of metal-metal and metal-ligand interactions

Accurate calculations on the bond length, dipole moment, and harmonic frequency of CO are presented, using large basis sets and high levels of electron correlation. The geometric structure, force constants and binding energies of Cr(CO)6 and Cr(CO)5 are computed using large basis sets and high levels of electron correlation. The molecule O4(+) is studied, using large basis sets and high levels of electron correlation, including the CASSCF, CASSI,and CASPT2 methods. Binding energies, geometries and frequencies are computed. Symmetry breaking is a particular problem for the antisymmetric stretch, which is addressed using the CASSI method. The symmetry breaking problem in O4(+) has also been studied using the Brueckner coupled-cluster method. This gives results in good agreement with CASSI. A multi-region numerical integration scheme is investigated for use in Density Functional Calculations. This scheme is found to give comparable results to a widely used scheme based on the Euler-Maclaurin technique.

Barnes, Leslie A.↗

Equation of state of boron carbide B 4 ⁢C

We present the results of recent experiments conducted on the Sandia Z machine and a new tabular equation of state for B 4 ⁢C. The equation of state was calibrated to a combination of density functional calculations reported here and fits to preexisting data. It was constructed partly to recover the effects of a shock-driven, polymorphic phase transition of unknown character beginning at particle velocities of just under 3 km/s (shock pressures of 95 GPa). Some of the Z experiments included sound speeds determined by the overtaking rarefaction method, from which we calculate the Grüneisen parameter and compare with previous experiments conducted at the OMEGA laser [Fratanduono et al ., Phys. Rev. B 94 , 184107 (2016)], our own first principles calculations, and another recent tabular equation of state [Zhang et al ., Phys. Rev. E 102 , 053203 (2020)]. We also compare our results with previous static compression, thermophysical, and melt studies, finding mixed consistency. We predict the onset and completion of shock melting at 225 and 265 GPa, respectively, and predict a melt curve that is largely flat to pressures of several hundred GPa.

36 MATERIALS SCIENCE↗

Monolayer TiS 2 Nanosheets on Au(111)–Structural Characterization and Effect of Edge Stability for Shape Control

Transition metal dichalcogenides are promising alternatives to noble metal catalysts, e.g., for (photo-)activation of greenhouse gases or hydrogenations. Herein, a direct synthetic route for 2D TiS 2 nanosheets on Au(111) by titanium deposition in the presence of a mild, organic, non-oxidizing sulfur source is presented. High-resolution scanning tunneling microscopy (STM) is used to gain atomic-level insights into the TiS 2 nanosheet morphology. In contrast to the literature, this protocol gains mostly hexagonal and truncated triangular nanosheets with an increased edge contrast in STM, analog to metallic edge states in MoS 2 . Synchrotron-based photoelectron spectroscopy allows insights into compositional details, specifically to distinguish different S sites on the TiS 2 sheets and other S species on the sample. Further, a minimum size is identified (9 S atoms side length), which underlines the importance of moiré reconstructions for stress relief. The TiS 2 sheets coexist with [Au]Ti 1 S 3 clusters, in which a single gold atom is alloyed into the surface and capped by three S atoms. Together with the finding of a critical sheet size, this points toward on-surface Ostwald ripening as a relevant process in the sheet formation. Ab-initio calculations (density functional theory) underscore that the chemical potential of S is an essential descriptor to maintain shape control.

2D nanomaterials↗

Stabilizing lithium superoxide formation in lithium-air batteries by Janus chalcogenide catalysts

Solid lithium peroxide (Li 2 O 2 ) is the major discharge product in Li-air batteries. However, the electronically insulating nature of Li 2 O 2 tends to affect the battery’s performance such as the polarization gap and cyclability. On the other hand, lithium superoxide (LiO 2 ), generated through a one-electron transfer process, offers greater electronic conductivity, lower charge transfer resistance, and thus reduced charge potential. Nevertheless, LiO 2 long-term stabilization as a final product remains a significant challenge. Here, in this study, we present the molybdenum (Mo)-based Janus chalcogenide family featuring asymmetric structures as a new generation of cathode catalysts for Li-air batteries. These catalysts demonstrate remarkable efficacy in stabilizing LiO 2 discharge products, even under high current densities of 5000 mA/g (corresponding to 0.5 mA/cm 2 ). Our density functional calculations provide an understanding of why the asymmetric Mo-Janus chalcogenides result in LiO 2 formation whereas the symmetric Mo-dichalcogenides produce Li 2 O 2 as the discharge product. These results pave the way to explore a new generation of advanced catalysts for superoxide-based Li-air batteries.

Chalcogenide↗

High-Throughput Screening for Boride Superconductors

A high-throughput screening using density functional calculations is performed to search for stable boride superconductors from the existing materials database. The workflow employs the fast frozen-phonon method as the descriptor to evaluate the superconducting properties quickly. Twenty-three stable candidates were identified during the screening. The superconductivity was obtained earlier experimentally or computationally for almost all found binary compounds. Previous studies on ternary borides are very limited. Here our extensive search among ternary systems confirmed superconductivity in known systems and found several new compounds. Among these discovered superconducting ternary borides, TaMo 2 B 2 shows the highest superconducting temperature of ∼12 K. Most predicted compounds were synthesized previously; therefore, our predictions can be examined experimentally. Our work also demonstrates that the boride systems can have diverse structural motifs that lead to superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Defect Stability in CdTe Based on Formation Energies and Migration Barriers

Native point defects are thought to play a key role in CdTe, either as compensation centers in intentionally doped material, as a source of conductivity in nominally undoped material, or as electron–hole recombination centers. Here, the discussion of their concentration and impact has often centered only on formation energies and transition levels. Using hybrid density functional calculations, including the effects of spin–orbit coupling (SOC), we discuss the stability of native point defects in CdTe based on their formation energies and migration barriers. We show that although Cd interstitials are the lowest energy donor defects, they are unstable at room temperature due to a low migration barrier. They are important for maintaining charge neutrality during growth or annealing at high temperatures, but once the material is brought to room temperature, they are not frozen in as often assumed and are expected to anneal out, leaving the other more stable defects to determine the conductivity. Taking this into account in the solution of the charge neutrality equation, we are able to predict the conductivity type and carrier concentrations that are in good agreement with experimental observations.

14 SOLAR ENERGY↗

Selective Adsorption of Thiol-Containing Molecules on Copper Sulfide Surfaces via Molecule–Surface Disulfide Bridges

Recent results in the fields of nanoenhanced agriculture and expanding interest in prebiotic chemistry have placed increased emphasis on understanding the chemically selective interaction of small molecules with the surfaces of metal sulfides. Here, we present an integrated experimental and computational study of the interaction of thiol-containing molecules with copper sulfide (covellite) surfaces in aqueous media. In situ Fourier-transform infrared (FTIR) measurements and ex situ X-ray photoelectron spectroscopy (XPS) measurements show that molecules bearing free thiol groups, including glutathione and cysteine, bind strongly to CuS (covellite) nanoparticles and to CuS (001) single crystals, while control studies show that similar molecules lacking the free thiol group exhibit much less binding. Additional experiments show that these thiol-containing molecules interact transiently with CuO nanoparticle surfaces but are readily removed by rinsing. The FTIR and XPS experiments demonstrate that adsorption of molecular thiols to CuS surfaces occurs in a chemically selective manner. Further experimental studies and density functional calculations show that the preferred mode of binding is through the surface S atoms, forming a Solid–S–S–Molecule disulfide linkage. While the role of disulfide linkages in controlling structure and function of proteins and other biomolecules is widely known, the formation of surface disulfide linkages as a motif for covalent molecular binding at surfaces has not been established previously.

36 MATERIALS SCIENCE↗

Untangling Sources of Error in the Density-Functional Many-Body Expansion

The many-body expansion provides a framework for data-driven applications of electronic structure theory, including parametrization of classical force fields and machine learning. In this article, we demonstrate that its use significantly amplifies quadrature grid errors when modern density-functional approximations are employed. Standard grids that work well in conventional density-functional calculations result in runaway error accumulation when used with the many-body expansion. At the same time, delocalization error is also exacerbated, leading to exaggerated estimates of nonadditive n-body interactions. This is illustrated for anion–water clusters using the SCAN, r2SCAN, ωB97X-V and ωB97M-V functionals. By employing dense quadrature grids, the inherent self-interaction error is exposed, which can then be mitigated using a variety of other strategies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗