“…The M06-L functional is used in this work since in a recent comparative study on the performance of nine different DFT functionals to calculate the interaction enthalpy of methane with a cluster model of Na-ETS-10, the M06-L approach was found to provide the best results [21].Besides,the M06-L approach together with a cluster model was found to provide very accurate energetic, vibrational and structural information for some other difficult systems such as i) adsorption of CO on MgO(001) [30], ii) the adsorption of CO and NO on Ni-doped MgO(001) [31], and iii) the interaction of water with open metal sites of CuBTC [32]. Very encouragingly, recent works concerning the interaction of CO with the (111) surfaces of Rh, Pt, Cu, Ag and Pd [33] or the interaction of benzene, pyridine, thymine and cytosine with Au(111) [34]showed that the M06-L method is suitable for treating dispersion dominated interactionsof adsorbateswith metallic surfaces [35].The M06-L approach was also effectively employed to predict reaction energies in organometallic systems [36][37][38] and the benchmark study on the interaction of hydrocarbons and smaller gases with several zeolite clusters byZhao and Truhlar [39]shows that the M06 family of functionals is well suited to study this kind of systems.Additionally, in our very recent work on the interaction of five atmospheric gases with Na-ETS-10 [21], calculated enthalpic data for methane, nitrogen and, at some extent, carbon dioxide, and calculated vibrational data for methane, nitrogen, hydrogen, carbon dioxide and water were found to be in very satisfactory agreement with experimental isosteric heats [40,41] and infrared results [42][43][44][45], respectively. The same 2Ti cluster model used in Ref.…”