Steering CO2 hydrogenation toward C-C coupling to hydrocarbons using porous organic polymer/metal interfaces.

The conversion of CO2 into fuels and chemicals is an attractive option for mitigating CO2 emissions. Controlling the selectivity of this process is beneficial to produce desirable liquid fuels, but C-C coupling is a limiting step in the reaction that requires high pressures. Here, we propose a strategy to favor C-C coupling on a supported Ru/TiO2 catalyst by encapsulating it within the polymer layers of an imine-based porous organic polymer that controls its selectivity. Such polymer confinement modifies the CO2 hydrogenation behavior of the Ru surface, significantly enhancing the C2+ production turnover frequency by 10-fold. We demonstrate that the polymer layers affect the adsorption of reactants and intermediates while being stable under the demanding reaction conditions. Our findings highlight the promising opportunity of using polymer/metal interfaces for the rational engineering of active sites and as a general tool for controlling selective transformations in supported catalyst systems.

Reduction of Cofed Carbon Dioxide Modifies the Local Coordination Environment of Zeolite-Supported, Atomically Dispersed Chromium to Promote Ethane Dehydrogenation.

The reduction of CO2 is known to promote increased alkene yields from alkane dehydrogenations when the reactions are cocatalyzed. The mechanism of this promotion is not understood in the context of catalyst active-site environments because CO2 is amphoteric, and even general aspects of the chemistry, including the significance of competing side reactions, differ significantly across catalysts. Atomically dispersed chromium cations stabilized in highly siliceous MFI zeolite are shown here to enable the study of the role of parallel CO2 reduction during ethylene-selective ethane dehydrogenation. Based on infrared spectroscopy and X-ray absorption spectroscopy data interpreted through calculations using density functional theory (DFT), the synthesized catalyst contains atomically dispersed Cr cations stabilized by silanol nests in micropores. Reactor studies show that cofeeding CO2 increases stable ethylene-selective ethane dehydrogenation rates over a wide range of partial pressures. Operando X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine-structure (EXAFS) spectra indicate that during reaction at 650 °C the Cr cations maintain a nominal 2+ charge and a total Cr-O coordination number of approximately 2. However, CO2 reduction induces a change, correlated with the CO2 partial pressure, in the population of two distinct Cr-O scattering paths. This indicates that the promotional effect of parallel CO2 reduction can be attributed to a subtle change in Cr-O bond lengths in the local coordination environment of the active site. These insights are made possible by simultaneously fitting multiple EXAFS spectra recorded in different reaction conditions; this novel procedure is expected to be generally applicable for interpreting operando catalysis EXAFS data.

Understanding and Harnessing Nanoscale Immiscibility in Ru-In Alloys for Selective CO2 Hydrogenation.

Bimetallic alloys made from immiscible elements are characterized by their tendency to segregate on the macroscopic scale, but their behavior is known to change at the nanoscale. Here, we demonstrate that in the Ru-In system, In atoms preferentially decorate the surface of 6 nm Ru nanoparticles, forming Ru-In superficial immiscible alloys. This surface decoration dramatically affects the catalytic performance of the system, even at small atomic fractions of In added to Ru. The interfaces between Ru and In enabled unexplored methanol productivity from CO2 hydrogenation, which outperformed not only the individual constituents but also ordered RuIn3 intermetallic alloys. Our work highlights that the formation of superficial immiscible alloys could offer new insights into the understanding and design of heterogeneous catalysts.

Structure and Site Evolution of Framework Ni Species in MIL-127 MOFs for Propylene Oligomerization Catalysis.

