Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “bulk density”

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 1,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Experimental and theoretical studies on the elasticity of tungsten to 13 GPa

Compressional (V_P) and shear wave (V_S) velocities of polycrystalline tungsten have been measured up to ~13 GPa at room temperature using ultrasonic interferometry in a multi-anvil apparatus. Using finite strain equation of state approaches, the elastic bulk and shear moduli and their pressure dependences are derived yielding $K_{S0}$= 325.9 ± 4.8 GPa, $G_0$ =164.1±2.5 GPa, $K^{'}_{S0}$ =3.65 ± 0.05 and $G^{'}_0$ =1.28 ± 0.02. On the basis of the current experimental data, the high-pressure behavior of Young’s modulus, Poisson’s ratio and ductility/brittleness for tungsten are also investigated. Complementary to the experimental data, the single crystal elastic constants, as well as the elastic anisotropy of tungsten are computed using density functional theory (DFT). The Voigt-Reuss-Hill average of the bulk and shear moduli calculated using the single crystal elastic constants from DFT are found comparable to the current experimental results within about 5%. Here, the present study offers a dataset for the elasticity of polycrystalline bcc tungsten to a maximum pressure more than 25-fold higher than other previous ultrasonic studies, which can further our understanding about the elastic, mechanical and electronic properties of tungsten under extreme conditions as well as thermodynamic modelling of its alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Structural, electronic, mechanical, and dynamical properties of graphene oxides: A first principles study

We report the results of a theoretical study on the structural, electronic, mechanical, and vibrational properties of some graphene oxide models (GDO, a-GMO, z-GMO, ep-GMO and mix-GMO) at ambient pressure. The calculations are based on the ab-initio plane-wave pseudo potential density functional theory, within the generalized gradient approximations for the exchange and correlation functional. The calculated values of lattice parameters, bulk modulus, and its first order pressure derivative are in good agreement with other reports. A linear response approach to the density functional theory is used to derive the phonon frequencies. We discuss the contribution of the phonons in the dynamical stability of graphene oxides and detailed analysis of zone centre phonon modes in all the above mentioned models. Our study demonstrates a wide range of energy gap available in the considered models of graphene oxide and hence the possibility of their use in nanodevices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
NASA NTRS1993

Porosity and the ecology of icy satellites

The case for a significant role for porosity in the structure and evolution of icy bodies in the Solar System has been difficult to establish. We present a relevant new data set and a series of structure models including a mechanical compression, not thermal creep, model for porosity that accounts satisfactorily for observed densities, moments of inertia, geologic activity, and sizes of tectonic features on icy satellites. Several types of observational data sets have been used to infer significant porosity, but until recently, alternative explanations have been preferred. Our first area of concern is the occurrence of cryovolcanism as a function of satellite radius; simple radiogenic heating models of icy satellites suggest minimum radii for melting and surface cryovolcanism to be 400 to 500 km, yet inferred melt deposits are seen on satellites half that size. One possible explanation is a deep, low conductivity regolith which lowers conductivity and raises internal temperatures, but other possibilities include tidal heating or crustal compositions of low conductivity. Our second area of concern is the occurrence and magnitude of tectonic strain; tectonic structures have been seen on icy satellites as small as Mimas and Proteus. The structures are almost exclusively extensional, with only a few possible compression Al features, and inferred global strains are on the order of 1 percent expansion. Expansions of this order in small bodies like Mimas and prevention of late compressional tectonics due to formation of ice mantles in larger bodies like Rhea are attained only in structure models including low-conductivity, and thus possibly high porosity, crusts. Thirdly, inferred moments of inertia less than 0.4 in Mimas and Tethys can be explained by high-porosity crusts, but also by differentiation of a high density core. Finally, the relatively low densities of smaller satellites like Mimas and Miranda relative to larger neighbors can be explained by deep porosity, but also by bulk compositional differences. Recent work has strengthened the case for significant porosity. Halley's nucleus was found to have a density near 0.6 g/cu cm, Janus and Epimethus were proposed to have densities near 0.7 g/cu cm, densities almost certainly due to high porosity. The irregular-spherical shape transition of icy satellites was quantitatively explained by low conductivity regoliths. A creative structure/thermal history model for Mimas simultaneously accounts quantitatively for Mimas' low density and moment of inertia by invoking initial high-porosity and subsequent compaction in the deep interior by thermal creep. The main problem with this promising model is that approximately 7 percent predicts a reduction in Mimas' radius, implying significant compressional failure and prevention of extensional tectonics, in contradiction to the observed extensional features and inferred 1 percent expansion in radius.

