Impact sensitivity, one of the most important screening factors for novel high energy density materials (HEDMs), was predicted by use of quantitative structure-property relationship (QSPR) based on the electrostatic potential (ESP) values calculated on the van der Waals molecular surface (MSEP). Among various 3D descriptors derived from MSEP, we utilized total and positive variance of MSEP, and devised a new QSPR equation by combining three other parameters. We employed 37 HEDMs bearing a benzene scaffold and nitro substituents, which were also utilized by Rice and Hare. All the molecular structures were optimized at the B3LYP/6-31G(d) level of theory and confirmed as minima by the frequency calculations. Our new QSPR equation provided a good result to predict the impact sensitivities of the molecules in the training set including zwitterionic molecules.