Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Vacancy Formation Energy”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The Defects Genome of Janus Transition Metal Dichalcogenides

2D Janus Transition Metal Dichalcogenides (TMDs) have attracted much interest due to their exciting quantum properties arising from their unique two-faced structure, broken-mirror symmetry, and consequent colossal polarization field within the monolayer. While efforts are made to achieve high-quality Janus monolayers, the existing methods rely on highly energetic processes that introduce unwanted grain-boundary and point defects with still unexplored effects on the material's structural and excitonic properties Through high-resolution scanning transmission electron microscopy (HRSTEM), density functional theory (DFT), and optical spectroscopy measurements; this work introduces the most encountered and energetically stable point defects. It establishes their impact on the material's optical properties. HRSTEM studies show that the most energetically stable point defects are single (V S and V Se ) and double chalcogen vacancy (V S −V Se ), interstitial defects (Mi), and metal impurities (MW) and establish their structural characteristics. DFT further establishes their formation energies and related localized bands within the forbidden band. Cryogenic excitonic studies on h-BN-encapsulated Janus monolayers offer a clear correlation between these structural defects and observed emission features, which closely align with the results of the theory. Finally, the overall results introduce the defect genome of Janus TMDs as an essential guideline for assessing their structural quality and device properties.

2D materials↗
DOE OSTI2024

Atomically precise engineering of spin–orbit polarons in a kagome magnetic Weyl semimetal

Atomically precise defect engineering is essential to manipulate the properties of emerging topological quantum materials for practical quantum applications. However, this remains challenging due to the obstacles in modifying the typically complex crystal lattice with atomic precision. Here, we report the atomically precise engineering of the vacancy-localized spin–orbit polarons in a kagome magnetic Weyl semimetal Co 3 Sn 2 S 2 , using scanning tunneling microscope. We achieve the step-by-step repair of the selected vacancies, leading to the formation of artificial sulfur vacancies with elaborate geometry. We find that that the bound states localized around these vacancies undergo a symmetry dependent energy shift towards Fermi level with increasing vacancy size. As the vacancy size increases, the localized magnetic moments of spin–orbit polarons become tunable and eventually become itinerantly negative due to spin–orbit coupling in the kagome flat band. These findings provide a platform for engineering atomic quantum states in topological quantum materials at the atomic scale.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Porosity in Sr 1-x Ca x FeO 3-δ oxygen carriers: The role of surface area and pretreatment on storage activity

Perovskite oxides have generated interest as robust, low-temperature oxygen carrier materials for a variety of clean energy applications, including chemical looping gasification and air separations. Methods to improve O 2 desorption kinetics are vital to allow these carriers to compete economically with traditional metal oxide carriers or cryogenic separations. In this report we investigated the cumulative roles that surface area, pretreatment, and elemental composition have on the oxygen storage properties of a state-of-the-art carrier system, Sr 1-x Ca x FeO 3-δ (x = 0.20, 0.25, 0.30) synthesized using multiple methods. Porous materials synthesized by the Pechini, or citrate, method had their surface area controlled using the synthesis temperature. The high surface area of the Sr 0.7 Ca 0.3 FeO 3-δ materials is most beneficial at low operating temperatures, such as 350 and 400 °C, as their reduction rates are twice as fast as those obtained with their bulk counterparts. These effects are observed at higher operating temperatures and within a single composition, but temperature tunability using variable Ca 2+ substitution in the Sr 1-x Ca x FeO 3-δ overshadows the improvements gained from higher surface areas. Additionally, we establish that pretreatment in N2 at an elevated temperature is necessary to enhance kinetics further. For maximum efficiency, pretreatment at the synthesis temperature is suggested for the Pechini method-synthesized systems, whereas 800 °C is adequate for bulk materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Experimental and theoretical studies of native deep-level defects in transition metal dichalcogenides

