计算化学公社

 找回密码 Forget password
 注册 Register
Views: 1942|回复 Reply: 8
打印 Print 上一主题 Last thread 下一主题 Next thread

[CASTEP/Dmol3/MS] doml3中过渡态频率搜索不收敛

[复制链接 Copy URL]

11

帖子

0

威望

67

eV
积分
78

Level 2 能力者

dmol3中过渡态结构优化顺利但是频率计算不收敛,在改变charge和spin后没太大影响,改变diis对收敛有帮助,但是diis只有我手动改input调到11才能收敛,请问大佬们有什么建议

(C9IW3EMK43ZYD4F7RSEJJK.png (128.38 KB, 下载次数 Times of downloads: 10)

(C9IW3EMK43ZYD4F7RSEJJK.png

[ZMNR75BDI7O($RX~)2KEIS.png (337.03 KB, 下载次数 Times of downloads: 6)

[ZMNR75BDI7O($RX~)2KEIS.png

1155

帖子

0

威望

4429

eV
积分
5584

Level 6 (一方通行)

2#
发表于 Post on 2023-2-27 10:08:06 | 只看该作者 Only view this author
你要优化什么过渡态?

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

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

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

11

帖子

0

威望

67

eV
积分
78

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2023-2-27 10:30:53 | 只看该作者 Only view this author
乐平 发表于 2023-2-27 10:08
你要优化什么过渡态?

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

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

1155

帖子

0

威望

4429

eV
积分
5584

Level 6 (一方通行)

4#
发表于 Post on 2023-2-27 10:46:19 | 只看该作者 Only view this author
zjf1240291644 发表于 2023-2-27 10:30
我是CH2OH的过渡态,我们这个计算应该是不固定表面原子的

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

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

11

帖子

0

威望

67

eV
积分
78

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2023-2-27 10:57:27 | 只看该作者 Only view this author
乐平 发表于 2023-2-27 10:46
不固定表面的原子? 那你想想有多少振动模式吧 ……

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

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

3809

帖子

3

威望

1万

eV
积分
20343

Level 6 (一方通行)

围观吃瓜群众

6#
发表于 Post on 2023-2-27 12:48:29 | 只看该作者 Only view this author
1、就改DIIS就行,是奏效的办法。
2、DMol3中计算频率固定原子应该不起什么作用,要用partial Hessian的方式来做,但此时没办法用来计算热化学性质。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

11

帖子

0

威望

67

eV
积分
78

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2023-2-27 13:45:23 | 只看该作者 Only view this author
卡开发发 发表于 2023-2-27 12:48
1、就改DIIS就行,是奏效的办法。
2、DMol3中计算频率固定原子应该不起什么作用,要用partial Hessian的方 ...

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

3809

帖子

3

威望

1万

eV
积分
20343

Level 6 (一方通行)

围观吃瓜群众

8#
发表于 Post on 2023-2-27 17:37:19 | 只看该作者 Only view this author
zjf1240291644 发表于 2023-2-27 13:45
大佬,请问calculation参数diis设置里面最大到10,我这改到11去发文章别人会认吗

没有问题,castep用的是20,VASP默认是40。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

11

帖子

0

威望

67

eV
积分
78

Level 2 能力者

9#
 楼主 Author| 发表于 Post on 2023-2-27 21:10:40 | 只看该作者 Only view this author
谢谢大佬的解答,十分感谢

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-21 07:27 , Processed in 0.174835 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list