Search NASA⌕ Search

SEARCH · Search NASA

Results for “hybrid functional”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Migration of Ga vacancies and interstitials in β – Ga 2 O 3

Pathways and energy barriers for the migration of Ga vacancies (V Ga ) and Ga interstitials (Ga i ) in β–Ga 2 O 3 are explored using hybrid functional calculations and the nudged elastic band method. Considering β–Ga 2 O 3 as primarily being an n-type semiconductor, we focus on defect charge states relevant under such conditions: $V^{3 –}_{Ga}$, Ga$^{3+}_{i}$, and Ga$^{+}_{i}$. Notably, we describe a mechanism by which V Ga can transform between its different split configurations. In all cases, the intermediate state consists of a vacancy split between three Ga sites—a three-split vacancy—which enables passage over a significantly lower energy barrier. This is because it avoids the unfavorable simple vacancy at the tetrahedral Ga site. The proposed mechanism lowers the overall barrier for $V^{3–}_{Ga}$ diffusion along the [001] crystal direction from 1.73 to 0.97 eV, whereas the 2.08 eV barrier for the [100] and [010] directions is unaffected. For Ga$^{3+}_{i}$, we obtain similar overall migration barriers of 0.72, 0.80, and 1.02 eV for the [010], [001], and [100] directions, respectively. Ga$^{+}_{i}$ exhibits a strong preference for diffusion within the large eight-sided channel; the overall migration barrier is 0.92 eV for the [010] direction, and 2.16 eV for the [001] and [100] directions. The limiting step for the two latter directions involves ionization of Ga$^{+}_{i}$ followed by a jump to an adjacent large eight-sided channel as Ga$^{3+}_{i}$. Finally, our results are discussed in light of experimental observations of thermally activated recovery processes in irradiated material.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Defect Properties, Anion Ordering, and Photochromic Mechanism in Yttrium Oxyhydride

Yttrium oxyhydride (YHO) undergoes a reversible photochromic transition when exposed to ultraviolet light. However, the mechanism for this transformation is not fully understood, and the structure and precise chemical composition of YHO remain under debate. Here, we use first-principles density functional theory calculations with a hybrid functional to study the structure, chemical stability, and point defect properties of YHO. As experiments have shown, we find that YHO prefers a cubic structure, with H and O anions present in equal concentrations and located on tetrahedral sites. Stoichiometric and ordered YHO is chemically stable, but it has a wide band gap of 5.01 eV, considerably larger than that measured in experiments (2.4–3.8 eV). On the other hand, Y4H10O has a smaller band gap of 2.97 eV and also has a region of chemical stability; thus, the actual material may include some fraction of this H-rich structure. The defect chemistry of YHO is dominated by anionic antisite species (H O and O H ), with hydrogen interstitials (H i ) and vacancies (V H ) also present in reasonably high concentrations. We show that antisite disorder lowers the band gap relative to the perfectly ordered structure, bringing the magnitude of the gap into closer agreement with experiment. Based on our calculations of defect migration and the positions of defect states relative to the band edges, we link the onset of photochromic behavior to the reaction H O – → V O 0 + H i – , which follows photoexcitation of a H O + defect. H i – can subsequently migrate away and be trapped by additional H O + defects, contributing to the persistence of the reaction, while the resultant oxygen vacancy, V O 0 , introduces an occupied defect state that leads to optical absorption at visible wavelengths. Our results can explain reported discrepancies between experimental and computational results for YHO, and they allow us to propose specific atomic-scale processes that can lead to photochromism. In conclusion, understanding these mechanisms is key for unlocking YHO’s application in devices ranging from smart windows and optoelectronics to electrochemical synapses for neural networks.

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

Ab initio calculations of the rate of carrier trapping and release at dopant sites in NaI: Tl beyond the harmonic approximation