Abstract Transition metal dichalcogenides (TMDs), especially in two-dimensional (2D) form, exhibit many properties desirable for device applications. However, device performance can be hindered by the presence of defects. Here, we combine state of the art experimental and computational approaches to determine formation energies and charge transition levels of defects in bulk and 2D MX 2 (M = Mo or W; X = S, Se, or Te). We perform deep level transient spectroscopy (DLTS) measurements of bulk TMDs. Simultaneously, we calculate formation energies and defect levels of all native point defects, which enable identification of levels observed in DLTS and extend our calculations to vacancies in 2D TMDs, for which DLTS is challenging. We find that reduction of dimensionality of TMDs to 2D has a significant impact on defect properties. This finding may explain differences in optical properties of 2D TMDs synthesized with different methods and lays foundation for future developments of more efficient TMD-based devices.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Determining the Adsorption Energetics of 2,3-Butanediol on RuO 2 (110): Coupling First-Principles Calculations With Global Optimizers

As climate change continues to pose a threat to the Earth due to the disrupted carbon cycles and fossil fuel resources remain finite, new sources of sustainable hydrocarbons must be explored. 2,3-butanediol is a potential source to produce butene because of its sustainability as a biomass-derived sugar. Butene is an attractive product because it can be used as a precursor to jet fuel, categorizing this work in the alcohol-to-jet pathway. While studies have explored the conversion of 2,3-butanediol to butene, little is understood about the fundamental reaction itself. We quantify the energetics for three pathways that were reported in the literature in the absence of a catalyst. One of these pathways forms a 1,3-butadiene intermediate, which is a highly exothermic process and thus is unlikely to occur since 2,3-butanediol likely gets thermodynamically trapped at this intermediate. We further determined the corresponding energetics of 2,3-butanediol adsorption on an ensemble of predetermined binding sites when it interacts with a defect-free stoichiometric RuO 2 (110) surface. Within this ensemble of adsorption sites, the most favorable site has 2,3-butanediol covering a Ru 5–coordinated cation. This approach is compared to that obtained using the global optimization algorithm as implemented in the Northwest Potential Energy Surface Search Engine. When using such a global optimization algorithm, we determined a more favorable ground-state structure that was missed during the manual adsorption site testing, with an adsorption energy of -2.61 eV as compared to -2.34 eV when using the ensemble-based approach. We hypothesize that the dehydration reaction requires a stronger chemical bond, which could necessitate the formation of oxygen vacancies. As such, this study has taken the first step toward the utilization of a global optimization algorithm for the rational design of Ru-based catalysts toward the formation of butene from sustainable resources.

09 BIOMASS FUELS↗
DOE OSTI2022

Valence Disproportionation of GeS in the PbS Matrix Forms Pb 5 Ge 5 S 12 Inclusions with Conduction Band Alignment Leading to High n-Type Thermoelectric Performance

Converting waste heat into useful electricity using solid-state thermoelectrics has a potential for enormous global energy savings. Lead chalcogenides are among the most prominent thermoelectric materials, whose performance decreases with an increase in chalcogen amounts (e.g., PbTe>PbSe>PbS). Herein, we demonstrate the simultaneous optimization of the electrical and thermal transport properties of PbS-based compounds by alloying with GeS. The addition of GeS triggers a complex cascade of beneficial events as follows: Ge 2+ substitution in Pb 2+ and discordant off-center behavior; formation of Pb 5 Ge 5 S 12 as stable second phase inclusions through valence disproportionation of Ge 2+ to Ge0 and Ge 4+ . PbS and Pb 5 Ge 5 S 12 exhibit good conduction band energy alignment that preserves the high electron mobility; the formation of Pb 5 Ge 5 S 12 increases the electron carrier concentration by introducing S vacancies. Sb doping as the electron donor produces a large power factor and low lattice thermal conductivity (κ lat ) of ~0.61 Wm -1 K -1 . The highest performance was obtained for the 14% GeS-alloyed samples, which exhibited an increased room temperature electron mobility of ~121 cm 2 V -1 s -1 for 3 × 10 19 cm -3 carrier density, and a ZT, of 1.32 at 923 K. This is ~ 55% greater that the corresponding Sb-doped PbS sample and is one of the highest reported for the n-type PbS system. Moreover, the average ZT (ZT avg ) of ~0.76 from 400 to 923 K is the highest for PbS-based systems.

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

