The structure and chemistry of CH3 and CH radicals adsorbed on Ni(111)

Plasma-catalysis for VOCs decomposition: A review on micro- and macroscopic modeling

Plasma-catalysis has been recognized as a promising method to decompose hazardous volatile organic compounds (VOCs) since many years ago. To understand the fundamental mechanisms of VOCs decomposition by plasma-catalysis systems, both experimental and modeling studies have been extensively carried out. However, literature on summarized modeling methodologies is still scarce. In this short review, we therefore present a comprehensive overview of modeling methodologies ranging from microscopic to macroscopic modeling in plasma-catalysis for VOCs decomposition. The modeling methods of VOCs decomposition by plasma and plasma-catalysis are classified and summarized. The roles of plasma and plasma-catalyst interactions in VOCs decomposition are also critically examined. Taking the current advances in understanding the decomposition mechanisms of VOCs into account, we finally provide our perspectives for future research directions. This short review aims to stimulate the further development of plasma-catalysis for VOCs decomposition in both fundamental studies and practical applications with advanced modeling methods.

Formation of Graphene on Gold-Nickel Surface Alloys

Nickel (Ni)-catalyzed growth of a single- or rotated-graphene layer is a well-established process above 800 K. In this report, a Au-catalyzed, low-temperature, and facile route at 500 K for graphene formation is described. The substantially lower temperature is enabled by the presence of a surface alloy of Au atoms embedded within Ni(111), which catalyzes the outward segregation of carbon atoms buried in the Ni bulk at temperatures as low as 400-450 K. The resulting surface-bound carbon in turn coalesces into graphene above 450-500 K. Control experiments on a Ni(111) surface show no evidence of carbon segregation or graphene formation at these temperatures. Graphene is identified by its out-of-plane optical phonon mode at 750 cm^-1 and its longitudinal/transverse optical phonon modes at 1470 cm^-1 while surface carbon is identified by its C-Ni stretch mode at 540 cm^-1, as probed by high-resolution electron energy-loss spectroscopy. Dispersion measurements of the phonon modes confirm the presence of graphene. Graphene formation is observed to be maximum at 0.4 ML Au coverage. The results of these systematic molecular-level investigations open the door to graphene synthesis at the low temperatures required for integration with complementary metal-oxide-semiconductor processes.

Automated exploitation of the big configuration space of large adsorbates on transition metals reveals chemistry feasibility

Mechanistic understanding of large molecule conversion and the discovery of suitable heterogeneous catalysts have been lagging due to the combinatorial inventory of intermediates and the inability of humans to enumerate all structures. Here, we introduce an automated framework to predict stable configurations on transition metal surfaces and demonstrate its validity for adsorbates with up to 6 carbon and oxygen atoms on 11 metals, enabling the exploration of ~10⁸ potential configurations. It combines a graph enumeration platform, force field, multi-fidelity DFT calculations, and first-principles trained machine learning. Clusters in the data reveal groups of catalysts stabilizing different structures and expose selective catalysts for showcase transformations, such as the ethylene epoxidation on Ag and Cu and the lack of C-C scission chemistry on Au. Deviations from the commonly assumed atom valency rule of small adsorbates are also manifested. This library can be leveraged to identify catalysts for converting large molecules computationally.

The discovery of heterogeneous catalysts for large molecule conversion has been lagging due to the combinatorial inventory of intermediates. Here, the author presents an automated framework to explore the chemical space of reaction intermediates.

Infrared Activities of Adsorbed Species on Metal Surfaces: The Puzzle of Adsorbed Methyl (CH3)

Reflection-absorption infrared spectroscopy (RAIRS) is widely used to identify molecular adsorbates on metals during surface chemical reactions, but the interpretation of RAIRS data can be difficult with experiment alone. Here, we reveal from first-principles calculations the origin of the contrasting RAIRS spectra of methyl adsorbed on Pt(111) and Ni(111). We find that the dynamic dipole associated with the symmetric C-H stretch vibration of CH₃ along surface normal is significant on Pt(111) but negligibly small on Ni(111), explaining the strong IR activity in the former and the absence of any RAIRS peaks in the latter. This difference is correlated to different charge transfer patterns between metals and the adsorbate, which are determined by the different preferred adsorption sites of methyl on the two surfaces. This work highlights the need of electronic structure calculations in interpreting RAIRS spectra of adsorbates on metal surfaces.

