-This work presents results of molecular simulations on adsorption of the four xylene isomers (o-, m-, p-xylene and ethylbenzene) in the porous zirconium terephthalate UiO-66. Grand-canonical Monte Carlo simulations (GCMC) are compared to multi-component adsorption equilibrium data obtained by breakthrough experiments. Additionally, the adsorption behavior of each isomer in order to better understand the multicomponent adsorption was evaluated. The simulations confirm that the experimentally observed ortho-selectivity is preferential in relation to the other isomers. Access to the small cages seems to be favored for o-xylene due to its more compact structure that provides more advantages in terms of rotation degrees of freedom. Furthermore, it was found that there is a competition between the other three isomers for adsorption in the UiO-66.