Search

Understanding the mechanism of CO2 reduction on iron is crucial for the design of more efficient and cheaper iron electrocatalyst for CO2 conversion. In the present study, we have employed spin-polarized density functional theory calculations within the generalized gradient approximation (DFT-GGA) to elucidate the mechanism of CO2 reduction into carbon monoxide and formic acid on the Fe (100) facet. We also sort to understand the transformations of the other isomers of adsorbed CO2 on iron as earlier mechanistic studies are centred on the transformations of the C2v geometry alone and not the other possible conformations i.e., flip-C2v and Cs modes. Two alternative reduction routes were considered i.e., the direct CO2 dissociation against the hydrogen-assisted CO2 transformation through formate and carboxylate into CO and formic acid. Our results show that CO2 in the C2v mode is the precursor to the formation of both products i.e., CO and formic acid. Both the formation and transformation of CO2 in the Cs and flip-C2v is challenging kinetically and thermodynamically compared to the C2v mode. The formic acid formation is favoured over CO via the reverse water gas shift reaction mechanism on Fe (100). Both formic acid formation and CO formation will proceed via the carboxylate intermediate since formate is a stable intermediate whose transformation into formic acid is challenging both kinetically and thermodynamically. Graphic abstract

Insight into the detailed mechanism of the Sabatier reaction on iron is essential for the design of cheap, environmentally benign, efficient and selective catalytic surfaces for CO2 reduction. Earlier attempts to unravel the mechanism of CO2 reduction on pure metals including inexpensive metals focused on Ni and Cu; however, the detailed mechanism of CO2 reduction on iron is not yet known. We have, thus, explored with spin-polarized density functional theory calculations the relative stabilities of intermediates and kinetic barriers associated with methanation of CO2 via the CO and non-CO pathways on the Fe (111) surface. Through the non-CO (formate) pathway, a dihydride CO2 species (H2CO2), which decomposes to aldehyde (CHO), is further hydrogenated into methoxy, methanol and then methane. Through the CO pathway, it is observed that the CO species formed from dihydroxycarbene is not favorably decomposed into carbide (both thermodynamically and kinetically challenging) but CO undergoes associative hydrogenation to form CH2OH which decomposes into CH2, leading to methane formation. Our results show that the transformation of CO2 to methane proceeds via the CO pathway, since the barriers leading to alkoxy transformation into methane are high via the non-CO pathway. Methanol formation is more favored via the non-CO pathway. Iron (111) shows selectivity towards CO methanation over CO2 methanation due to differences in the rate-determining steps, i.e., 91.6 kJ mol−1 and 146.2 kJ mol−1, respectively.

Iron and nickel are known active sites in the enzyme carbon monoxide dehydrogenases (CODH) which catalyzes CO2 to CO reversibly. The presence of nickel impurities in the earth abundant iron surface could provide a more efficient catalyst for CO2 degradation into CO, which is a feedstock for hydrocarbon fuel production. In the present study, we have employed spin-polarized dispersion-corrected density functional theory calculations within the generalized gradient approximation to elucidate the active sites on Fe (100)-Ni bimetals. We sort to ascertain the mechanism of CO2 dissociation to carbon monoxide on Ni deposited and alloyed surfaces at 0.25, 0.50 and 1 monolayer (ML) impurity concentrations. CO2 and (CO + O) bind exothermically i.e., -0.87 eV and − 1.51 eV respectively to the bare Fe (100) surface with a decomposition barrier of 0.53 eV. The presence of nickel generally lowers the amount of charge transferred to CO2 moiety. Generally, the binding strengths of CO2 were reduced on the modified surfaces and the extent of its activation was lowered. The barriers for CO2 dissociation increased mainly upon introduction of Ni impurities which is undesired. However, the 0.5 ML deposited (FeNi0.5(A)) surface is promising for CO2 decomposition, providing a lower energy barrier (of 0.32 eV) than the pristine Fe (100) surface. This active 1-dimensional defective FeNi0.5(A) surface provides a stepped surface and Ni-Ni bridge binding site for CO2 on Fe (100). Ni-Ni bridge site on Fe (100) is more effective for both CO2 binding or sequestration and dissociation compared to the stepped surface providing the Fe-Ni bridge binding site.

