Quantitative structure–property relationship (QSPR) modeling is a powerful approach for predicting environmental behavior of organic pollutants with their structure descriptors. This study reports an optimal QSPR model for estimating logarithmic soil sorption coefficients (log K OC ) of polycyclic aromatic hydrocarbons (PAHs). Quantum chemical descriptors computed using density functional theory at the B3LYP/6-31G(d) level and partial least squares (PLS) analysis with an optimizing procedure were used to generate QSPR models for log K OC of PAHs. The correlation coefficient of the optimal model was 0.993, and the results of a cross-validation test ([Formula: see text]) showed this optimal model had high fitting precision and good predicting ability. The log K OC values predicted by the optimal model are very close to those observed. The PLS analysis indicated that PAHs with larger electronic spatial extent tend to more easily adsorb and accumulate in soils and sediments, whereas those with higher molecular total energy and larger energy gap between the lowest unoccupied and the highest occupied molecular orbital adsorb and accumulate in soils and sediments less readily.