We present ab initio calculations of electron capture and release coefficients as a function of temperature for Tl dopants in the widely used scintillator NaI, which is a soft ionic crystal. The modelled capture and release events occur by transitions mediated by multiple phonon absorption and emission between states with significantly different local geometries around the trapping site. We demonstrate that such transitions are not well-described by the normal harmonic approximation to the nuclear dynamics. We avoid the harmonic approximation by numerically solving the vibrational Schrödinger equation for the motion along/of a single effective phonon coordinate. The localized trapped state is not correctly described by semilocal density functionals, so we employ a hybrid functional tuned to reproduce the band gap and lattice constant of the host material. Our calculations, which combine plane-wave and cluster electronic structure methods in a novel way, are in reasonable agreement with available experimental data for detrapping and predict an unusual temperature dependence for trapping rates.

Prange, Micah P.↗
DOE OSTI2023

Energy Frontier Research Centers: Center for the Computational Design of Functional Layered Materials (CCDM) August 1, 2014 - July 31, 2018; Center for Complex Materials from First Principles (CCM) August 1, 2018 - July 31, 2021 (Final Report)

The mission of the DOE Energy Frontier Research Centers CCDM (2014-2018) and CCM (2018-2021) was to theoretically develop, computationally apply, and experimentally validate electronic structure methods for all materials, with a focus on the complex materials, especially layered and two-dimensional materials, strongly-correlated materials, and liquid water. This was achieved by over 200 published journal articles authored by about 17 senior investigators from physics and chemistry and from theory, computation, and experiment, plus their collaborators. In particular, the Centers confirmed the predictive power of the SCAN (strongly constrained and appropriately normed) density functional, which was constructed to satisfy 17 known exact constraints and several appropriate norms. Without being fitted to real bonded systems, and at a modest computational cost, SCAN correctly predicted covalent, ionic, metallic, hydrogen, and van der Waals bonds in many challenging materials. SCAN gave an improved description of defects in semiconductors, surface properties of metals, seven phases of ice, liquid water, liquid and supercooled silicon, subtle structural distortions in ferroelectrics, formation energies and structural predictions for solids, and critical pressures for structural phase transitions. Perhaps most remarkably, SCAN correctly described some strongly-correlated materials that were previously believed to be beyond the reach of density-functional approximations. SCAN is the only density functional that correctly predicts the band gap closing under chemical doping of the cuprate high-temperature superconducting materials. SCAN also predicts a landscape of competing stripe and magnetic phases in the cuprates. For some materials with some codes, SCAN has convergence problems that are greatly reduced by the CCM-developed r 2 SCAN, without loss of accuracy or rigor. SCAN and r 2 SCAN still make some self-interaction error, which is greatly reduced by the CCDM/ CCM-developed local orbital scaling correction (LOSC). These Centers further proved that the fundamental energy gaps of a solid from an orbital energy difference and from total energy differences are the same for a large class of generalized Kohn-Sham (GKS) functionals, including SCAN and standard hybrid functionals, and that symmetry breaking arises when a dynamic density fluctuation drops to zero frequency. The Centers identified new mechanisms for catalysis in layered materials with ions intercalated between the layers, investigated charge density waves both in model systems and in real layered materials, studied changes of band gap with the number of layers, and explored topological ultrathin films, bent nanoribbons, and defects.

08 HYDROGEN↗
NASA NTRS1995

The Structure of Li7 and K7

The self-consistent-field (SCF) approach and density functional theory, using the B3LYP hybrid functional, yield three low-lying structures for Li7(-). The relative separations differ for the SCF and B3LYP approaches, however the B3LYP results are in good agreement with the coupled cluster results. For K7(-), only an octahedron with one face capped is found to be a minimum; this the second most stable structure for Li7(-). A comparison of the computed separations between the low-lying states of K7 and the photoelectron detachment spectra does not allow an unambiguous assignment of the structure of K7(-).

Bauschlicher. Charles W., Jr.↗
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 OSTI2021

57Fe Mössbauer isomer shift of pure iron and iron oxides at high pressure—An experimental and theoretical study