Mechanistic insight into carbon-carbon bond formation on cobalt under simulated Fischer-Tropsch synthesis conditions

Facile C-C bond formation is essential to the formation of long hydrocarbon chains in Fischer-Tropsch synthesis. Various chain growth mechanisms have been proposed previously, but spectroscopic identification of surface intermediates involved in C-C bond formation is scarce. We here show that the high CO coverage typical of Fischer-Tropsch synthesis affects the reaction pathways of C₂H_x adsorbates on a Co(0001) model catalyst and promote C-C bond formation. In-situ high resolution x-ray photoelectron spectroscopy shows that a high CO coverage promotes transformation of C₂H_x adsorbates into the ethylidyne form, which subsequently dimerizes to 2-butyne. The observed reaction sequence provides a mechanistic explanation for CO-induced ethylene dimerization on supported cobalt catalysts. For Fischer-Tropsch synthesis we propose that C-C bond formation on the close-packed terraces of a cobalt nanoparticle occurs via methylidyne (CH) insertion into long chain alkylidyne intermediates, the latter being stabilized by the high surface coverage under reaction conditions.

The mechanism by which C-C bonds form during Fischer-Tropsch synthesis remains debated while spectroscopic identification of reaction intermediates remains scarce. Here, the authors identify alkylidynes as reactive intermediates for C-C bond formation on cobalt terrace sites and moreover show that these intermediates are stabilized by the high surface coverage typical for Fischer-Tropsch synthesis.

DFT studies of hydrocarbon combustion on metal surfaces

Catalytic combustion of hydrocarbons is an important technology to produce energy. Compared to conventional flame combustion, the catalyst enables this process to operate at lower temperatures; hence, reducing the energy required for efficient combustion. The reaction and activation energies of direct combustion of hydrocarbons (CH → C + H) on a series of metal surfaces were investigated using density functional theory (DFT). The data obtained for the Ag, Au, Al, Cu, Rh, Pt, and Pd surfaces were used to investigate the validity of the Brønsted-Evans-Polanyi (BEP) and transition state scaling (TSS) relations for this reaction on these surfaces. These relations were found to be valid (R² = 0.94 for the BEP correlation and R² = 1.0 for the TSS correlation) and were therefore used to estimate the energetics of the combustion reaction on Ni, Co, and Fe surfaces. It was found that the estimated transition state and activation energies (E_TS = -69.70 eV and E_a = 1.20 eV for Ni, E_TS = -87.93 eV and E_a = 1.08 eV for Co and E_TS = -92.45 eV and E_a = 0.83 eV for Fe) are in agreement with those obtained by DFT calculations (E_TS = -69.98 eV and E_a = 1.23 eV for Ni, E_TS = -87.88 eV and E_a = 1.08 eV for Co and E_TS = -92.57 eV and E_a = 0.79 eV for Fe). Therefore, these relations can be used to predict energetics of this reaction on these surfaces without doing the time consuming transition state calculations. Also, the calculations show that the activation barrier for CH dissociation decreases in the order Ag ˃ Au ˃ Al ˃ Cu ˃ Pt ˃ Pd ˃ Ni > Co > Rh > Fe.

Understanding FTS selectivity: the crucial role of surface hydrogen

Monomeric forms of carbon play a central role in the synthesis of long chain hydrocarbons via the Fischer-Tropsch synthesis (FTS). We explored the chemistry of C₁H_xad species on the close-packed surface of cobalt. Our findings on this simple model catalyst highlight the important role of surface hydrogen and vacant sites for product selectivity. We furthermore find that CO_ad affects hydrogen in multiple ways. It limits the adsorption capacity for H_ad, lowers its adsorption energy and inhibits dissociative H₂ adsorption. We discuss how these findings, extrapolated to pressures and temperatures used in applied FTS, can provide insights into the correlation between partial pressure of reactants and product selectivity. By combining the C₁H_x stability differences found in the present work with literature reports of the reactivity of C₁H_x species measured by steady state isotope transient kinetic analysis, we aim to shed light on the nature of the atomic carbon reservoir found in these studies.

