First‐ and second‐dimension retention times for a series of alkyl phosphates were predicted for multiple column combinations in GC×GC. This was accomplished through the use of a three‐parameter thermodynamic model where the analytes’ interactions with the stationary phases in both dimensions are known. Ionic liquid columns were employed to impart unique selectivity for alkyl phosphates, and it was determined that for alkyl phosphate compounds, ionic liquid columns are best used in the primary dimension. Retention coordinates for unknown phosphates are predicted from the thermodynamic parameters of a set standard alkyl phosphates. Additionally, we present changing retention properties of alkyl phosphates on some ionic liquid columns, due to suspected reaction between the analyte and column. This makes it difficult to accurately predict their retention properties, and in general poses a problem for ionic liquid columns with these types of analytes.