zyj19831206 发表于 2015-11-6 19:29
这个办法好，看来小卒那个用dmol3预先优化是对的，问题是dmol3优化出来的虚频太多了，不知道该怎么办？

至于dmol3能优化出来这个反应的过渡态我还是有一定质疑的。

如果以它为初始结构用高斯能找到过渡态，一般只有一个虚频，如果虚频不止一个，再用int=ultrafine结合tight的opt收敛限，还不行再用calcall。

sobereva 发表于 2015-11-6 19:24
如果dmol3确实能找到过渡态，拿它的过渡态作为初猜用opt(ts,calcfc,noeigen)来找
能用TS解决的就不建议用Q ...

这个办法好，看来小卒那个用dmol3预先优化是对的，问题是dmol3优化出来的虚频太多了，不知道该怎么办？

如果dmol3确实能找到过渡态，拿它的过渡态作为初猜用opt(ts,calcfc,noeigen)来找
能用TS解决的就不建议用QST2

sobereva 发表于 2015-11-6 18:09
确信反应是这样进行的么？看结构，这样的反应不像是有过渡态的样子...

需要酸催化剂，根据实验也就是估计后期加入一些氢原子，但是暂时我想先这么算算，但是dmol 貌似能算出来啊？为啥同样的gaussian 反而还不行？泛函多应该可选性强，更容易收敛啊？

确信反应是这样进行的么？看结构，这样的反应不像是有过渡态的样子...

sobereva 发表于 2015-11-6 15:32
这是多聚甲醛？输入文件上看不是啊。

抱歉，可能输入文件建的有问题，已经重新上传。

这是多聚甲醛？输入文件上看不是啊。