A QSRR study was performed to develop a predictive model correlating the observed kovats retention indices of three congeneric aromatic series with their molecular structures in gas chromatography. The congeneric aromatic series included the substituted benzene, benzaldehyde and acethophenone compounds, which had been studied, previously, on six OV stationary phases with different phenyl percentages. At first, a model was generated for six columns separately, using only calculated descriptors and MLR technique. Then a combined model, added a polarity term of stationary phase (M), was also developed for all these columns, and the result was apparently satisfactory (R 2 ¼ 0.991, F ¼ 1009.828, SE ¼ 23.98). Since the intercept had a high value in this model, the neural network back propagation algorithm was applied for comparison, and it was found that the neural network could exceed the level of the multiple regression method. The stability and validity of both models were tested by cross-validation technique and by predicted response values for the prediction set and test set. The results of the study indicated that a seven parameter equation can be utilized for prediction of retention indices of compounds on different OV stationary phases simultaneously using ANN technique, for which there are no empirical RI values.