Croft, Steven K.↗
DOE OSTI2015

Density and x-ray emission profile relationships in highly ionized high-Z laser-produced plasmas

We present a benchmark measurement of the electron density profile in the region where the electron density is 10{sup 19 }cm{sup –3} and where the bulk of extreme ultraviolet (EUV) emission occurs from isotropically expanding spherical high-Z gadolinium plasmas. It was found that, due to opacity effects, the observed EUV emission is mostly produced from an underdense region. We have analyzed time-resolved emission spectra with the aid of atomic structure calculations and find the multiple ion charge states around 18+ during the laser pulse irradiation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Geodyn Material Library: Pseudocap models for dry porous tocks

This report describes the second edition of the Pseudocap Strength models for porous rocks implemented in GEODYN material library. The first model was developed in 2007 and calibrated for concrete. Then, the model parameters were calibrated based on triaxial tests reported for limestones and sandstones of various porosities. In these models some key parameters were chosen as functions of the reference porosity,Φ. Since then multiple modifications were implemented in the model, therefore, it has been recalibrated for some common porous materials such as limestones, sandstones, alluvium, tuffs and granite. Two types of models are described in this repot. The first type (called Pseudocap Model or PM) is for rocks from a specific location. Parameters were calibrated for several specific geologic materials. The second type (called Generic Pseudocap Model or GPM) is useful for the sites where only basic information (rock type, porosity) is available. Generic models include built-in correlations between porosities and other mechanical properties observed for certain rock types. The models of both types were validated by comparing not only to quasi-static triaxial tests for these materials but also to shock Hugoniot data and spherical explosion data for some materials. All models were derived in the frame of isotropic plasticity. They are designed to be used in explicit finite element/difference codes. Tangent stiffness tensor is not provided but can be calculated numerically for the model to be used in implicit finite element codes. For an isotropic material the stress can be decomposed into volumetric and deviatoric parts. The volumetric part is modeled using an Equation of state (EOS) which calculates the pressure and the bulk sound speed as functions of the internal specific energy and density. Here a simple, Mie-Gruneisen EOS is presented, but tabulated EOS (LEOS) provided by the library can be used as well. The stress is limited by the yield surface which depends on three invariants of the stress tensor and specific internal energy. In addition, to capture the strain-rate dependence a simple multiplier is used for the yield surface which depends on the equivalent plastic strain rate. The failure surface (the ultimate yield, Y f , defined later) is chosen in the Hoek-Brown form, commonly used in rock mechanics. It includes measurable parameters such as Unconfined Compressive Strength (UCS) as well as scale parameters characterizing the quality of the rock such as GSI (Geologic Strength Index). Thus, even though the model is calibrated for small samples it offers a way to extrapolate the strength to the field scale using geological characterization of the rock mass. The model captures effects of brittle-ductile transition in rocks by introducing a cap multiplier to the yield function. The rate of dilatancy (bulking) is proportional to the slope of the yield surface affected by the cap. Therefore, it takes place only at low confinements when the pressure is less than the brittle-ductile transition pressure, P BD . On the contrary, the porous compaction takes place at pressures higher than P BD . The cap moves as the porosity is compacted or new porosity is generated due to dilatancy. The porous compaction is modeled using an evolution equation which includes deviatoric stress so that the onset of compaction corresponds to the cap surface. The model captures effects of shear-enhanced compaction which is an important for porous rocks. Section 2 describes the modeling framework and Section 3 presents the model calibration procedure. Section 4 compares experimental data for various rocks versus model predictions. The model parameters used for this comparison are given in Appendix. The files with material constants are available with the latest GEODYN material library distribution.

58 GEOSCIENCES↗
DOE OSTI2017

Theoretical Insights to Bulk Activity Towards Oxygen Evolution in Oxyhydroxides

