Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hybrid density functional 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

The Fermi–Löwdin self-interaction correction for ionization energies of organic molecules

(Semi)-local density functional approximations (DFAs) suffer from self-interaction error (SIE). When the first ionization energy (IE) is computed as the negative of the highest occupied orbital (HO) eigenvalue, DFAs notoriously underestimate them compared to quasi-particle calculations. Here, the inaccuracy for the HO is attributed to SIE inherent in DFAs. We assessed the IE based on Perdew-Zunger self-interaction correction on 14 small to moderate-sized organic molecules relevant in organic electronics and polymer donor materials. Though self-interaction corrected DFAs were found to significantly improve the IE relative to the uncorrected DFAs, they overestimate. However, when the self-interaction correction is interiorly scaled using a function of the iso-orbital indicator zs , only the regions where SIE is significant get a correction. We discuss these approaches and show how these methods significantly improve the description of the HO eigenvalue for the organic molecules.

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

Evaluating the stability and activity of dilute Cu-based alloys for electrochemical CO 2 reduction

Cu-based catalysts currently offer the most promising route to actively and selectively produce value-added chemicals via electrochemical reduction of CO 2 (eCO 2 R); yet further improvements are required for their wide-scale deployment in carbon mitigation efforts. Here, we systematically investigate a family of dilute Cu-based alloys to explore their viability as active and selective catalysts for eCO 2 R through a combined theoretical–experimental approach. Using a quantum–classical modeling approach that accounts for dynamic solvation effects, we assess the stability and activity of model single-atom catalysts under eCO 2 R conditions. Our calculations identify that the presence of eCO2R intermediates, such as CO*, H*, and OH*, may dynamically influence the local catalyst surface composition. Additionally, we identify through binding energy descriptors of the CO*, CHO*, and OCCO* dimer intermediates that certain elements, such as group 13 elements (B, Al, and Ga), enhance the selectivity of C 2+ species relative to pure Cu by facilitating CO dimerization. The theoretical work is corroborated by preliminary testing of eCO 2 R activity and selectivity of candidate dilute Cu-based alloy catalyst films prepared by electron beam evaporation in a zero-gap gas diffusion electrode-based reactor. Of all studied alloys, dilute CuAl was found to be the most active and selective toward C 2+ products like ethylene, consistent with the theoretical predictions. Finally, we attribute the improved performance of dilute CuAl alloys to more favorable dimerization reaction energetics of bound CO species relative to that on pure Cu. In a broader context, the results presented here demonstrate the power of our simulation framework in terms of rational catalyst design.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

Structural, Optical, and Electronic Properties of Two Quaternary Chalcogenide Semiconductors: Ag 2 SrSiS 4 and Ag 2 SrGeS 4

Quaternary chalcogenide materials have long been a source of semiconductors for optoelectronic applications. Recent studies on the I2-II-IV-X4 (I = Ag, Cu, Li; II = Ba, Sr, Eu, Pb; IV = Si, Ge, Sn; X = S, Se) materials have shown particular versatility and promise among these compounds. These semiconductors take advantage of a diverse bonding scheme and chemical differences among cations to target a degree of anti-site defect resistance. Within this set of compounds, the materials containing both Ag and Sr have not been experimentally studied and leave a gap in the full understanding of the family. Here, we have synthesized powders and single crystals of two Ag- and Sr-containing compounds, Ag 2 SrSiS 4 and Ag 2 SrGeS 4 , each found to form in the tetragonal I4¯2m structure of Ag 2 BaGeS 4 . During the synthesis targeting the title compounds, two additional materials, Ag 2 Sr 3 Si 2 S 8 and Ag 2 Sr 3 Ge 2 S 8 , have also been identified. Furthermore, these cubic compounds represent impurity phases during the synthesis of Ag 2 SrSiS 4 and Ag 2 SrGeS 4 . We show through hybrid density functional theory calculations that Ag 2 SrSiS 4 and Ag 2 SrGeS 4 have highly dispersive band edge states and indirect band gaps, experimentally measured as 2.08(1) and 1.73(2) eV, respectively. Second-harmonic generation measurements on Ag 2 SrSiS 4 and Ag 2 SrGeS 4 powders show frequency doubling capabilities in the near-infrared range.

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