A mixed-valence oxotrimer metal-organic framework (MOF), Ni-MIL-127, with a fully coordinated nickel atom and two iron atoms in the inorganic node, generates a missing linker defect upon thermal treatment in helium (>473 K) to engender an open coordination site on nickel which catalyzes propylene oligomerization devoid of any cocatalysts or initiators. This catalyst is stable for ∼20 h on stream at 500 kPa and 473 K, unprecedented for this chemistry. The number of missing linkers on synthesized and activated Ni-MIL-127 MOFs is quantified using temperature-programmed oxidation, 1H nuclear magnetic resonance spectroscopy, and X-ray absorption spectroscopy to be ∼0.7 missing linkers per nickel; thus, a majority of Ni species in the MOF framework catalyze propylene oligomerization. In situ NO titrations under reaction conditions enumerate ∼62% of the nickel atoms as catalytically relevant to validate the defect density upon thermal treatment. Propylene oligomerization rates on Ni-MIL-127 measured at steady state have activation energies of 55-67 kJ mol-1 from 448 to 493 K and are first-order in propylene pressures from 5 to 550 kPa. Density functional theory calculations on cluster models of Ni-MIL-127 are employed to validate the plausibility of the missing linker defect and the Cossee-Arlman mechanism for propylene oligomerization through comparisons between apparent activation energies from steady-state kinetics and computation. This study illustrates how MOF precatalysts engender defective Ni species which exhibit reactivity and stability characteristics that are distinct and can be engineered to improve catalytic activity for olefin oligomerization.

Dynamic Evolution of Palladium Single Atoms on Anatase Titania Support Determines the Reverse Water-Gas Shift Activity.

Research interest in single-atom catalysts (SACs) has been continuously increasing. However, the lack of understanding of the dynamic behaviors of SACs during applications hinders catalyst development and mechanistic understanding. Herein, we report on the evolution of active sites over Pd/TiO2-anatase SAC (Pd1/TiO2) in the reverse water-gas shift (rWGS) reaction. Combining kinetics, in situ characterization, and theory, we show that at T ≥ 350 °C, the reduction of TiO2 by H2 alters the coordination environment of Pd, creating Pd sites with partially cleaved Pd-O interfacial bonds and a unique electronic structure that exhibit high intrinsic rWGS activity through the carboxyl pathway. The activation by H2 is accompanied by the partial sintering of single Pd atoms (Pd1) into disordered, flat, ∼1 nm diameter clusters (Pdn). The highly active Pd sites in the new coordination environment under H2 are eliminated by oxidation, which, when performed at a high temperature, also redisperses Pdn and facilitates the reduction of TiO2. In contrast, Pd1 sinters into crystalline, ∼5 nm particles (PdNP) during CO treatment, deactivating Pd1/TiO2. During the rWGS reaction, the two Pd evolution pathways coexist. The activation by H2 dominates, leading to the increasing rate with time-on-stream, and steady-state Pd active sites similar to the ones formed under H2. This work demonstrates how the coordination environment and nuclearity of metal sites on a SAC evolve during catalysis and pretreatments and how their activity is modulated by these behaviors. These insights on SAC dynamics and the structure-function relationship are valuable to mechanistic understanding and catalyst design.

Genesis of Active Pt/CeO2 Catalyst for Dry Reforming of Methane by Reduction and Aggregation of Isolated Platinum Atoms into Clusters.

Atomically dispersed metal catalysts offer the advantages of efficient metal utilization and high selectivities for reactions of technological importance. Such catalysts have been suggested to be strong candidates for dry reforming of methane (DRM), offering prospects of high selectivity for synthesis gas without coke formation, which requires ensembles of metal sites and is a challenge to overcome in DRM catalysis. However, investigations of the structures of isolated metal sites on metal oxide supports under DRM conditions are lacking, and the catalytically active sites remain undetermined. Data characterizing the DRM reaction-driven structural evolution of a cerium oxide-supported catalyst, initially incorporating atomically dispersed platinum, and the corresponding changes in catalyst performance are reported. X-ray absorption and infrared spectra show that the reduction and agglomeration of isolated cationic platinum atoms to form small platinum clusters/nanoparticles are necessary for DRM activity. Density functional theory calculations of the energy barriers for methane dissociation on atomically dispersed platinum and on platinum clusters support these observations. The results emphasize the need for in-operando experiments to assess the active sites in such catalysts. The inferences about the catalytically active species are suggested to pertain to a broad class of catalytic conversions involving the rate-limiting dissociation of light alkanes.

CatMass: software for calculating optimal sample masses for X-ray absorption spectroscopy experiments involving complex sample compositions.

This paper presents software for calculating the optimal mass of samples with complex compositions (e.g. supported metal catalysts) for X-ray absorption spectroscopy (XAS) and scattering measurements. The ability to calculate the sample mass and other relevant parameters needed for an XAS measurement allows experimentalists to be better prepared in terms of detector selection, energy range of scan and overall time needed to complete the measurement, thus increasing efficiency. CatMass builds on existing sample mass calculators allowing users to determine the optimum sample preparation, collection geometry, usable energy range for a scan and approximate edge step of the absorption event. Visualization tools present the absorption calculation results in a format familiar to XAS experimentalists, with the added ability to save calculations and plots for future reference or recalculation. CatMass is a program broadly applicable in catalysis and is helpful for users with complex samples due to composition/stoichiometry or multiple competing elements.

