Atomically Precise Single-Site Catalysts via Exsolution in a Polyoxometalate-Metal-Organic-Framework Architecture.

Single-site catalysts (SSCs) achieve a high catalytic performance through atomically dispersed active sites. A challenge facing the development of SSCs is aggregation of active catalytic species. Reducing the loading of these sites to very low levels is a common strategy to mitigate aggregation and sintering; however, this limits the tools that can be used to characterize the SSCs. Here we report a sintering-resistant SSC with high loading that is achieved by incorporating Anderson-Evans polyoxometalate clusters (POMs, MMo6O24, M = Rh/Pt) within NU-1000, a Zr-based metal-organic framework (MOF). The dual confinement provided by isolating the active site within the POM, then isolating the POMs within the MOF, facilitates the formation of isolated noble metal sites with low coordination numbers via exsolution from the POM during activation. The high loading (up to 3.2 wt %) that can be achieved without sintering allowed the local structure transformation in the POM cluster and the surrounding MOF to be evaluated using in situ X-ray scattering with pair distribution function (PDF) analysis. Notably, the Rh/Pt···Mo distance in the active catalyst is shorter than the M···M bond lengths in the respective bulk metals. Models of the active cluster structure were identified based on the PDF data with complementary computation and X-ray absorption spectroscopy analysis.

Investigation of Rare Earth Element Binding to a Surface-Bound Affinity Peptide Derived from EF-Hand Loop I of Lanmodulin.

Bioinspired strategies have been given extensive attention for the recovery of rare earth elements (REEs) from waste streams because of their high selectivity, regeneration potential, and sustainability as well as low cost. Lanmodulin protein is an emerging biotechnology that is highly selective for REE binding. Mimicking lanmodulin with shorter peptides is advantageous because they are simpler and potentially easier to manipulate and optimize. Lanmodulin-derived peptides have been found to bind REEs, but their properties have not been explored when immobilized on solid substrates, which is required for many advanced separation technologies. Here, two peptides, LanM1 and scrambled LanM1, are designed from the EF-hand loop 1 of lanmodulin and investigated for their binding affinity toward different REEs when surface-bound. First, the ability of LanM1 to bind REEs was confirmed and characterized in solution using circular dichroism (CD), nuclear magnetic resonance (NMR), and molecular dynamics (MD) simulations for Ce(III) ions. Isothermal titration calorimetry (ITC) was used to further analyze the binding of the LanM1 to Ce(III), Nd(III), Eu(III), and Y(III) ions and in low-pH conditions. The performance of the immobilized peptides on a model gold surface was examined using a quartz crystal microbalance with dissipation (QCM-D). The studies show that the LanM1 peptide has a stronger REE binding affinity than that of scrambled LanM1 when in solution and when immobilized on a gold surface. QCM-D data were fit to the Langmuir adsorption model to estimate the surface-bound dissociation constant (Kd) of LanM1 with Ce(III) and Nd(III). The results indicate that LanM1 peptides maintain a high affinity for REEs when immobilized, and surface-bound LanM1 has no affinity for potential competitor calcium and copper ions. The utility of surface-bound LanM1 peptides was further demonstrated by immobilizing them to gold nanoparticles (GNPs) and capturing REEs from solution in experiments utilizing an Arsenazo III-based colorimetric dye displacement assay and ultraviolet-visible (UV-vis) spectrophotometry. The saturated adsorption capacity of GNPs was estimated to be around 3.5 µmol REE/g for Ce(III), Nd(III), Eu(III), and Y(III) ions, with no binding of non-REE Ca(II) ions observed.

Poisoning of Pt/γ-Al2O3 Aqueous Phase Reforming Catalysts by Ketone and Diketone-Derived Surface Species.

