Search NASA⌕ Search

SEARCH · Search NASA

Results for “optimized energy density functionals”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Atomic Data and Spectral Line Intensities for Ne III

A number of satellites and rockets have been launched to observe radiation from the Sun and other astrophysical objects. Line radiation is emitted when the electron impact excited levels decay to the lower levels by photon emission. From this radiation, the physical parameters such as electron temperature and density of the astrophysical plasma, elemental abundance, and opacity can be inferred. Ne III lines have been observed in H II regions, Ne-rich filaments in supernovae, and planetary nebulae. The allowed line at 489.50 Angstroms due to the transition 2s(sup 2) 2p(sup 5) (sup 3) P2 (goes to) 2s(sup 2)2p(sup 4)(sup 3)P2 has been identified in the solar spectrum by Vernazza and Reeves using Skylab observations. Other Ne III lines in the solar EUV spectrum have been reported by Thomas and Neupert based on observations from the Solar EUV Rocket Telescope and Spectrograph (SERTS) instrument. Atomic data for Ne III have been calculated by using a set of programs developed at, University College, London. The Superstructure and Distorted Wave (DW) programs have been updated over the years. In the Superstructure program, configuration interaction can be taken into account and radial functions are calculated in a modified Thomas-Fermi-Amaldi potential. This is a statistical potential and depends on parameters lambda 1 which are determined by optimizing the weighted sum of term energies. They are found to be lambda(sub 0)=1.2467, lambda(sub 1)=1.1617, and lambda(sub 2)=1.0663. The relativistic corrections are included by using the Breit-Pauli Hamiltonian as a perturbation to the nonrelativistic Hamiltonian. The same potential is used to calculate reactance matrices in the DW approximation in LS coupling. Collision strengths in intermediate coupling are obtained by using term coupling coefficients obtained from the Superstructure program. In this calculation, the configurations used are 2s(sup 2)2p(sup 4), 2s2p(sup 5), 2s(sup 2)2p(sup 3)3s, 2s(sup 2)p(sup 3)3d giving rise to 57 fine-structure levels in intermediate coupling.

Bhatia, A. K.↗
DOE OSTI2021

Computing molecular excited states on a D-Wave quantum annealer

Abstract The possibility of using quantum computers for electronic structure calculations has opened up a promising avenue for computational chemistry. Towards this direction, numerous algorithmic advances have been made in the last five years. The potential of quantum annealers, which are the prototypes of adiabatic quantum computers, is yet to be fully explored. In this work, we demonstrate the use of a D-Wave quantum annealer for the calculation of excited electronic states of molecular systems. These simulations play an important role in a number of areas, such as photovoltaics, semiconductor technology and nanoscience. The excited states are treated using two methods, time-dependent Hartree–Fock (TDHF) and time-dependent density-functional theory (TDDFT), both within a commonly used Tamm–Dancoff approximation (TDA). The resulting TDA eigenvalue equations are solved on a D-Wave quantum annealer using the Quantum Annealer Eigensolver (QAE), developed previously. The method is shown to reproduce a typical basis set convergence on the example $$\hbox {H}_2$$ H 2 molecule and is also applied to several other molecular species. Characteristic properties such as transition dipole moments and oscillator strengths are computed as well. Three potential energy profiles for excited states are computed for $$\hbox {NH}_3$$ NH 3 as a function of the molecular geometry. Similar to previous studies, the accuracy of the method is dependent on the accuracy of the intermediate meta-heuristic software called qbsolv.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Optimizing accuracy and efficacy in data-driven materials discovery for the solar production of hydrogen

The production of hydrogen fuels, via water splitting, is of practical relevance for meeting global energy needs and mitigating the environmental consequences of fossil-fuel-based transportation. Water photoelectrolysis has been proposed as a viable approach for generating hydrogen, provided that stable and inexpensive photocatalysts with conversion efficiencies over 10% can be discovered, synthesized at scale, and successfully deployed. While a number of first-principles studies have focused on the data-driven discovery of photocatalysts, in the absence of systematic experimental validation, the success rate of these predictions may be limited. We address this problem by developing a screening procedure with co-validation between experiment and theory to expedite the synthesis, characterization, and testing of the computationally predicted, most desirable materials. Here, starting with 70 150 compounds in the Materials Project database, the proposed protocol yielded 71 candidate photocatalysts, 11 of which were synthesized as single-phase materials. Experiments confirmed hydrogen generation and favorable band alignment for 6 of the 11 compounds, with the most promising ones belonging to the families of alkali and alkaline-earth indates and orthoplumbates. This study shows the accuracy of a nonempirical, Hubbard-corrected density-functional theory method to predict band gaps and band offsets at a fraction of the computational cost of hybrid functionals, and outlines an effective strategy to identify photocatalysts for solar hydrogen generation.

