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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2024

Spectroscopic and Computational Evidence of Uranium Dihydrogen Complexes

Dihydrogen complexation, a phenomenon with robust precedent in the transition metal series, is spectroscopically detected for a uranium(III) complex and thereby extended for the first time to the 5f series. The vacant coordination site and low valence of (C 5 H 4 SiMe 3 ) 3 U prove to be key to the reversible formation of (C 5 H 4 SiMe 3 ) 3 U–H 2 (complex 1), and the paramagnetism of the f 3 center facilitates the detection of complex 1 by NMR spectroscopy. Density functional theory calculations reveal that the delocalization of the 5f electron density from (C 5 H 4 SiMe 3 ) 3 U onto the side-on dihydrogen ligand is crucial to complex formation, an unusual bonding situation for an actinide acid–base complex. Here, the spectroscopic and computational results are compared to those reported for lanthanide metallocenes to yield insight into the nature of–and future possibilities for–f-element dihydrogen complexation.

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

Structures and localized vibrational states of defects in graphite by tight-binding calculations

The structural and vibrational properties of pristine graphite and point defects in graphite are studied by tight-binding (TB) calculations using a three-center TB potential model. We showed that the three-center TB potential without ‘‘ad hoc’’ van der Waals interaction corrections can accurately describe the inter-layer distance of graphite and the lowest-energy structures and stabilities of typical point defects in graphite. The results from our TB calculations are in good agreement with those from density-functional theory calculations with van der Waals interaction corrections. We also investigated the vibrational properties to gain better understanding on the localization of vibrational states induced by the point defects. Our calculation results show that although localized or quasi-localized vibrational modes can be found in all defected graphite, the localization induced by Frenkel pair, dual-vacancy, and dual-interstitial defects is much stronger. Atomic displacements associated with the localized vibrational modes induced by these three point defects are also analyzed

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

Magnetic properties of the quasi-one-dimensional S = 1 spin chain antiferromagnet BaNi Te 2 O 7

We report a quasi-one-dimensional S=1 spin chain compound BaNiTe 2 O 7 . This magnetic system has been investigated by magnetic susceptibility, specific heat, and neutron powder diffraction. These results indicate that BaNiTe 2 O 7 develops a short-range magnetic correlation around T ~ 22K. With further cooling, an antiferromagnetic phase transition is observed at T N ~ 5.4K. Neutron powder diffraction revealed antiferromagnetic noncollinear order with a commensurate propagation vector k = (1/2,1,0). The refined magnetic moment size of Ni 2+ at 1.5 K is 1.84 μB , and its noncollinear spin texture is confirmed by first-principles calculations. In conclusion, inelastic neutron-scattering results and density functional theory calculations confirmed the quasi-one-dimensional nature of the spin systems.

1-dimensional spin chains↗
DOE OSTI2021

Analysis of defects in In 2 O 3 :H synthesized in presence of water vapor and hydrogen gas mixture

Using hybrid functional-based density functional theory calculations, we analyze the structure and kinetics of defects formed in two competing synthesis routes to prepare hydrogen-doped In 2 O 3 films, using a hydrogen and oxygen gas mixture vs synthesis in the presence of water vapor. For both of these synthesis routes, we find that H + is the dominant defect species: when the Fermi level is close to the conduction band, H + has a lower formation energy than other intrinsic or extrinsic defects. Our results also suggest that water molecules spontaneously split into H + (which occupies octahedral voids) and OH – interstitials (which occupies vacant oxygen lattice sites or oxygen vacancies). From the analysis of the binding energies between these different defects, we conclude that these defects do not cluster and are most likely to stay spatially distributed throughout the films. In addition, the sum of formation energies of an oxygen (i.e., O$^{2–}_i$) and a H + interstitial is close to the formation energy of a OH – interstitial, meaning that water molecules are completely split into 2H + and O$^{2–}_i$ at the synthesis conditions. Further, in the presence of H 2 + O 2 gas mixture, oxygen interstitials occupy oxygen vacancies while hydrogen interstitials occupy vacant oxygen lattice sites and form bonds with lattice oxygens. Our analysis of the defect equilibria suggests that the hydrogen content in films synthesized in the presence of water vapor is higher than films synthesized in the presence of a hydrogen gas mixture. At high dopant concentrations, a hydrogen bond network is formed in the system and this leads to large distortions in the lattice.

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

Energy barriers for Dy and H penetrating graphene on 6 H -SiC(0001) and freestanding bilayer graphene from first-principles calculations