Cubic Crystal Structure Formation and Optical Properties within the Ag–B II –M IV –X (B II = Sr, Pb; M IV = Si, Ge, Sn; X = S, Se) Family of Semiconductors

Quaternary chalcogenide semiconductors are promising materials for energy conversion and nonlinear optical applications, with properties tunable primarily by varying the elemental composition and crystal structure. Here, we first analyze the connections among several cubic crystal structure types, as well as the orthorhombic Ag 2 PbGeS 4 -type structure, reported for select members within the Ag–B II –M IV –X (B II = Sr, Pb; M IV = Si, Ge, Sn; X = S, Se) compositional space. Focusing on the Ag–Pb–Si–S and Ag–Sr–Sn–S systems, we show that one structure type, with the formulas Ag 2 Pb 3 Si 2 S 8 and Ag 2 Sr 3 Sn 2 S 8 , is favored. We have prepared powder and single-crystal samples of Ag 2 Pb 3 Si 2 S 8 and Ag 2 Sr 3 Sn 2 S 8 , showing that each takes on the noncentrosymmetric cubic space group I$\bar{43}$d and is isostructural to the previously reported compound Ag 2 Sr 3 Ge 2 Se 8 . Through hybrid density functional theory calculations, these cubic compounds are demonstrated to be (quasi-)direct band gap semiconductors with high densities of states at the band maxima. Furthermore, the band-gap energies are measured by reflectance spectroscopy as 1.95(3) and 2.66(4) eV for Ag 2 Pb 3 Si 2 S 8 and Ag 2 Sr 3 Sn 2 S 8 , respectively. We further measure the optical properties and show the electronic band structures of three other isostructural A I –B II –M IV –X-type materials, i.e., Ag 2 Sr 3 Si 2 S 8 , Ag 2 Sr 3 Ge 2 S 8 , and Ag 2 Sr 3 Ge 2 Se 8 , showing that the band gaps can be predictably tuned by element substitution. Detailed visual analyses of the different structures and of their relationships with other members of the Ag–B II –M IV –X compositional family provide a basis for a broader understanding of the structure formation and optoelectronic properties within the quaternary chalcogenide semiconductor family.

14 SOLAR ENERGY↗
DOE OSTI2024

Hybrid Density Functional Theory Comparison of Oxygen Release and Solvent Decomposition Kinetics on Li x NiO 2 Surfaces

High-nickel-content layered oxides are among the most promising electric vehicle battery cathode materials. However, their interfacial reactivity with electrolytes and tendency toward oxygen release (possibly yielding reactive 1 O 2 ) remain degradation concerns. Elucidating the most relevant (i.e., fastest) interfacial degradation mechanism will facilitate future mitigation strategies. We apply screened hybrid density functional (HSE06) calculations to compare the reaction kinetics of Li x NiO 2 surfaces with ethylene carbonate (EC) with those of O 2 release. Further, on both the (001) and (104) facets, EC oxidative decomposition exhibits lower activation energies than O 2 release. Our calculations, coupled with previously computed liquid-phase reaction rates of 1 O 2 with EC, strongly question the role of “reactive 1 O 2 ” species in electrolyte oxidative degradation. The possible role of other oxygen species is discussed. To deal with the challenges of modeling Li x NiO 2 surface reactivity, we emphasize a “local structure” approach instead of pursuing the global energy minimum.

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

Factors Affecting the Electron Conductivity in Single Crystal Li 7 La 3 Zr 2 O 12 and Li 7 P 3 S 11

