Ab initio study of O2 precursor states on the Pd(111) surface

Accelerating the prediction of inorganic surfaces with machine learning interatomic potentials

The surface properties of solid-state materials often dictate their functionality, especially for applications where nanoscale effects become important. The relevant surface(s) and their properties are determined, in large part, by the material's synthesis or operating conditions. These conditions dictate thermodynamic driving forces and kinetic rates responsible for yielding the observed surface structure and morphology. Computational surface science methods have long been applied to connect thermochemical conditions to surface phase stability, particularly in the heterogeneous catalysis and thin film growth communities. This review provides a brief introduction to first-principles approaches to compute surface phase diagrams before introducing emerging data-driven approaches. The remainder of the review focuses on the application of machine learning, predominantly in the form of learned interatomic potentials, to study complex surfaces. As machine learning algorithms and large datasets on which to train them become more commonplace in materials science, computational methods are poised to become even more predictive and powerful for modeling the complexities of inorganic surfaces at the nanoscale.

A Vacuum Ultraviolet-Enhanced Oxidation Mechanism for Pd: Near-Surface Oxidation for Atomic Layer Etching

Density functional theory (DFT) is used to better understand the oxidation of Pd metal using vacuum ultraviolet (VUV) light co-exposed with O₂, which is known to produce O and O₃. The oxidation of Pd metal arising from O, O₂, and O₃ is assessed on bare Pd, Pd with a 0.25 monolayer of adsorbed atomic O, and Pd with increasing O incorporation into the substrate. DFT calculations are complemented experimentally by co-exposing 20 nm Pd films to 1 Torr of O₂ and VUV photons (6.5 < hν < 11.3 eV) from a D₂ lamp at temperatures ranging from 50 to 200 °C and times from 30 s to 40 min. Oxidation of Pd is characterized using in situ X-ray photoelectron spectroscopy. Co-exposures at 50 °C and 1 Torr O₂ are performed with the Pd illuminated by the VUV light and shadowed from the VUV light in attempting to select for the oxidant that impinges on the Pd surface and causes oxidation. Results suggest that atomic O incident from the gas phase is responsible for oxidation of Pd, as no PdO_x formation is observed for the same time period with the sample shadowed. Growth of PdO_x via O diffusion is studied with the nudged elastic band method. Atomic O diffusion through Pd has an activation energy barrier of ∼2.87 eV with respect to a surface O. This decreases to ∼1.80 eV once the 0.25 monolayer of O occupies the surface. The extent of Pd oxidation is limited to the near-surface Pd region for all times and temperatures investigated. PdO_x formation does not appear to exceed one to two atomic layers of Pd for conditions explored herein.

Superoxide Decay Pathways in Oxygen Reduction Reaction on Carbon-Based Catalysts Evidenced by Theoretical Calculations

Superoxide decay behavior is evaluated by contrastive potential energy surface analyses from theoretical calculations. The presence of carbon surfaces enhances the superoxide decay through disproportionation instead of electroreduction, which is a non-electrochemical reaction that impedes the increase in energy efficiency for relevant energy conversion applications. Nitrogen doping of carbon effectively retards the disproportionation by influencing the oxygen stretching vibration and changing the proton-coupled electron transfer trend.

Highly efficient Pd-based core–shell nanowire catalysts for O2 dissociation

With the consideration of the stability and cost, we found that the Fe, Co, Ni, Cu, Ru, Ir atoms have lower price than the Pd and favor at the core even with O adatom at the surface. The formed M@Pd core–shell nanowires are active for O₂ dissociation with activation barriers no larger than 0.25 eV.

Surface structural characteristics and tunable electronic properties of wet-chemically prepared Pd nanoparticles

