The purpose of the present work was to systematically study the chromatographic behaviour of different aromatic stationary phases in a subcritical fluid mobile phase. We attempted to assess the chemical origin of the differences in retention characteristics between the different columns. Various types of aromatic stationary phases, all commercially available, were investigated. The effect of the nature of the aromatic bonding on interactions between solute and stationary phases and between solute and carbon dioxide-methanol mobile phase was studied by the use of a linear solvation energy relationship (LSER): the solvation parameter model. This study was performed to provide a greater knowledge of the properties of these phases in subcritical fluid chromatography, and to allow a more rapid and efficient choice of aromatic stationary phase in regard of the chemical nature of the solutes to be separated. Charge transfer interactions naturally contribute to the retention on all these stationary phases but are completed by various other types of interactions, depending on the nature of the aromatic group. The solvation vectors were used to compare the different phase properties. In particular, the similarities in the chromatographic behaviour of porous graphitic carbon (PGC), polystyrene-divinylbenzene (PS-DVB) and aromatic-bonded silica stationary phases are evidenced.