Strong adsorption of ketone and diketone byproducts and their fragmentation products during the aqueous phase reforming of biomass derived oxygenates is believed to be responsible for the deactivation of supported Pt catalysts. This study involves a combined experimental and theoretical approach to demonstrate the interactions of several model di/ketone poisons with Pt/γ-Al2O3 catalysts. Particular di/ketones were selected to reveal the effects of hydroxyl groups (acetone, hydroxyacetone), conjugation with C=C bonds (mesityl oxide), intramolecular distance between carbonyls in diketones (2,3-butanedione, 2,4-pentanedione), and length of terminal alkyl chains (3,4-hexanedione). The formation of adsorbed carbon monoxide (1900-2100 cm-1) as a decarbonylation product was probed using infrared spectroscopy and to calculate the extent of poisoning during subsequent methanol dehydrogenation based on the reduction of the ν(C≡O) band integral relative to experiments in which only methanol was dosed. Small Pt particles appeared less active in decarbonylation and were perhaps poisoned by strongly adsorbed di/ketones on undercoordinated metal sites and bulky conjugated species formed on the γ-Al2O3 support from aldol self-condensation. Larger Pt particles were more resistant to di/ketone poisoning due to higher decarbonylation activity yet still fell short of the expected yield of adsorbed CO from subsequent methanol activity. Vibrational spectra acquired using inelastic neutron scattering showed evidence for strongly binding methyl and acyl groups resulting from di/ketone decarbonylation on a Pt sponge at 250 °C. Adsorption energies and molecular configurations were obtained for di/ketones on a Pt(111) slab using density functional theory, revealing potential descriptors for predicting decarbonylation activity on highly coordinated metal sites. Calculated reaction energies suggest it is energetically favorable to reform surface methyl groups into adsorbed CO and H. However, the rate of this surface reaction is limited by a high activation barrier indicating that either improved APR catalyst designs or regeneration procedures may be necessary.

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.

Synthetic Access to a Framework-Stabilized and Fully Sulfided Analogue of an Anderson Polyoxometalate that is Catalytically Competent for Reduction Reactions.

Polyoxometalates (POMs) featuring 7, 12, 18, or more redox-accessible transition metal ions are ubiquitous as selective catalysts, especially for oxidation reactions. The corresponding synthetic and catalytic chemistry of stable, discrete, capping-ligand-free polythiometalates (PTMs), which could be especially attractive for reduction reactions, is much less well developed. Among the challenges are the propensity of PTMs to agglomerate and the tendency for agglomeration to block reactant access of catalyst active sites. Nevertheless, the pervasive presence of transition metal sulfur clusters metalloenzymes or cofactors that catalyze reduction reactions and the justifiable proliferation of studies of two-dimensional (2D) metal-chalcogenides as reduction catalysts point to the promise of well-defined and controllable PTMs as reduction catalysts. Here, we report the fabrication of agglomeration-immune, reactant-accessible, capping-ligand-free CoIIMo6IVS24n- clusters as periodic arrays in a water-stable, hierarchically porous Zr-metal-organic framework (MOF; NU1K) by first installing a disk-like Anderson polyoxometalate, CoIIIMo6VIO24m-, in size-matched micropores where the siting is established via difference electron density (DED) X-ray diffraction (XRD) experiments. Flowing H2S, while heating, reduces molybdenum(VI) ions to Mo(IV) and quantitatively replaces oxygen anions with sulfur anions (S2-, HS-, S22-). DED maps show that MOF-templated POM-to-PTM conversion leaves clusters individually isolated in open-channel-connected micropores. The structure of the immobilized cluster as determined, in part, by X-ray photoelectron spectroscopy (XPS), X-ray absorption fine structure (XAFS) analysis, and pair distribution function (PDF) analysis of total X-ray scattering agrees well with the theoretically simulated structure. PTM@MOF displays both electrocatalytic and photocatalytic competency for hydrogen evolution. Nevertheless, the initially installed PTM appears to be a precatalyst, gaining competency only after the loss of ∼3 to 6 sulfurs and exposure to hydride-forming metal ions.

Differences in solvation thermodynamics of oxygenates at Pt/Al2O3 perimeter versus Pt(111) terrace sites.

A prominent role of water in aqueous-phase heterogeneous catalysis is to modify free energies; however, intuition about how is based largely on pure metal surfaces or even homogeneous solutions. Using multiscale modeling with explicit liquid water molecules, we show that the influence of water on the free energies of adsorbates at metal/support interfaces is different than that on pure metal surfaces. We specifically compute free energies of solvation for methanol and its constituents on a Pt/Al2O3 catalyst and compare the results to analogous values calculated on a pure Pt catalyst. We find that the more hydrophilic Pt/Al2O3 interface leads to smaller (more positive) free energies of solvation due to an increased entropy penalty resulting from the additional work necessary to disrupt the interfacial water structure and accommodate the interfacial species. The results will be of interest in other fields, including adsorption and proteins.