The ligand substitution reaction, Pd L(3,2,1)-edge and S K-edge x-ray absorption fine structure (XAFS), XAFS simulations, and valence-band and core-level x-ray photoelectron spectroscopy (XPS) have been used to systematically study the surface chemical and electronic properties of wet-chemically prepared Pd nanoparticles of varied size, molecular capping, and metal composition. It was found that the replacement of weakly interacting capping molecules (amine and tetra-alkylphosphonium bromide) with strongly binding thiols caused a considerable change in the surface bonding of Pd nanoparticles. However, the Pd d-electron counts (number of d electrons) remained almost unchanged before and after ligand substitution, which is unexpected since Pd atoms normally lose electrons to the more electronegative S atoms. XAFS results and simulations provided useful insights into the surface structural characteristics of Pd nanoparticles and satisfactorily accounted for the unexpected d-electron behavior involved in the ligand substitution process. XPS valence and core-level spectra further revealed a size-dependent d-band narrowing and presented complementary information to XAFS about the surface electronic properties of Pd atoms. The small weakly bound Pd nanoparticles seem inevitably to have a net d-electron depletion due to the influence of the surface effect (chemical adsorption by oxygen), which is more significant than the d-electron enriching nanosize effect. However, it was demonstrated that by forming Pd-Ag alloy nanoparticles, a net increase of the Pd d-electron counts can be realized. Therefore, it is illustrated that by manipulating the surface, size, and alloying effects, the electronic properties of Pd nanoparticles can be possibly tuned.

Adsorption-state-dependent subpicosecond photoinduced desorption dynamics

Femtosecond laser excitation has been used to initiate desorption of molecular oxygen from the (111) surface of Pd and to study the adsorption-state dependence of the substrate-adsorbate coupling. The relative populations of the two chemical states, peroxo (O2(2-)) and superoxo (O2-), were varied by changing the total coverage. Two-pulse correlation measurements exhibit a dominant 400 fs response and a slower 10 ps decay that are relatively independent of the initial O2 coverage. In contrast, the photodesorption yield and the nonlinearity of the fluence dependence show a systematic coverage dependence. The coverage-independent subpicosecond response indicates that the photoinduced desorption from the two states is driven primarily by the same electron-mediated mechanism, while the coverage dependence of the yield indicates that the desorption efficiency from the superoxo state is greater than that from the peroxo state. These results are discussed in the context of the electron-phonon two-temperature model with an empirical adsorbate-electron frictional coupling that depends on both the electronic temperature and the activation energy for desorption. With a coupling strength that decreases as the activation energy decreases, the trends with varying coverage, absorbed fluence, and time delay can all be reproduced. The model is consistent with a transition from a resonantly enhanced (diabatic) regime to an adiabatic regime as the system relaxes, accounting for the biexponential correlation behavior.

In Situ Extended X-ray Absorption Fine Structure Study during Selective Alcohol Oxidation over Pd/Al2O3 in Supercritical Carbon Dioxide

High-pressure in situ X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) data are reported during the selective oxidation of benzyl alcohol to benzaldehyde in supercritical carbon dioxide over a Pd/Al(2)O(3) catalyst (shell impregnated). For this purpose, a continuous-flow system with a spectroscopic cell suitable for in situ X-ray absorption studies on heterogeneous catalysts up to 200 degrees C and 200 bar has been developed. Due to the high contribution of the dense fluid to the overall X-ray absorption, high stability of the process pressure is mandatory, particularly when recording EXAFS spectra. According to EXAFS and XANES results, the palladium particles were fully reduced after exposure to benzyl alcohol in scCO(2). In contrast to Pd-catalyzed liquid-phase oxidation, a higher oxygen tolerance of the catalyst was observed. Palladium was partially oxidized on the surface under typical reaction conditions (0.9 mol % benzyl alcohol/0.5 mol % O(2) in carbon dioxide), which gradually increased when the concentration of oxygen in the feed was raised. Both XANES and EXAFS data uncovered that palladium is mainly oxidized on the surface or within the outermost layers. These results are in accordance with simulations of the XANES data using the FEFF8.20 code (program for ab initio calculations on multiple scattering XAS) and EXAFS data fitting/simulation.

C–O bond scission of methoxide on Pd nanoparticles: A density functional study

C-O bond scission of methoxide species adsorbed at the surface of Pd nanoparticle was studied by DF calculations for the example of cuboctahedral Pd(79). To investigate different locations of adsorbed intermediates as well as the transition state of C-O bond scission, a substrate model was used, which allows one to consider adsorbates without any local geometry restrictions. In contrast to reaction sites on the flat Pd(111) surface and on extended facets, scission of the C-O bond of methoxide at cluster edges is exothermic by approximately 40 kJ mol(-1) and the decomposition product CH(3) is found to be stabilized there. However, the high calculated activation barrier, approximately 140 mol(-1), implies only a very slow reaction compared to dehydrogenation of CH(3)O.

Relationship between Kinetic and Thermodynamic Characteristics of Oxygen Dissociative Adsorption on Close-Packed Metal Surfaces

