计算化学公社

标题: doml3中过渡态频率搜索不收敛 [打印本页]

作者
Author: zjf1240291644 时间: 2023-2-27 09:43
标题: doml3中过渡态频率搜索不收敛
dmol3中过渡态结构优化顺利但是频率计算不收敛，在改变charge和spin后没太大影响，改变diis对收敛有帮助，但是diis只有我手动改input调到11才能收敛，请问大佬们有什么建议

作者
Author: 乐平 时间: 2023-2-27 10:08
你要优化什么过渡态？

吸附的 CO 和 CH2 之间形成 C-O 的过渡态？

如果是，显然从你截图的结构来看，并不是过渡态的样子，C 和 O 太远了。

另外，频率分析的时候，催化剂表面的原子都固定了吗？

作者
Author: zjf1240291644 时间: 2023-2-27 10:30

乐平发表于 2023-2-27 10:08
你要优化什么过渡态？

吸附的 CO 和 CH2 之间形成 C-O 的过渡态？

我是CH2OH的过渡态，我们这个计算应该是不固定表面原子的

作者
Author: 乐平 时间: 2023-2-27 10:46

zjf1240291644 发表于 2023-2-27 10:30
我是CH2OH的过渡态，我们这个计算应该是不固定表面原子的

不固定表面的原子？那你想想有多少振动模式吧 ……

3N - 6，其中 N 是原子个数。容易收敛吗？

作者
Author: zjf1240291644 时间: 2023-2-27 10:57

乐平发表于 2023-2-27 10:46
不固定表面的原子？那你想想有多少振动模式吧 ……

3N - 6，其中 N 是原子个数。容易收敛吗？

但是催化剂就是表面发生作用，所以我只固定了下面两层的原子

作者
Author: 卡开发发 时间: 2023-2-27 12:48
1、就改DIIS就行，是奏效的办法。
2、DMol3中计算频率固定原子应该不起什么作用，要用partial Hessian的方式来做，但此时没办法用来计算热化学性质。

作者
Author: zjf1240291644 时间: 2023-2-27 13:45

卡开发发发表于 2023-2-27 12:48
1、就改DIIS就行，是奏效的办法。
2、DMol3中计算频率固定原子应该不起什么作用，要用partial Hessian的方 ...

大佬，请问calculation参数diis设置里面最大到10，我这改到11去发文章别人会认吗

作者
Author: 卡开发发 时间: 2023-2-27 17:37

zjf1240291644 发表于 2023-2-27 13:45
大佬，请问calculation参数diis设置里面最大到10，我这改到11去发文章别人会认吗

没有问题，castep用的是20，VASP默认是40。

作者
Author: zjf1240291644 时间: 2023-2-27 21:10
谢谢大佬的解答，十分感谢

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