Computational and Experimental Characterization of the Ligand Environment of a Ni-Oxo Catalyst Supported in the Metal-Organic Framework NU-1000.

Heterogeneous catalysts exhibit significant changes in composition due to the influence of operating conditions, and these compositional changes can have dramatic effects on catalytic performance. For traditional bulk metal heterogeneous catalysts, relationships between composition and catalytic operating conditions are well documented. However, the influence of operating conditions on the compositions of single-site heterogeneous catalysts remains largely unresolved. To address this, we report a combined computational and experimental characterization of a Ni oxo catalyst under catalytic hydrogenation conditions. Specifically, pair distribution function (PDF) analysis is combined with ab initio thermodynamic modeling to investigate ligand environments present on a Ni oxo cluster supported in the metal-organic framework NU-1000. Comparisons of the experimentally observed and simulated Ni-O coordination numbers and Ni-O, Ni···Ni, and Ni···Zr distances provide insight into the Ni ligand environment under H2 (g). These comparisons suggest significant OH and H2O content and, further, that different Ni ions within the cluster and/or NU-1000 structure may comprise subtly different numbers of these ligands. Further, the observation of significant H2O content under H2 (g) suggests that the NU-1000 support supplies H2O to the cluster. Examples of ligand environments that could lead to the observed PDFs are provided. The combination of simulations and experiments provides new insights into the ligand environment for Ni-NU-1000 catalysts that will be useful for understanding the ligand environments of other single-site Ni catalysts as well.

Effect of manganese substitution of ferrite nanoparticles on particle grain structure.

To investigate the influence of manganese substitution on the saturation magnetization of manganese ferrite nanoparticles, samples with various compositions (Mn x Fe3-x O4, x = 0, 0.25, 0.5, 0.75, and 1) were synthesized and characterized. The saturation magnetization of such materials was both calculated using density functional theory and measured via vibrating sample magnetometry. A discrepancy was found; the computational data demonstrated a positive correlation between manganese content and saturation magnetization, while the experimental data exhibited an inverse correlation. X-ray diffraction (XRD) and magnetometry results indicated that the crystallite diameter and the magnetic diameter decrease when adding more manganese, which could explain the loss of magnetization of the particles. For 20 nm nanoparticles, with increasing manganese substitution level, the crystallite size decreases from 10.9 nm to 6.3 nm and the magnetic diameter decreases from 15.1 nm to 3.5 nm. Further high resolution transmission electron microscopy (HRTEM) analysis confirmed the manganese substitution induced defects in the crystal lattice, which encourages us to find ways of eliminating crystalline defects to make more reliable ferrite nanoparticles.

Influence of an electrified interface on the entropy and energy of solvation of methanol oxidation intermediates on platinum(111) under explicit solvation.

Liquid water and electric fields play significant roles in phenomena occurring at catalytic and electrocatalytic interfaces; however, how their interplay influences interfacial energetics remains uncertain. Electric fields control the orientations of water molecules, which we hypothesized would influence the solvation thermodynamics of surface species. To explore this hypothesis, we used multiscale simulations involving density functional theory and classical molecular dynamics. We computed the energies and entropies of solvation of surface species on Pt(111), specifically, adsorbed CH3OH, COH, and CO, which are intermediates in the pathway of methanol oxidation, in the presence of electric fields spanning -0.5 to +0.5 V Å-1. We found that both the energy and entropy of solvation depend on the strength and direction of the field, with the entropy of solvation being significantly impacted. Both the energy and entropy dependence on the field can be ascribed to water molecule orientations. Specifically, more positive fields orient water molecules so that they can more effectively hydrogen bond with surface species, which strengthens the energies of solvation. However, at more negative fields, competition with the surface species causes interfacial water molecules to reorient, which leads to disorder in the water structure and hence increased entropy.