The nature of the electrochemical water splitting activity of layered pure and Fe-doped NiOOH is investigated using density functional theory calculations. We find similar thermodynamics for the oxygen evolution reaction (OER) intermediates between the layers of oxyhydroxides, that is, in the bulk of the materials as on the (001) surface. The effect of interlayer spacing on adsorption energy is affected by both the crystal structure and the level of hydrogenation of the active sites. For the Fe-doped NiOOH, we observe general weakening of binding between the different OER intermediates and the catalyst material. The calculated OER activity depends both on doping and interlayer spacing, and our results are generally congruent with available experimental data. In conclusion, these results suggest that such interlayer “bulk” sites may contribute to measured OER activity for both the pure and Fe-doped NiOOH catalysts.

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

THE COS-HALOS SURVEY: RATIONALE, DESIGN, AND A CENSUS OF CIRCUMGALACTIC NEUTRAL HYDROGEN

We present the design and methods of the COS-Halos survey, a systematic investigation of the gaseous halos of 44 z = 0.15-0.35 galaxies using background QSOs observed with the Cosmic Origins Spectrograph aboard the Hubble Space Telescope. This survey has yielded 39 spectra of z{sub em} ≅ 0.5 QSOs with S/N ∼10-15 per resolution element. The QSO sightlines pass within 150 physical kpc of the galaxies, which span early and late types over stellar mass log M{sub *}/M{sub ☉} = 9.5-11.5. We find that the circumgalactic medium exhibits strong H I, averaging ≅ 1 Å in Lyα equivalent width out to 150 kpc, with 100% covering fraction for star-forming galaxies and 75% covering for passive galaxies. We find good agreement in column densities between this survey and previous studies over similar range of impact parameter. There is weak evidence for a difference between early- and late-type galaxies in the strength and distribution of H I. Kinematics indicate that the detected material is bound to the host galaxy, such that ∼> 90% of the detected column density is confined within ±200 km s{sup –1} of the galaxies. This material generally exists well below the halo virial temperatures at T ∼< 10{sup 5} K. We evaluate a number of possible origin scenarios for the detected material, and in the end favor a simple model in which the bulk of the detected H I arises in a bound, cool, low-density photoionized diffuse medium that is generic to all L* galaxies and may harbor a total gaseous mass comparable to galactic stellar masses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Polder maps: Improving OMIT maps by excluding bulk solvent

The crystallographic maps that are routinely used during the structure-solution workflow are almost always model-biased because model information is used for their calculation. As these maps are also used to validate the atomic models that result from model building and refinement, this constitutes an immediate problem: anything added to the model will manifest itself in the map and thus hinder the validation. OMIT maps are a common tool to verify the presence of atoms in the model. The simplest way to compute an OMIT map is to exclude the atoms in question from the structure, update the corresponding structure factors and compute a residual map. It is then expected that if these atoms are present in the crystal structure, the electron density for the omitted atoms will be seen as positive features in this map. This, however, is complicated by the flat bulk-solvent model which is almost universally used in modern crystallographic refinement programs. This model postulates constant electron density at any voxel of the unit-cell volume that is not occupied by the atomic model. Consequently, if the density arising from the omitted atoms is weak then the bulk-solvent model may obscure it further. A possible solution to this problem is to prevent bulk solvent from entering the selected OMIT regions, which may improve the interpretative power of residual maps. This approach is called a polder (OMIT) map. Polder OMIT maps can be particularly useful for displaying weak densities of ligands, solvent molecules, side chains, alternative conformations and residues both in terminal regions and in loops. As a result, the tools described in this manuscript have been implemented and are available in PHENIX .

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

Bulk vs Intrinsic Activity of NiFeO x Electrocatalysts in the Oxygen Evolution Reaction: The Influence of Catalyst Loading, Morphology, and Support Material