The 57Fe isomer shift (IS) of pure iron has been measured up to 100 GPa using synchrotron Mössbauer spectroscopy in the time domain. Apart from the expected discontinuity due to the α→ε structural and spin transitions, the IS decreases monotonically with increasing pressure. The absolute shifts were reproduced without semi-empirical calibrations by periodic density functional calculations employing extensive localized basis sets with several common density functionals. However, the best numerical agreement is obtained with the B1WC hybrid functional. Extension of the calculations to 350 GPa, a pressure corresponding to the Earth’s inner core, predicted the IS range of 0.00 to -0.85 mm/s, covering the span from Fe(0) to Fe(VI) compounds measured at ambient pressure. The calculations also reproduced the pressure trend from polymorphs of prototypical iron oxide minerals, FeO and Fe2O3. Analysis of the electronic structure shows a strong donation of electrons from oxygen to iron at high pressure. The assignment of formal oxidation to the Fe atom becomes ambiguous under this condition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1996

Improved Damage Resistant Composite Materials Incorporating Shape Memory Alloys

Metallic shape memory alloys (SMA) such as nitinol have unique shape recovery behavior and mechanical properties associated with a material phase change that have been used in a variety of sensing and actuation applications. Recent studies have shown that integrating nitinol-SMA actuators into composite materials increases the composite material's functionality. Hybrid composites of conventional graphite/epoxy or glass/epoxy and nitinol-SMA elements can perform functions in applications where monolithic composites perform inadequately. One such application is the use of hybrid composites to function both in load bearing and armor capacities. While monolithic composites with high strength-to-weight ratios function efficiently as loadbearing structures, because of their brittle nature, impact loading can cause significant catastrophic damage. Initial composite failure modes such as delamination and matrix cracking dissipate some impact energy, but when stress exceeds the composite's ultimate strength, fiber fracture and material perforation become dominant. One of the few methods that has been developed to reduce material perforation is hybridizing polymer matrix composites with tough kevlar or high modulus polyethynylene plies. The tough fibers increase the impact resistance and the stiffer and stronger graphite fibers carry the majority of the load. Similarly, by adding nitinol-SMA elements that absorb impact energy through the stress-induced martensitic phase transformation, the composites' impact perforation resistance can be greatly enhanced. The results of drop-weight and high velocity gas-gun impact testing of various composite materials will be presented. The results demonstrate that hybridizing composites with nitinol-SMA elements significantly increases perforation resistance compared to other traditional toughening elements. Inspection of the composite specimens at various stages of perforation by optical microscope illustrates the mechanisms by which perforation is initiated. Results suggest that the out-of-plane transverse shear properties of the composite and nitinol elements have a significant effect on the perforation resistance. Applications that can utilize the hybrid composites effectively will also be presented with the experimental studies.

Paine, Jeffrey S. N.↗
NASA NTRS1995

Using Density Functional Theory (DFT) for the Calculation of Atomization Energies

The calculation of atomization energies using density functional theory (DFT), using the B3LYP hybrid functional, is reported. The sensitivity of the atomization energy to basis set is studied and compared with the coupled cluster singles and doubles approach with a perturbational estimate of the triples (CCSD(T)). Merging the B3LYP results with the G2(MP2) approach is also considered. It is found that replacing the geometry optimization and calculation of the zero-point energy by the analogous quantities computed using the B3LYP approach reduces the maximum error in the G2(MP2) approach. In addition to the 55 G2 atomization energies, some results for transition metal containing systems will also be presented.

Bauschlicher, Charles W., Jr.↗
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 OSTI2026

Mott vs Kondo: Influence of various density functional based methods on the Ce isostructural phase transition mechanism