08 HYDROGEN↗
DOE OSTI2022

Resolving atomistic structure and oxygen evolution activity in nickel antimonates

The oxygen evolution reaction (OER) requires electrodes that are not only catalytically active, but also stable under harsh electrochemical environments to enable efficient, durable technologies. Our recent report of a stable amorphous Ni 0.5 Sb 0.5 O z OER photoanode established Ni–Sb–O as an important system for computational understanding of both the structural and catalytic behavior of these complex oxides. In the present work we show that Ni x Sb 1-x O z with x > 0.33 crystallizes into a previously unknown phase. Guided by experimental X-ray diffraction, here, we use density functional theory calculations to perform a prototype phase search to identify a broad family of stable and metastable mixed rutile and hexagonal-like phases for x = 0.33, 0.50, and 0.66 compositions. For the identified phases, we predict favorable oxygen vacancy formation energies for Ni-rich compositions under the reducing synthesis conditions which match measured Ni K-edge X-ray absorption spectra. The calculated overpotential for the most active site decreases with increasing Ni content, from 0.91 V (x = 0.33) to 0.49 V (x = 0.66), which captures the experimentally observed trend. We find the active site changes from the Ni–O–Sb bridge to a Ni–O–Ni bridge at increasing Ni concentrations, rather than the commonly studied singly under-coordinated sites. Finally, detailed Pourbaix analysis of the identified phases show excellent electrochemical stability, consistent with experimentally measured low metal ion concentrations in the electrolyte of photoelectrochemical cells. Collectively, our consideration of an ensemble of structures enables identification of the most catalytically prolific structural motifs, aiding the understanding of crystalline and amorphous catalysts and elucidating the co-optimization of activity and durability in nickel antimonates.

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

Insights into Supported Subnanometer Catalysts Exposed to CO via Machine-Learning-Enabled Multiscale Modeling

Subnanometer catalysts offer high noble metal utilization and superior performance for several reactions. However, understanding their structures and properties on an atomic scale under working conditions is challenging due to the large configurational space. Here, we introduce an efficient multiscale framework to predict their stability exposed to an adsorbate. The framework integrates a comprehensive toolset including density functional theory (DFT) calculations, cluster expansion, machine learning, and structure optimization. The end-to-end machine-learning workflow guides DFT data generation and enables significant computational acceleration. We demonstrate the approach for CO-adsorbed Pdn (n = 1–55) clusters on CeO 2 (111). Simulation results reveal that CO can facilitate restructuring by stabilizing smaller planar structures and bilayer structures of specific intermediate sizes, consistent with experimental reports. Metal–support interactions, preferential CO adsorption, and metal nuclearity and structure control catalyst stability. As a result, the framework allows automatic discovery of stable catalyst structures and a systematic strategy to exploit properties in the subnanometer scale.

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

First-Principles Modeling of Sodium Ion and Water Intercalation into Titanium Disulfide Interlayers for Water Desalination

Recent experiments revealed the possibility of using titanium disulfide (TiS 2 ) as the cathode material in capacitive deionization (CDI) devices for water desalination. Although it performed stably up to 70 cycles with a salt removal capacity of 14 mg/g (corresponding to a sodium ion removal capacity of 35.8 mg/g) at high molar concentration (600 mM NaCl), the maximum capacity of TiS 2 as a CDI electrode was much lower (~70 mAh/g) than as a supercapacitor (239 mAh/g). Understanding why the ion capacities of these two configurations of the same material differ will entail elucidating detailed charge/discharge mechanisms at the atomic scale. Here, we present density functional theory simulations of sodium intercalation into TiS 2 interlayers in order to gain such understanding. We systematically investigated TiS 2 stacking patterns and ion intercalation sites, energetics of the intercalated compounds, and phase transformation during sodium intercalation. The calculated structural evolution and capacitance–voltage curve agree quite well with previous measurements. We conclude that the different maximum capacities of TiS 2 measured in aqueous and dry environments originate from weaker interlayer interactions with respect to shear strain after ~33% intercalation of Na + –H 2 O pairs, which is detrimental to the mechanical stability of TiS 2 . In conclusion, this study sheds light on the underlying mechanisms of ion intercalation into layered materials and contributes to understanding requirements for future design and optimization of CDI electrode materials for water desalination.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Modern quantum chemistry with [Open]Molcas

MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree–Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functions.

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

Analysis of core asymmetries in inertial confinement fusion implosions using three-dimensional hot-spot reconstruction

Three-dimensional effects play a crucial role during the hot-spot formation in inertial confinement fusion (ICF) implosions. A data analysis technique for 3D hot-spot reconstruction from experimental observables has been developed to characterize the effects of low modes on 3D hot-spot formations. In nuclear measurements, the effective flow direction, governed by the maximum eigenvalue in the velocity variance of apparent ion temperatures, has been found to agree with the measured hot-spot flows for implosions dominated by mode ℓ = 1. Asymmetries in areal-density (ρR) measurements were found to be characterized by a unique cosine variation along the hot-spot flow axis. In x-ray images, a 3D hot-spot x-ray emission tomography method was developed to reconstruct the 3D hot-spot plasma emissivity using a generalized spherical-harmonic Gaussian function. The gradient-descent algorithm was used to optimize the mapping between the projections from the 3D hot-spot emission model and the measured x-ray images along multiple views. Furthermore, this work establishes a platform to analyze 3D low-mode core asymmetries in ICF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effect of the Precursor Chemistry on the Crystallization of Triple Cation Mixed Halide Perovskites

Triple cation, mixed halide perovskite compositions have been reported to be more thermally stable, exhibit fewer phase impurities, and show higher power conversion efficiency and better reproducibility than single cation perovskites. In this work, we explain the formation of Cs0.05FA0.81MA0.14Pb(I0.85Br0.15)3 via a multimodal in situ study combining structural information from synchrotron grazing-incidence wide-angle X-ray scattering (GIWAXS) and optical properties from photoluminescence (PL) spectroscopy with density functional theory calculations (DFT). The focus here is on the effects of the solvent and antisolvent during crystallization. The predominantly used solvents N,N-dimethylformamide (DMF), dimethyl sulfoxide (DMSO), and the antisolvent chlorobenzene (CB) as well as the solvent-antisolvent-precursor interactions are investigated. Given the high elemental complexity and mutual interdependencies between solvent, antisolvent, and perovskite precursors, we found significant differences in the crystallization pathways. DMF-pure precursors show the formation of the DMF-containing intermediate phase and the nucleation of compositionally distinct perovskite phases, while when DMSO is added, only crystalline α- and δ-phases were found. In addition, the presence of DMSO helps the formation of α-perovskite. Coordination energy and bond order (BO) calculations support our experimental findings. Dripping of CB induces nucleation at room temperature, slows the α-phase formation rate, and appears to reduce the nucleation radius. These findings provide novel insights into solvent, antisolvent, and perovskite precursor interactions and their formation pathways. The complexity of interactions between solvents and reagents highlights the importance of understanding these effects to further improve the reproducibility and optimize processing conditions.

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

Overcoming Metastable CO 2 Adsorption in a Bulky Diamine-Appended Metal-Organic Framework