Improved catalytic activity of rhodium monolayer modified nickel (110) surface for the methane dehydrogenation reaction: a first-principles study

The catalytic activity of pure Ni (110) and single Rh layer deposited Ni (110) surface for the complete dehydrogenation of methane is theoretically investigated by means of gradient-corrected periodic density functional theory. A detailed kinetic study, based on the analysis of the optimal reaction pathway for the transformation of CH4 to C and H through four elementary steps (CH4 → CH3 + H; CH3 → CH2 + H; CH2 → CH + H; CH → C + H) is presented for pure Ni (110) and Rh/Ni (110) surfaces and compared with pure Rh (110) surface. Through systematic examination of adsorbed geometries and transition states, we show that single layer deposition of Rh on Ni (110) surface has a striking influence on lowering the activation energy barrier of the dehydrogenation reaction. Moreover, it is found that a pure Ni (110) surface has a tendency for carbon deposition on the catalytic surface during the methane dissociation reaction which decreases the stability of the catalyst. However, the deposition of carbon is largely suppressed by the addition of a Rh overlayer on the pure Ni (110) surface. The physical origin of stronger chemisorption of carbon on Ni (110) relative to Rh/Ni (110) has been elucidated by getting insight into the electronic structures and d-band model of the catalytic surfaces. Considering the balance in both the catalytic activity as well as the catalyst stability, we propose that the Rh/Ni (110) surface possesses much improved catalytic property compared to pure Ni (110) and pure Rh (110) surfaces.

Temperature effects on adsorption and diffusion dynamics of CH3CH2(ads) and H3C-C≡C(ads) on Ag(111) surface and their self-coupling reactions: ab initio molecular dynamics approach

Density functional theory (DFT)-based molecular dynamics (DFTMD) simulations in combination with a Fourier transform of dipole moment autocorrelation function are performed to investigate the adsorption dynamics and the reaction mechanisms of self-coupling reactions of both acetylide (H3C-C(β)≡C(α) (ads)) and ethyl (H3C(β)-C(α)H2(ads)) with I(ads) coadsorbed on the Ag(111) surface at various temperatures. In addition, the calculated infrared spectra of H3C-C(β)≡C(α)(ads) and I coadsorbed on the Ag(111) surface indicate that the active peaks of -C(β)≡C(α)- stretching are gradually merged into one peak as a result of the dominant motion of the stand-up -C-C(β)≡C(α)- axis as the temperature increases from 200 K to 400 K. However, the calculated infrared spectra of H3C(β)-C(α)H2(ads) and I coadsorbed on the Ag(111) surface indicate that all the active peaks are not altered as the temperature increases from 100 K to 150 K because only one orientation of H3C(β)-C(α)H2(ads) adsorbed on the Ag(111) surface has been observed. These calculated IR spectra are in a good agreement with experimental reflection absorption infrared spectroscopy results. Furthermore, the dynamics behaviors of H3C-C(β)≡C(α)(ads) and I coadsorbed on the Ag(111) surface point out the less diffusive ability of H3C-C(β)≡C(α)(ads) due to the increasing s-character of Cα leading to the stronger Ag-Cα bond in comparison with that of H3C(β)-C(α)H2(ads) and I coadsorbed on the same surface. Finally, these DFTMD simulation results allow us to predict the energetically more favourable reaction pathways for self-coupling of both H3C-C(β)≡C(α)(ads) and H3C(β)-C(α)H2(ads) adsorbed on the Ag(111) surface to form 2,4-hexadiyne (H3C-C≡C-C≡C-CH3(g)) and butane (CH3-CH2-CH2-CH3(g)), respectively. The calculated reaction energy barriers for both H3C-C≡C-C≡C-CH3(g) (1.34 eV) and CH3-CH2-CH2-CH3(g) (0.60 eV) are further employed with the Redhead analysis to estimate the desorption temperatures approximately at 510 K and 230 K, respectively, which are in a good agreement with the experimental low-coverage temperature programmed reaction spectroscopy measurements.