The oxygen adsorption and subsequent reduction on the {100} and {110} surfaces of 25% Ba-doped LaMnO3 (LBM25) have been studied at the density functional theory (DFT) with Hubbard correction and the results compared with adsorption on 25% Ca-doped LaMnO3 (LCM25) and Sr-doped LaMnO3 (LSM25). The trend in the reduction energies at the Mn cation sites are predicted to be in the order LSM25 2 precursors at the Mn cation sites of LCM25, LSM25 and LBM25 are thermodynamically stable, when compared directly with the adsorption energies (Eads = − 0.56 to − 1.67 eV) reported for the stable molecular O2 precursors on the Pt, Ni, Pd, Cu and Ir {111} surfaces. The predicted Gibbs energies as a function of temperature (T = 500–1100 °C) and pressures (p = 0.2 atm) for the adsorption and dissociation on the surfaces were negative, an indication of the feasibility of oxygen reduction reaction on the {100} and {110} surfaces at typical operating temperatures reported in this work.

Graphene despite its high surface area has very limited activity towards the oxygen reduction reaction (ORR), demonstrating selectivity towards the unfavorable two-electron mechanism. We have employed the spin polarized density functional theory (DFT) method to investigate the effect of the heteroatom p-type beryllium (Be) dopant on the oxygen reduction reaction activity of graphene. The preferred doping sites, active sites and reaction mechanism available on the doped graphene surfaces were investigated with increasing Be concentrations of 0.03 ML, 0.06 ML and 0.09 ML. Our results reveal that oxygen does not bind strongly to bare graphene, and Be at the lattice sites provides site for the oxygen adsorption and ORR. Oxygen binds dissociatively on the doped surfaces preferentially in the order 0.06 ML > 0.09 ML > 0.03 ML. From this studies introduction of Be impurities in a single honeycomb ring of graphene has significant impact on the binding and adsorbate-adsorbent interactions which leads to dissociative adsorption of oxygen, favouring the 4e- ORR pathway. The reaction is kinetically favoured most on the surface with a stronger O* binding and weaker OH* binding. Overall, the 0.03 ML and 0.06 ML doped surfaces are the most active facets for the ORR showing exergonic reaction energies at 0 V.

Despite series of experimental studies that reveal unique activities of late 3d transition metals and their role in microorganisms known for CO2 conversion, these surfaces are not industrially viable yet. An insight into the elementary steps of surface catalytic processes is crucial for effective surface modification and design. The mechanisms of CO2 transformation into CO, through the reverse water gas shift and methane reforming, are being studied. Mechanisms of CO2 methanation is also being explored by the Sabatier reaction into methane. This review covers both experimental and theoretical studies into the mechanisms of CO2 reduction into CO and methane, on single metals and bimetals of late 3d transition metals, i.e. Fe, Co, Ni, Cu and Zn. This paper highlights progress and gaps still existing in our knowledge of the reaction mechanisms. These mechanistic studies reveal CO2 activation and reduction mechanisms are specific to both composition and surface facet. Surfaces with least CO2 binding potential are seen to favour CO and O binding and provide higher barriers to dissociation. No direct correlation has been seen between binding strength of CO2 and its degree of activation. Hydrogen-assisted dissociation is seen to be generally favoured kinetically on Cu and Ni surfaces over direct dissociation except on the Ni (211) surface. Methane production on Cu and Ni surfaces is seen to occur via the non-formate pathway. Hydrogenation reactions have focused on Cu and Ni, and more needs to be done on other surfaces, i.e. Co, Fe and Zn.

Understanding the mechanism of the catalytic upgrade of bio-oils via the process of hydrodeoxygenation (HDO) is desirable to produce targeted oxygen-deficient bio-fuels. We have used calculations based on the density functional theory to investigate the reaction mechanism of HDO of guaiacol over Cu (111) surface in the presence of H2, leading to the formation of catechol and anisole. Our analysis of the thermodynamics and kinetics involved in the reaction process shows that catechol is produced via direct demethylation, followed by dehydrogenation of –OH and re-hydrogenation of catecholate in a concerted fashion. The de-methylation step is found to be the rate-limiting step for catechol production with a barrier of 1.97 eV. Formation of anisole will also proceed via the direct dehydroxylation of guaiacol followed by hydrogenation. Here, the rate-limiting step is the dehydroxylation step with an energy barrier of 2.07 eV. Thermodynamically, catechol formation is favored while anisole formation is not favored due to the weaker interaction seen between anisole and the Cu (111) surface, where the binding energies of guaiacol, catechol, and anisole are -1.90 eV, −2.18 eV, and −0.72 eV, respectively. The stepwise barriers also show that the Cu (111) surface favors catechol formation over anisole as the rate-limiting barrier is higher for anisole production. For catechol, the overall reaction is downhill, implying that this reaction path is thermodynamically and kinetically preferred and that anisole, if formed, will more easily transform.