One of the serious challenges in all solid-state Li ion batteries is neutral Li intrusion into the solid-state electrolyte that can ultimately cause catastrophic failure. One possibility for this is due to n-type electron conductivity that induces the reaction Li + + e – → Li 0 at sites where the potential is less than the Li + /Li potential. This paper reports hybrid density functional theory calculations of the electronic conductivity in two prototype single crystalline solid-state electrolytes, cubic Li 7 La 3 Zr 2 O 12 (c-LLZO) and Li 7 P 3 S 11 (LPS). The formation energies of important point defects that can affect electron conductivity are determined, and we find that the mechanism of n-type electron conductivity for both solid-state electrolytes is via “small” electron polaron hopping, where the quotes signify that substantial Li ion rearrangement is associated with the polaron formation and its migration. In both electrolytes, the formation energies for the small polarons at the Fermi energy are too high to generate measurable electron conductivity at room temperature. For c-LLZO, the concentration of electron polarons necessary to ensure charge neutrality from positively charged oxygen vacancies formed in synthesis can be significantly higher. Hence, the electron conductivity could be significant when measured with ion-blocking metal electrodes, and we discuss how the synthesis conditions could affect this magnitude. However, in the solid-state battery, these polarons are replaced by negatively charged Li vacancies so that the electron conductivity should remain minimal. For LPS single crystals, the inherent minimal electron conductivity is independent of synthesis conditions. Here we also show that the cost of forming Li 0 in bulk c-LLZO is enormous due to strain effects so that it could only potentially form at voids, grain boundaries, or around vacancy defects which relax the lattice strain.

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

Insights into the Reactivity of Brookite TiO 2 Nanorods in Liquid Water from Ab Initio Molecular Dynamics Simulations

Brookite TiO 2 , a rare natural polymorph of TiO 2 , has been reported to be an excellent photocatalyst for the production of hydrogen from water and aqueous alcohol solutions, especially when it is reduced and synthesized in the form of nanorods. Here, we investigate the reactivity of stoichiometric and reduced brookite nanorods in liquid water using ab initio molecular dynamics and hybrid density functional theory calculations. Our simulations show a much higher water dissociation fraction on reduced nanorods than on stoichiometric ones, with an accumulation of the resulting bridging hydroxyls (O br H) and terminal hydroxyls (Ti –OH) on different facets of the nanorod. ObrH groups accumulate preferentially on low-energy (210) facets, where they are stabilized by adjacent reduced Ti (Ti 3+ ) sites, while Ti –OH groups prefer to form at the four-fold coordinated Ti atoms on high-energy (010) facets. This hydroxylation pattern also favors the spatial localization of excited holes on the (010) facets. This coupling between water-induced surface chemistry and charge separation underpins the enhanced photocatalytic activity of brookite nanorods, providing useful information for the design of more efficient TiO2-based nanostructures for solar-driven hydrogen evolution.

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

Transition-Metal-Related Quantum Emitters in Wurtzite AlN and GaN

Transition-metal centers exhibit a paramagnetic ground state in wide-bandgap semiconductors and are promising for nanophotonics and quantum information processing. Specifically, there is a growing interest in discovering prominent paramagnetic spin defects that can be manipulated using optical methods. Here, we investigate the electronic structure and magneto-optical properties of Cr and Mn substitutional centers in wurtzite AlN and GaN. We use state-of-the-art hybrid density functional theory calculations to determine level structure, stability, optical signatures, and magnetic properties of these centers. The excitation energies are calculated using the constrained occupation approach and rigorously verified with the complete active space configuration interaction approach. Our simulations of the photoluminescence spectra indicate that $Cr$$^{1+}_{Al}$ in AlN and $Cr$$^{1+}_{Ga}$ in GaN are responsible for the observed narrow quantum emission near 1.2 eV. We compute the zero-field splitting (ZFS) parameters and outline an optical spin polarization protocol for $Cr$$^{1+}_{Al}$ and $Cr$$^{1+}_{Ga}$. Our results demonstrate that these centers are promising candidates for spin qubits.

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