Currently, intercalation of foreign guest atoms into two-dimensional (2D) layered van der Waals materials is an active research area motivated in part by the development of next-generation energy-storage technologies and optoelectronic devices. One such extensively studied 2D material is the graphene-on-SiC system. To realize and control the desired intercalated structures, it is fundamentally important to understand the kinetic process of intercalation. For the intercalation of a guest atom into graphene layers on SiC substrate, a critical kinetic parameter is the energy barrier of a guest atom penetrating the perfect graphene top layer into the gallery under it. However, accurate theoretical calculations for such penetration barriers are unavailable in literature. From our first-principles density functional theory calculations, we obtain the global energy barriers of 3.47 and 1.80 eV for single Dy and H atoms penetrating the graphene top layer on a graphene buffer layer supported by a Si-terminated 6H-SiC(0001) substrate, respectively. Furthermore, for comparison as well as for examining the lateral strain effects, we also obtain the global barriers of 5.05 and 1.50 eV for single Dy and H atoms penetrating freestanding bilayer graphene with a tensile strain of about 8.8% to match our model for supported graphene, as well as the global barriers of 7.21 and 4.18 eV for penetrating unstrained freestanding bilayer graphene, respectively. From corresponding minimum energy paths with multiple energy minima and saddle points, we can also obtain various local energy barriers and the global backward barrier from the graphene gallery back to the top surface.

2D materials↗
DOE OSTI2023

Integrated Quantum-Classical Protocol for the Realistic Description of Solvated Multinuclear Mixed-Valence Transition-Metal Complexes and Their Solvatochromic Properties

Linear cyanide-bridged polymetallic complexes, which undergo photoinduced metal-to-metal charge transfer, represent prototypical systems for studying long-range electron-transfer reactions and understanding the role played by specific solute–solvent interactions in modulating the excited-state dynamics. Here, to tackle this problem, while achieving a statistically meaningful description of the solvent and of its relaxation, one needs a computational approach capable of handling large polynuclear transition-metal complexes, both in their ground and excited states, as well as the ability to follow their dynamics in several environments up to nanosecond time scales. Here, we present a mixed quantum classical approach, which combines large-scale molecular dynamics (MD) simulations based on an accurate quantum mechanically derived force field (QMD-FF) and self-consistent QMD polarized point charges, with IR and UV–vis spectral calculations to model the solvation dynamics and optical properties of a cyano-bridged trinuclear mixed-valence compound (trans-[(NC) 5 Fe III (μ-CN)Ru II (pyridine) 4 (μ-NC)Fe III (CN) 5 ] 4– ). We demonstrate the reliability of the QMD-FF/MD approach in sampling the solute conformational space and capturing the local solute–solvent interactions by comparing the results with higher-level quantum mechanics/molecular mechanics (QM/MM) MD reference data. The IR spectra calculated along the classical MD trajectories in different solvents correctly predict the red shift of the CN stretching band in the aprotic medium (acetonitrile) and the subtle differences measured in water and methanol, respectively. By explicitly including the solvent molecules around the cyanide ligands and calculating the thermal averaged absorption spectra using time-dependent density functional theory calculations within the Tamm–Dancoff approximation, the experimental solvatochromic shift is quantitatively reproduced going from water to methanol, while it is overestimated for acetonitrile. This discrepancy can likely be traced back to the lack of important dispersion interactions between the solvent cyano groups and the pyridine substituents in our micro solvation model. The proposed protocol is applied to the ground state in water, methanol, and acetonitrile and can be flexibly generalized to study excited-state nonequilibrium solvation dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2023

Structure determination of ζ-N 2 from single-crystal X-ray diffraction and theoretical suggestion for the formation of amorphous nitrogen

The allotropy of solid molecular nitrogen is the consequence of a complex interplay between fundamental intermolecular as well as intramolecular interactions. Understanding the underlying physical mechanisms hinges on knowledge of the crystal structures of these molecular phases. That is especially true for ζ-N 2 , key to shed light on nitrogen’s polymerization. Here, we perform single-crystal X-ray diffraction on laser-heated N 2 samples at 54, 63, 70 and 86 GPa and solve and refine the hitherto unknown structure of ζ-N 2 . In its monoclinic unit cell (space group C 2/ c ), 16 N 2 molecules are arranged in a configuration similar to that of ε-N 2 . The structure model provides an explanation for the previously identified Raman and infrared lattice and vibrational modes of ζ-N 2 . Density functional theory calculations give an insight into the gradual delocalization of electronic density from intramolecular bonds to intermolecular space and suggest a possible pathway towards nitrogen’s polymerization.

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

The role of water vapor during the synthesis of hydrogen doped In 2 O 3

Hydrogen-doped In 2 O 3 synthesized using magnetron sputtering or atomic layer deposition in the presence of water vapor has high transparency in the near infra-red region and mobility values that are almost three to four times of those in commercially used In 2 O 3 :Sn at a much lower carrier concentration. However, simple questions like what happens when water molecules enter In 2 O 3 or why the dominant charge carrier changes from an oxygen vacancy in In 2 O 3 :Sn to H + in In 2 O 3 are not clear. Using hybrid functional based density functional theory calculations, in this study we show that water molecules spontaneously split into H + and OH – . The H + interstitial bonds with an oxygen in In 2 O 3 , and the OH – interstitial occupies a vacant lattice site or is trapped by an oxygen vacancy. Thus, in agreement with experimental findings, our results suggest that the concentration of oxygen vacancies decreases as the partial pressure of water increases. The charge transition level, ϵ ( – / + ) , of a hydrogen interstitial and the Fermi level at which the formation energies of OH – and H + are equal lie above the conduction band minimum, suggesting that the H + interstitial has a lower formation energy. In addition, water molecules do not exhibit any propensity to form clusters, meaning that H + ions can diffuse and remain spatially separated from each other.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Crystal structure of the $τ_{11}$ Al 4 Fe 1.7 Si phase from neutron diffraction and ab initio calculations

