Coverage-Dependent CO Adsorption and Dissociation Mechanisms on Iron Surfaces from DFT Computations

A DFT study towards dynamic structures of iron and iron carbide and their effects on the activity of the Fischer-Tropsch process

The Fe-based Fischer-Tropsch synthesis (FTS) catalyst shows a rich phase chemistry under pre-treatment and FTS conditions. The exact structural composition of the active site, whether iron or iron carbide (FeCx), is still controversial. Aiming to obtain an insight into the active sites and their role in affecting FTS activity, the swarm intelligence algorithm is implemented to search for the most stable Fe(100), Fe(110), Fe(210) surfaces with different carbon ratios. Then, ab initio atomistic thermodynamics and Wulffman construction were employed to evaluate the stability of these surfaces at different chemical potentials of carbon. Their FTS reactivity and selectivity were later assessed by semi-quantitative micro-kinetic equations. The results show that stability, reactivity, and selectivity of the iron are all affected by the carbonization process when the carbon ratio increases. Formation of the carbide, a rather natural process under experimental conditions, would moderately increase the turnover frequency (TOF), but both iron and iron carbide are active to the reaction.

H2O Derivatives Mediate CO Activation in Fischer-Tropsch Synthesis: A Review

The process of Fischer-Tropsch synthesis is commonly described as a series of reactions in which CO and H₂ are dissociated and adsorbed on the metals and then rearranged to produce hydrocarbons and H₂O. However, CO dissociation adsorption is regarded as the initial stage of Fischer-Tropsch synthesis and an essential factor in the control of catalytic activity. Several pathways have been proposed to activate CO, namely direct CO dissociation, activation hydrogenation, and activation by insertion into growing chains. In addition, H₂O is considered an important by-product of Fischer-Tropsch synthesis reactions and has been shown to play a key role in regulating the distribution of Fischer-Tropsch synthesis products. The presence of H₂O may influence the reaction rate, the product distribution, and the deactivation rate. Focus on H₂O molecules and H₂O-derivatives (H*, OH* and O*) can assist CO activation hydrogenation on Fe- and Co-based catalysts. In this work, the intermediates (C*, O*, HCO*, COH*, COH*, CH*, etc.) and reaction pathways were analyzed, and the H₂O and H₂O derivatives (H*, OH* and O*) on Fe- and Co-based catalysts and their role in the Fischer-Tropsch synthesis reaction process were reviewed.

Initial stage of carbonization of iron during hydrocarbons dissociation: a molecular dynamics study

The carbonization of iron is a very important early phenomenon in the field of heterogeneous catalysis and the petrochemical industry, but the mechanism is still controversial. In this work, the carbonization mechanism and carbonization structure of iron nanoparticles by different carbon sources (CH₄, C₂H₆, C₂H₄, C₂H₂) were systematically investigated using the reactive molecular dynamics method. The results show that saturated alkanes are dehydrogenated while adsorbed, but unsaturated olefins and alkynes undergo bond-breaking while adsorbed. The C-H bond is more likely to break than the C-C bond. Hydrocarbons with high carbon content have a strong ability to carbonize Fe nanoparticles under the same conditions. For C₂H₄ and C₂H₂, the C atoms generated from dissociation form a large number of long carbon chains intertwined with branched chains and multiple carbon rings. The C2 species formed by C₂H₂ after complete dehydrogenation diffuse rapidly to the interior of the nanoparticles, releasing the surface active sites and accelerating the carbonization process. Carbon-rich iron carbides (FeC_x) with different Fe/C ratios were obtained by carbonization with different carbon sources. In addition, the Fe(110) surface exhibits the strongest carburizing ability. These findings provide systematic insights into the initial stages of metal Fe carburization.

Water Gas Shift Reaction Activity on Fe (110): A DFT Study

Metal Fe is one of the phases existing on iron-based catalysts for a high-temperature water gas shift reaction (WGSR), but research on the activity of metal Fe in WGSR is almost not reported. In this work, the density functional theory (DFT) method was used to systematically study the reaction activity and mechanisms of WGSR on metal Fe (110), including the dissociation of H2O, the transformation of CO and the formation of H2, as well as the analysis of surface electronic properties. The results show that (1) the direct dissociation of H2O occurs easily on Fe (110) and the energy barrier is less than 0.9 eV; (2) the generation of CO2 is difficult and its energy barrier is above 1.8 eV; (3) H migrates easily on the Fe surface and the formation of H2 also occurs with an energy barrier of 1.47 eV. Combined with the results of Fe3O4, it can be concluded that the active phase should be Fe3O4 with O vacancy defects, and the iron-rich region plays an important role in promoting the formation of H2 in WGSR.