The activation energies and rate constants characterizing dissociative oxygen adsorption on transition and noble metal surfaces are calculated using a previously developed analytical formalism. Three approaches are considered: classical and quantum nonadiabatic approaches as well as a classical adiabatic approach. The results show a nonlinear relationship between the activation energies and the corresponding reaction energies. The effects due to the quantum description of the O-O bond and due to the Condon approximation are discussed. Calculated activation energies are compared with published experimental data.

Transient and steady state CO oxidation kinetics on nanolithographically prepared supported Pd model catalysts: Experiments and simulations

Applying molecular-beam methods to a nanolithographically prepared planar PdSiO2 model catalyst, we have performed a detailed study of the kinetics of CO oxidation. The model catalyst was prepared by electron-beam lithography, allowing individual control of particle size and position. The sample was structurally characterized by atomic force microscopy and scanning electron microscopy before and after reaction. In the kinetic experiments, the O-rich and CO-rich regimes were investigated systematically with respect to their transient and steady-state behaviors, both under bistable and monostable reaction conditions. Separate molecular beams were used in order to supply the reactants, allowing individual control over the reactant fluxes. The desorbing CO2 was detected by both angle-resolved and angle-integrated mass spectrometries. The experimental results were analyzed using different types of microkinetic models, including a detailed reaction-diffusion model, which takes into account the structural parameters of the catalyst as well as scattering of the reactants from the support. The model quantitatively reproduces the results as a function of the reactant fluxes and the surface temperature. Various kinetic effects observed are discussed in detail on the basis of the model. Specifically, it is shown that under conditions of limited oxygen mobility, the switching behavior between the kinetic regimes is largely driven by the surface mobility of CO.

Local reaction rates and surface diffusion on nanolithographically prepared model catalysts: Experiments and simulations

Combining molecular beam methods and angular resolved mass spectrometry, we have studied the angular distribution of desorbing products during CO oxidation on a planar Pd/silica supported model catalyst. The model catalyst was prepared by means of electron beam lithography, allowing individual control of particle size, position, and aspect ratio, and was characterized by atomic force microscopy and scanning electron microscopy before and after reaction. In the experiment, both oxygen and CO rich regimes were investigated using separate molecular beams for the two reactants. This allows exploration of diffusion effects of reactants on the particles and of shadowing and backscattering phenomena. A reaction-diffusion model was developed in order to extract information about local reaction rates on the surface of the catalyst nanoparticles. The model takes into account the structural parameters of the catalyst as well as the backscattering of the reactants and products from the support. It allows a quantitative description of the experimental data and provides a detailed understanding of temperature and reactant flux dependent effects. Moreover, information on the surface mobility of oxygen under steady-state reaction conditions could be obtained by comparison with the experimental results.

CH3O decomposition on PdZn(111), Pd(111), and Cu(111). A theoretical study

Methanol steam re-forming, catalyzed by Pd/ZnO, is a potential hydrogen source for fuel cells, in particular in pollution-free vehicles. To contribute to the understanding of pertinent reaction mechanisms, density functional slab model studies on two competing decomposition pathways of adsorbed methoxide (CH(3)O) have been carried out, namely, dehydrogenation to formaldehyde and C-O bond breaking to methyl. For the (111) surfaces of Pd, Cu, and 1:1 Pd-Zn alloy, adsorption complexes of various reactants, intermediates, transition states, and products relevant for the decomposition processes were computationally characterized. On the surface of Pd-Zn alloy, H and all studied C-bound species were found to prefer sites with a majority of Pd atoms, whereas O-bound congeners tend to be located on sites with a majority of Zn atoms. Compared to Pd(111), the adsorption energy of O-bound species was calculated to be larger on PdZn(111), whereas C-bound moieties were less strongly adsorbed. C-H scission of CH(3)O on various substrates under study was demonstrated to proceed easier than C-O bond breaking. The energy barrier for the dehydrogenation of CH(3)O on PdZn(111) (113 kJ mol(-)(1)) and Cu(111) (112 kJ mol(-)(1)) is about 4 times as high as that on Pd(111), due to the fact that CH(3)O interacts more weakly with Pd than with PdZn and Cu surfaces. Calculated results showed that the decomposition of methoxide to formaldehyde is thermodynamically favored on Pd(111), but it is an endothermic process on PdZn(111) and Cu(111) surfaces.