We report carbon capture at fossil fuel-fired power plants is a critical strategy to mitigate anthropogenic contributions to global warming, but widespread deployment of this technology is hindered by a lack of energy-efficient materials that can be optimized for CO 2 capture from a specific flue gas. As a result of their tunable, step-shaped CO 2 adsorption profiles, diamine-functionalized metal-organic frameworks (MOFs) of the form diamine-Mg 2 (dobpdc) (dobpdc 4- = 4,4'-dioxidobiphenyl-3,3'-dicarboxylate) are among the most promising materials for carbon capture applications. Here, we present a detailed investigation of dmen-Mg 2 (dobpdc) (dmen = 1,2-diamino-2-methylpropane), one of only two MOFs with an adsorption step near the optimal pressure for CO 2 capture from coal flue gas. While prior characterization suggested that this material only adsorbs CO 2 to half capacity (0.5 CO 2 per diamine) at 1 bar, we show that the half-capacity state is actually a metastable intermediate. Under appropriate conditions, the MOF adsorbs CO 2 to full capacity, but conversion from the half-capacity structure happens on a very slow time scale, rendering it inaccessible in traditional adsorption measurements. Data from solid-state magic angle spinning nuclear magnetic resonance spectroscopy, coupled with van der Waals-corrected density functional theory, indicate that ammonium carbamate chains formed at half capacity and full capacity adopt opposing configurations, and the need to convert between these states likely dictates the sluggish post-half-capacity uptake. By use of the more symmetric parent framework Mg 2 (pc-dobpdc) (pc-dobpdc 4- = 3,3'-dioxidobiphenyl-4,4'-dicarboxylate), the metastable trap can be avoided and the full CO 2 capacity of dmen-Mg 2 (pc-dobpdc) accessed under conditions relevant for carbon capture from coal-fired power plants.

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

Multipole Expansion of Atomic Electron Density Fluctuation Interactions in the Density-Functional Tight-Binding Method

The accuracy of the density-functional tight-binding (DFTB) method in describing noncovalent interactions is limited due to its reliance on monopole-based spherical charge densities. In this study, we present a multipole-extended second-order DFTB (mDFTB2) method that takes into account atomic dipole and quadrupole interactions. Furthermore, we combine the multipole expansion with the monopole-based third-order contribution, resulting in the mDFTB3 method. To assess the accuracy of mDFTB2 and mDFTB3, we evaluate their performance in describing noncovalent interactions, proton transfer barriers, and dipole moments. Here, our benchmark results show promising improvements even when using the existing electronic parameters optimized for the original DFTB3 model. Both mDFTB2 and mDFTB3 outperform their monopole-based counterparts, DFTB2 and DFTB3, in terms of accuracy. While mDFTB2 and mDFTB3 perform comparably for neutral and positively charged systems, mDFTB3 exhibits superior performance over mDFTB2 when dealing with negatively charged systems and proton transfers. Overall, the incorporation of the multipole expansion significantly enhances the accuracy of the DFTB method in describing noncovalent interactions and proton transfers.

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

Study of Material Densification of In718 in the Higher Throughput Parameter Regime

Selective Laser Melting (SLM) is a powder bed fusion additive manufacturing process used increasingly in the aerospace industry to reduce the cost, weight, and fabrication time for complex propulsion components. Previous optimization studies for SLM using the Concept Laser M1 and M2 machines at NASA Marshall Space Flight Center have centered on machine default parameters. The objective of this project is to characterize how heat treatment affects density and porosity from a microscopic point of view. This is performs using higher throughput parameters (a previously unexplored region of the manufacturing operating envelope for this application) on material consolidation. Density blocks were analyzed to explore the relationship between build parameters (laser power, scan speed, and hatch spacing) and material consolidation (assessed in terms of density and porosity). The study also considers the impact of post-processing, specifically hot isostatic pressing and heat treatment, as well as deposition pattern on material consolidation in the higher energy parameter regime. Metallurgical evaluation of specimens will also be presented. This work will contribute to creating a knowledge base (understanding material behavior in all ranges of the AM equipment operating envelope) that is critical to transitioning AM from the custom low rate production sphere it currently occupies to the world of mass high rate production, where parts are fabricated at a rapid rate with confidence that they will meet or exceed all stringent functional requirements for spaceflight hardware. These studies will also provide important data on the sensitivity of material consolidation to process parameters that will inform the design and development of future flight articles using SLM.

Cordner, Samuel↗
DOE OSTI2023

Evaluation of Thermal Neutron Scattering Cross Section of Uranium Silicide with Ab Initio Lattice Dynamics