We used a combination of ultrahigh vacuum surface science techniques, X-ray spectroscopy, electrochemistry, and density functional theory (DFT), to characterize the influence of catalyst morphology, loading/coverage, and substrate material on the bulk (all atoms) and intrinsic (electrochemically accessible atoms) activity of NiFeO x electrocatalysts in the oxygen evolution reaction (OER). NiFeO x catalysts were grown on both Au(111) and highly oriented pyrolytic graphite (HOPG) electrodes. DFT predicted Fe edge-site atoms at the NiFeO x /Au(111) interface to be the most thermodynamically favorable reaction center, and X-ray absorption spectroscopy data indicated small NiFeO x catalyst particles on Au(111) contained a high population of OER active Fe edge-site atoms. However, restructuring of the Au(111) surface due to repeated oxidation and reduction cycles of the OER CV measurements encapsulated small NiFeO x nanoparticles at catalyst loadings below ~1.5 nmol metal /cm 2 , passivated catalyst edges and reduced bulk OER activity of Au-supported NiFeO x compared with HOPG-supported ones. Analysis of intrinsic activity revealed that the Au(111) support strongly benefited electrochemically accessible NiFeO x atoms, and we observed a 2–3 fold activity enhancement compared with HOPG-supported catalysts for loadings above ~1 nmol metal /cm 2 . Overall, evaluating bulk vs intrinsic activity and identifying loading/coverage-dependent support effects is important for accurately probing fundamental interfacial chemistry, choosing suitable catalyst loadings and supports, and optimizing system parameters to maximize the performance of electrocatalyst systems.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Effects of temperature and gas–liquid mass transfer on the operation of small electrochemical cells for the quantitative evaluation of CO 2 reduction electrocatalysts

In the last few years, there has been increased interest in electrochemical CO 2 reduction (CO2R). Many experimental studies employ a membrane separated, electrochemical cell with a mini H-cell geometry to characterize CO2R catalysts in aqueous solution. This type of electrochemical cell is a mini-chemical reactor and it is important to monitor the reaction conditions within the reactor to ensure that they are constant throughout the study. Here we show that operating cells with high catalyst surface area to electrolyte volume ratios (S/V) at high current densities can have subtle consequences due to the complexity of the physical phenomena taking place on electrode surfaces during CO2R, particularly as they relate to the cell temperature and bulk electrolyte CO 2 concentration. Both effects were evaluated quantitatively in high S/V cells using Cu electrodes and a bicarbonate buffer electrolyte. Electrolyte temperature is a function of the current/total voltage passed through the cell and the cell geometry. Even at a very high current density, 20 mA cm -2 , the temperature increase was less than 4 °C and a decrease of < 10% in the dissolved CO 2 concentration is predicted. In contrast, limits on the CO 2 gas-liquid mass transfer into the cells produce much larger effects. By using the pH in the cell to measure the CO 2 concentration, significant undersaturation of CO 2 is observed in the bulk electrolyte, even at more modest current densities of 10 mA cm -2 . Undersaturation of CO 2 produces large changes in the faradaic efficiency observed on Cu electrodes, with H 2 production becoming increasingly favored. Finally, we show that the size of the CO 2 bubbles being introduced into the cell is critical for maintaining the equilibrium CO 2 concentration in the electrolyte, and we have designed a high S/V cell that is able to maintain the near-equilibrium CO 2 concentration at current densities up to 15 mA cm -2 .

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

First-principles study of structural, elastic, electronic, vibrational and thermodynamic properties of uranium aluminides

The structural, elastic, electronic, vibrational and thermodynamic properties of uranium aluminides UAl x (x = 2, 3, 4) were studied by density functional theory calculations. The single crystal elastic constants of UAl x were predicted using the stress-strain method, which were further used to calculate the polycrystalline aggregate properties of UAl x , including bulk modulus, shear modulus, Young’s modulus, and Poisson’s ratio, etc. The calculated electronic density of states confirm that UAl x compounds are metallic phases with majority states at the Fermi level contributed by U 5 f electrons. The phonon dispersion relations and density of states show that the low frequency acoustic phonon modes of UAl x are dominated by the lattice vibration of uranium atoms while high frequency optical phonon modes are from the vibration of aluminum atoms. Using quasiharmonic approximation, thermodynamic properties of UAl x , including Gibbs free energy, entropy, heat capacity, and linear thermal expansion coefficient, were predicted by including both lattice vibrational and thermal electronic contributions. The thermal electronic energy was found to be crucial for the description of the temperature dependence of the thermodynamic properties. In conclusion, the derived Gibbs energy functions of UAl x are expected to be useful to the thermodynamic modeling of the ternary U-Mo-Al system.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Energetics of hydrogen and helium-vacancy complexes in bulk and near surfaces of tungsten: First-principles study