Reactivity and Passivation of Fe Nanoclusters on h-BN/Rh(111)

The reactivity of iron nanocluster arrays on h‐BN/Rh(111) was studied using in situ high‐resolution X‐ray photoelectron spectroscopy. The morphology and reactivity of the iron nanoclusters (Fe‐NCs) were investigated by CO adsorption. On‐top and hollow/edge sites were determined to be the available adsorption sites on the as‐prepared Fe‐NCs and CO dissociation was observed at 300 K. C‐ and O‐precovered Fe‐NCs showed no catalytic activity towards CO dissociation because the hollow/edge sites were blocked by the C and O atoms. Therefore, these adsorption sites were identified to be the most active sites of the Fe‐NCs.

Ab initio investigation of the role of the d-states occupation on the adsorption properties of H2, CO, CH4 and CH3OH on the Fe13, Co13, Ni13 and Cu13 clusters

Here, we report a theoretical investigation, based on density functional theory calculations, into the role of the occupation d-states on the adsorption properties of CH₄, CO, H₂ and CH₃OH on 3d 13-atom transition-metal (TM₁₃) clusters (TM = Fe, Co, Ni, Cu). Except for Cu₁₃, a gradual increase in the occupation of the d-states, i.e., from Fe₁₃ to Ni₁₃, increases the magnitude of the adsorption energy almost linearly for the H₂/TM₁₃ and CO/TM₁₃ systems, which can be explained by the enhancement of the sp-d hybridization due to the shift of the d-states towards the highest occupied molecular orbital (HOMO). For Cu₁₃, the d-states are located well below the HOMO, which reduces the sp-d hybridization, and hence, a smaller adsorption energy is obtained. However, this picture does not hold for CH₄/TM₁₃ and CH₃OH/TM₁₃, where the adsorption energy has nearly the same value for all TM₁₃ clusters, which can be explained by electrostatic effects such as local polarization of the molecules and nearby TM atoms, and hence, the basic features of physisorption systems. Based on the electron density difference, the polarization effects are slightly larger for systems with empty d-states.

Exploring direct and hydrogen-assisted CO activation on iridium surfaces – surface dependent activity

To understand CO activation on iridium surfaces, direct dissociation, H-assisted activation and hydrogenation to methanol were computed on the flat Ir(111) and Ir(100), corrugated Ir(110) and Ir(210), and stepped Ir(311) and Ir(221) surfaces.

Iron assisted formation of CO2 over condensed CO and its relevance to interstellar chemistry

Catalytic conversion of CO to CO₂ assisted by neutral iron atoms has been investigated in ultrahigh vacuum (UHV) under cryogenic conditions (10 K).

Altering CO binding on gold cluster cations by Pd-doping

The introduction of dopant atoms into metal nanoparticles is an effective way to control the interaction with adsorbate molecules and is important in many catalytic processes. In this work, experimental and theoretical evidence of the influence of Pd doping on the bonding between small cationic AuN+ clusters and CO is presented. The CO adsorption is studied by combining low-pressure collision cell reactivity and infrared multiple photon dissociation spectroscopy experiments with density functional theory calculations. Measured dissociation rates of cluster-CO complexes (N ≤ 21) allow the estimation of cluster-CO binding energies, showing that Pd doping increases the CO adsorption energy to an extent that is size-dependent. These trends are reproduced by theoretical calculations up to N = 13. In agreement with theory, measurements of the C-O vibrational frequency suggest that for the doped PdAuN-1+ (N = 3-5, 11) clusters, CO adsorbs on an Au atom, while for N = 6-10 and N = 12-14, CO interacts directly with the Pd dopant. A pronounced red-shifting of the C-O vibrational frequency is observed when CO interacts directly with the Pd dopant, indicating a significant back-donation of electron charge from Pd to CO. In contrast, the blue-shifted frequencies, observed when CO interacts with an Au atom, indicate that σ-donation dominates the Au-CO interaction. Studying such systems at the sub-nanometre scale enables a fundamental comprehension of the interactions between adsorbates, dopants and the host (Au) species at the atomic level.

Adsorption mechanism of CO molecule on Al(111) surface: periodic DFT investigation