The cerium iso-structural phase transition (γ–α) is dominated by f-electron localization changes that results in a magnetic ordering change and a volume collapse. Generally, these physics are difficult to capture with ab initio and first principles methods. However, previous works have shown various methods to be successful in predicting at least some of the physics of the γ to α phase transition. Therefore, here, we perform a broad survey of density functional based methods across three levels of theory and types of functions [GGA (generalized gradient approximation), MetaGGA, and hybrid functionals] and compare the results, focusing on hydrostatic compression across the phase boundary at zero kelvin. For the methods that best reproduce experimental results, we directly probe the predicted mechanisms and frame the results in the Mott/Kondo debate, assessing how the underlying methods and assumptions of different functionals can assess the physical drivers in the phase transition, providing insight into the governing processes of this unique phase transition.

36 MATERIALS SCIENCE↗
NASA NTRS1996

The Successive OH Binding Energies of Sc(OH)n+ for n=1-3

The geometries of Sc(OH)n+, for n = 1-3, have been optimized using density functional theory, in conjunction with the B3LYP hybrid functional. The zero-point energies are computed at the same level of theory. The successive OH bond energies have been computed at the CCSD(T) level for ScOH+ and Sc(OH)2+. The computed result for ScOD+ is in excellent agreement with the recent experiment of Armentrout and co-workers. There is a dramatic drop for the third OH, because Sc+ has only two valence electrons and therefore the bonding changes when the third OH is added. The difference between the B3LYP and CCSD(T) OH binding energies for the first two OH groups is discussed.

Bauschlicher, Charles W., Jr.↗
NASA NTRS1997

Quantum Chemistry Study of Cycloaddition Pathways for the Reaction of o-Benzyne with Fullerenes and Carbon Nanotubes

Functionalization of fullerenes via the [2+2] cycloaddition reaction with o-benzyne has been demonstrated in the laboratory. In contrast, [2+4) cycloaddition products are formed when benzyne reacts with planar polycyclic aromatic hydrocarbons. Using density functional theory (DFT) calculations with Becke's hybrid functional and small contracted gaussian basis sets, we are able to reproduce these product preferences. The objective of this work is to explore the functionalization of carbon nanotubes. We have studied o-benzyne cycloaddition products with a [14,0] single-walled nanotube. We find both the [2+2] and [2+4] adducts to be stable, with the latter product being somewhat favored.

Jaffe, Richard↗
DOE OSTI2020

Self-trapped hole and impurity-related broad luminescence in $β$-Ga 2 O 3

Here we explore the luminescence properties of self-trapped holes and impurity-related acceptors using one-dimensional configuration coordinate diagrams derived from hybrid functional calculations. The photoluminescence spectrum of as-grown $β$-Ga 2 O 3 typically consists of a broad band in the wavelength region from ultraviolet to green and is often dominated by an impurity independent ultraviolet band that is commonly attributed to self-trapped holes. Here, we use the self-trapped hole as a benchmark to evaluate the accuracy of the theoretical defect luminescence spectra and estimate the optical properties of Mg Ga , Be Ga , Ca Ga , Cd Ga , Zn Ga , Li Ga , and N O acceptor impurities, as well as their complexes with hydrogen donors. We also explore $V$ Ga acceptors complexed with hydrogen and Si Ga donor impurities. The results show that these defects can give rise to broad luminescence bands peaking in the infrared to visible part of the spectrum, making them potential candidates for the defect origin of broad luminescence bands in $β$-Ga 2 O 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2022

Nonunique fraction of Fock exchange for defects in two-dimensional materials

By investigating the vacancy and substitutional defects in monolayer WS 2 with hybrid functionals, we find that there is no unique amount of Fock exchange that concurrently satisfies the generalized Koopmans' condition and reproduces the band-gap and band-edge positions. Fixing the mixing parameter of Fock exchange based on the band gap can lead to qualitatively incorrect defect physics in two-dimensional materials. Instead, excellent agreement is achieved with both experiment and many-body perturbation theory within the GW approximation once the mixing parameters are tuned individually for the defect species and the band edges. Here, we show the departure from a unique optimized mixing parameter is inherent to two-dimensional systems as the band edges experience a reduced screening while the localized defects are subject to bulklike screening.