Templated encapsulation of platinum-based catalysts promotes high-temperature stability to 1,100 °C.

Stable catalysts are essential to address energy and environmental challenges, especially for applications in harsh environments (for example, high temperature, oxidizing atmosphere and steam). In such conditions, supported metal catalysts deactivate due to sintering-a process where initially small nanoparticles grow into larger ones with reduced active surface area-but strategies to stabilize them can lead to decreased performance. Here we report stable catalysts prepared through the encapsulation of platinum nanoparticles inside an alumina framework, which was formed by depositing an alumina precursor within a separately prepared porous organic framework impregnated with platinum nanoparticles. These catalysts do not sinter at 800 °C in the presence of oxygen and steam, conditions in which conventional catalysts sinter to a large extent, while showing similar reaction rates. Extending this approach to Pd-Pt bimetallic catalysts led to the small particle size being maintained at temperatures as high as 1,100 °C in air and 10% steam. This strategy can be broadly applied to other metal and metal oxides for applications where sintering is a major cause of material deactivation.

Aqueous Structure of Lanthanide-EDTA Coordination Complexes Determined by a Combined DFT/EXAFS Approach.

Sustainable production of rare earth elements (REEs) is critical for technologies needed for climate change mitigation, including wind turbines and electric vehicles. However, separation technologies currently used in REE production have large environmental footprints, necessitating more sustainable strategies. Aqueous, affinity-based separations are examples of such strategies. To make these technologies feasible, it is imperative to connect aqueous ligand structure to ligand selectivity for individual REEs. As a step toward this goal, we analyzed the extended X-ray absorption fine structure (EXAFS) of four lanthanides (La, Ce, Pr, and Nd) complexed by a common REE chelator, ethylenediaminetetraacetic acid (EDTA) to determine the aqueous-phase structure. Reference structures from density functional theory (DFT) were used to help fit the EXAFS spectra. We found that all four Ln-EDTA coordination complexes formed 9-coordinate structures with 6 coordinating atoms from EDTA (4 carboxyl oxygen atoms and 2 nitrogen atoms) and 3 oxygen atoms from water molecules. All EXAFS fits were of high quality (R-factor < 0.02) and showed decreasing average first-shell coordination distance across the series (2.62-2.57 Å from La-Nd), in agreement with DFT (2.65-2.56 Å from La-Nd). The insights determined herein will be useful in the development of ligands for sustainable rare earth elements (REE) separation technologies.

Size-Dependent Dispersion of Rhodium Clusters into Isolated Single Atoms at Low Temperature and the Consequences for CO Oxidation Activity.

Understanding the dynamic structural evolution of supported metal clusters under reaction conditions is crucial to develop structure reactivity relations. Here, we followed the structure of different size Rh clusters supported on Al2 O3 using in situ/operando spectroscopy and ex situ aberration-corrected electron microscopy. We report a dynamic evolution of rhodium clusters into thermally stable isolated single atoms upon exposure to oxygen and during CO oxidation. Rh clusters partially disperse into single atoms at room temperature and the extent of dispersion increases as the Rh size decreases and as the reaction temperature increases. A strong correlation is found between the extent of dispersion and the CO oxidation kinetics. More importantly, dispersing Rh clusters into single atoms increases the activity at room temperature by more than two orders of magnitude due to the much lower activation energy on single atoms (40 vs. 130 kJ/mol). This work demonstrates that the structure and reactivity of small Rh clusters are very sensitive to the reaction environment.

Catalytic Performance and Near-Surface X-ray Characterization of Titanium Hydride Electrodes for the Electrochemical Nitrate Reduction Reaction.

