计算化学公社

标题: 求助怎么可以扫出非最稳定态的势能面 [打印本页]

作者
Author: Dwppppp 时间: 2017-6-16 01:21
标题: 求助怎么可以扫出非最稳定态的势能面
小弟刚刚入门计算，正在计算一个分子NMP（N-methylpiperidine) 的势能面（在不同twist的程度下）。
这个分子有船式（boat）和椅式（chair）结构，我在扫基态和离子态（正一价）的时候，都发生了在势能面某处出现了尖点和突变为另一条线。
我用gaussview观察得知是分子由椅式结构突变为船式结构或者反过来突变所造成的结果。
那么我就想得到两种形式的完整的势能面，怎么能让gaussian不去自动优化为能量最低的结构，而是沿着船式或者椅式结构一直算下去呢？

谢谢大家！求教！

作者
Author: sobereva 时间: 2017-6-16 01:48
opt=modredundant，恰当在坐标末尾后面增加对某些几何变量的冻结设定，使得船式和椅式之间没法变换

作者
Author: Dwppppp 时间: 2017-6-16 01:57
谢谢sob老师的回复！
我可以这样做，但是我并不知道在固定条件下那个几何变量的最稳定值是多少，也就意味着我通过这样方式找到的势能面，可能和真正relax下的势能面有一些微小的偏差。这个是不是在gaussian里没办法避免了？
谢谢！

作者
Author: wbn 时间: 2017-6-16 02:14

Dwppppp 发表于 2017-6-16 01:57
谢谢sob老师的回复！
我可以这样做，但是我并不知道在固定条件下那个几何变量的最稳定值是多少，也就意味 ...

用之前得到的结构松开固定的变量，只固定扫描的变量，调低步长再优化一下就行了

作者
Author: Dwppppp 时间: 2017-6-16 02:25

wbn 发表于 2017-6-16 02:14
用之前得到的结构松开固定的变量，只固定扫描的变量，调低步长再优化一下就行了

谢谢老师
所以说调低步长确实会让系统找相对接近的结构来进行计算，我试一下

作者
Author: liyuanhe211 时间: 2017-6-16 02:28
分段解决，自某点翻成椅式就先解决那个点自己的问题(适当初猜，小步长，限制性优化这个点)，解决之后以它为初猜做其后的点，最后把曲线拼接起来就可以了。如果这个点怎么也无法找到想要的构象，那它可能不是这个子空间的极小点，而是过渡态之类的，那就无法得到一个柔性扫描曲线都是船式
。

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