Too big, too small, or just right? A benchmark assessment of density functional theory for predicting the spatial extent of the electron density of small chemical systems

Multipole moments are the first-order responses of the energy to spatial derivatives of the electric field strength. The quality of density functional theory prediction of molecular multipole moments thus characterizes errors in modeling the electron density itself, as well as the performance in describing molecules interacting with external electric fields. However, only the lowest non-zero moment is translationally invariant, making the higher-order moments origin-dependent. Therefore, instead of using the 3 × 3 quadrupole moment matrix, we utilize the translationally invariant 3 × 3 matrix of second cumulants (or spatial variances) of the electron density as the quantity of interest (denoted by K). The principal components of K are the square of the spatial extent of the electron density along each axis. A benchmark dataset of the principal components of K for 100 small molecules at the coupled cluster singles and doubles with perturbative triples at the complete basis set limit is developed, resulting in 213 independent K components. The performance of 47 popular and recent density functionals is assessed against this Var213 dataset. Several functionals, especially double hybrids, and also SCAN and SCAN0 predict reliable second cumulants, although some modern, empirically parameterized functionals yield more disappointing performance. The H, Li, and Be atoms, in particular, are challenging for nearly all methods, indicating that future functional development could benefit from the inclusion of their density information in training or testing protocols.

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

Hydration structure of flat and stepped MgO surfaces

We investigate the solvation structure of flat and stepped MgO(001) in neutral liquid water using ab initio molecular dynamics based on a hybrid density functional with dispersion corrections. Our simulations show that the MgO surface is covered by a densely packed layer of mixed intact and dissociated adsorbed water molecules in a planar arrangement with strong intermolecular H-bonds. The water dissociation fractions in this layer are >20% and >30% on the flat and stepped surfaces, respectively. Slightly above the first water layer, we observe metastable OH groups perpendicular to the interface, similar to those reported in low temperature studies of water monolayers on MgO. Here, these species receive hydrogen bonds from four nearby water molecules in the first layer and have their hydrophobic H end directed toward bulk water, while their associated protons are bound to surface oxygens. The formation of these OH species is attributed to the strong basicity of the MgO surface and can be relevant for understanding various phenomena from morphology evolution and growth of (nano)crystalline MgO particles to heterogeneous catalysis.

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

Computational discovery of stable and metastable ternary oxynitrides

Materials design from first principles enables exploration of uncharted chemical spaces. Extensive computational searches have been performed for mixed-cation ternary compounds, but mixed-anion systems are gaining increased interest as well. Central to computational discovery is the crystal structure prediction, where the trade-off between reliance on prototype structures and size limitations of unconstrained sampling has to be navigated. We approach this challenge by letting two complementary structure sampling approaches compete. We use the kinetically limited minimization approach for high-throughput unconstrained crystal structure prediction in smaller cells up to 21 atoms. On the other hand, ternary—and, more generally, multinary—systems often assume structures formed by atomic ordering on a lattice derived from a binary parent structure. Thus, we additionally sample atomic configurations on prototype lattices with cells up to 56 atoms. Using this approach, we searched 65 different charge-balanced oxide–nitride stoichiometries, including six known systems as the control sample. The convex hull analysis is performed both for the thermodynamic limit and for the case of synthesis with activated nitrogen sources. We identified 34 phases that are either on the convex hull or within a viable energy window for potentially metastable phases. We further performed structure sampling for “missing” binary nitrides whose energies are needed for the convex hull analysis. Among these, we discovered metastable Ce3N4 as a nitride analog of the tetravalent cerium oxide, which becomes stable under slightly activated nitrogen condition µN > +0.07 eV. Given the outsize role of CeO2 in research and application, Ce3N4 is a potentially important discovery.

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

Structural and dynamic properties of solvated hydroxide and hydronium ions in water from ab initio modeling

