计算化学公社

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

[ORCA] 做AIMD底物分子的初始相对位置和结构选择的问题

[复制链接 Copy URL]

146

帖子

0

威望

435

eV
积分
581

Level 4 (黑子)

各位老师好,

      我在看社长的文章《使用ORCA做从头算动力学(AIMD)的简单例子》, 在初始的氨气分子中的氮原子和水分子上的氢相对位置的选择上是不是需要靠得相对近一些?如果太远的话,我感觉这整个质子转移的发生过程可能需要很久的模拟时间才能实现,甚至于可能都观察不到质子的转移(距离太远的话)。而且如果距离太远的话,像这种孤立体系的话,两个分子之间是不是会一直远离?

1.png (57.82 KB, 下载次数 Times of downloads: 4)

1.png

1万

帖子

0

威望

7400

eV
积分
18164

Level 6 (一方通行)

2#
发表于 Post on 2022-9-11 15:12:38 | 只看该作者 Only view this author
你是实际跑过了发现质子确实转移不了,还是只是理论推测如此?
这个体系[Al(OH2)6]3+带电荷很多,而NH3又有比较大的偶极矩,所以两者之间有很强的吸引作用,并不像两个普通分子放在一起那样容易远离
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

977

帖子

0

威望

1929

eV
积分
2906

Level 5 (御坂)

偽語言學家+註冊摸魚師

3#
发表于 Post on 2022-9-11 15:20:00 | 只看该作者 Only view this author
不会一直远离,可能要久一点。还有如果初速度较大可能远离一阵,你的热浴温度有没有升得太快?
ORCA大法好!算得快,还免费,我爷爷的爷爷都说不赖!
山河自落蒼梧月,風雨猶驅草木兵。
san huɔᵊ ki lɔk tʰɔŋ ŋ ᵑgut,fʊŋ ʝi ʝiu kʰui tʰɔu ᵐbʊk pɪŋ。

146

帖子

0

威望

435

eV
积分
581

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2022-9-11 15:50:25 | 只看该作者 Only view this author
chands 发表于 2022-9-11 15:20
不会一直远离,可能要久一点。还有如果初速度较大可能远离一阵,你的热浴温度有没有升得太快?

我只是猜想一下,感觉可能会远离

146

帖子

0

威望

435

eV
积分
581

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2022-9-11 15:55:09 | 只看该作者 Only view this author
wzkchem5 发表于 2022-9-11 15:12
你是实际跑过了发现质子确实转移不了,还是只是理论推测如此?
这个体系[Al(OH2)6]3+带电荷很多,而NH3又 ...

我在做自己的体系(也是一个质子转移的过程)的时候,就是发现如果两个分子初始位置相距远一点,AIMD的模拟基本就观察不到两个相互靠近发生质子转移,就只有相互远离。所以在考虑这种模拟的时候初始位置的选择是不是也十分重要。

或者在我的体系中是不是通过增加一个正交的限制势,然后增加模拟的时间,也能观察质子转移的过程?

1万

帖子

0

威望

7400

eV
积分
18164

Level 6 (一方通行)

6#
发表于 Post on 2022-9-11 16:56:54 | 只看该作者 Only view this author
wwp 发表于 2022-9-11 08:55
我在做自己的体系(也是一个质子转移的过程)的时候,就是发现如果两个分子初始位置相距远一点,AIMD的模 ...

是不是因为你的体系里面质子给体不带电荷,或者带电荷不多?
一般情况下建议加限制势,不过如果要定量计算反应速率的话,得注意限制势对反应速率的影响
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

146

帖子

0

威望

435

eV
积分
581

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2022-9-11 18:52:45 | 只看该作者 Only view this author
本帖最后由 wwp 于 2022-9-11 18:58 编辑
wzkchem5 发表于 2022-9-11 16:56
是不是因为你的体系里面质子给体不带电荷,或者带电荷不多?
一般情况下建议加限制势,不过如果要定量计 ...

对的,质子给体是中性的。 我不用计算反应速率,应该不用考虑这个。谢谢老师

188

帖子

2

威望

3206

eV
积分
3434

Level 5 (御坂)

8#
发表于 Post on 2023-10-22 13:55:19 | 只看该作者 Only view this author
wzkchem5 发表于 2022-9-11 16:56
是不是因为你的体系里面质子给体不带电荷,或者带电荷不多?
一般情况下建议加限制势,不过如果要定量计 ...

王老师,请问这个限制势ORCA里好加吗?我的两个中性的也老是跑开。

1万

帖子

0

威望

7400

eV
积分
18164

Level 6 (一方通行)

9#
发表于 Post on 2023-10-22 16:22:34 | 只看该作者 Only view this author
量化小菜鸡 发表于 2023-10-22 06:55
王老师,请问这个限制势ORCA里好加吗?我的两个中性的也老是跑开。

在%md里面加cell关键字,具体看说明书
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

188

帖子

2

威望

3206

eV
积分
3434

Level 5 (御坂)

10#
发表于 Post on 2023-10-22 16:45:53 | 只看该作者 Only view this author
wzkchem5 发表于 2023-10-22 16:22
在%md里面加cell关键字,具体看说明书

好的,谢谢王老师。

本版积分规则 Credits rule

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

GMT+8, 2024-11-26 07:28 , Processed in 0.168531 second(s), 23 queries , Gzip On.

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