Site Densities, Rates, and Mechanism of Stable Ni/UiO-66 Ethylene Oligomerization Catalysts.

Nickel-functionalized UiO-66 metal organic frameworks (MOFs) oligomerize ethylene in the absence of cocatalysts or initiators after undergoing ethylene-pressure-dependent transients and maintain stable oligomerization rates for >15 days on stream. Higher ethylene pressures shorten induction periods and engender more active sites for ethylene oligomerization; these sites exhibit invariant selectivity-conversion characteristics to justify that only one type of catalytic center is relevant for oligomerization. The number of active sites is estimated using in situ NO titration to disambiguate the effect of increased reaction rates upon exposure to increasing ethylene pressures. After accounting for augmented site densities with increasing ethylene pressures, ethylene oligomerization is first order in ethylene pressure from 100 to 1800 kPa with an activation energy of 81 kJ mol-1 at temperatures from 443-503 K on Ni/UiO-66. A representative Ni/UiO-66 cluster model that mimics high ethylene pressure process conditions is validated with ab initio thermodynamic analysis, and the Cossee-Arlman mechanism is posited based on comparisons between experimental and computed activation enthalpies from density functional theory calculations on these cluster models of Ni/UiO-66. The insights gained from experiment and theory help rationalize evolution in structure and stability for ethylene oligomerization Ni/UiO-66 MOF catalysts.

A Method for Obtaining Liquid-Solid Adsorption Rates from Molecular Dynamics Simulations: Applied to Methanol on Pt(111) in H2O.

Adsorption is an important step in heterogeneous catalysis as it predetermines how many reactant molecules can participate in a surface reaction per unit time. While the rate of adsorption processes is well studied in gas-solid adsorption in both theory and experiment, such rates are still not well studied for liquid-solid adsorption. This is partly because the ever-changing configurations of liquid-phase solvent molecules impede the ability to study a molecule approaching a surface from a liquid phase by either experiment or theory. In this work, we develop a method using molecular dynamics (MD) simulations to study the rate of adsorption in liquid-solid adsorption processes. Specifically, we use MD to model the diffusion of a methanol molecule in aqueous solvent and its adsorption to a Pt(111) surface. We find that by approximating the solute motion as following the same displacement rates as a random walk model, the adsorbed and non-adsorbed states of the methanol molecule near the Pt(111) surface can be discerned and quantified. In particular, this methodology enables extracting a sticking coefficient and a macroscopically relatable adsorption rate. This method can be applied to arbitrary types of reactants and surfaces, as well as different liquid environments, thus providing a general tool for predicting quantitative adsorption rates of liquid-solid adsorption systems.

Insights into how the aqueous environment influences the kinetics and mechanisms of heterogeneously-catalyzed COH* and CH3OH* dehydrogenation reactions on Pt(111).

Water influences catalytic reactions in multiple ways, including energetic and mechanistic effects. While simulations have provided significant insight into the roles that H2O molecules play in aqueous-phase heterogeneous catalysis, questions still remain as to the extent to which H2O structures influence catalytic mechanisms. Specifically, influences of the configurational variability in the water structures at the catalyst interface are yet to be understood. Configurational variability is challenging to capture, as it requires multiscale approaches. Herein, we apply a multiscale sampling approach to calculate reaction thermodynamics and kinetics for COH* dehydrogenation to CO* and CH3OH* dehydrogenation to CH2OH* on Pt(111) catalysts under liquid H2O. We explore various pathways for these dehydrogenation reactions that could influence the overall mechanism of methanol decomposition by including participation of H2O structures both energetically and mechanistically. We find that the liquid H2O environment significantly influences the mechanism of COH* dehydrogenation to CO* but leaves the mechanism of CH3OH* dehydrogenation to CH2OH* largely unaltered.

Multiscale Sampling of a Heterogeneous Water/Metal Catalyst Interface using Density Functional Theory and Force-Field Molecular Dynamics.