Predicting the asymmetric structure and dynamics of solvated hydroxide and hydronium in water from ab initio molecular dynamics (AIMD) has been a challenging task. The difficulty mainly comes from a lack of accurate and efficient exchange–correlation functional in elucidating the amphiphilic nature and the ubiquitous proton transfer behaviors of the two ions. By adopting the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation functional in AIMD simulations, we systematically examine the amphiphilic properties, the solvation structures, the electronic structures, and the dynamic properties of the two water ions. In particular, we compare these results to those predicted by the PBE0-TS functional, which is an accurate yet computationally more expensive exchange–correlation functional. We demonstrate that the general-purpose SCAN functional provides a reliable choice for describing the two water ions. Specifically, in the SCAN picture of water ions, the appearance of the fourth and fifth hydrogen bonds near hydroxide stabilizes the pot-like shape solvation structure and suppresses the structural diffusion, while the hydronium stably donates three hydrogen bonds to its neighbors. We apply a detailed analysis of the proton transfer mechanism of the two ions and find the two ions exhibit substantially different proton transfer patterns. In conclusion, our AIMD simulations indicate that hydroxide diffuses more slowly than hydronium in water, which is consistent with the experimental results.

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

Reinforcement learning based hybrid bond-order coarse-grained interatomic potentials for exploring mesoscale aggregation in liquid–liquid mixtures

Exploring mesoscopic physical phenomena has always been a challenge for brute-force all-atom molecular dynamics simulations. Although recent advances in computing hardware have improved the accessible length scales, reaching mesoscopic timescales is still a significant bottleneck. Coarse-graining of all-atom models allows robust investigation of mesoscale physics with a reduced spatial and temporal resolution but preserves desired structural features of molecules, unlike continuum-based methods. Here, we present a hybrid bond-order coarse-grained forcefield (HyCG) for modeling mesoscale aggregation phenomena in liquid–liquid mixtures. The intuitive hybrid functional form of the potential offers interpretability to our model, unlike many machine learning based interatomic potentials. We parameterize the potential with the continuous action Monte Carlo Tree Search (cMCTS) algorithm, a reinforcement learning (RL) based global optimizing scheme, using training data from all-atom simulations. The resulting RL-HyCG correctly describes mesoscale critical fluctuations in binary liquid–liquid extraction systems. cMCTS, the RL algorithm, accurately captures the mean behavior of various geometrical properties of the molecule of interest, which were excluded from the training set. The developed potential model along with the RL-based training workflow could be applied to explore a variety of other mesoscale physical phenomena that are typically inaccessible to all-atom molecular dynamics simulations.

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

Analytical harmonic vibrational frequencies with VV10-containing density functionals: Theory, efficient implementation, and benchmark assessments

VV10 is a powerful nonlocal density functional for long-range correlation that is used to include dispersion effects in many modern density functionals, such as the meta-generalized gradient approximation (mGGA), B97M-V, the hybrid GGA, ωB97X-V, and the hybrid mGGA, ωB97M-V. While energies and analytical gradients for VV10 are already widely available, this study reports the first derivation and efficient implementation of the analytical second derivatives of the VV10 energy. The additional compute cost of the VV10 contributions to analytical frequencies is shown to be small in all but the smallest basis sets for recommended grid sizes. Here, this study also reports the assessment of VV10-containing functionals for predicting harmonic frequencies using the analytical second derivative code. The contribution of VV10 to simulating harmonic frequencies is shown to be small for small molecules but important for systems where weak interactions are important, such as water clusters. In the latter cases, B97M-V, ωB97M-V, and ωB97X-V perform very well. The convergence of frequencies with respect to the grid size and atomic orbital basis set size is studied, and recommendations are reported. Finally, scaling factors to allow comparison of scaled harmonic frequencies with experimental fundamental frequencies and to predict zero-point vibrational energy are presented for some recently developed functionals (including r2SCAN, B97M-V, ωB97X-V, M06-SX, and ωB97M-V).

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

Deep chalcogen donors and electron localization in aluminum nitride

