“…Notice that diffuse functions were only considered for the atoms in the hydrocarbon molecule and in the central region of the Na-ETS-10 cluster model, i.e., in the central Na cation (Na2), and in the framework oxygen atoms neighboring the Na2 cation (atoms labeled O1, O2, O5, O6, O7 and Na2 in Figure 1a). 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.…”