The intermetallic τ 11 Al 4 Fe 1.7 Si phase is of interest for high-temperature structural application due to its combination of low density and high strength. We determine the crystal structure of the τ 11 phase through a combination of powder neutron diffraction and density functional theory calculations. Using Pawley and Rietveld refinements of the neutron diffraction data provides an initial crystal structure model. Since Al and Si have nearly identical neutron scattering lengths, we use density-functional calculations to determine their preferred site occupations. The τ 11 phase exhibits a hexagonal crystal structure with space group P6 3 /mmc and lattice parameters of a = 7.478 Å and c = 7.472 Å. The structure comprises five Wyckoff positions; Al occupies the 6h and 12k sites, Fe the 2a and 6h sites, and Si the 2a sites. Here we observe site disorder and partial occupancies on all sites with a large fraction of 80% Fe vacancies on the 2d sites, indicating an entropic stabilization of the τ 11 phase at high temperature.

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

99 TcO 4 - removal from legacy defense nuclear waste by an alkaline-stable 2D cationic metal organic framework

Removal of 99 TcO 4 - from legacy defense nuclear tank waste at Savannah River Site is highly desirable for the purpose of nuclear safety and environmental protection, but currently not achievable given the extreme conditions including high alkalinity, high ionic strength, and strong radiation field. Herein, we present a potential solution to this long-term issue by developing a two-dimensional cationic metal organic framework SCU-103, showing ultrahigh stability in alkaline aqueous media and great resistance to both β and γ radiation. More importantly, it is very effective for 99 TcO 4 - separation from aqueous media as demonstrated by fast exchange kinetics, high sorption capacity, and superior selectivity, leading to the successful removal of 99 TcO 4 - from actual Savannah River Site high level tank waste for the first time, to the best of our knowledge. In addition, the uptake mechanism is comprehensively elucidated by molecular dynamics simulation and density functional theory calculation, showing a unique chemical recognition of anions with low charge density.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2020

Radiative capture rates at deep defects from electronic structure calculations

We present a methodology to calculate radiative carrier capture coefficients at deep defects in semiconductors and insulators from first principles. Electronic structure and lattice relaxations are accurately described with hybrid density functional theory. Calculations of capture coefficients provide an additional validation of the accuracy of these functionals in dealing with localized defect states. We also discuss the validity of the Condon approximation, showing that even in the event of large lattice relaxations the approximation is accurate. We test the method on GaAs:VGa-TeAs and GaN:CN, for which reliable experiments are available, and demonstrate very good agreement with measured capture coefficients.

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

First-Principles Quantum Transport Modeling of Spin-Transfer and Spin-Orbit Torques in Magnetic Multilayers

A unified approach for computing (i) spin-transfer torque in magnetic trilayers like spin valve and magnetic tunnel junction, where injected charge current flows perpendicularly to interfaces, and (ii) spin-orbit torque in magnetic bilayers of the type ferromagnet/spin-orbit-coupled material, where injected charge current flows parallel to the interface, was reviewed. The experimentally explored and technologically relevant spin-orbit-coupled materials include 5d heavy metals, topological insulators, Weyl semimetals, and transition metal dichalcogenides. This approach requires to construct the torque operator for a given Hamiltonian of the device and the steady-state nonequilibrium density matrix, where the latter is expressed in terms of the nonequilibrium Green’s functions and split into three contributions. Tracing these contributions with the torque operator automatically yields field-like and damping-like components of spin-transfer torque or spin-orbit torque vector, which is particularly advantageous for spin-orbit torque where the direction of these components depends on the unknown-in-advance orientation of the current-driven nonequilibrium spin density in the presence of spin-orbit coupling. Illustrative examples are provided by computing spin-transfer torque in a one-dimensional toy model of a magnetic tunnel junction and realistic Co/Cu/Co spin valve, both of which are described by first-principles Hamiltonians obtained from noncollinear density functional theory calculations, as well as by computing spin-orbit torque in a ferromagnetic layer described by a tight-binding Hamiltonian which includes spin-orbit proximity effect within ferromagnetic monolayers assumed to be generated by the adjacent monolayer transition metal dichalcogenide. In addition, it is shown here how spin-orbit proximity effect, quantified by computing (via first-principles retarded Green’s function) spectral functions and spin textures on monolayers of realistic ferromagnetic material like Co in contact with heavy metal or monolayer transition metal dichalcogenide, can be tailored to enhance the magnitude of spin-orbit torque. Furthermore, errors made in the calculation of spin-transfer torque are quantified when using Hamiltonian from collinear density functional theory, with rigidly rotated magnetic moments to create noncollinear magnetization configurations, instead of proper (but computationally more expensive) self-consistent Hamiltonian obtained from noncollinear density functional theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