刚接触计算化学计算过的一个分子，镁卟啉：
我分别用B3LYP,BHandHLYP，PBE0优化基态和激发态结构
对优化基态就没什么说的，很好的D4h点群结构。
优化激发态，D2h对称，有文献佐证：在turbomole里分别用这三种泛函优化，定义激发态分子对称性为D2h，只有BHandHLYP优化出来，算频率没有虚频，B3LYP，PBE0优化出来算频率有虚频。B3LYP很诡异，去年算的那会没有虚频，现在算又有虚频，对照control文件没有任何差别。这点我始终不明白，
用高斯优化，我取基态优化好的坐标，用CIS，TD方法优化，这样出来的结构高斯判断是C2h，并且都没有虚频。但我放到turbomole下判断是D2h对称
接着，我就取turbomole里面B3LYP，PBE0优化好的激发态结构（turbomole算有虚频），放到高斯里去做频率，没虚频！！！ Darn it !
，对照control文件没有任何差别。这点我始终不明白，
Darn it ! 
回复此楼
+1/673
