Multiscale simulation of liquid chromatography: Effective diffusion in macro–mesoporous beds and the B-term of the plate height equation

We performed multiscale simulations of analyte sorption and diffusion in hierarchical porosity models of monolithic silica columns for reversed-phase liquid chromatography to investigate how the mean mesopore size of the chromatographic bed and the analyte-specific interaction with the chromatographic interface influence the analyte diffusivity at various length scales. The reproduced experimental conditions comprised the retention of six analyte compounds of low to moderate solute polarity on a silica-based, endcapped, C₁₈ stationary phase with water‒acetonitrile and water–methanol mobile phases whose elution strength was varied via the volumetric solvent ratio. Detailed information about the analyte-specific interfacial dynamics received from molecular dynamics simulations was incorporated through appropriate linker schemes into Brownian dynamics diffusion simulations in three hierarchical porosity models received from physical reconstructions of silica monoliths with a mean macropore size of 1.23 µm and mean mesopore sizes of 12.3, 21.3, or 25.7 nm. The mean mesopore size was found to have a similar influence on the effective mesopore diffusivity as the analyte polarity and the mobile-phase elution strength, which together determine the analyte residence time on a column. A smaller mesopore size attenuated the increase of the effective mesopore diffusivity with increasing mobile-phase elution strength significantly. The effective bed diffusivity was limited by the analyte residence time rather than by morphological details of the mesopore space. The stronger an analyte was retained by the chromatographic interface inside the mesopores, the slower became the mass transfer between the pore space hierarchies and the lower was the effective bed diffusivity. The B-terms of the plate height equation were finally generated with the bed diffusivities and phase-based retention factors derived from the hierarchical porosity models using additional information about the stationary-phase limit obtained from the analysis of analyte–bonded phase contacts. The B-terms highlight analyte- and mobile phase-specific behavior relevant to isocratic and gradient elution conditions in chromatographic practice. In particular, U-shaped B-term curves are observed due to the dominating contribution of the retention factor and the bed diffusivity to the B-term at low and high elution strength of the mobile phase, respectively.