Geometrical optimization and electrostatic potential calculations have been performed at the HF/6-31G* level of theory for investigated persistent organic pollutants (POPs). A number of statistically based parameters have been obtained. Relationship between soot-water partition coefficients (log K-SC) of POPs and the structural descriptors has been established by the multiple linear regression method. The result shows that the quantities derived from electrostatic potential <(V-s(-))over bar> and V-s,V-max, together with molecular surface area (A(S)) and the energy of the highest occupied molecular orbital (E-HOMO) can be well used to express the quantitative relationship between structure and log K-SC (QSPR) of POPs. Predictive capability of the model has been demonstrated by leave-one-out cross-validation with the cross-validated correlation coefficient of 0.9797. Furthermore, the predictive power of this model was further examined for the external test set with the correlation coefficient of 0.9811 between observed and predicted log K-SC, validating the robustness and good predictive ability of our model. Furthermore, in order to further investigate the applicability of these parameters derived from electrostatic potential in prediction of soot-water partition coefficient for organic pollutants, eleven polycyclic aromatic hydrocarbons (PAHs), eleven polychlorinated biphenyls (PCBs) and nine phenyl urea herbicides (PUHs) from other source have also been studied. The QSPR models established may provide a new powerful method for predicting soot-water partition coefficients (log K-SC) of organic pollutants. (C) 2012 Elsevier Inc. All rights reserved.