The electrochemical nitrate reduction reaction (NO3RR) on titanium introduces significant surface reconstruction and forms titanium hydride (TiHx, 0 < x ≤ 2). With ex situ grazing-incidence X-ray diffraction (GIXRD) and X-ray absorption spectroscopy (XAS), we demonstrated near-surface TiH2 enrichment with increasing NO3RR applied potential and duration. This quantitative relationship facilitated electrochemical treatment of Ti to form TiH2/Ti electrodes for use in NO3RR, thereby decoupling hydride formation from NO3RR performance. A wide range of NO3RR activity and selectivity on TiH2/Ti electrodes between -0.4 and -1.0 VRHE was observed and analyzed with density functional theory (DFT) calculations on TiH2(111). This work underscores the importance of relating NO3RR performance with near-surface electrode structure to advance catalyst design and operation.

Propane Dehydrogenation Catalyzed by Isolated Pt Atoms in ≡SiOZn-OH Nests in Dealuminated Zeolite Beta.

Atomically dispersed noble metal catalysts have drawn wide attention as candidates to replace supported metal clusters and metal nanoparticles. Atomic dispersion can offer unique chemical properties as well as maximum utilization of the expensive metals. Addition of a second metal has been found to help reduce the size of Pt ensembles in bimetallic clusters; however, the stabilization of isolated Pt atoms in small nests of nonprecious metal atoms remains challenging. We now report a novel strategy for the design, synthesis, and characterization of a zeolite-supported propane dehydrogenation catalyst that incorporates predominantly isolated Pt atoms stably bonded within nests of Zn atoms located within the nanoscale pores of dealuminated zeolite Beta. The catalyst is stable in long-term operation and exhibits high activity and high selectivity to propene. Atomic resolution images, bolstered by X-ray absorption spectra, demonstrate predominantly atomic dispersion of the Pt in the nests and, with complementary infrared and nuclear magnetic resonance spectra, determine a structural model of the nested Pt.

Beyond Radical Rebound: Methane Oxidation to Methanol Catalyzed by Iron Species in Metal-Organic Framework Nodes.

Recent work has exploited the ability of metal-organic frameworks (MOFs) to isolate Fe sites that mimic the structures of sites in enzymes that catalyze selective oxidations at low temperatures, opening new pathways for the valorization of underutilized feedstocks such as methane. Questions remain as to whether the radical-rebound mechanism commonly invoked in enzymatic and homogeneous systems also applies in these rigid-framework materials, in which resisting the overoxidation of desired products is a major challenge. We demonstrate that MOFs bearing Fe(II) sites within Fe3-µ3-oxo nodes active for conversion of CH4 + N2O mixtures (368-408 K) require steps beyond the radical-rebound mechanism to protect the desired CH3OH product. Infrared spectra and density functional theory show that CH3OH(g) is stabilized as Fe(III)-OCH3 groups on the MOF via hydrogen atom transfer with Fe(III)-OH groups, eliminating water. Consequently, upon addition of a protonic zeolite in inter- and intrapellet mixtures with the MOF, we observed increases in CH3OH selectivity with increasing ratio and proximity of zeolitic H+ to MOF-based Fe(II) sites, as methanol is protected within the zeolite. We infer from the data that CH3OH(g) is formed via the radical-rebound mechanism on Fe(II) sites but that subsequent transport and dehydration steps are required to protect CH3OH(g) from overoxidation. The results demonstrate that the radical-rebound mechanism commonly invoked in this chemistry is insufficient to explain the reactivity of these systems, that the selectivity-controlling steps involve both chemical and physical rate phenomena, as well as offering a strategy to mitigate overoxidation in these and similar systems.

Mechanistic and Electronic Insights into a Working NiAu Single-Atom Alloy Ethanol Dehydrogenation Catalyst.

Elucidation of reaction mechanisms and the geometric and electronic structure of the active sites themselves is a challenging, yet essential task in the design of new heterogeneous catalysts. Such investigations are best implemented via a multipronged approach that comprises ambient pressure catalysis, surface science, and theory. Herein, we employ this strategy to understand the workings of NiAu single-atom alloy (SAA) catalysts for the selective nonoxidative dehydrogenation of ethanol to acetaldehyde and hydrogen. The atomic dispersion of Ni is paramount for selective ethanol to acetaldehyde conversion, and we show that even the presence of small Ni ensembles in the Au surface results in the formation of undesirable byproducts via C-C scission. Spectroscopic, kinetic, and theoretical investigations of the reaction mechanism reveal that both C-H and O-H bond cleavage steps are kinetically relevant and single Ni atoms are confirmed as the active sites. X-ray absorption spectroscopy studies allow us to follow the charge of the Ni atoms in the Au host before, under, and after a reaction cycle. Specifically, in the pristine state the Ni atoms carry a partial positive charge that increases upon coordination to the electronegative oxygen in ethanol and decreases upon desorption. This type of oxidation state cycling during reaction is similar to the behavior of single-site homogeneous catalysts. Given the unique electronic structure of many single-site catalysts, such a combined approach in which the atomic-scale catalyst structure and charge state of the single atom dopant can be monitored as a function of its reactive environment is a key step toward developing structure-function relationships that inform the design of new catalysts.

