by combining sol-gel synthesis with a proper surfactant templating agent. The latter was not a single solvated organic or inorganic species, but a self-assembled surfactant array. Three main mesophases are typically produced, namely, the cubic MCM-48, the lamellar MCM-50, and the hexagonal MCM-41. The latter is characterized by a regular 2D hexagonal array (P6mm space group symmetry) of highly uniform parallel non-interconnected channels, with the pore diameter depending on the specific templating agent. Compared to the other members of the M41S mesoporous family, MCM-41 has a higher thermal stability, and it is easier to prepare, thus attracting greater attention. [1,12] The ordered pore architecture of MCM-41 can be clearly observed through transmission electron microscopy (TEM) [2,13] and assessed through low-angle X-ray diffraction (LA-XRD) techniques or small-angle x-ray scattering (SAXS). [2] Overall, porosity can be as high as 80%, with the pore walls composed of amorphous silica. [1] The analysis of the N 2 adsorption isotherms according to the Brunauer-Emmett-Teller (BET) theory [14] is widely employed to estimate the Specific Surface Area (SSA), which is typically >700 m 2 g -1 . [11] Analysis of the same adsorption isotherm through the Barrett-Joyner-Halenda (BJH) model [15] is the typical choice to determine the pore size distribution (PSD). However, other approaches are available to retrieve these textural details, which are more recent and have a more rigorous physical base, such as the Density Functional Theory (DFT), allowing the determination of both the SSA, the pore volume, and the PSD. [16] The former models and theories have been recently questioned. [17] Computer calculations and simulations are able to unveil the microscopic details and the energetic aspects of the interaction between the guest molecules and the solid surface. Adsorption capacity, selectivity, diffusivity, and so on, can be investigated at the microscopic level under different pressure, temperature and humidity, so to obtain fundamental information to tailor material synthesis and functionalization to meet specific goals. Computational approaches can be used to complement the experiments, but also to possibly predict the system behavior, and to save time and resources on the laboratory bench. To this aim, the preparation of a 3D model that matches the characteristics of the real sample and that is able to reproduce the macroscopic behavior, as seen by the experiments, is crucial. Mesoporous materials cannot be modeled as simple flat surfaces or slits. At Among mesoporous ordered silica-based materials, MCM-41 is widely employed in both experimental and computational investigations of gas adsorption. Specific surface area and pore size distribution are typically obtained by analyzing the N 2 adsorption isotherm. To this aim, different models and theories are available, and which one is more accurate is under debate. On the computational side, the in silico model ought to match the characteristics of real samples and reproduc...