The mechanisms of LMnO3 (L = O−, Cl, NPH3, CH3, and Cp)-catalyzed oxidation of ethyne has been studied on the singlet and triplet hypersurfaces at the M06/6-311G(d) level of theory. For the first step, the [3 + 2] pathways to the formation of the metalla-2,5-dioxol-3-ene intermediate are kinetically and thermodynamically the most favored pathways in all the complexes studied; it is favored over the [2 + 2] addition pathways to the metallaoxetene intermediate. The formation of the oxirene precursor that could give the oxirene the reported key intermediates in the ozonolysis of alkynes would most likely result from the oxidation of ethyne by MnO3Cl on the triplet potential energy surface (PES).

Metal clusters of both iron (Fe) and nickel (Ni) have been found in nature as active electro-catalytic sites, for example in the enzyme carbon mono-oxide dehydrogenase found in autotrophic organisms. Thus, surface modification of iron with nickel could improve the surface work function to enhance catalytic applications. The effects of surface modifications of iron by nickel on the structural and electronic properties have been studied using spin-polarised density functional theory calculations within the generalised gradient approximation. The thermodynamically preferred sites for Ni adsorption on the Fe (100), (110) and (111) surfaces have been studied at varying monolayer coverages (including 0.25 ML and 1 ML). The work function of the bare Fe surfaces is found to be of the order (100) ∼ (111) < (110) i.e. 3.80 eV ∼ 3.84 eV < 4.76 eV, which is consistent with earlier studies. The adsorption energies show that monolayer Ni deposition is thermodynamically favoured on the (100) and (111) surfaces, but not on the (110) surface. Expansion of the first interlayer spacing (d12) of all three Fe surfaces is observed upon Ni deposition with the extent of expansion decreasing in the order (111) > (110) > (100), i.e. 6.78% > 5.76% > 1.99%. The extent of relaxation is magnified on the stepped (111) surface (by 1.09% to 30.88%), where the Ni coordination number is highest at 7 compared to 5 on the (100) facet and 4 on the (110) facet. The Ni deposition changes the work functions of the various surfaces due to charge reordering illustrated by charge density plots, where the work function is reduced only on the (110) surface by 0.04 eV, 0.16 eV and 0.17 eV at 1 ML, 0.5 ML and 0.25 ML respectively, with a concomitant increase in the surface dipole (polarity). This result implies enhanced electron activity and electrochemical reactivity on the most stable and therefore frequently occurring Ni-doped (110) facet compared to the clean (110) facet, which has implications for the development of improved Fe electro-catalysts (for example for CO2 activation and reduction). These findings improve our understanding of the role of surface topology and stability on the extent of Ni interactions with Fe surfaces and the extent to which the Fe surface structures and properties are altered by the Ni deposition.

The mechanisms of the oxidation of tetramethylethylene (TME) by permanganyl chloride (MnO3Cl) have been explored on the singlet and triplet potential energy surfaces at the B3LYP LANL2DZ/6-31G (d) level of theory. The results show that the pathway leading to the formation of the five-membered dioxylate through concerted [3 + 2] addition is favored kinetically and thermodynamically over the three other possible pathways, namely the [2 + 2] addition via the transient metallaoxetane intermediate, epoxidation, and hydrogen transfer pathways. The epoxide precursor that on hydrolysis would yield the epoxide product will most likely arise from a stepwise path through the intermediacy of an organometallic intermediate. This pathway affords the product that is more stable (thermodynamically favorable). However, kinetically, both the stepwise and the concerted [2 + 1] addition pathways leading to the epoxide precursors are very competitive (activation barrier difference of