Search for: All records

Atomic-scale imaging reveals highly enantiospecific restructuring that is dependent on both the enantiomer of tartaric acid and surface chirality. 

Many computational studies of catalytic surface reaction kinetics have demonstrated the existence of linear scaling relationships between physical descriptors of catalysts and reaction barriers on their surfaces. In this work, the relationship between catalyst activity, electronic structure, and alloy composition was investigated experimentally using a AgxPd1−x Composition Spread Alloy Film (CSAF) and a multichannel reactor array that allows measurement of steady-state reaction kinetics at 100 alloy compositions simultaneously. Steady-state H2 −D2 exchange kinetics were measured at atmospheric pressure on AgxPd1−x catalysts over a temperature range of 333−593 K and a range of inlet H2 and D2 partial pressures. X-ray photoelectron spectroscopy (XPS) was used to characterize the CSAF by determining the local surface compositions and the valence band electronic structure at each composition. The valence band photoemission spectra showed that the average energy of the valence band, ε̅v, shifts linearly with composition from −6.2 eV for pure Ag to −3.4 eV for pure Pd. At all reaction conditions, the H2 −D2 exchange activity was found to be highest on pure Pd and gradually decreased as the alloy was diluted with Ag until no activity was observed for compositions with xPd < 0.58. Measured H2 −D2 exchange rates across the CSAF were fit using the Dual Subsurface Hydrogen (2H′) mechanism to extract estimates for the activation energy barriers to dissociative adsorption, ΔEads ‡ , associative desorption, ΔEdes ‡ , and the surface-to-subsurface diffusion energy, ΔEss, as a function of alloy composition, xPd. The 2H′ mechanism predicts ΔEads ‡ = 0−10 kJ/mol, ΔEdes ‡ = 30−65 kJ/mol, and ΔEss = 20−30 kJ/mol for all alloy compositions with xPd ≥ 0.64, including for the pure Pd catalyst (i.e., xPd = 1). For these Pd-rich catalysts, ΔEdes ‡ and ΔEss appeared to increase by ∼5 kJ/mol with decreasing xPd. However, due to the coupling of kinetic parameters in the 2H′ mechanism, we are unable to exclude the possibility that the kinetic parameters predicted when xPd ≥ 0.64 are identical to those predicted for pure Pd. This suggests that H2 −D2 exchange occurs only on bulk-like Pd domains, presumably due to the strong interactions between H2 and Pd. In this case, the decrease in catalytic activity with decreasing xPd can be explained by a reduction in the availability of surface Pd at high Ag compositions. 

Abstract Surface segregation, whereby the surface composition of an alloy differs systematically from the bulk, has historically been hard to study, because it requires experimental and modeling methods that span alloy composition space. In this work, we study surface segregation in catalytically relevant noble and platinum‐group metal alloys with a focus on three ternary systems: AgAuCu, AuCuPd, and CuPdPt. We develop a data set of 2478 fcc slabs with those compositions including all three low‐index crystallographic orientations relaxed with Density Functional Theory using the PBEsol functional with D3 dispersion corrections. We fine‐tune a machine learning model on this data and use the model in a series of 1800 Monte Carlo simulations spanning ternary composition space for each surface orientation and ternary chemical system. The results of these simulations are validated against prior experimental surface segregation data collected using composition spread alloy films for AgAuCu and AuCuPd. Our findings reveal that simulations conducted using the (110) orientation most closely match experimentally observed surface segregation trends, and while predicted trends qualitatively match observation, biases in the PBEsol functional limit numeric accuracy. This study advances understanding of surface segregation and the utility of computational studies and highlights the need for further improvements in simulation accuracy. 

