We managed to develop isocratic method for simultaneous analysis of triethylene glycol and cetylpyridinium chloride on Obelisc N mixed-mode HILIC column (4.6x150 mm). Glycol is retained by HILIC mechanism, cetylpyridinium by cation-exchange mechanism. Mobile phase 90% ACN with 5 mmol ammonium acetate pH 5, detection ELSD/LC/MS (no UV for glycols). Retention time 3.2 min for glycol and 5.5 min for cetylpyridinium chloride. I don't have ability to post chromatigrams here, but can send it by email.
The task was to analyze both compounds on one column with one set of conditions, there is no problem to retain neither benzalkonium chloride (C8-C18 chains) nor cetylpiridinium chloride (C18 quat). There are tons of methods for this type of compounds. Main problem is peak tailing with symmetry around 0.5-0.6. If you do gradient you can improve peak shape, but it is still going to tail due to residual silanol interaction (reference above has tailing even with gradient). If you shield silanols with basic group (how it is done on Primesep D column) you will have a perfect peak shape even with isocratic conditions, but none of RP columns (including Primesep D), will provide SIGNIFICANT retention for glycols.
http://www.sielc.com/application_151.html (you can use other buffers with this column, but peak shape will be the same)
So you need either two methods with GC for glyclos and Primesep D for benzalkonium, or mixed-mode HILIC on Obelisc N column. ELSD might not be accurate for glycols due to volatility (we used 30*C and linearity is not that great)
Other choices include derivatization of glycols with UV-active derivatizing agent which will provide additional hydrophobicity and UV activity. In this case you can use Primesep D for perfect peak shape.