A series of quantitative structure-retention relationship models were developed to predict gas chromatographic relative retention times （GC-RRTs） for 209 polybrominated diphenyl ether （PBDE） congeners on 10 stationary phases. A genetic algorithm with twofold leave-multiple-out cross validation （LMOCV） was used to select optimal subsets from large-size molecular descriptors. Overall multiple-linear regression fitting correlation coefficients （R2） are greater than 0.988, except for the CP-Sil 19 colunm, in which Q^uocv （correlation coefficient of LMOCV）, Q^oocv （correlation coefficient of leave-one-out cross validation, LOOCV）, and Rp2re （correlation coefficients of prediction set） are larger than 0.98. The excellent statistical parameters reveal that the models are robust and have high internal and external predictive capability. According to the descriptors for constructing the models, the GC-RRTs in various stationary phases are highly dependent on distances among atoms, branches of molecules, and molecular properties. PBDE congeners with 1, 9, and 10 bromines are major outliers based on the results of the application domain.
