QSRR is performed to accelerate the screening of five different stationary phases in HILIC method development. QSRR is established based on DFT computed molecular descriptors utilizing a GA-PLS for variable selection. GA descriptor sets are optimized to improve prediction accuracy. QSRR is used to obtain some insight into the HILIC mechanism.