Abstract We show with hybrid density functional theory calculations that chalcogen donors other than oxygen (i.e. S N , Se N , and Te N ) give rise to deep donor states in aluminum nitride. These donors trap a localized electron in their neutral charge state, leading to deep ( + / 0 ) donor levels that are 0.45 eV or more from the conduction-band edge. As such, this behavior is distinct from the DX behavior leads to deep (+/−) levels which affects other donors such as O N and Si Al . We highlight how these results hint at the formation of small electron polarons in AlN, which are found to be unstable in the bulk, but metastable when bound to donor dopants like Si Al and the chalocogens, with activation energies on the order of 0.2–0.3 eV. These results indicate that S, Se, and Te are not shallow donor dopants in aluminum nitride and identify origins of the experimentally observed ∼200–300 meV activation energies for dopant activation in donor-doped samples.

aluminum nitride↗
DOE OSTI2021

Photoexcitation of bulk polarons in rutile TiO 2

The excitation of surface-localized polaronic states has recently been discussed as an additional photocatalytic channel to band gap excitation for rutile TiO 2 . A contribution from photoexcitation of bulk polarons could, in principle, provide a greater contribution because of their higher number and their protection from oxidation. However, determining such a contribution to the photoyield is challenging and has not been achieved thus far. In this work, we use two photon photoemission spectroscopy measurements to separate bulk and surface polaron photoexcitation. We find that bulk polarons are less bound by 0.2 eV compared with polarons at the surface, consistent with our results of hybrid density functional theory calculations. Because the excited state is also shifted to higher energy, bulk polarons have the same photoexcitation resonance energy as at the surface (3.6 eV) with a threshold at 3.1 eV. This is degenerate with the band gap, suggesting that bulk polarons could also provide an additional contribution to the photoyield.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Electronic structure and small-hole polarons in YTiO 3

As a prototypical Mott insulator with ferromagnetic ordering, YTiO 3 (YTO) is of great interest in the study of strong electron correlation effects and orbital ordering. Here we report the first molecular beam epitaxy (MBE) growth of YTO films, combined with theoretical and experimental characterizations of the electronic structure and charge transport properties. The obstacles of YTO MBE growth are discussed and potential routes to overcome them are proposed. DC transport and Seebeck measurements on thin films and bulk single crystals identify p-type Arrhenius transport behavior with an activation energy of ~0.17 eV in thin films, consistent with the energy barrier for small hole polaron migration from hybrid density functional theory calculations. Hard x-ray photoelectron spectroscopy measurements show the lower Hubbard band at 1.1 eV below the Fermi level, whereas a Mott-Hubbard band gap of ~1.5 eV is determined from photoluminescence measurements. These findings provide critical insight into the electronic band structure of YTO and related materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Hybrid density functional study of band gap engineering of SrTi O 3 photocatalyst via doping for water splitting

Perovskite SrTiO 3 (STO) is an attractive photocatalyst for solar water splitting but suffers from a limited photoresponse in the visible spectral range due to its wide band gap. By means of hybrid density functional theory calculations, we systematically study engineering its band gap via doping 4$\textit{d}$ and 5$\textit{d}$ transition-metals $\textit{M}$ (M = Zr, Nb, Mo, Tc, Ru, Rh, Pd, Hf, Ta, W, Re, Os, Ir, and Pt) and chalcogen elements $\textit{Y}$ ($\textit{Y}$ = S and Se).We find that transition-metal dopant M either has no effect on the STO band gap or introduces detrimental midgap states except for Pd and Pt that are able to reduce the STO band gap. In contrast, doping S and Se significantly reduces STO’s direct band gap, thus, leading to appreciable optical absorption transitions in the visible spectral range. In this work, our findings provide that Pd-, S-, and Se-doped STO are potential promising photocatalysts for water splitting under visible light irradiation thereby providing insightful theoretical guides for experiments to improve the photocatalytic activity of STO.

36 MATERIALS SCIENCE↗