A computational study of the role of cobalt in thiophene adsorption on small Mo and MoCo clusters as site models for the HDS process

Non periodic density functional theory calculations are used to investigate the role of cobalt atoms in the adsorption of thiophene on small Mo and MoCo clusters. Metallic aggregates play the role of those active sites found in the true catalysts. Two interaction modes between thiophene and metallic sites are considered, namely, the S-mode, in which the organosulfur molecule interacts through the S atom, and the R-mode, in which the interaction takes place through the thiophene ring. A large number of sites, in which thiophene effectively adsorbs, was found, both in the monometallic case and in the bimetallic one. Considerably larger adsorption energies were found when thiophene interacts via the R-mode than when adsorption occurs through the S-mode. The activation of C-S bonds is also more important for R-mode cases than for S-mode ones. Further analysis made on some selected systems and based on density of states and molecular orbital overlap population-projected density of states reveals that thiophene and metallic clusters interact in an energy range around −6.0 eV with respect to the Fermi energy. Bands observed at energies below −6.0 eV correspond to thiophene states that become shifted with respect to the values obtained for isolated thiophene depending on the strength of the interaction. Bands above -6.0 eV describe how C and S atoms interact with Co and Mo ones, providing both bonding and antibonding patterns that helps to understand the overall interaction. Most important is the finding that cobalt atoms seem to play no relevant role during the adsorption of thiophene on metallic sites. Thus, present results obtained using non periodic GGA density functional theory seem to point to cobalt taking part in another step of the overall HDS process, hydrogen adsorption or hydrogen attack to C-S bonds, for instance.