Uranium silicide (U 3 Si 2 ) is a candidate material for the high-density nuclear fuel in commercial light water reactors [1], [2]. Its higher uranium density, 11.3 g-U/cm3, compared to that of uranium dioxide (UO 2 ), 9.7 g-U/cm3, can improve the performance of a nuclear reactor while using low enriched uranium (LEU) and diversify the choice of cladding materials [1]–[3]. It also has a higher thermal conductivity than UO 2 , which can reduce the thermal stress on the material caused by a temperature gradient across the fuel pellet and provide a larger margin for some postulated accidents [1], [2], [4]. Furthermore, compared to U3Si, another high-density fuel candidate, it has better resistance to in-pile swelling due to less irradiation-induced rapid amorphization [1], [3]. Corresponding to its importance in nuclear engineering, many previous studies have reported the properties of U3Si2. Experiments showed that U 3 Si 2 is a paramagnetic (PM) metal, where a slight linear increase in magnetic susceptibility was measured with increasing temperature [5], [6]. In addition, thermodynamic quantities such as thermal expansion coefficient, heat capacity, and thermal conductivity were experimentally determined over a wide temperature range [1], [7], [8]. In several computational studies, ab initio atomistic simulations based on density functional theory (DFT) were performed to calculate various properties including elastic constants, electronic density of states (DOS), and phonon dispersion curves [9]–[12]. Nevertheless, thermal neutron scattering cross sections, which are critical to the prediction of the parameters in reactor physics that are ultimately related to reactor criticality, have not yet been evaluated for U3Si2. The scattering cross section can be calculated from the phonon DOS, or the energy spectrum of lattice vibrations, of the crystalline system [13], [14]. However, there is also no experimental data available for the phonon DOS of U 3 Si 2 . While some computational studies reported the phonon DOS and/or dispersion curves from ab initio simulations [9]–[12], the accuracy cannot be guaranteed because it is unclear whether the spin-polarization behavior of PM U 3 Si 2 was properly described. In the present study, the thermal neutron scattering cross section for U 3 Si 2 is evaluated for the first time by calculating the phonon DOS for U3Si2 from ab initio lattice dynamics (AILD) simulations based on DFT. First, U 3 Si 2 is modeled based on the experimental structure, and AILD simulations are performed on the modeled U3Si2 to optimize the structure. Next, AILD simulations are performed for supercells with atomic displacement to calculate Hellmann-Feynman forces. Based on the calculated forces, partial phonon DOSs for U and Si are obtained, and the thermal neutron scattering law (TSL) for U 3 Si 2 is finally evaluated. To verify the accuracy of the calculations in the present study, the calculation results are compared with experimental data on the structure and heat capacity of U3Si2 [1], [7], [8], [15].

Geometry Optimization↗
DOE OSTI2023

Active sampling for neural network potentials: Accelerated simulations of shear-induced deformation in Cu–Ni multilayers

Neural network potentials (NNPs) can greatly accelerate atomistic simulations relative to ab initio methods, allowing one to sample a broader range of structural outcomes and transformation pathways. In this work, we demonstrate an active sampling algorithm that trains an NNP that is able to produce microstructural evolutions with accuracy comparable to those obtained by density functional theory, exemplified during structure optimizations for a model Cu–Ni multilayer system. We then use the NNP, in conjunction with a perturbation scheme, to stochastically sample structural and energetic changes caused by shear-induced deformation, demonstrating the range of possible intermixing and vacancy migration pathways that can be obtained as a result of the speedups provided by the NNP. The code to implement our active learning strategy and NNP-driven stochastic shear simulations is openly available at https://github.com/pnnl/Active-Sampling-for-Atomistic-Potentials .

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

Linear magnetoresistance with a universal energy scale in the strong-coupling superconductor Mo8Ga41 without quantum criticality

The recent discovery of a nonsaturating linear magnetoresistance in several correlated electron systems near a quantum critical point has revealed an interesting interplay between the linear magnetoresistance and the zero-field linear-in-temperature resistivity. These studies suggest a possible role of quantum criticality on the observed linear magnetoresistance. Here we report our discovery of a nonsaturating, linear magnetoresistance in Mo 8 Ga 41 , a nearly isotropic strong electron-phonon coupling superconductor with a linear-in-temperature resistivity from the transition temperature to ~ 55 K. The growth of the resistivity in field is comparable to that in temperature, provided that both quantities are measured in the energy unit. Overall, our data sets are remarkably similar to magnetoresistance data of the optimally doped La 2 - x Sr x CuO 4 , despite the clearly different crystal and electronic structures, and the apparent absence of quantum critical physics in Mo 8 Ga 41 . A new empirical scaling formula is developed, which is able to capture the key features of the low-temperature magnetoresistance data of Mo 8 Ga 41 , as well as the data of La 2 - x Sr x CuO 4 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2014

Novel Catalysts and Processing Technologies for Production of Aerospace Fuels from Non-Petroleum Raw Materials