The adsorption mechanism of the CO molecule on Al(111) surface has been investigated systematically at the atom-molecule level by the method of periodic density functional theory. The adsorption energies, adsorption structures, charge transfer, and density of states have been calculated in a wide range of coverage. It is found that the hcp-hollow site is the energetically favorable site. A significant positive correlation has been found between the adsorption energy (Eads) and coverage. The adsorbed CO molecules are almost perpendicular on the surface with the C atom facing the surface. There is an obvious charge transfer from Al atoms to the C atom; the Al atoms that have interaction with the C atom offer the most charge. The 4σ, 1π, and 5σ molecular orbitals of CO are found to contribute to bonding with the Al. The charges filling in the 2π molecular orbital contribute to C–O bond activation. In conclusion, the passivation of aluminum surface and the activation of CO molecule occur simultaneously in the adsorption of CO on Al surface.

Interaction between H2O, N2, CO, NO, NO2 and N2O molecules and a defective WSe2 monolayer

In this study, the interaction between gas molecules, including H₂O, N₂, CO, NO, NO₂ and N₂O, and a WSe₂ monolayer containing an Se vacancy (denoted as V_Se) has been theoretically studied. Theoretical results show that H₂O and N₂ molecules are highly prone to be physisorbed on the V_Se surface. The presence of the Se vacancy can significantly enhance the sensing ability of the WSe₂ monolayer toward H₂O and N₂ molecules. In contrast, CO and NO molecules highly prefer to be molecularly chemisorbed on the V_Se surface with the non-oxygen atom occupying the Se vacancy site. Furthermore, the exposed O atoms of the molecularly chemisorbed CO or NO can react with additional CO or NO molecules, to produce C-doped or N-doped WSe₂ monolayers. The calculated energies suggest that the filling of the CO or NO molecule and the removal of the exposed O atom are both energetically and dynamically favorable. Electronic structure calculations show that the WSe₂ monolayers are p-doped by the CO and NO molecules, as well as the C and N atoms. However, only the NO molecule and N atom doped WSe₂ monolayers exhibit significantly improved electronic structures compared with V_Se. The NO₂ and N₂O molecules will dissociate directly to form an O-doped WSe₂ monolayer, for which the defect levels due to the Se vacancy can be completely removed. The calculated energies suggest that although the dissociation processes for NO₂ and N₂O molecules are highly exothermic, the N₂O dissociation may need to operate at an elevated temperature compared with room temperature, due to its large energy barrier of ∼1 eV.

Suppression by Pt of CO adsorption and dissociation and methane formation on Fe5C2(100) surfaces

To understand the chemical origin of platinum promotion effects on iron based Fischer–Tropsch synthesis catalysts, the effects of Pt on CO adsorption and dissociation as well as surface carbon hydrogenation on the Fe₅C₂(100) facet with different surface C* contents have been studied using the spin-polarized density functional theory method.

Ethanol synthesis from syngas over Cu(Pd)-doped Fe(100): a systematic theoretical investigation

The selectivity of ethanol formation are in the order of Fe₃Cu₆/Fe(100) > Fe₃Pd₆/Fe(100) > Fe₉/Fe(100) > Cu₉/Fe(100).

Adsorption of atmospheric gases on cementite 010 surfaces

We study the adsorption of a series of small molecules on the nonstoichiometric {010} surface of cementite (θ-Fe3C) by means of first-principles calculations. We find that CO, N2, H2O, and CH4 prefer to adsorb over iron atoms in an atop configuration. O2, CO2, and the OH radical prefer a configuration bridging two iron atoms and CH2O adsorbs in a configuration bridging a surface iron atom and a surface carbon atom. Adsorption energies are small for H2, CO2, and CH4, indicating a physisorption process, while those for CO, CH2O and especially for O2 and the OH radical are large, indicating a strong chemisorption process. H2O and N2 display adsorption energies between these two extremes, indicating moderate chemisorption. The dissociation of H2, CH2O, the OH radical, and O2 is favoured on this surface. Comparison with adsorption on Fe{100} surfaces indicates that most of these gases have similar adsorption energies on both surfaces, with the exception of CO and the OH radical. In addition, we find similarities between the reactivities of cementite and Mo2C surfaces, due to the similar covalent character of both carbides.

Mechanism of coverage dependent CO adsorption and dissociation on the Mo(100) surface

The mechanism of coverage dependent CO adsorption and dissociation on the Mo(100) surface was investigated using periodic density functional theory. Structure optimization and frequency calculation were carried out using the GGA-PBE method and a p(3 × 3) supercell model. Energetic data have been obtained using the revised PBE method and the PBE optimized structures. CO adsorption prefers tilted adsorption configuration at the 4-fold hollow sites at low coverage and tilted and atop configurations at high coverage. The computed C-O stretching frequencies of the tilted and atop adsorbed CO molecules agree very well with the experimental results. Starting from the saturation coverage, five binding states have been found: two for molecular (α) CO adsorption and three for dissociative (β) CO adsorption, which are in agreement with the temperature-programmed desorption experiments. In addition, CO prefers dissociation with very low barriers in all coverages as long as free sites are available and is coverage independent; this nicely explains the observed CO dissociation at very low temperatures. All such agreements validate our computational methods and provide the basis of further studies.