Defect graph neural networks for materials discovery in high-temperature clean-energy applications

We present a graph neural network approach that fully automates the prediction of defect formation enthalpies for any crystallographic site from the ideal crystal structure, without the need to create defected atomic structure models as input. Here we used density functional theory reference data for vacancy defects in oxides, to train a defect graph neural network (dGNN) model that replaces the density functional theory supercell relaxations otherwise required for each symmetrically unique crystal site. Interfaced with thermodynamic calculations of reduction entropies and associated free energies, the dGNN model is applied to the screening of oxides in the Materials Project database, connecting the zero-kelvin defect enthalpies to high-temperature process conditions relevant for solar thermochemical hydrogen production and other energy applications. The dGNN approach is applicable to arbitrary structures with an accuracy limited principally by the amount and diversity of the training data, and it is generalizable to other defect types and advanced graph convolution architectures. In conclusion, it will help to tackle future materials discovery problems in clean energy and beyond.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Dielectric screening and vacancy formation for large neutral and charged Si n H m ( n > 1500 ) nanocrystals using real-space pseudopotentials

A commonly used procedure for computing the properties of defects in crystalline materials is to consider a large supercell that includes the defect of interest. This is a straightforward technique as standard energy band codes can be used for such computations. For neutral defects, the only impediment of such an approach is to avoid defect-defect interactions between adjoining cells. However, this procedure can be complex if the defect of interest is charged as the system at large contains Coulombic divergences. Moreover, some have recently argued that the conventional definition of formation energies for charged defects cannot be reconciled with statistical mechanics. Here, we focus on an alternative approach. We consider large nanocrystals wherein a charged defect can be placed. Since the system is confined, a charged defect within the nanocrystal does not result in a Coulombic divergence. The chief impediment is computational, i.e., while no defect-defect or Coulombic divergences are present, the nanocrystal must be sufficiently large to allow the system to properly replicate a bulklike configuration. With the development of new algorithms and hardware advances, computations for systems of sufficient size to address this issue are feasible. In particular, we solve the Kohn-Sham equation in real space using pseudopotential-density-functional theory for large silicon nanocrystals, which contain thousands of atoms. Further, we focus on (i) the screening of a point charge and (ii) the formation of a charged vacancy in hydrogen-terminated silicon nanocrystals. This approach allows us to examine the role of quantum confinement in addition to exploring the bulk limit. Comparisons to other methods confirm the viability of this approach.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Vacancy-mediated complex phase selection in high entropy alloys

Phase selection in Ti-Zr-Hf-Al high-entropy alloys was investigated by in-situ high-energy x-ray diffraction, single-crystal x-ray diffraction, and KKR-CPA electronic-structure methods that address disorder and vacancies, predicting formation enthalpy and chemical short-range order (SRO). Samples with varying Al content were synthesized, characterized, and computationally assessed to ascertain the composition-dependent phase selection, as increased Al content often acts as a stabilizer of a body-centered-cubic structure. Equiatomic TiZrHfAl was especially interesting due to its observed bcc superstructure – a variant of γ-brass with 4 vacancies per cell (not 2 as in γ-brass). We highlight how vacancy ordering mediates selection of this variant of γ-brass, which is driven by vacancy-atom SRO that dramatically suppress all atomic SRO. As vacancies are inherent in processing refractory systems, we expect that similar discoveries await in other high entropy alloys or in revisiting older experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2020

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

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

36 MATERIALS SCIENCE↗
DOE OSTI2021

Developing Accelerated Test Protocols and Tuning Microstructures of the Common Materials to Improve Robustness, Reliability, and Endurance of SOFC Cells

