On the Capabilities of Transition Metal Carbides for Carbon Capture and Utilization Technologies.

The search for cheap and active materials for the capture and activation of CO2 has led to many efforts aimed at developing new catalysts. In this context, earth-abundant transition metal carbides (TMCs) have emerged as promising candidates, garnering increased attention in recent decades due to their exceptional refractory properties and resistance to sintering, coking, and sulfur poisoning. In this work, we assess the use of Group 5 TMCs (VC, NbC, and TaC) as potential materials for carbon capture and sequestration/utilization technologies by combining experimental characterization techniques, first-principles-based multiscale modeling, vibrational analysis, and catalytic experiments. Our findings reveal that the stoichiometric phase of VC exhibits weak interactions with CO2, displaying an inability to adsorb or dissociate it. However, VC often exhibits the presence of surface carbon vacancies, leading to significant activation of CO2 at room temperature and facilitating its catalytic hydrogenation. In contrast, stoichiometric NbC and TaC phases exhibit stronger interactions with CO2, capable of adsorbing and even breaking of CO2 at low temperatures, particularly notable in the case of TaC. Nevertheless, NbC and TaC demonstrate poor catalytic performance for CO2 hydrogenation. This work suggests Group 5 TMCs as potential materials for CO2 abatement, emphasizes the importance of surface vacancies in enhancing catalytic activity and adsorption capability, and provides a reference for using the infrared spectra as a unique identifier to detect oxy-carbide phases or surface C vacancies within Group 5 TMCs.

Comprehensive Density Functional and Kinetic Monte Carlo Study of CO2 Hydrogenation on a Well-Defined Ni/CeO2 Model Catalyst: Role of Eley-Rideal Reactions.

A detailed multiscale study of the mechanism of CO2 hydrogenation on a well-defined Ni/CeO2 model catalyst is reported that couples periodic density functional theory (DFT) calculations with kinetic Monte Carlo (kMC) simulations. The study includes an analysis of the role of Eley-Rideal elementary steps for the water formation step, which are usually neglected on the overall picture of the mechanism, catalytic activity, and selectivity. The DFT calculations for the chosen model consisting of a Ni4 cluster supported on CeO2 (111) show large enough adsorption energies along with low energy barriers that suggest this catalyst to be a good option for high selective CO2 methanation. The kMC simulations results show a synergic effect between the two 3-fold hollow sites of the supported Ni4 cluster with some elementary reactions dominant in one site, while other reactions prefer the another, nearly equivalent site. This effect is even more evident for the simulations explicitly including Eley-Rideal steps. The kMC simulations reveal that CO is formed via the dissociative pathway of the reverse water-gas shift reaction, while methane is formed via a CO2 → CO → HCO → CH → CH2 → CH3 → CH4 mechanism. Overall, our results show the importance of including the Eley-Rideal reactions and point to small Ni clusters supported on the CeO2 (111) surface as potential good catalysts for high selective CO2 methanation under mild conditions, while very active and selective toward CO formation at higher temperatures.

Structural, electronic, and magnetic properties of Ni nanoparticles supported on the TiC(001) surface.

Metals supported on transition metal carbides are known to exhibit good catalytic activity and selectivity, which is interpreted in terms of electron polarization induced by the support. In the present work we go one step further and investigate the effect that a titanium carbide (TiC) support has on the structural, electronic, and magnetic properties of a series of Ni nanoparticles of increasing size exhibiting a two- or three-dimensional morphology. The obtained results show that three-dimensional nanoparticles are more stable and easier to form than their homologous two-dimensional counterparts. Also, comparison to previous results indicates that, when used as the support, transition metal carbides have a marked different chemical activity with respect to oxides. The analysis of the magnetic moments of the supported nanoparticles evidences a considerable quenching of the magnetic moment that affects mainly the Ni atoms in close contact with the TiC substrate indicating that these atoms are likely to be responsible for the catalytic activity reported for these systems. The analysis of the electronic structure reveals the existence of chemical interactions between the Ni nanoparticles and the TiC support, even if the net charge transfer between both systems is negligible.

Assessing the usefulness of transition metal carbides for hydrogenation reactions.

Transition Metal Carbides (TMCs) are proposed as replacements for and expensive late Transition Metals (TMs) as heterogeneous catalysts, often implying hydrogenation reactions or steps. Present density functional theory based calculations support using group IV TMCs and δ-MoC as viable TM alternatives, given the moderate exoergicity and affordable reaction step energy barriers.

Unexpectedly large impact of van der Waals interactions on the description of heterogeneously catalyzed reactions: the water gas shift reaction on Cu(321) as a case example.