CO2 reforming of CH4 on Ni(111): a density functional theory calculation

CO(2) reforming of CH(4) on Ni(111) was investigated by using density functional theory. On the basis of thermodynamic analyses, the first step is CH(4) sequential dissociation into surface CH (CH(4) --> CH(3) --> CH(2) --> CH) and hydrogen, and CO(2) dissociation into surface CO and O (CO(2) --> CO + O). The second step is CH oxygenation into CHO (CH + O --> CHO), which is more favored than dissociation into C and hydrogen (CH --> C + H). The third step is the dissociation of CHO into surface CO and H (CHO --> CO + H). This can explain the enhanced selectivity toward the formation of CO and H(2) on Ni catalysts. It is found that surface carbon formation by the Bouduard back reaction (2CO = C((ads)) + CO(2)) is more favored than by CH(4) sequential dehydrogenation. The major problem of CO(2) reforming of CH(4) is the very strong CO adsorption on Ni(111), which results in the accumulation of CO on the surface and hinders the subsequent reactions and promotes carbon deposition. Therefore, promoting CO desorption should maintain the reactivity and stability of Ni catalysts. The computed energy barriers of the most favorable elementary reaction identify the CH(4) activation into CH(3) and H as the rate-determining step of CO(2) reforming of CH(4) on Ni(111), in agreement with the isotopic experimental results.

A detailed analysis of vibrational excitations in x-ray photoelectron spectra of adsorbed small hydrocarbons

The vibrational fine structure of x-ray photoelectron (XP) spectra of a number of different small hydrocarbon molecules and reaction intermediates adsorbed on Pt(111) and Ni(111) has been investigated in detail. The data for methyl, methylidyne, acetylene, and ethylene can consistently be analyzed within the linear coupling model. The S factor, i.e., the intensity ratio of the first vibrationally excited to the adiabatic transition, is obtained to be 0.17+/-0.02 per C-H bond; for the deuterated species a value of 0.23+/-0.02 is obtained. Therefore, the vibrational fine structure can be used for fingerprinting in the analysis of XP spectra and for identifying unknown reaction intermediates. From the data, Deltar, the change of the minimum in the potential energy curve upon core ionization, is calculated within the linear coupling model using a first order correction. For all adsorbates, including the deuterated ones, a value of Deltar=0.060+/-0.004 A is obtained. Furthermore, from the binding energy of the adiabatic peak and from the energy of the vibrational excitation in the ionic final state some information on the adsorbate/substrate bond and the adsorption site can be derived.

The relationship between adsorption energies of methyl on metals and the metallic electronic properties: A first-principles DFT study

A theoretical study of CH3 adsorbed on the (111) surface of some transition and noble metal surfaces M (M = Cu, Ni, Rh, Pt, Pd, Ag, Au) and on the Fe(100) is presented. We find that the hollow site is preferred more than the top one for Fe, Ni, Rh, and Cu, but it is the other way for Pt, Pd, Au, and Ag. In addition, a good linear relationship was observed between the chemisorption energy and d-band center for Group VIII metals or the square of the coupling matrix element for Group IB metals at the hollow site. Interestingly, with a detailed comparison of the adsorption energies at the top and hollow sites, we find that the adsorption energies among each group are very similar on the top site, which supports the theoretical model of Hammer and Norskov that the coupling between the HOMO of adsorbate and sp states of the metal is dominant and almost equal, and that the second coupling to the d-band contributes less but reflects the change of the adsorption energy. It confirms that the coupling to the d band comprises two opposite factors, that is, the d-band center was attractive and the square of the coupling matrix element was repulsive, such that the contributions from the two factors can counteract each other at the top site.

Density Functional Theory Study of Hydrogen Adsorption on Fe5C2(001), Fe5C2(110), and Fe5C2(100)

