用Gaussian算这个属于瞎折腾，既不适合计算，而且计算那么大的簇也非常昂贵。直接用免费又非常高效的CP2K做周期性体系的几何优化、振动分析、分子动力学才是首选。

通过北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）可以一次性学明白，而且会发现CP2K结合Multiwfn和gview使用做这些计算相当简单。

换第一性原理软件比如CP2K。