This work is built on our prior work where we found that phase transformation in praseodymium nickelates, e.g. Pr 2 NiO 4 (PNO) and (Pr 1-x Nd x ) 2 NiO 4 (PNNO), can be electrochemically driven, and is substantially faster when compared to thermal annealing studies. The first task aims at an attempt to further accelerate the phase transformation in the oxygen electrode by alternating the current input in the cells, which lead to the development of accelerated test protocols (ATPs). ATPs showed up to 60x faster phase transformation and up to 10x faster performance degradation in (Pr 0.50 Nd 0.50 ) 2 NiO 4 electrodes, when compared to long-term operation under constant current density. Furthermore, the phase stable Nd 2 NiO 4 and (La 0.6 Sr 0.4 )(Co 0.8 Fe 0.2 )O 3 (LSCF6482) electrodes were tested in full cells under ATPs, and showed up to 10x faster performance degradation within 1,100 hours in a comparison with long-term thermal annealing studies and electrochemical operation under constant current density. The second task aims at the quantification of the contributions of cell components to the total impedance of a solid oxide fuel cell (SOFC) using electrochemical impedance spectroscopy (EIS) and distribution of relaxation times (DRT). Specifically, the role of gas composition at both anode and cathode was systematically studied, aiming at deconvoluting, identifying and quantifying the contributions of different electrode processes. This was achieved by first tuning the partial pressure of H 2 at the anode and subsequently varying the partial pressure of O 2 at the cathode. The results suggest that, while DRT offers a viable way of deconvoluting different times distributions, additional attention is needed before assigning a peak to a specific electrode process due to the significant overlap of the contributions from the cathode and the anode. Density function theory studies show that both Pr-vacancies and O-defects play a key role on the activity and stability for nickelates towards oxygen reduction reaction. The resident O-interstitials and oxygen ions in the PrO layer form peroxide (O 2 2- ) nearby Pr vacancies. The O 2 2- limits oxygen-ion transport due to the required additional energy to break its O-O bond. We further calculated the formation and segregation energies for different Ln ions (La, Pr, Nd, Pm, Sm, Gd, Tb, Dy, and Ho) in PNO and CeO 2 (111) surfaces. In addition to Nd, Pm and La are suggested as potential dopants in PNO to enhance it stability without decomposition due to their more negative formation energies, lower diffusion energies, and positive separation energies.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

Impact of Sr-Containing Secondary Phases on Oxide Conductivity in Solid-Oxide Electrolyzer Cells

Solid-oxide electrolyzer cells (SOECs) based on a yttria-stabilized zirconia (YSZ) oxide electrolyte produce hydrogen from water with the assistance of excess thermal energy; however, Sr diffusion within the Gd-doped CeO 2 (GDC) barrier layer during processing or operation can lead to the formation of unwanted secondary phases such as SrO and SrZrO 3 . Here, to establish and compare the degree of impact of these phases on SOEC performance, we conduct first-principles calculations to study their bulk oxide conductivities and compare them to that of the YSZ electrolyte. We find that SrO has a low conductivity arising from the poor mobility and low concentration of mobile oxygen vacancies, and its presence in SOECs should therefore be avoided. SrZrO 3 also has a lower oxide conductivity than YSZ; however, this discrepancy is primarily due to lower vacancy concentrations rather than low mobility. We find that sufficient levels of Y-doping on the Zr site can increase oxygen vacancy concentrations in SrZrO 3 to achieve an oxide ionic conductivity on par with that of YSZ, thereby mitigating any potential deleterious effect on transport performance. Energy-dispersive X-ray spectroscopy confirms that Y is the most common minority element present in SrZrO 3 forming near the GDC–YSZ interface, alleviating concerns regarding the impact of SrZrO 3 on device performance. These results from our combined computational–experimental analysis can inform future engineering strategies designed to limit the detrimental effects of Sr-induced secondary phase formation on SOEC performance.

08 HYDROGEN↗
DOE OSTI2026

Tensile creep properties of NbTaTi and NbTaTiV refractory complex concentrated alloys

