较大体系多个分子结构优化和频率分析的基组选择问题

tym1999

想用Gaussian计算多个分子，如下图所示，（gjf文件如下，大约一共有170个原子，分子元素为C、H、O)的opt+fre,之后用于计算结合能和静电势的分析,使用的超算平台，核数为256，尝试过用b3lyp/6-311g(d,p)，但是耗时很长，而且算出来会有非常多的小虚频，很难消除，所以想问问这种体系可以用Gaussian算吗，如果可以，应该选用什么基组呢？有没有其他更好的软件计算这种体系呢？

wzkchem5

可以考虑用ORCA的B97-3c等方法，开RI以后计算很快。此外注意构象搜索问题

sobereva

不要把freq写成fre

先用GFN-xTB或者PM6-D3H4做预优化，然后b3lyp-D3(BJ)/6-31G*优化，当前计算条件算这体系根本没难度
用B3LYP不带D3优化这个是开玩笑，好好看
谈谈“计算时是否需要加DFT-D3色散校正？”
http://sobereva.com/413（http://bbs.keinsci.com/thread-9772-1-1.html）
谈谈量子化学研究中什么时候用B3LYP泛函优化几何结构是适当的
http://sobereva.com/557（http://bbs.keinsci.com/thread-17899-1-1.html）

虚频该怎么消怎么消。本来弱相互作用体系就容易出虚频，我推荐的这个级别已经是属于最不容易出虚频的级别范畴里的了，换其它的也不会有任何好处。
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278（http://bbs.keinsci.com/thread-633-1-1.html）


原则上来说需要结合molclus（http://www.keinsci.com/research/molclus.html）等程序做构型/构象搜索。至于实际是否需要，看具体研究目的，如果构象搜索做不做对研究的问题影响不大可以忽略此问题