A Theory-Guided X-ray Absorption Spectroscopy Approach for Identifying Active Sites in Atomically Dispersed Transition-Metal Catalysts.

Atomically dispersed supported metal catalysts offer new properties and the benefits of maximized metal accessibility and utilization. The characterization of these materials, however, remains challenging. Using atomically dispersed platinum supported on crystalline MgO (chosen for its well-defined bonding sites) as a prototypical example, we demonstrate how systematic density functional theory calculations for assessing all the potentially stable platinum sites, combined with automated analysis of extended X-ray absorption fine structure (EXAFS) spectra, leads to unbiased identification of isolated, surface-enveloped platinum cations as the catalytic species for CO oxidation. The catalyst has been characterized by atomic-resolution imaging and EXAFS and high-energy resolution fluorescence detection X-ray absorption near edge spectroscopy. The proposed platinum sites are in agreement with experiment. This theory-guided workflow leads to rigorously determined structural models and provides a more detailed picture of the structure of the catalytically active site than what is currently possible with conventional EXAFS analyses. As this approach is efficient and agnostic to the metal, support, and catalytic reaction, we posit that it will be of broad interest to the materials characterization and catalysis communities.

Tailoring the Local Environment of Platinum in Single-Atom Pt1 /CeO2 Catalysts for Robust Low-Temperature CO Oxidation.

A single-atom Pt1 /CeO2 catalyst formed by atom trapping (AT, 800 °C in air) shows excellent thermal stability but is inactive for CO oxidation at low temperatures owing to over-stabilization of Pt2+ in a highly symmetric square-planar Pt1 O4 coordination environment. Reductive activation to form Pt nanoparticles (NPs) results in enhanced activity; however, the NPs are easily oxidized, leading to drastic activity loss. Herein we show that tailoring the local environment of isolated Pt2+ by thermal-shock (TS) synthesis leads to a highly active and thermally stable Pt1 /CeO2 catalyst. Ultrafast shockwaves (>1200 °C) in an inert atmosphere induced surface reconstruction of CeO2 to generate Pt single atoms in an asymmetric Pt1 O4 configuration. Owing to this unique coordination, Pt1 δ+ in a partially reduced state dynamically evolves during CO oxidation, resulting in exceptional low-temperature performance. CO oxidation reactivity on the Pt1 /CeO2 _TS catalyst was retained under oxidizing conditions.

Selective Methanol Carbonylation to Acetic Acid on Heterogeneous Atomically Dispersed ReO4/SiO2 Catalysts.

Methanol carbonylation to acetic acid (AA) is a large-scale commodity chemical production process that requires homogeneous liquid-phase organometallic catalysts with corrosive halide-based cocatalysts to achieve high selectivity and activity. Here, we demonstrate a heterogeneous catalyst based on atomically dispersed rhenium (ReO4) active sites on an inert support (SiO2) for the halide-free, gas phase carbonylation of methanol to AA. Atomically dispersed ReO4 species and nanometer sized ReOx clusters were deposited on a high surface area (700 m2/g) inert SiO2 using triethanolamine as a dispersion promoter and characterized using aberration corrected scanning transmission electron microscopy (AC-STEM), UV-vis spectroscopy, and X-ray absorption spectroscopy (XAS). Reactivity measurements at atmospheric pressure with 30 mbar of methanol and CO (1:1 molar ratio) showed that bulk Re2O7 and ReOx clusters on SiO2 (formed at >10 wt %) were selective for dimethyl ether formation, while atomically dispersed ReO4 on SiO2 (formed at <10 wt %) exhibited stable (for 60 h) > 93% selectivity to AA with single pass conversion >60%. Kinetic analysis, in situ FTIR, and in situ XAS measurements suggest that the AA formation mechanism involves methanol activation on ReO4, followed by CO insertion into the terminal methyl species. Further, the introduction of ∼0.2 wt % of atomically dispersed Rh to 10 wt % atomically dispersed ReO4 on SiO2 resulted in >96% selectivity toward AA production at volumetric reaction rates comparable to homogeneous processes. This work introduces a new class of promising heterogeneous catalysts based on atomically dispersed ReO4 on inert supports for alcohol carbonylation.