Refractory complex concentrated alloys (RCCAs) have emerged as promising candidates for high-temperature structural applications due to their high melting points and potential for exceptional strength, whereas their creep properties remain relatively underexplored. Here, we investigate the tensile creep behavior of NbTaTi and NbTaTiV in high vacuum to elucidate the influence of alloying on creep mechanisms. NbTaTiV exhibits higher creep resistance than NbTaTi, which is attributed to V-induced edge-dislocation-glide-controlled deformation. Microstructural analysis reveals the formation of Ti- and interstitial impurity-rich secondary phases during creep, which enhances the creep resistance without significant embrittlement. Stress exponents (∼3) and activation energies (∼those for vacancy diffusion) align with thermally activated dislocation glide-controlled creep. Despite NbTaTiV surpassing Ni-based superalloys in yield strength at elevated temperatures, its single-phase BCC structure limits creep resistance under engineering-relevant conditions, highlighting the necessity of deliberate secondary phase design to achieve competitive high-temperature performance.

Complex concentrated alloys↗
DOE OSTI2021

Multisublattice cluster expansion study of short-range ordering in iron-substituted strontium titanate

Owing to the challenges in obtaining realistic atomic configurations in large chemical phase spaces, it is not straightforward to describe structure–property relations in materials exhibiting configurational disorder. One example is iron-substituted strontium titanate (SrTi 1–x Fe x O 3–d , STF), a promising perovskite-derivative cathode material in solid oxide fuel cells that exhibits full solid solubility 0 ≤ x ≤ 1 and a tendency to exhibit short-range order. Here we demonstrate a multisublattice cluster expansion (CE) framework and apply it to STF across the full composition range. The CE approach is distinct from more traditional CE formulations in that clusters are defined explicitly by the chemical species distributed among multiple sublattices, rather than via cluster functions of occupation variables with decoration. The modified CE approach makes it easy to distinguish meaningful chemical interactions that are harder to extract from conventional CE, since for the latter chemical identity in a cluster is expressed as a product of site occupations. The least absolute shrinkage and selection operator (LASSO) is implemented as a regression analysis tool to select key clusters and avoid overfitting. We demonstrate this formulation on STF, and show that it can accurately predict configurational energies in comparison to conventional CE. From the key clusters, we identify that short-range ordering between substitutional Fe and oxygen vacancies (V O ) results in the formation of Fe–VO strings. In addition, we consider the stability of STF through CE-based Monte Carlo (MC) simulations and confirm the presence of superstructures that were previously observed in transmission electron microscopy. In this work, analysis of atomic configurations from MC samples reveals variations in the oxidation state of Fe atoms, which can be explained by the ordering tendency of Fe and V O . The cluster description and selection formalism described here may be applied to other disordered multisublattice systems for accurate and efficient material modeling.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Assessing Carrier Mobility, Dopability, and Defect Tolerance in the Chalcogenide Perovskite Ba Zr S 3

The chalcogenide perovskite Ba Zr S 3 has attracted much attention as a promising solar absorber for thin-film photovoltaics. Here we use first-principles calculations to evaluate its carrier transport and defect properties. We find that Ba Zr S 3 has a phonon-limited electron mobility of 37 cm 2 / V s, which is comparable to that in halide perovskites, but lower hole mobility of 11 cm 2 / V s. The defect computations indicate that Ba Zr S 3 is intrinsically n -type due to shallow sulfur vacancies, but that strong compensation by sulfur vacancies will prevent attempts to make it p -type. We also establish that Ba Zr S 3 shows some degree of defect tolerance, presenting only few low formation energy, deep intrinsic defects. Among the deep defects, sulfur interstitials are the dominant nonradiative recombination centers but exhibit a moderate capture coefficient. Our work highlights the material’s intrinsic limitations in carrier mobility and p -type doping, and suggests focusing on suppressing the formation of sulfur interstitials to achieve longer carrier lifetime. Published by the American Physical Society 2024

Yuan, Zhenkun (ORCID:0000000272005973)↗
DOE OSTI2024

Thermodynamic Modeling of Point Defects in Triple Conducting Perovskite Ba 0.95 La 0.05 FeO 3- $_δ$ with Incorporation of the Hydride Defect Formation Reaction for Solid Oxide Cells

