The adsorption of short linear alkanes has been explored in the highly flexible MIL-53(Cr) porous metal-organic framework by means of molecular simulations based on configurational bias grand canonical Monte Carlo. The unusual shape of the adsorption isotherms with the existence of steps has been successfully modelled by creating a (narrow pore, large pore) phase mixture domain, the composition of which varies with pressure. A further step consisted of combining our computational approach with several experimental tools including microcalorimetry, gravimetry and in situ X-ray diffraction, to fully characterize the adsorption behaviour of the isostructural MIL-47(V) rigid MOF, i.e. the preferential arrangement of each type of alkane inside the pores and the resulting interaction energy. Finally, relationships are established between the adsorption enthalpies and both alkyl chain length and polarisability of the alkanes that can be further utilised to predict the energetics of the adsorption process for longer alkane chains.

International audienceThe adsorption properties of linear long chain alkanes (from n-pentane to n-nonane) within the rigid MOF MIL-47 (V) have been explored by combining gravimetry measurements and molecular simulations. Both experimental absolute isotherms and enthalpies of adsorption for all n-alkanes were compared with those obtained by configurational bias grand canonical Monte Carlo simulations (CB-GCMC) based on two different force fields. From a fair agreement between experimental and simulated data, a further step consisted of investigating the microscopic adsorption mechanism in play to shed some light onto the preferential orientations and conformations of all investigated n-alkanes. Whereas the trans conformation is predominantly observed for all n-alkanes, the proportion of the n-alkane conformations lying parallel to the direction of the tunnel significantly increases with the chain length, emphasizing that the confinement effect is stronger for the longer chain n-alkanes. Finally, molecular dynamics simulations allowed us to emphasize that all n-alkanes follow a pathway along the direction of the tunnel, leading to a 1D type diffusion mechanism, the motions being mainly centered around the middle of the pores at low loading, whereas they are significantly shifted toward the pore wall when the alkane concentration increases