计算化学公社

标题: 何种情况下过渡态两个虚频投稿可以接受呢？ [打印本页]

作者
Author: airyang 时间: 2017-3-24 23:47
标题: 何种情况下过渡态两个虚频投稿可以接受呢？
我计算的自交换TS应该有个C2对称性：
保证C2对称性的情况下总是无法消除掉一个很小的虚频（-8.15，对应了两个片段之间的扭动，对应反应历程的虚频是-1076.4）
如果破坏掉对称性可以得到单虚频（-1147.34），C1结构相对自由能也较C2稍低1.3kcal。两个结构电子能差别不大（0.3kcal），1.3kcal的自由能差别主要来源于熵差1.4kcal。

印象看到过Truhlar有篇计算，在SI中说明了有几个TS有第两个小的虚频，也就那么用了，不过印象只是定性计算。

想请问，我的情况可以通过声明第二个虚频很小，而用C2结构么？在意最终的活化能数值。

作者
Author: sobereva 时间: 2017-3-25 00:56
如果从原理上确定无疑就是应该有C2对称性，则如果是DFT计算用int=ultrafine再试，如果还有小虚频再calcall，如果还不行把基组再进一步提升（到def2-TZVP或cc-pVTZ这个档为止）。
即便是小虚频也应当尽量消掉，要不然审稿人总会吐槽一下

作者
Author: beefly 时间: 2017-3-25 07:34
非阿贝尔群的高对称体系，是可以出现高阶过渡态的。比如D5h对称性的C5H5+就是二阶过渡态，有两个相等的虚频，其中一个导致极小点，另一个导致一阶过渡态。也就是Jahn-Teller效应

作者
Author: airyang 时间: 2017-3-25 12:11

beefly 发表于 2017-3-25 07:34
非阿贝尔群的高对称体系，是可以出现高阶过渡态的。比如D5h对称性的C5H5+就是二阶过渡态，有两个相等的虚频 ...

谢谢老师~

作者
Author: airyang 时间: 2017-3-25 12:20

sobereva 发表于 2017-3-25 00:56
如果从原理上确定无疑就是应该有C2对称性，则如果是DFT计算用int=ultrafine再试，如果还有小虚频再calcall ...

目前的优化用的是ultrafine，M062X-D3/6-31+Gdp, 有机时了试下calcall或者提升基组

不过，目前那个不想要的小虚频对应着破坏C2对称性的摆动，calcall会类似手动displacement直接掉到C1对称性么。。。

谢谢sob老师~

作者
Author: sobereva 时间: 2017-3-25 23:45

airyang 发表于 2017-3-25 12:20
目前的优化用的是ultrafine，M062X-D3/6-31+Gdp, 有机时了试下calcall或者提升基组

不过，目前那个不 ...

calcall通常没有此效果。破坏对称性需要自己按照虚频来移动

作者
Author: airyang 时间: 2017-3-26 03:33
soga，算起来~

作者
Author: zjxitcc 时间: 2017-3-26 16:10

beefly 发表于 2017-3-25 07:34
非阿贝尔群的高对称体系，是可以出现高阶过渡态的。比如D5h对称性的C5H5+就是二阶过渡态，有两个相等的虚频 ...

长知识了，请问有没这方面（高对称体系，高阶鞍点）的理论文献或书籍？

欢迎光临计算化学公社 (http://bbs.keinsci.com/)