Nitrile rubber (NBR, 39:61 wt% of acrylonitrile:butadiene) and styrene butadiene rubber (SBR, 50:50 wt% of styrene: butadiene) matrices have been equilibrated by molecular dynamics (MD) simulations. Transition-state approach is used to calculate the diffusion and solubility coefficients of small penetrants in these matrices, indicating quite low values in NBR and reasonable agreement with experimental results. MD simulations have been performed to analyze water diffusion in these matrices. Aggregation of water molecules is observed in the hydrophobic matrix SBR. MD simulations with fictitious nonpolar water molecules inhibit aggregation and lead to enhanced diffusion in SBR. In NBR there is a slight increase in diffusion for fictitious water molecules. The lower diffusion constants in NBR result from slower local relaxation of the matrix due to tighter intermolecular packing and higher cohesive energy density. The free volume distribution that affects solubility coefficients is not a major determining factor for the diffusion coefficients in these matrices. (c) 2006 Elsevier Ltd. All rights reserved.