The role of the average pore diameter (APD) of 1.7μm AtlantisTM Premier BEHTM Particles derivatized with a zwitterionic group (propylsulfobetaine) on the efficiency of their 2.1 × 50 mm hydrophilic interaction liquid chromatography (HILIC) packed columns is investigated experimentally. Van Deemter plots for toluene (neutral, hydrophobic), cytosine (neutral, polar), tosylate (negatively charged), bretylium and atenolol (positively charged) were measured on three HILIC columns packed with BEH Z-HILIC Particles having APDs of 95 Å, 130 Å, and 300 Å. The intraparticle diffusivities of the analytes across these three BEH Z-HILIC Particles were measured by the peak parking method.The experimental data reveal that the slope of the C-branch of the van Deemter plots can be reduced by factors of about 15 for toluene, 2.5 for cytosine, 6 for atenolol, 5 for tosylate, and 14 for bretylium with increasing the APD from 95 Å to 300 Å. This observation is explained by: (1) the reduced amount of the highly viscous water diffuse layer and subsequent increase of the amount of acetonitrile-rich eluent in the mesopores, (2) the localized electrostatic adsorption of the retained analytes onto the zwitterion-bonded BEH Particles, and (3) depletion/excess of the analytes into the water diffuse layer. A general model of intraparticle diffusivity was then proposed to account for the impact of the APD of Z-HILIC Particles on the solid-to-liquid mass transfer resistance of small molecules. The model highlights the relevance of the thickness of the water diffuse layer, the access of the bulk eluent into the mesopore, the localized electrostatic adsorption, and the partitioning constant of the retained analyte between the bulk eluent and the water diffuse layer.