Uniformity Is Key in Defining Structure-Function Relationships for Atomically Dispersed Metal Catalysts: The Case of Pt/CeO2.

Catalysts consisting of atomically dispersed Pt (Ptiso) species on CeO2 supports have received recent interest due to their potential for efficient metal utilization in catalytic convertors. However, discrepancies exist between the behavior (reducibility, interaction strength with adsorbates) of high surface area Ptiso/CeO2 systems and of well-defined surface science and computational model systems, suggesting differences in Pt local coordination in the two classes of materials. Here, we reconcile these differences by demonstrating that high surface area Ptiso/CeO2 synthesized at low Pt loadings (<0.1% weight) exhibit resistance to reduction and sintering up to 500 °C in 0.05 bar H2 and minimal interactions with CO-properties previously seen only for model system studies. Alternatively, Pt loadings >0.1 weight % produce a distribution of sub-nanometer Pt structures, which are difficult to distinguish using common characterization techniques, and exhibit strong interactions with CO and weak resistance to sintering, even in 0.05 bar H2 at 50 °C-properties previously seen for high surface area materials. This work demonstrates that low metal loadings can be used to selectively populate the most thermodynamically stable adsorption sites on high surface area supports with atomically dispersed metals. Further, the site uniformity afforded by this synthetic approach is critical for the development of relationships between atomic scale local coordination and functional properties. Comparisons to recent studies of Ptiso/TiO2 suggest a general compromise between the stability of atomically dispersed metal catalysts and their ability to interact with and activate molecular species.

Structural evolution of atomically dispersed Pt catalysts dictates reactivity.

The use of oxide-supported isolated Pt-group metal atoms as catalytic active sites is of interest due to their unique reactivity and efficient metal utilization. However, relationships between the structure of these active sites, their dynamic response to environments and catalytic functionality have proved difficult to experimentally establish. Here, sinter-resistant catalysts where Pt was deposited uniformly as isolated atoms in well-defined locations on anatase TiO2 nanoparticle supports were used to develop such relationships. Through a combination of in situ atomic-resolution microscopy- and spectroscopy-based characterization supported by first-principles calculations it was demonstrated that isolated Pt species can adopt a range of local coordination environments and oxidation states, which evolve in response to varied environmental conditions. The variation in local coordination showed a strong influence on the chemical reactivity and could be exploited to control the catalytic performance.

Insights into Copper Sulfide Formation from Cu and S K edge XAS and DFT studies.

An understanding of the fundamentals of the reaction between CuO with trace amounts of H2S to form CuS products is critical for the optimal utilization of this process in sulfur removal applications. Unfortunately, CuS is a complex material, featuring various Cu2-xS compounds (with 0 ≤ x ≤ 1), distorted crystal phases, and varying electronic structures and coordination environments of Cu and S ions. In this work, we combine ex situ and in situ X-ray absorption spectroscopy (XAS) at S and Cu K edges, fixed bed sorption experiments, DFT simulations, and other characterization techniques to speciate the CuS products formed at different temperatures (298-383 K) and from CuO sorbents with different crystallite sizes (2.8-40 nm). The results of our analysis identify the formation of a distorted CuS layer at the surface of CuO crystals with disulfide groups with shorter Cu-S bonds and higher delocalization of the positive charge of the Cu center into (S1-)2. This distorted CuS layer dominates the XAS signal at lower temperatures (298-323 K) and at the initial stages of sulfidation at higher temperatures (353 and 383 K) where conversion is low (<40%). First-principles atomistic simulations confirm the thermodynamic favorability of the formation of surface (S1-)2 on both CuO (111) and (1̅11) surfaces, providing further support for our experimental observations. Furthermore, these simulations reveal that the presence of disulfide bonds stabilized surface hydroxyl groups, leading to lower Gibbs Free Energies of their surface migration.