计算化学公社

标题: OPT优化过程的问题 [打印本页]

作者
Author: 小书童 时间: 2016-8-19 11:49
标题: OPT优化过程的问题
做一个结构优化，初始结构和最终优化结构相差巨大，初始结构的力已经接近0了，但是优化过程中又不断远离0，最终经过120多步才走到收敛的机构，总觉得我的优化过程是有问题的，但是又说不上来，麻烦各位老师给看看。此外，对于力在很小范围内振动的情况怎么解决？麻烦老师指点一下，

作者
Author: sobereva 时间: 2016-8-19 15:40
步长上限设3，calcall再试试

PS：即便能找到和初猜结构很相近的极小点，当前情况也说明这个极小点很不稳定，实际中跑到这个结构后立即就会跑到你之前优化出来的结构去。当前情况也可能是势能面二分所导致的。

作者
Author: 小书童 时间: 2016-8-23 22:11

sobereva 发表于 2016-8-19 15:40
步长上限设3，calcall再试试

PS：即便能找到和初猜结构很相近的极小点，当前情况也说明这个极小点很不稳 ...

opt=(calcall,maxstep=3,notrust）,优化出来的结构仍然是不是IRC连接的一端，仍然偏离。走IRC的时候我又采用了一个特别大步长的方法走，似乎是可以走到现在OPT的结构，但是现在就又矛盾了，我的TS对应的振动方式并没有指向这个变化方向。现在茫然不知道该如何做了，麻烦sob老师再指点一下，多谢老师！

作者
Author: sobereva 时间: 2016-8-24 10:29

小书童发表于 2016-8-23 22:11
opt=(calcall,maxstep=3,notrust）,优化出来的结构仍然是不是IRC连接的一端，仍然偏离。走IRC的时候我又 ...

估计就是势能面分叉
http://bbs.keinsci.com/forum.php?mod=viewthread&tid=3066

作者
Author: liyuanhe211 时间: 2016-8-24 11:08

sobereva 发表于 2016-8-24 10:29
估计就是势能面分叉
http://bbs.keinsci.com/forum.php?mod=viewthread&tid=3066

是如何判断出这种情况是势能面二分造成的呢？

作者
Author: 小书童 时间: 2016-8-24 14:44
分子表面不同的静电势，有负值，也有正值，静电势为正是说明该处的原子带正电荷吗？

作者
Author: sobereva 时间: 2016-8-24 17:36

小书童发表于 2016-8-24 14:44
分子表面不同的静电势，有负值，也有正值，静电势为正是说明该处的原子带正电荷吗？

没法这么说。只能说正值地方是核电荷的贡献超过电子的贡献，负值则是反过来。
而且原子电荷很多情况根本连定性正确反映分子静电势也做不到，比如卤键涉及的卤原子，原子电荷一般为负，但sigma-hole位置静电势为正。

作者
Author: 小书童 时间: 2016-8-24 20:01

sobereva 发表于 2016-8-24 17:36
没法这么说。只能说正值地方是核电荷的贡献超过电子的贡献，负值则是反过来。
而且原子电荷很多情况根本 ...

哦哦，明白了，谢谢老师指点！

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