Pi electron stacking, how does it work?
You going along the right track. The stacking attraction occurs when two molecules with $\pi$ orbitals come face to face with one another, typically their separation is 0.34 nm. However, for the interaction to be attractive the two molecules have to be displaced slightly, this prevents direct electron repulsion between $\pi$ electrons on one molecule with those on the other. The attraction is in general between the $\pi$ electrons on one molecule and the $\sigma$ framework on the other$^1$ and vice - versa. 

The model assumes that a $\pi$ orbital can be split into three charged parts, separated by a distance $\delta$. The lobes of the orbital have charge -1/2 each and at the position of the atom's nuclei is placed a charge of +1. Thus there are also three interactions $\sigma - \sigma , \pi - \sigma $ and $ \pi - \pi$ and three positions of charge, making 9 interactions in all for each atom. The interaction is taken to be electrostatic, i.e. Coulomb and each term has the form $E_{ij} = q_iq_j/r_{ij} $ for atoms i and j and is summed over all atoms. 

In addition a van-der-Waals term is added to the calculation of the total interaction, this has the form $E_{vdw} = \Sigma _{ij}(A_{ij}exp(-\alpha_{ij}R_{ij})-C_{ij}/r_{ij}^r))$ where A, $\alpha$ and C have standard values$^2$. In solution the value of this energy is much reduced vs. its value in the crystal making the $\pi-\pi$ interaction the dominant one.

The charges are shown schematically below for carbon atoms. For other atoms, e.g. nitrogen the charge is given as 1.5 with 1.5 electrons associated with them. (All images are from ref 1 )



charge distribution



optimal geometry for Zn porphyrin $\pi-\pi$ interaction. The molecules are displaced (in x and y) relative to one another at a separation of 0.34 nm. 



The contour plot for $\pi-\pi$ interaction and van-der-waals (vdw) energy as a function of their centre to centre offset as shown above (energies in $\pu{ kJ mol^{-1}}$). The black squares are geometries from $\pi-\pi $ stacking in crystals. The shape of the contour plot with or without the vdw interaction is the same as with just the $\pi-\pi$ interaction alone, i.e this interaction determines geometry but  the vdw generally lowers the energy and more clearly show regions of low energy as negative values.

$^1$ The model was developed by Hunter &amp; Sanders J.Am. Chem. Soc. 1990, v112, p 5525-5534

$^2$ Van der waals parameters Caillet &amp; Claverie, Acta Crystallogr. sect A1975, v31, p448