如果只考虑一个阴-阳离子复合物构型，严格来说需要用molclus（http://www.keinsci.com/research/molclus.html）之类程序做构型搜索，找到结合后能量最低结构。
更简单的做法是先计算阳离子，并用Multiwfn计算阳离子表面静电势（参考下文），在静电势最正的位点附近合适位置放一个阴离子，然后优化，也大概率能得到能量最低结构。如果几何优化的级别就够，HOMO-LUMO gap直接读就完了。
使用Multiwfn的定量分子表面分析功能预测反应位点、分析分子间相互作用
http://sobereva.com/159
使用Multiwfn+VMD快速地绘制静电势着色的分子范德华表面图和分子间穿透图（含视频演示）
http://sobereva.com/443（http://bbs.keinsci.com/thread-11080-1-1.html）