A significant number of heterogeneously-catalyzed chemical processes occur under liquid conditions, but simulating catalyst function under such conditions is challenging when it is necessary to include the solvent molecules. The bond breaking and forming processes modeled in these systems necessitate the use of quantum chemical methods. Since molecules in the liquid phase are under constant thermal motion, simulations must also include configurational sampling. This means that multiple configurations of liquid molecules must be simulated for each catalytic species of interest. The goal of the protocol presented here is to generate and sample trajectories of configurations of liquid water molecules around catalytic species on flat transition metal surfaces in a way that balances chemical accuracy with computational expense. Specifically, force field molecular dynamics (FFMD) simulations are used to generate configurations of liquid molecules that can subsequently be used in quantum mechanics-based methods such as density functional theory or ab initio molecular dynamics. To illustrate this, in this manuscript, the protocol is used for catalytic intermediates that could be involved in the pathway for the decomposition of glycerol (C3H8O3). The structures that are generated using FFMD are modeled in DFT in order to estimate the enthalpies of solvation of the catalytic species and to identify how H2O molecules participate in catalytic decompositions.

Free Energies of Catalytic Species Adsorbed to Pt(111) Surfaces under Liquid Solvent Calculated Using Classical and Quantum Approaches.

Solvent plays an important role in liquid phase heterogeneous catalysis; however, methods for calculating the free energies of catalytic phenomena at the solid-liquid interface are not well-established. For example, solvent molecules alter the energies of catalytic species and participate in catalytic reactions and can thus significantly influence catalytic performance. In this work, we begin to establish methods for calculating the free energies of such phenomena, specifically, by employing an explicit solvation method using a multiscale sampling (MSS) approach. This MSS approach combines classical molecular dynamics with density functional theory. We use it to calculate the free energies of solvation of catalytic species, specifically adsorbed NH*, NH2*, CO*, COH*, CH2OH*, and C3H7O3* on Pt(111) surfaces under aqueous phase and under a mixed H2O/CH3OH solvent. We compare our calculated values with analogous values from implicit solvation for validation and to identify situations where implicit solvation is sufficient versus where explicit solvent is needed to compute adsorbate free energies. Our results indicate that explicit quantum-based methods are needed when adsorbates form chemical bonds and/or strong hydrogen bonds with H2O solvent. Using MSS, we further separate the calculated free energies into energetic and entropic contributions in order to understand how each influences the free energy. We find that adsorbates that exhibit strong energies also exhibit strong and negative entropies, and we attribute this relationship to hydrogen bonding between the adsorbates and the solvent molecules, which provides a large energetic contribution but reduces the overall mobility of the solvent.

Sinter-Resistant Platinum Catalyst Supported by Metal-Organic Framework.

Single atoms and few-atom clusters of platinum are uniformly installed on the zirconia nodes of a metal-organic framework (MOF) NU-1000 via targeted vapor-phase synthesis. The catalytic Pt clusters, site-isolated by organic linkers, are shown to exhibit high catalytic activity for ethylene hydrogenation while exhibiting resistance to sintering up to 200 °C. In situ IR spectroscopy reveals the presence of both single atoms and few-atom clusters that depend upon synthesis conditions. Operando X-ray absorption spectroscopy and X-ray pair distribution analyses reveal unique changes in chemical bonding environment and cluster size stability while on stream. Density functional theory calculations elucidate a favorable reaction pathway for ethylene hydrogenation with the novel catalyst. These results provide evidence that atomic layer deposition (ALD) in MOFs is a versatile approach to the rational synthesis of size-selected clusters, including noble metals, on a high surface area support.

Optimizing Open Iron Sites in Metal-Organic Frameworks for Ethane Oxidation: A First-Principles Study.

Activation of the C-H bonds in ethane to form ethanol is a highly desirable, yet challenging, reaction. Metal-organic frameworks (MOFs) with open Fe sites are promising candidates for catalyzing this reaction. One advantage of MOFs is their modular construction from inorganic nodes and organic linkers, allowing for flexible design and detailed control of properties. In this work, we studied a series of single-metal atom Fe model systems with ligands that are commonly used as MOF linkers and tried to understand how one can design an optimal Fe catalyst. We found linear relationships between the binding enthalpy of oxygen to the Fe sites and common descriptors for catalytic reactions, such as the Fe 3d energy levels in different reaction intermediates. We further analyzed the three highest-barrier steps in the ethane oxidation cycle (including desorption of the product) with the Fe 3d energy levels. Volcano relationships are revealed with peaks toward higher Fe 3d energy and stronger electron-donating group functionalization of linkers. Furthermore, we found that the Fe 3d energy levels positively correlate with the electron-donating strength of functional groups on the linkers. Finally, we validated our hypotheses on larger models of MOF-74 iron sites. Compared with MOF-74, functionalizing the MOF-74 linkers with NH2 groups lowers the enthalpic barrier for the most endothermic step in the reaction cycle. Our findings provide insight for catalyst optimization and point out directions for future experimental efforts.