Understanding the interaction between hydrogen (H) and helium-vacancy (He-V) complexes in tungsten (W) is important for the development of plasma-facing materials in fusion reactors. H trapping by He x V y complexes in bulk W, as well as the H solution behavior and H trapping by He x V complexes near W(100), W(111), and W(110) surfaces, has been investigated by first-principles computer simulations using density function theory. The results show that the sequential H binding energies to He x V complexes in bulk W decrease with the increasing number of H and He. For the He x V 2 complexes in bulk W, H prefers to trap at interstitial sites near the junction of the di-vacancy, where the H can minimize the isosurface of optimal charge density. The most stable interstitial sites for H below W surfaces are dependent on the surface orientation. Our calculations indicate that H atoms tend to prefer a depth of 0.3 nm below the W(100) and W(111) surfaces due to the surface reconstruction. Here, the binding energy of H to a HeV complex near W surfaces has the most significant orientation dependence below the W(111) surface, followed by the W(100) and W(110) surfaces. Compared with the bulk value, the largest difference in the average binding energy of H to the stable He x V complexes at the three W surfaces is about 0.2 eV. Furthermore, the effect of surfaces on the H binding energy to He x V complexes can be ignored for depths greater than 0.65 nm.

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

Density functional calculations for structural, electronic, and magnetic properties of gadolinium-oxide clusters

Gadolinium-oxide clusters in various sizes and stoichiometries have been systematically studied by employing the density functional theory with the generalized gradient approximation. The clusters in bulk stoichiometry are relatively more stable and their binding energies increase with the increasing size. Stoichiometric (Gd{sub 2}O{sub 3}){sub n} clusters of n = 1–3 prefer cage-like structures, whereas the clusters of n = 4–30 prefer compact structures layered by wedge-like units and exhibit a rough feature toward the bulk-like arrangement with small disorders of atomic positions. The polyhedral-cages analogous to carbon-fullerenes are stable isomers yet not the minimum energy configurations. Their stabilities can be improved by embedding one oxygen atom or a suitable cage to form core-shell configurations. The mostly favored antiferromagnetic couplings between adjacent Gd atoms are nearly degenerated in energy with their ferromagnetic couplings, resulting in super-paramagnetic characters of gadolinium-oxide clusters. The Ruderman-Kittel-Kasuya-Yosida (RKKY)-type mechanism together with the superexchange-type mechanism plays cooperation role for the magnetic interactions in clusters. We present, as a function of n, calculated binding energies, ionization potential, electron affinity, and electronic dipole moment.

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

Describing the Diverse Geometries of Gold from Nanoclusters to Bulk—A First-Principles-Based Hybrid Bond-Order Potential

Molecular dynamics simulations using empirical force fields (EFFs) are crucial for gaining fundamental insights into atomic structure and long time scale dynamics of Au nanoclusters with far-reaching applications in energy and devices. This approach is thwarted by the failure of currently available EFFs in describing the size-dependent dimensionality and diverse geometries exhibited by Au clusters (e.g., planar structures, hollow cages, tubes, pyramids, space-filled structures). Here, we mitigate this issue by introducing a new hybrid bond-order potential (HyBOP), which accounts for (a) short-range interactions via Tersoff-type BOP terms that accurately treat bond directionality and (b) long-range dispersion effects by a scaled Lennard–Jones term whose contribution depends on the local atomic density. We optimized the independent parameters for our HyBOP using a global optimization scheme driven by genetic algorithms. Moreover, to ensure good transferability of these parameters across different length scales, we used an extensive training data set that encompasses structural and energetic properties of 1000 13-atom Au clusters, surface energies, as well as bulk polymorphs, obtained from density functional theory (DFT) calculations. Our newly developed HyBOP has been found to accurately describe (a) global minimum energy configurations at different cluster sizes as well as order of stability of various cluster configurations at any size, (b) critical size of transition from planar to globular clusters, (c) evolution of structural motifs with cluster size, and (c) thermodynamics, structure, elastic properties, and energetic ordering of bulk condensed phases as well as surfaces, in excellent agreement with DFT calculations and spectroscopic experiments. This makes our newly developed HyBOP a valuable, computationally robust but inexpensive tool for investigating a wide range of materials phenomena occurring in Au at the atomistic level.