The 2D Ising model is well-formulated to address problems in adsorption thermodynamics. It is particularly well-suited to describing the adsorption isotherms predicting the surface enantiomeric excess, ees, observed during competitive co-adsorption of enantiomers onto achiral surfaces. Herein, we make the direct one-to-one correspondence between the 2D Ising model Hamiltonian and the Hamiltonian used to describe competitive enantiomer adsorption on achiral surfaces. We then demonstrate that adsorption from racemic mixtures of enantiomers and adsorption of prochiral molecules are directly analogous to the Ising model with no applied magnetic field, i.e., the enantiomeric excess on chiral surfaces can be predicted using Onsager’s solution to the 2D Ising model. The implication is that enantiomeric purity on the surface can be achieved during equilibrium exposure of prochiral compounds or racemic mixtures of enantiomers to achiral surfaces. 

Comprehensive mapping of enantiospecific surface reactivity versus the crystallographic orientation of Cu( hkl ) surfaces vicinal to Cu(111) has been conducted using a spherically shaped single crystal on which the surface normal vectors, [ hkl ], span all possible orientations lying with 14° of the [111] direction. This has allowed direct measurement on 169 different Cu( hkl ) surfaces of the two rate constants, k (hkl)i and k (hkl)e, that determine the kinetics of the vacancy-mediated, explosive decomposition of tartaric acid (TA). The initiation rate constant, k (hkl)i, quantifies the kinetics of an initiation step that creates vacancies in the adsorbed TA monolayer. The explosion rate constant, k (hkl)e, quantifies the kinetics of a vacancy-mediated explosion step that results in TA decomposition and product desorption. Enantiospecificity is revealed by the dependence of TA decomposition kinetics on the chirality of the local surface orientation. Diastereomerism is demonstrated by the fact that d -TA is more reactive than l -TA on S surfaces while l -TA is more reactive on R surfaces. The time to reach half coverage, t (hkl)1/2, during isothermal TA decomposition at 433 K allowed determination of the most enantiospecific surface orientation; Cu(754). The ideal Cu(754) surface structure consists of (111) terraces separated by monoatomic steps formed by the (100) and (110) microfacets. 

Kinetic parameters have been estimated for the H2– D2 exchange reaction on a thin film Pd catalyst by fitting reaction data from T = 333 to 593 K over a range of inlet partial pressures, Pin H2 and Pin D2 . A rigorous approach to estimating the 95% confidence regions of the kinetic parameters reveals some of the issues and complexities that are not routinely considered in the estimation of kinetic parameter uncertainty from catalytic data. Three different mechanistic models were used to assess the influence of subsurface hydrogen, H′: the traditional Langmuir–Hinshelwood (LH) mechanism, the Single Subsurface Hydrogen (1H′) mechanism, and the Dual Subsurface Hydrogen (2H′) mechanism. The fitting was performed by fixing the preexponential factors for all Arrhenius rate constants and equilibrium constants to their transition state theory values. The diffusion of H and D atoms from the surface into the subsurface was constrained to be endothermic (i.e. ΔE ss > 0) and represented as an equilibrium process. Performance of the fitting routine was evaluated on a noiseless simulated dataset (created using ΔE‡ ads = 0, ΔE‡ des = 43, and ΔE ss = 25 kJ/mol) and the same simulated dataset with the inclusion of 3% Gaussian noise. In both cases, the solver was able to return the chosen values of ΔE‡ ads , ΔE‡ des , and ΔE ss . Mapping of the behavior of the residual sum of squared errors, 2 , about its global minimum within 3D ( ads , des , ss ) parameter space allowed quantification and visualization of the 95% confidence regions using 2D error ellipses for each pair of fitting parameters. For the experimental dataset on the Pd catalyst, fitting to the LH model predicted that H2– D2 exchange is adsorption rate limited, with ΔE‡ ads = 51.1 ± 0.6 kJ/mol with 95% confidence. On the other hand, fitting to both the 1H′ and 2H′ models led to predictions of ΔE‡ ads = 0, consistent with the current understanding that the barrier to H2 dissociation on Pd is low. Thus, the results detailed herein provide supporting evidence for a non-LH mechanism for H2– D2 exchange on Pd while also illustrating the issues associated with quantification of uncertainty in kinetic parameter estimation.