Theoretical insights into the effect of terrace width and step edge coverage on CO adsorption and dissociation over stepped Ni surfaces

We rationalized Ni(211) as a representative model for stepped surfaces and explored the effect of coverage on CO activation.

Reaction of CO, H2O, H2 and CO2 on the clean as well as O, OH and H precovered Fe(100) and Fe(111) surfaces

The adsorption and activation of CO, H₂O, CO₂ and H₂ on the clean as well as O, OH and H precovered Fe(100) surface at 0.25 ML coverage and Fe(111) surface at 0.33 ML coverage were computed (GGA-PBE) to investigate the catalytic activity of metallic iron.

Influence of surface vacancy defects on the carburisation of Fe 110 surface by carbon monoxide

Adsorption and dissociation of gaseous carbon monoxide (CO) on metal surfaces is one of the most frequently occurring processes of carburisation, known as primary initiator of metal dusting corrosion. Among the various factors that can significantly influence the carburisation process are the intrinsic surface defects such as single surface vacancies occurring at high concentrations due to their low formation energy. Intuitively, adsorption and dissociation barriers of CO are expected to be lowered in the vicinity of a surface vacancy, due to the strong attractive interaction between the vacancy and the C atom. Here the adsorption energies and dissociation pathways of CO on clean and defective Fe 110 surface are explored by means of density functional theory. Interestingly, we find that the O adatom, resulting from the CO dissociation, is unstable in the electron-deficit neighbourhood of the vacancy due to its large electron affinity, and raises the barrier of the carburisation pathway. Still, a full comparative study between the clean surface and the vacancy-defected surface reveals that the complete process of carburisation, starting from adsorption to subsurface diffusion of C, is more favourable in the vicinity of a vacancy defect.

Adsorption and dissociation of H2O on the (001) surface of uranium mononitride: energetics and mechanism from first-principles investigation

The interfacial interaction of uranium mononitride (UN) with water from the environment unavoidably leads to corrosion of nuclear fuels, which affects a lot of processes in the nuclear fuel cycle. In this work, the microscopic adsorption behaviors of water on the UN(001) surface as well as water dissociation and accompanying H2 formation mechanisms have been investigated on the basis of DFT+U calculations and ab initio atomistic thermodynamics. For adsorption of one H2O monomer, the predicted adsorption energies are -0.88, -2.07, and -2.07 eV for the most stable molecular, partially dissociative, and completely dissociative adsorption, respectively. According to our calculations, a water molecule dissociates into OH and H species via three pathways with small energy barriers of 0.78, 0.72, and 0.85 eV, respectively. With the aid of the neighboring H atom, H2 formation through the reaction of H* + OH* can easily occur via two pathways with energy barriers of 0.61 and 0.36 eV, respectively. The molecular adsorption of water shows a slight coverage dependence on the surface while this dependence becomes obvious for partially dissociative adsorption as the water coverage increases from 1/4 to 1 ML. In addition, based on the "ab initio atomistic thermodynamic" simulations, increasing H2O partial pressure will enhance the stability of the adsorbed system and water coverage, while increasing temperature will decrease the H2O coverage. We found that the UN(001) surface reacts easily with H2O at room temperature, leading to dissolution and corrosion of the UN fuel materials.

Theoretical study about Mo2C(101)-catalyzed hydrodeoxygenation of butyric acid to butane for biomass conversion

To explore the conversion mechanism of fatty acids to long-chain alkanes using molybdenum carbide catalysts, the full potential energy surface of the hydrogenation of butyric acid to butane on the H-pre-covered hexagonal Mo₂C(101) surface has been systematically computed.

Structures of seven molybdenum surfaces and their coverage dependent hydrogen adsorption

Surface morphology of metallic molybdenum and coverage dependent hydrogen adsorption.

Insight into both coverage and surface structure dependent CO adsorption and activation on different Ni surfaces from DFT and atomistic thermodynamics

CO adsorption and activation from low to high coverage on Ni catalyst.

Direct versus hydrogen-assisted CO dissociation over stepped Ni and Ni3Fe surfaces: a computational investigation

