Molecular simulation of the adsorption and structure of benzene confined in mesoporous silicas

Grand Canonical Monte Carlo simulations are used to study the adsorption of benzene at 298 K in an atomistic cylindrical silica nanopore of a diameter 3.6 nm. The adsorption involves a transition from a partially filled pore (a two layers thick film at the pore surface) to a completely filled pore configuration. Strong layering of the benzene molecules at the pore surface is observed. It is found that the layering decays as the distance to the pore surface increases. The position of the peaks for the density of the C, H atoms and the center of mass of the molecules shows that benzene molecules prefer an orientation in which their ring is perpendicular to the pore surface. This result is corroborated by calculating orientational order parameters and examining the distribution of the distances between the H and C atoms of the benzene molecules and the H and O atoms of the silica substrate.