The molecular mechanisms of the water gas shift reaction on Cu(321) have been chosen to investigate the effect of dispersion terms on the description of the energy profile and reaction rates. The present study based on periodic DFT calculations shows that including dispersion terms does not change the qualitative picture of the overall reaction, maintaining the rate determining step and the predominant route. However, the effect of dispersion is different for different adsorbates - reactants, intermediates or products - with a clear net effect and with no compensation of errors. Thus, in the OH + OH → H2O + O process the dispersion effects imply up to three orders of magnitude in the calculated reaction rates; the formation of carboxyl is highly disfavoured when dispersion terms are explicitly included and finally, the reaction rate for CO2 production (at 463 K) through cis-COOH dissociation is enhanced by three orders of magnitude by including dispersion terms in the calculation of the energy barrier. Consequently, the inclusion of dispersion terms largely affects the overall potential energy profile and produces tremendous changes in the predicted reaction rates. Therefore, dispersion terms must be included when aiming at obtaining information from macroscopic simulations employing for instance microkinetic or kinetic Monte Carlo approaches, where these effects should be clearly shown.

ReaxFF molecular dynamics simulations of CO collisions on an O-preadsorbed silica surface.

A quasiclassical trajectory dynamics study was performed for carbon monoxide collisions over an oxygen preadsorbed ß-cristobalite (001) surface. A reactive molecular force field (ReaxFF) was used to model the potential energy surface. The collisions were performed fixing several initial conditions: CO rovibrational states (v = 0-5 and j = 0, 20, 35), collision energies (0.05 ≤ E(col) ≤ 2.5 eV), incident angles (θ(v) = 0°, 45°) and surface temperatures (T(surf) = 300 K, 900 K). The principal elementary processes were the molecular reflection and the non-dissociative molecular adsorption. CO2 molecules were also formed in minor extension via an Eley-Rideal reaction although some of them were finally retained on the surface. The scattered CO molecules tend to be translationally colder and internally hotter (rotationally and vibrationally) than the initial ones. The present study supports that CO + O(ad) reaction should be less important than O + O(ad) reaction over silica for similar initial conditions of reactants, in agreement with experimental data.

Theoretical study of the dynamics and kinetics of the O + CS → CO + S chemical laser reaction, where CO shows a very high vibrational excitation.

The dynamics and kinetics of the O((3)P) + CS(X(1)Σ(+)) → CO(X(1)Σ(+)) + S((3)P) chemical laser reaction was studied theoretically in detail for the first time, as a function of collision energy (0.0388-2.0 eV) and rovibrational excitation of CS. This was made using the quasi-classical trajectory (QCT) method and employing the best ab initio analytical ground potential energy surface (1(3)A' PES) available. A broad set of properties was determined, including scalar and vector properties, and the reaction mode. The behaviors observed and the considerable formation of OCS collision complexes were interpreted from some characteristics of the PES (early barrier, shallow minimum in the exit channel, and high exoergicity (mainly channeled into CO vibration; up to ∼81% of the available energy)) and the kinematics. The QCT vibrational and rotational CO populations and the vector properties show a quite good agreement with experiments, but the QCT rate constants disagree. To better account for the kinetics, we performed CASPT2/aug-cc-pVTZ ab initio calculations on the stationary points along the minimum energy path of the ground and first excited (1(3)A'') PESs. The transition state theory, which can be satisfactorily applied here, leads to rate constants (100-2000 K) that are quite close to the measured ones, where comparison is possible (150-300 K). We expect that these results will encourage further theoretical and experimental developments.

Recombination and chemical energy accommodation coefficients from chemical dynamics simulations: O/O2 mixtures reacting over a ß-cristobalite (001) surface.

A microkinetic model is developed to study the reactivity of an O/O(2) gas mixture over a ß-cristobalite (001) surface. The thermal rate constants for the relevant elementary processes are either inferred from quasiclassical trajectory calculations or using some statistical approaches, resting on a recently developed interpolated multidimensional potential energy surface based on density functional theory. The kinetic model predicts a large molecular coverage at temperatures lower than 1000 K, in contrary to a large atomic coverage at higher temperatures. The computed atomic oxygen recombination coefficient, mainly involving atomic adsorption and Eley-Rideal recombination, is small and increases with temperature in the 700-1700 K range (0.01 < γ(O) < 0.02) in good agreement with experiments. In the same temperature range, the estimated chemical energy accommodation coefficient, the main contribution to which is the atomic adsorption process is almost constant and differs from unity (0.75 < ß(O) < 0.80).