Density functional theory calculations have been carried out for hydrogen adsorption on the (001), (110), and (100) surfaces of Fe5C2. At 1/3 and 2/3 monolyer (ML) on (001), the most stable hydrocarbon species is CsH, while CsH and CsH3 can coexist at 1 ML. On (110), only dissociated hydrogen is found at 2/5 ML, while CsH is the most stable hydrogen carbon species at 4/5 ML, and CsH and CH3 coexist at 6/5 ML. On (001) and (110) surfaces, CsH2 is less stable and can dissociate into CsH or convert into CsH3, respectively. These results are in agreement with the experimental observations. On the metallic Fe5C2(100) surface which lacks surface carbon atoms on the surface monolayer, dissociated hydrogen is found at 1/2 ML, while both dissociated hydrogen and activated H2 are found at 1 ML.

Preference for Vibrational over Translational Energy in a Gas-Surface Reaction

State-resolved gas-surface reactivity measurements revealed that vibrational excitation of nu3 (the antisymmetric C-H stretch) activates methane dissociation more efficiently than does translational energy. Methane molecules in the vibrational ground state require 45 kilojoules per mole (kJ/mol) of translational energy to attain the same reactivity enhancement provided by 36 kJ/mol of nu3 excitation. This result contradicts a key assumption underlying statistical theories of gas-surface reactivity and provides direct experimental evidence of the central role that vibrational energy can play in activating gas-surface reactions.

Chemisorption and reaction characteristics of methyl radicals on Cu(110)

Methyl radicals are generated by pyrolysis of azomethane, and the condition for achieving neat adsorption on Cu(110) is described for studying their chemisorption and reaction characteristics. The radical-surface system is examined by X-ray photoemission spectroscopy, ultraviolet photoemission spectroscopy, temperature-programmed desorption, low-energy electron diffraction (LEED), and high-resolution electron energy loss spectroscopy under ultrahigh vacuum conditions. It is observed that a small fraction of impinging CH3 radicals decompose into methylene possibly on surface defect sites. This type of CH2 radical has no apparent effect on CH3(ads) surface chemistry initiated by dehydrogenation to form active CH2(ads) followed by chain reactions to yield high-mass alkyl products. All thermal desorption products, such as H2, CH4, C2H4, C2H6, and C3H6, are detected with a single desorption peak near 475 K. The product yields increase with surface coverage until saturation corresponding to 0.50 monolayer of CH3(ads). The mass distribution is, however, invariant with initial CH3(ads) coverage, and all desorbed species exhibit first-order reaction kinetics. LEED measurement reveals a c(2 x 2) adsorbate structure independent of the amount of gaseous exposure. This strongly suggests that the radicals aggregate into close-packed two-dimensional islands at any exposure. The islanding behavior can be correlated with the reaction kinetics and is deemed to be essential for the chain propagation reactions. Some relevant aspects of the CH3/Cu(111) system are also presented. The new results are compared with those of prior studies employing methyl halides as radical sources. Major differences are found in the product distribution and desorption kinetics, and these are attributed to the influence of surface halogen atoms present in those earlier investigations.

Methane activation by nickel cluster cations, Ni[sub n][sup +] (n=2–16): Reaction mechanisms and thermochemistry of cluster-CH[sub x] (x=0–3) complexes

The kinetic energy dependences of the reactions of Ni+(n) (n=2-16) with CD(4) are studied in a guided ion beam tandem mass spectrometer over the energy range of 0-10 eV. The main products are hydride formation Ni(n)D+, dehydrogenation to form Ni(n)CD+(2), and double dehydrogenation yielding Ni(n)C+. These primary products decompose at higher energies to form Ni(n)CD+, Ni(n-1)D+, Ni(n-1)C+, Ni(n-1)CD+, and Ni(n-1)CD+(2). Ni(n)CD(2) (+) (n=5-9) and Ni(n-1)CD(2) (+) (n > or =4) are not observed. In general, the efficiencies of the single and double dehydrogenation processes increase with cluster size. All reactions exhibit thresholds, and cross sections for the various primary and secondary reactions are analyzed to yield reaction thresholds from which bond energies for nickel cluster cations to C, CD, CD(2), and CD(3) are determined. The relative magnitudes of these bond energies are consistent with simple bond order considerations. Bond energies for larger clusters rapidly reach relatively constant values, which are used to estimate the chemisorption energies of the C, CD, CD(2), and CD(3) molecular fragments to nickel surfaces.