Using degrees of rate control to improve selective n-butane oxidation over model MOF-encapsulated catalysts: sterically-constrained Ag3Pd(111).

Metal nanoparticles encapsulated within metal organic frameworks (MOFs) offer steric restrictions near the catalytic metal that can improve selectivity, much like in enzymes. A microkinetic model is developed for the regio-selective oxidation of n-butane to 1-butanol with O2 over a model for MOF-encapsulated bimetallic nanoparticles. The model consists of a Ag3Pd(111) surface decorated with a 2-atom-thick ring of (immobile) helium atoms which creates an artificial pore of similar size to that in common MOFs, which sterically constrains the adsorbed reaction intermediates. The kinetic parameters are based on energies calculated using density functional theory (DFT). The microkinetic model was analysed at 423 K to determine the dominant pathways and which species (adsorbed intermediates and transition states in the reaction mechanism) have energies that most sensitively affect the reaction rates to the different products, using degree-of-rate-control (DRC) analysis. This analysis revealed that activation of the C-H bond is assisted by adsorbed oxygen atoms, O*. Unfortunately, O* also abstracts H from adsorbed 1-butanol and butoxy as well, leading to butanal as the only significant product. This suggested to (1) add water to produce more OH*, thus inhibiting these undesired steps which produce OH*, and (2) eliminate most of the O2 pressure to reduce the O* coverage, thus also inhibiting these steps. Combined with increasing butane pressure, this dramatically improved the 1-butanol selectivity (from 0 to 95%) and the rate (to 2 molecules per site per s). Moreover, 40% less O2 was consumed per oxygen atom in the products. Under these conditions, a terminal H in butane is directly eliminated to the Pd site, and the resulting adsorbed butyl combines with OH* to give the desired 1-butanol. These results demonstrate that DRC analysis provides a powerful approach for optimizing catalytic process conditions, and that highly selectivity oxidation can sometimes be achieved by using a mixture of O2 and H2O as the oxidant. This was further demonstrated by DRC analysis of a second microkinetic model based on a related but hypothetical catalyst, where the activation energies for two of the steps were modified.

A modelling approach for MOF-encapsulated metal catalysts and application to n-butane oxidation.

Metal nanoparticles (NP) encapsulated by metal-organic frameworks (MOFs) are novel composite materials that have shown promise as regioselective catalysts. The regioselectivity in these materials arises from steric constraints imposed by the porous MOF structure, which limit the way molecules approach and interact with the metal surface. Here we introduce a conceptually simple DFT approach to model reactions under such steric constraints. This approach is computationally efficient and accounts for the steric constraints imposed by a MOF pore in a general way. The adsorption of reactants, intermediates, and products associated with oxidation of n-butane to 1-butanol (and 2-butanol) on clean and oxygen-covered palladium surfaces is investigated with (and without) the constraints of a pore. Reaction energies are calculated, and we find that the thermodynamic favorability of the intermediate reactions is affected by the presence of steric constraints, oxygen coverage, and the exposed crystal surface of the metal. Based on these results, the Pd(111) surface with 0.25 ML oxygen coverage and steric constraints (which could be provided by a suitable MOF) seems promising to favor the desired sequence of reactions that would lead to the conversion of n-butane to 1-butanol.

Stepwise adsorption in a mesoporous metal-organic framework: experimental and computational analysis.

Stepwise adsorption in a metal-organic framework with both micro- and meso-pores is caused by adsorbates first filling the micropores, then adsorbing along the mesopore walls, and finally filling the mesopores.