Transportation fuels production (including aerospace propellants) from non-traditional sources (gases, waste materials, and biomass) has been an active area of research and development for decades. Reducing terrestrial waste streams simultaneous with energy conversion, plentiful biomass, new low-cost methane sources, and/or extra-terrestrial resource harvesting and utilization present significant technological and business opportunities being realized by a new generation of visionary entrepreneurs. We examine several new approaches to catalyst fabrication and new processing technologies to enable utilization of these non-traditional raw materials. Two basic processing architectures are considered: a single-stage pyrolysis approach that seeks to basically re-cycle hydrocarbons with minimal net chemistry or a two-step paradigm that involves production of supply or synthesis gas (mainly carbon oxides and hydrogen) followed by production of fuel(s) via Sabatier or methanation reactions and/or Fischer-Tropsch synthesis. Optimizing the fraction of product stream relevant to targeted aerospace (and other transportation) fuels via modeling, catalyst fabrication and novel reactor design are described. Energy utilization is a concern for production of fuels for either terrestrial or space operations; renewable sources based on solar energy and/or energy efficient processes may be mission enabling. Another important issue is minimizing impurities in the product stream(s), especially those potentially posing risks to personnel or operations through (catalyst) poisoning or (equipment) damage. Technologies being developed to remove (and/or recycle) heteroatom impurities are briefly discussed as well as the development of chemically robust catalysts whose activity are not diminished during operation. The potential impacts on future missions by such new approaches as well as balance of system issues are addressed.

propellants↗
DOE OSTI2024

Insights into the Electrochemical Oxidation and Reduction of Nickel Oxide Surfaces

Surface oxidation/reduction processes, driven by varying electrochemical potentials, can substantially impact catalyst effectiveness and, consequently, electrolyzer performance. Here, this study combines theoretical and experimental approaches to explore the surface redox behavior of nickel oxides, which are cost-effective and efficient catalysts for many electrochemical reactions. Surface Pourbaix diagrams for three different phases of nickel oxides, i.e., nickel hydroxide (Ni(OH) 2 ), nickel oxyhydroxide (NiOOH), and nickel dioxide (NiO 2 ), were constructed using density functional theory-based simulations. Various experimental methods, including cyclic voltammetry, in situ Raman spectroscopy, and electrochemical titration, were employed to probe the surface redox processes of nickel oxide thin films. Our findings indicate that the ABAB stacking sequence of Ni(OH) 2 lacks stability under oxidizing conditions to host the surface oxidation (deprotonation) events, while the AABBCC stacking sequence of NiOOH is energetically favorable due to the presence of interlayer hydrogen bonding. Rapid charge transfer facilitated by interlayer hydrogen bonding accounts for the higher reactivity of partially oxidized/reduced NiOOH (001) surfaces compared to Ni(OH) 2 (001) and NiO 2 (001) surfaces with the same stoichiometry, where interlayer hydrogen bonding is absent. Insights presented in this work can offer guidelines for optimizing operational conditions and tailoring the surface structures and oxidation states of nickel oxides to enhance performance in applications such as electrocatalysis and supercapacitors.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Polarons and electrical leakage in BaZrO 3 and BaCeO 3

The proton conductors barium zirconate (BaZrO 3 , or BZO), barium cerate (BaCeO 3 , or BCO), and their alloys have attracted considerable interest as solid-state electrolytes in solid-oxide fuel and electrolysis cells. However, the reasons for their non-negligible electrical conductivity, which can limit their performance, are not fully understood. To address that question, we use first-principles calculations based on density functional theory to study the properties of hole and electron polarons in BZO and BCO. We confirm that hole polarons form in both materials, slightly more favorably in BCO. Electron polarons, on the other hand, are stable only in BCO and in Ce-containing alloys; in pure BZO, electron polaron states are unfavorable relative to free carriers, though they become accessible when Ce impurities are introduced. In general, doped BZO and BCO will have negligible electron concentrations, but larger concentrations of holes and hole polarons may be present. Avoiding extreme O-rich conditions and limiting dopant concentrations are key strategies for reducing significant p -type electrical leakage. Furthermore, our results provide physical insights into the different electronic behaviors of BZO and BCO, which can be used to optimize their performance as pure ionic conductors.

08 HYDROGEN↗