The adsorption and dissociation of CO over stepped Ni and Ni3Fe surfaces were systematically studied using density functional theory slab calculations. Both (211)-like surface structure terminations (NiNi step and NiFe step, denoted as Ni3Fe(211)-AA and Ni3Fe(211)-AB) are considered for Ni3Fe. Direct scission of the C-O bond in CO is identified as the least likely one among the three proposed dissociation pathways and CO dissociation via a CHO intermediate appears to be most feasible at low CO coverage on pure and alloyed Ni(211) surfaces. The priority of H-assisted CO dissociation might originate from the more activated C-O bond in COH and CHO. Compared to Ni(211), the Ni3Fe(211)-AB surface could facilitate CO activation especially for the most possible CHO intermediate mechanism, whose rate-limiting step is found to be altered. The d-band center theory and Mulliken charge analysis are also employed to explain the activity difference between Ni3Fe(211)-AB and Ni3Fe(211)-AA. The significant structural sensitivity of CO dissociation highlights the importance of Fe locating in the step edge and the high reactivity of Ni3Fe(211)-AB is largely ascribed to the synergistic effect between Ni and Fe at the step edge.

Quantifying the origin of inter-adsorbate interactions on reactive surfaces for catalyst screening and design

The adsorption energy of reactant molecules and reaction intermediates is one of the key descriptors of catalytic activity of surfaces and is commonly used as a metric in screening materials for design of heterogeneous catalysts. The efficacy of such screening schemes depends on the accuracy of calculated adsorption energies under reaction conditions. These adsorption energies can depend strongly on interactions between adsorbed molecules in the adlayer. However, these interactions are typically not accounted for in screening procedures that use DFT-based zero-coverage adsorption energies. Identifying the physical mechanisms behind these interactions is essential to model realistic catalyst surfaces under reaction conditions and to understand the dependence of adsorption energies on reaction parameters like surface strain and composition. This article describes a method to quantitatively resolve the observed inter-adsorbate interactions into various direct adsorbate-adsorbate interactions (i.e. Coulombic and steric) and surface-mediated interactions (i.e. adsorbate-induced surface relaxation and change in electronic structure) by combining density functional theory and cluster-expansion calculations of coverage-dependent adsorption energies. The approach is implemented on a model catalyst surface of FeS2(100) reacting with H2S molecules. We find that the adsorption energy of H2S molecules can be affected by over 0.55 eV by the repulsive inter-adsorbate interactions caused primarily by the adsorbate-induced changes to the electronic structure of the FeS2 surface. These interactions also show a strong monotonic dependence on surface strain, being three times stronger on compressively strained surfaces than on surfaces under tensile strain. The large magnitude of inter-adsorbate interactions as well as their strong dependence on lattice strain demonstrate the need for using coverage-dependent adsorption energies for more accurate screening, for example for strained catalytic systems like core-shell and overlayer structures.

Coverage dependent water dissociative adsorption on Fe(110) from DFT computation

Using density functional theory calculations and ab initio atomistic thermodynamics, H2O adsorption and dissociation on the Fe(110) p(4 × 4) surface at different coverages have been computed. At the lowest coverage, the adsorbed H2O, OH, O and H species can migrate easily on the surface. For (H2O)n adsorption, H2O molecules donating H atoms for H-bonding adsorb more strongly than those accepting H atoms for H-bonding. Monomeric H2O dissociation is favored both thermodynamically and kinetically. On nO pre-covered Fe(110) surfaces (n = 1-8), H2O dissociation is accessible for nO + H2O (n = 1-7) both kinetically and thermodynamically, while H2O desorption instead of dissociation occurs at n = 8. With the increased number of surface O atoms, H2 dissociative adsorption energies vary in a narrow range for n = 1-4 and decrease for n = 5-7, while at n = 8, the surface does not adsorb H2. At low OH coverage (n = 2, 4), OH groups are perpendicularly adsorbed without H-bonding, while for n≥ 6, adsorbed OH groups are linearly arranged and stabilized by H-bonding. The maximal OH coverage (n = 12) is 0.75 ML and the reasonable O coverage (n = 7) is 0.44 ML, in line with the experiment. The calculated desorption temperatures of H2O and H2 agree well with the available experimental data. These results provide fundamental insights into water-involved reactions catalyzed by iron and interaction mechanisms of water interaction with metal surfaces.

Coverage dependent adsorption and co-adsorption of CO and H2 on the CdI2-antitype metallic Mo2C(001) surface

The adsorption and co-adsorption of CO and H₂ at different coverage on the CdI₂-antitype metallic Mo₂C(001) surface termination have been systematically computed at the level of periodic density functional theory.

High coverage adsorption and co-adsorption of CO and H2 on Ru(0001) from DFT and thermodynamics

The adsorption and co-adsorption of CO and H₂ at different coverages on p(4 × 4) Ru(0001) have been computed using periodic density functional theory (GGA-RPBE) and atomistic thermodynamics.