36 MATERIALS SCIENCE↗
DOE OSTI2011

Elasticity of Diopside to 8 GPa and 1073K and Implications for the Upper Mantle

Simultaneous measurements of elastic wave velocities and density have been conducted on diopside by a combined ultrasonic interferometry and X-ray diffraction methods at pressure and temperature conditions relevant to the Earth's upper mantle. The current study not only yielded the first direct measurement of the shear modulus at high pressures, but also enabled a simultaneous determination of the bulk and shear properties and their pressure and temperature derivatives from the measured density and velocities. Finite strain analysis of the experimental data results in K{sub S0} = 116.4(7) GPa, K'S0 = 4.9(1), ({partial_derivative}K{sub S}/{partial_derivative}T){sub P} = -0.012(1) GPa GPa, G'{sub 0} = 1.6(1) and ({partial_derivative}G/{partial_derivative}T){sub P} = -0.011(1) GPa/K. With these results and the same finite strain equations, the P and S wave velocities of diopsidic mantle clinopyroxene were calculated along a 1600 K adiabatic geotherm. In comparison with other mantle minerals, the seismic velocities of diopsidic clinopyroxene at upper mantle depths are 1-3% higher than those of orthopyroxene, 1-2% and 6-8% lower than those of olivine and majoritic garnet, respectively.

58 GEOSCIENCES↗
DOE OSTI2020

Utilizing ink composition to tune bulk-electrode gas transport, performance, and operational robustness for a Fe-N-C catalyst in polymer electrolyte fuel cell

With lower site density and turnover frequency, platinum group metal (PGM)-free catalysts based electrodes are often greater than 50 mu m thick in order to increase performance across the fuel cell operating range. Consequently, PGM-free electrodes have an additional bulk electrode transport resistance beyond the local or aggregate level transport in thin platinum-based electrodes. In parallel to the development of more active and durable PGM-free catalysts, advancements in understanding the interplay between PGM-free electrode fabrication, bulk-electrode transport, proton conductivity and performance are needed. Here, the relationship between ionic and gas phase transport through the electrode thickness is modified by adjusting electrocatalyst and ionomer flocculation/interaction at the ink level. The influence of the ink composition (water/n-propanol content) is examined via various in-situ electrochemical and ex-situ characterization techniques and the resulting electrode structure/performance relationship contrasted with electrode performance robustness across a range of relative humidity (RH). For the electrocatalyst examined here, a water-rich (82 wt% H2O) ink formulation was favorable for operation at high RH due to improved molecular diffusion through larger electrode pores. In contrast, the improved interactions between ionomer and electrocatalyst enabled a more robust electrode and higher performance during low RH operation for the 50 wt% H2O content ink.

08 HYDROGEN↗
DOE OSTI2018

New Bonding Model of Radical Adsorbate on Lattice Oxygen of Perovskites

In this paper, a new model of bonding between radical adsorbates and lattice oxygens is proposed that considers both the adsorbate–oxygen bonding and the weakening of the metal–oxygen bonds. Density functional calculations of SrMO 3 perovskites for M being 3d, 4d, and 5d transition metals are used to correlate the bulk electronic structure with the surface-oxygen reactivity. Occupation of the metal–oxygen antibonding states, examined via the crystal orbital Hamilton population (COHP), is found to be a useful bulk descriptor that correlates with the vacancy formation energy of the lattice oxygen and its hydrogen adsorption energy. Analysis of density-of-states and COHP indicates that H adsorption energy is a combined result of formation of the O–H bond and the weakening of the surface metal–oxygen bond due to occupation of the metal–oxygen antibonding states by the electron from H. Lastly, this insight will be useful in understanding the trends in surface reactivity of perovskites and transition-metal oxides in general.

36 MATERIALS SCIENCE↗