Optimal binding site of a methane molecule on the silanol covered (010) surface of silicalite-1: ONIOM calculations

The binding energies and the corresponding structures of a methane molecule on the silanol covered (010) surface of silicalite-1 have been investigated using ab initio methods. Different levels of calculations, HF/6-31G(d), MP2/6-31G(d) and ONIOM (MP2/6-31G(d):HF/6-31G(d)) including the correction of an error due to an unbalance of the basis set, known as basis set super position error (BSSE), as well as the size of the cluster representing the silicalite-1 surface, were systematically examined to validate the model used. The ONIOM method with the BSSE correction was found to be a compromise between accuracy and computer time required. The optimal binding site on the silicalite-1 surface was observed at the configuration where the methane molecule points one H atom toward the O atom of the silanol group. The corresponding binding energy is -1.71 kJ/mol. This value is significantly higher than that of -5.65 kJ/mol when the methane molecule approaches the center of the straight channel. At this configuration, the C atom of methane was observed to locate exactly at the center of the channel. This leads to the conclusion that the methane molecule will relatively seldom be adsorbed on the silanol covered (010) surface of silicalite-1. Instead, the adsorption process will take place directly at the center of the straight channel.