Distinct from the proton defect, the hydride defect species may be present in certain perovskite materials in reducing environments such as in fuel electrodes of solid oxide cells (SOCs) or on the reducing side of ceramic membranes. A generalized defect thermodynamic model was developed for the triple-conducting perovskites (La,Ba)Fe 1-x M x O 3-δ (M = Y and Zr) to allow inclusion of the hydride defect formation reaction in addition to the other three main defect reactions, namely, the oxygen vacancy formation, hydration, and charge disproportionation reactions. This comprehensive defect model also allows the incorporation of polynomial functional forms of oxygen nonstoichiometry δ to describe the defect reaction energies and entropies and to enable refinements of the defect reaction equilibrium constants in the defect thermodynamic analysis. As a first step, the developed model is applied to the Ba 0.95 La 0.05 FeO 3-δ material as an illustrative system to obtain its Brouwer diagrams with both the proton and hydride defects in relevant SOC conditions, particularly for more reducing environments. In conclusion, the results provide direct guidance on the influence of electronic and ionic defect concentrations upon thermodynamic properties and ultimately on the performance of Ba 0.95 La 0.05 FeO 3-δ and potentially other (La,Ba)Fe 1-x M x O 3-δ perovskite materials involved in SOC applications.

25 ENERGY STORAGE↗
DOE OSTI2025

Suppressing Polaronic Defect–Photocarrier Interaction in Halide Perovskites by Pre-distorting Its Lattice

In halide perovskites, photocarriers can have strong polaronic interactions with point defects. For iodide-deficient MAPbI 3 , we found that the Fermi level can shift significantly by 0.6–0.7 eV upon light illumination. This energy level shift is accompanied by the formation of deep electron traps. These experimental observations are consistent with the formation of a Pb–Pb dimer when photoexcited electrons are trapped at an iodide vacancy. Interestingly, we found that this polaronic interaction is suppressed when a portion of MA + cations is replaced by smaller Cs + ions. Density functional theory calculations reveal that Cs-doping can reduce the distance between two Pb atoms across an iodide vacancy, even without electron trapping. The predistortion of the lattice induced by cation replacement resembles the Pb–Pb dimer formed by electron trapping at the defect site, which explains the suppression of light-induced effects observed in the experiment. Furthermore, our finding unveils a counterintuitive strategy to enhance the photostability of halide perovskites by preintroducing distortions into its lattice.

14 SOLAR ENERGY↗
DOE OSTI2022

Ammonia synthesis on BaTiO 2.5 H 0.5 : computational insights into the role of hydrides

Perovskite oxyhydrides such as BaTiO 2.5 H 0.5 have been found to be able to catalyze NH3 synthesis, but the mechanism and the role of the catalyst's lattice hydrides in the catalytic reaction remain unknown. Here we employ first principles density functional theory to investigate the mechanism of ammonia synthesis and the role of lattice hydrides on a prototypical perovskite oxyhydride, BaTiO 2.5 H 0.5 (BTOH). Two mechanistic hypotheses, the distal and alternating pathways, have been tested on the Ti2O 2 termination of the BTOH (210) surface, previously determined to be the most stable surface termination under the reaction conditions considered. In the distal pathway, H atoms hydrogenate N 2 to form the *N–NH x key intermediates, followed by N–N bond breaking. In the alternating pathway, H atoms hydrogenate N 2 in an alternating fashion to form the *NH x –NH y intermediates before N–N bond breaking and formation of co-adsorbed *NH x /*NH y on the surface. We find that the subsurface hydride vacancy formed after reaction of *N 2 with the lattice hydride is key to the distal pathway, leading to surface nitride formation after breaking the *N–NH 3 bond, while the neighboring surface Ti sites are key to bridging and stabilizing the *NNH intermediate in the alternating pathway. In both pathways, desorption of NH 3 is the most uphill in energy. Our results provide important insights into the role of hydrides and surface vacancies in hydrogenation reactions over BTOH, which will be useful to guide future spectroscopic experiments such as operando IR and inelastic neutron scattering to verify the key intermediates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