A practical mixed basis set was developed to facilitate accurate calculations of potential energy surfaces for π-stacking interactions. Correlation consistent basis sets (cc-PVXZ) were augmented by p-type Gaussian functions placed above and below the planes of C6 moieties. Møller-Plesset (MP2, SCS-MP2) and coupled cluster [CCSD(T)] calculations show that such generated basis sets provide an accurate description of π-stacking systems with favorable computation times compared to the standard augmented basis sets. The addition of these off-center functions eliminates the linear dependence of the augmented basis sets, which is one of the most encountered numerical problems during calculation of the oligomers of polyaromatic hydrocarbons (PAH). In this work, we present a comparative study of the general characteristics of the potential energy surfaces for the parallel stacked and T-shape conformations of benzene and planar C6 clusters, using a combination of cc-PVXZ and our optimized functions. We discuss properties, such as the depth and curvature of the potential functions, short and long distance behavior, and the frictional forces between two model monomers.