计算化学公社

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

[新手求助] 求助两个离子的反应的irc和过渡态都走通了,但是在模拟ads的时候两个离子越来越远

[复制链接 Copy URL]

99

帖子

0

威望

239

eV
积分
338

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 han- 于 2022-11-23 10:05 编辑

各位老师 ,我在对两个离子的反应进行模拟的时候,irc 和过渡态都走通了,但是在模拟ads的时候 两个离子越来越远一直不收敛 报错该怎么设置呢?

com1.log

17.12 KB, 下载次数 Times of downloads: 4

com1.gjf

1.39 KB, 下载次数 Times of downloads: 8

1016

帖子

0

威望

2096

eV
积分
3112

Level 5 (御坂)

天上地下三山亖海五湖六合八荒九州之泛粤大典編輯委員會編外委員

14#
发表于 Post on 2022-11-24 15:05:40 | 只看该作者 Only view this author
han- 发表于 2022-11-24 14:55
老师 这样的话 我之前算出来的过渡态和irc的能量是不是也是不对的?

肯定有偏差,搞不好定性也不对。
ORCA大法好!CP2K大法好!
嶺南粤音 泛粤典 https://jyutjam.org/

山河自落蒼梧月,風雨猶驅草木兵。
san huɔɞ ki lɔk tʰɔŋ ŋ ᵑgut,fʊŋ ʝi ʝiu kʰui tʰɔu ᵐbʊk pɪŋ

99

帖子

0

威望

239

eV
积分
338

Level 3 能力者

13#
 楼主 Author| 发表于 Post on 2022-11-24 14:55:33 | 只看该作者 Only view this author
chands 发表于 2022-11-24 12:06
不行。因为这个体系里没有单独的离子,每个离子与介质(离子液体)的作用太强,不能不考虑。

老师 这样的话 我之前算出来的过渡态和irc的能量是不是也是不对的?

1016

帖子

0

威望

2096

eV
积分
3112

Level 5 (御坂)

天上地下三山亖海五湖六合八荒九州之泛粤大典編輯委員會編外委員

12#
发表于 Post on 2022-11-24 12:06:32 | 只看该作者 Only view this author
han- 发表于 2022-11-24 11:24
老师 如果直接用过渡态的能量减去单独优化好的两个离子的能量,是不是不行?

不行。因为这个体系里没有单独的离子,每个离子与介质(离子液体)的作用太强,不能不考虑。
ORCA大法好!CP2K大法好!
嶺南粤音 泛粤典 https://jyutjam.org/

山河自落蒼梧月,風雨猶驅草木兵。
san huɔɞ ki lɔk tʰɔŋ ŋ ᵑgut,fʊŋ ʝi ʝiu kʰui tʰɔu ᵐbʊk pɪŋ

99

帖子

0

威望

239

eV
积分
338

Level 3 能力者

11#
 楼主 Author| 发表于 Post on 2022-11-24 11:24:24 | 只看该作者 Only view this author
chands 发表于 2022-11-23 10:19
我昨天跑了你的输出文件,也是convergence failure,可能跟抗衡离子和背景电荷设置不对有关。

离子液 ...

老师 如果直接用过渡态的能量减去单独优化好的两个离子的能量,是不是不行?

1016

帖子

0

威望

2096

eV
积分
3112

Level 5 (御坂)

天上地下三山亖海五湖六合八荒九州之泛粤大典編輯委員會編外委員

10#
发表于 Post on 2022-11-23 10:19:23 | 只看该作者 Only view this author
han- 发表于 2022-11-23 10:07
老师 输出已上传 但是报错了
老师,这个阴离子模拟一定要从团簇里抠吗?
可以只加背景电荷不加抗衡离子 ...

我昨天跑了你的输出文件,也是convergence failure,可能跟抗衡离子和背景电荷设置不对有关。

离子液体的反应是比较复杂。最严格的方法是跑MD然后观测反应历程,用背景电荷是简化。但是背景电荷的设置估计也得跑个MD来决定,“从团簇里抠”的意思是只考虑相邻的区域,这样体系不用太大,计算耗时太长。

抗衡离子的意义是显式表示阴离子和周围阳离子的作用。这个作用用背景电荷描述太粗糙,结果搞不好连定性都不准。
ORCA大法好!CP2K大法好!
嶺南粤音 泛粤典 https://jyutjam.org/

山河自落蒼梧月,風雨猶驅草木兵。
san huɔɞ ki lɔk tʰɔŋ ŋ ᵑgut,fʊŋ ʝi ʝiu kʰui tʰɔu ᵐbʊk pɪŋ

99

帖子

0

威望

239

eV
积分
338

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2022-11-23 10:07:07 | 只看该作者 Only view this author
chands 发表于 2022-11-22 22:38
你把输出文件传上来看看。

老师 输出已上传 但是报错了
老师,这个阴离子模拟一定要从团簇里抠吗?
可以只加背景电荷不加抗衡离子吗?

1016

帖子

0

威望

2096

eV
积分
3112

Level 5 (御坂)

