In this study,DFT-B3LYP level of theory with the 6-311G~(**)basis set was used to calculate a set of molecular descriptom of 36 nitrobenzenes.Quantitative structure-activity relationship(QSAR)models of the toxicity(plGC_(50))of 36 nitrobenzenes were established using some of the following calculated descriptors:HOMO energy(E_(HOMO)),LUMO energy(E_(LUMO)),the difference between LUMO and HOMO energy values(△E=E_(LUMO)-E_(HOMO)),average polarizability (P),dipole moment(μ),molecular volume(ν)and the charge of nitro group(Q_(-NO2)).The QSAR models obtained by employing stepwise multiple regression techniques are aimed at correlating the structures with their reported experimental toxicity values.Among the models presented in this study,statistically the most significant one is a two-parameter linear equation with the correlation coefficient (R~2) of 0.896 and cross-validated correlation coefficient(R~2_(cv)) of 0.878.To further validate the predictive ability of the resulting model,external validation was carried out with R~2_ (ext) and Q~2_(ext) values of 0.894,0.900,0.909,0.874,0.885 and 0.903,respectively.This study also reconstructed the obtained models using AMl-based calculated descriptors to demonstrate the superiority of DFT over AMl for quantum calculations of mechanical descriptors.The results were discussed in the light of the main factors that influence the toxicity of nitrobenznenes.