An attempt is presented to infer the stacking interaction energy from the atom in molecules (AIM) properties computed in all critical points of electron density found between interacting molecules. In this study, benzene and pyridine homodimers in coaxial and parallel displaced configurations were used as model structures. We classify the structures and briefly discuss their peculiarities in the AIM theory framework. The maximum absolute error of 2 kcal/mol achieved for 300 potential energy surface point set indicates that a more sophisticated algorithm for estimating atom–atom interactions is needed still. © 2012 Wiley Periodicals, Inc.