天上地下三山亖海五湖六合八荒九州之泛粤大典編輯委員會編外委員

8#
发表于 Post on 2022-11-22 22:38:14 | 只看该作者 Only view this author
han- 发表于 2022-11-22 22:19
老师 麻烦您看一下我这个输入问题出在哪里了?我看相关背景电荷和抗衡离子的帖子都是关于晶体或团簇的;
...

你把输出文件传上来看看。
ORCA大法好!CP2K大法好!
嶺南粤音 泛粤典 https://jyutjam.org/

山河自落蒼梧月,風雨猶驅草木兵。
san huɔɞ ki lɔk tʰɔŋ ŋ ᵑgut,fʊŋ ʝi ʝiu kʰui tʰɔu ᵐbʊk pɪŋ

99

帖子

0

威望

239

eV
积分
338

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2022-11-22 22:19:52 | 只看该作者 Only view this author
chands 发表于 2022-11-17 11:17
加了抗衡离子之后,电荷数要按实际情况决定,比如原来体系电荷是-3,加了三个K+,当然电荷数现在是0。
...

老师 麻烦您看一下我这个输入问题出在哪里了?我看相关背景电荷和抗衡离子的帖子都是关于晶体或团簇的;
还有一个问题老师,背景电荷要保证和我的体系电性相同吗?(比如这是两个阴离子反应,所以背景电荷需要是-2吗?)

1016

帖子

0

威望

2096

eV
积分
3112

Level 5 (御坂)

天上地下三山亖海五湖六合八荒九州之泛粤大典編輯委員會編外委員

6#
发表于 Post on 2022-11-17 11:17:56 | 只看该作者 Only view this author
han- 发表于 2022-11-16 21:02
老师 我在周围添加了一些钾离子 选择了SMD 这样对吗?
还有就是老师 我加了抗衡离子之后 我需要把电荷 ...

加了抗衡离子之后,电荷数要按实际情况决定,比如原来体系电荷是-3,加了三个K+,当然电荷数现在是0。
用隐式溶剂模型要根据体系指定介电常数等参数,见
http://sobereva.com/431
ORCA大法好!CP2K大法好!
嶺南粤音 泛粤典 https://jyutjam.org/

山河自落蒼梧月,風雨猶驅草木兵。
san huɔɞ ki lɔk tʰɔŋ ŋ ᵑgut,fʊŋ ʝi ʝiu kʰui tʰɔu ᵐbʊk pɪŋ

99

帖子

0

威望

239

eV
积分
338

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2022-11-16 21:02:16 | 只看该作者 Only view this author
chands 发表于 2022-11-16 18:27
啊,这个要复杂点。实际体系中,由于抗衡离子和背景电荷的作用,两个阴离子之间的排斥没这么大。所以体系 ...

老师 我在周围添加了一些钾离子 选择了SMD 这样对吗?
还有就是老师 我加了抗衡离子之后 我需要把电荷由之前的-3 改为0 吗?

1016

帖子

0

威望

2096

eV
积分
3112

Level 5 (御坂)

天上地下三山亖海五湖六合八荒九州之泛粤大典編輯委員會編外委員

4#
发表于 Post on 2022-11-16 18:27:01 | 只看该作者 Only view this author
本帖最后由 chands 于 2022-11-16 18:31 编辑
han- 发表于 2022-11-16 17:36
老师 是固体反应处于熔融态

啊,这个要复杂点。实际体系中,由于抗衡离子和背景电荷的作用,两个阴离子之间的排斥没这么大。所以体系中要要加抗衡离子,还要体现背景电荷(比如用隐式溶剂模型)。
ORCA大法好!CP2K大法好!
嶺南粤音 泛粤典 https://jyutjam.org/

山河自落蒼梧月,風雨猶驅草木兵。
san huɔɞ ki lɔk tʰɔŋ ŋ ᵑgut,fʊŋ ʝi ʝiu kʰui tʰɔu ᵐbʊk pɪŋ

99

帖子

0

威望

239

eV
积分
338

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2022-11-16 17:36:33 | 只看该作者 Only view this author
chands 发表于 2022-11-16 16:21
你的体系是气相?溶液?气相中两个阴离子会相互排斥。

老师 是固体反应处于熔融态

1016

帖子

0

威望

2096

eV
积分
3112

Level 5 (御坂)

天上地下三山亖海五湖六合八荒九州之泛粤大典編輯委員會編外委員

2#
发表于 Post on 2022-11-16 16:21:28 | 只看该作者 Only view this author
你的体系是气相?溶液?气相中两个阴离子会相互排斥。
ORCA大法好!CP2K大法好!
嶺南粤音 泛粤典 https://jyutjam.org/

山河自落蒼梧月,風雨猶驅草木兵。
san huɔɞ ki lɔk tʰɔŋ ŋ ᵑgut,fʊŋ ʝi ʝiu kʰui tʰɔu ᵐbʊk pɪŋ

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

GMT+8, 2025-8-18 04:22 , Processed in 0.179747 second(s), 24 queries , Gzip On.

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