Diffusion of p-xylene in silicalite is simulated by using the Buckingham exp-6 model. The best agreement with experimentally determined activation energies for the diffusion (25–30 kJ/mol) is obtained when the framework atoms are relaxed (24±2 kJ/mol). Study of the interaction energy versus the rotation angle of p-xylene along its long molecular axis (010 direction) shows the existence of a single minimum corresponding to a stable orientation of the molecule at the channel intersection (site I). The importance of including the Si-H interactions between the sorbent and the sorbate is emphasized : neglecting these interactions (the widely used Lennard-Jones model, where “ the crucial contact is the O…H distance”) yields overestimated sorption energies and two possible orientations for the p-xylene molecule at site I. A significant simplification of the computer simulations might be achieved by considering only the Si⋯hydrocarbon interactons : in that case the calculated energies have to be multiplied by ≈2.5. Accordingly, as already shown by 29Si CP-MAS nmr, the crucial contact seems rather to be Si…H.