2-dimensional systems↗
DOE OSTI2026

Optimization of Random Phase Approximation Calculations for Improved Energies of Molecules, Solids, and Surfaces

We present an optimized random phase approximation method (optRPA26) that significantly improves upon conventional RPA through an optimized choice of reference orbitals and energy components, rather than a modification of the RPA correlation functional itself. The method employs an empirically constructed hybrid functional to generate DFT orbitals to evaluate the RPA correlation energy, which is then scaled by a constant. Comprehensive benchmarks across molecules, bulk solids, and surface systems demonstrate that optRPA26 consistently achieves high accuracy, with mean absolute errors of 0.05 eV for W4-11-RE reaction energies, 0.07 eV for cohesive energies, 0.09 eV for metal oxide formation energies, 0.11–0.12 eV for adsorption of small molecules on metals, and 0.06 eV for adsorption on oxides. In addition, optRPA26 correctly captures phase stability in metal oxides and magnetic metals. The optRPA26 approach can be run using standard RPA implementations, highlighting its potential as a general-purpose reference method that can accurately capture covalent, ionic, metallic, and van der Waals bonding in molecules, solids, and interfaces.

Adsorption↗
DOE OSTI2021

Electronic Structure of Metallophthalocyanines, MPc (M = Fe, Co, Ni, Cu, Zn, Mg) and Fluorinated MPc

In this work, we compute the electronic structure and optical excitation energies of metal-free and transition-metal phthalocyanines (H2Pc and MPc for M = Fe, Co, Ni, Cu, Zn, Mg) using density functional theory with optimally tuned range-separated hybrid functionals (OT-RSH). We show that the OT-RSH approach provides photoemission spectra in quantitative agreement with experiments as well as optical band gaps within 10% of their experimental values, capturing the interplay of localized d-states and delocalized pi-pi* states for these organometallic compounds. We examine the tunability of MPcs and H 2 Pc through fluorination, resulting in quasi-rigid shifts of the molecular orbital energies by up to 0.7 eV. Our comprehensive data set provides a new computational benchmark for gas-phase phthalocyanines, significantly improving upon other density-functional-theory-based approaches.

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

Accurate parameterization of the kinetic energy functional for calculations using exact-exchange

In this report, electronic structure calculations based on Kohn-Sham density functional theory (KSDFT) that incorporate exact exchange or hybrid functionals are associated with large computational expense, a consequence of the inherent cubic scaling bottleneck and large associated prefactor, which limits the length and time scales that can be accessed. Though orbital-free density functional theory (OFDFT) calculations scale linearly with system size and are associated with significantly smaller prefactor, they are limited by the absence of accurate density-dependent kinetic energy functionals. Therefore, the development of accurate density-dependent kinetic energy functionals is important for OFDFT calculations of large realistic systems. To this end, we propose a method to train kinetic energy functional models at the exact-exchange level of theory by using a dictionary of physically relevant terms that have been proposed in the literature in conjunction with linear or nonlinear regression methods to obtain the fitting coefficients. For our dictionary, we use gradient expansion of the kinetic energy, nonlocal models proposed in the literature and their nonlinear combinations, such as a model that incorporates spatial correlations between higher order derivatives of electron density at two points. The predictive capabilities of these models are assessed by using a variety of model one-dimensional systems that exhibit diverse bonding characteristics, such as a chain of eight hydrogens, LiF, LiH, C 4 H 2 , C 4 N 2 and C 3 O 2 . We show that by using data from model one-dimensional KSDFT calculations performed using the exact-exchange functional for only a few neutral structures, it is possible to generate models with high accuracy for charged systems as well as electron and kinetic energy densities during self-consistent field iterations. In addition, we show that it possible to learn both the orbital dependent terms, i.e., the kinetic energy and the exact-exchange energy, and models that incorporate additional nonlinearities in spatial correlations, such as a quadratic model, are needed to capture subtle features of the kinetic energy density that are present in exact-exchange-based KSDFT calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