计算化学公社

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

[ORCA] 求助orca AIMD模拟有机小分子加热分解过程的文献或资料

[复制链接 Copy URL]

29

帖子

0

威望

2088

eV
积分
2117

Level 5 (御坂)

各位前辈,如题,1. 想要AIMD模拟有机小分子真空状态下的加热分解温度和分解位点,想知道orca是否能够模拟?
2. 如果orca做的不好或者不能模拟,那么AIMD模拟我应该用什么软件更为合适?
2. 网上并没有找到类似模拟的文献或资料,想请教各位前辈是否能够提供文献和资料,我可以自学。
先谢谢各位!


生命不息,折腾不止!

325

帖子

0

威望

2926

eV
积分
3251

Level 5 (御坂)

计算化学路人甲

2#
发表于 Post on 2021-6-23 14:37:08 | 只看该作者 Only view this author
关注,楼主弄明白了欢迎分享经验。

2407

帖子

1

威望

5960

eV
积分
8387

Level 6 (一方通行)

3#
发表于 Post on 2021-6-23 14:48:57 | 只看该作者 Only view this author
AIMD, ab initio molecular dynamics?
Check out this link, https://brehm-research.de/orcamd,
as suggested by the ORCA manual.
High-Performance Computing for You
为您专属定制的高性能计算解决方案

更多讯息,请访问:
https://labitc.top
http://tophpc.top:8080
电邮: ask@hpc4you.top

1万

帖子

0

威望

9004

eV
积分
20764

Level 6 (一方通行)

4#
发表于 Post on 2021-6-23 15:48:57 | 只看该作者 Only view this author
AIMD只能跑ps级的时间尺度,如果用半经验也只能跑到ns级。
如果你的分子在几ns之内就有明显分解,那可以用AIMD算。但是一般热分解没有这么快,至少得是毫秒级的,如果是那样就算不动了,而且不止orca,任何支持AIMD的软件都算不动。
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

29

帖子

0

威望

2088

eV
积分
2117

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2021-6-23 16:02:15 | 只看该作者 Only view this author
wzkchem5 发表于 2021-6-23 15:48
AIMD只能跑ps级的时间尺度,如果用半经验也只能跑到ns级。
如果你的分子在几ns之内就有明显分解,那可以用 ...

前辈,
1.AIMD跑单个分子(100个原子)的分解也很难实现吗?
2.既然AIMD没办法模拟,那是不是只能用反应力场(reaxff)动力学来模拟?
3.但是反应力场普适性太差了,还有别的更好更合适的模拟手段吗?
生命不息,折腾不止!

6万

帖子

99

威望

5万

eV
积分
120201

管理员

公社社长

6#
发表于 Post on 2021-6-23 16:39:51 | 只看该作者 Only view this author
1 可以
2 xtb用GFN-xTB可以跑明显长得多得多的尺度
即便是被认为不稳定的分子,在一般温度下有限的模拟时间内也很难跑出分解现象,往往都需要把温度设得很高使得有限时间内分解现象得以出现

北京科音量子化学高级培训班里的例子:








北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

1万

帖子

0

威望

9004

eV
积分
20764

Level 6 (一方通行)

7#
发表于 Post on 2021-6-23 19:13:52 | 只看该作者 Only view this author
胖鱼加加 发表于 2021-6-23 09:02
前辈,
1.AIMD跑单个分子(100个原子)的分解也很难实现吗?
2.既然AIMD没办法模拟,那是不是只能用反 ...

如果这个分子的分解时间在ps级,那可以。但是既然你研究的是分解温度,按我理解在分解温度下,分解时间至少也是毫秒级吧?如果某个温度下几ps就分解了,那就意味着在低得多的温度下就能观察到分解,那这个温度就不是分解温度,而是远高于分解温度。这个和sob老师举的例子还不太一样,有的研究只关心分解途径,这时温度可以设得比实际温度高很多,把分解时间缩短到ps量级。
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

224

帖子

5

威望

4578

eV
积分
4902

Level 6 (一方通行)

8#
发表于 Post on 2021-6-23 20:26:13 | 只看该作者 Only view this author
wzkchem5 发表于 2021-6-23 19:13
如果这个分子的分解时间在ps级,那可以。但是既然你研究的是分解温度,按我理解在分解温度下,分解时间至 ...

主要是高温下的分解途径不一定对应于低温下的最优路径吧,感觉这种做法似乎需要经过比较严格的考察才成
我需要一些假日,但我不希望每天都是假日。因为我没有承担痛苦,因为那不是真正的自由。

1万

帖子

0

威望

9004

eV
积分
20764

Level 6 (一方通行)

9#
发表于 Post on 2021-6-23 21:48:42 | 只看该作者 Only view this author
ene 发表于 2021-6-23 13:26
主要是高温下的分解途径不一定对应于低温下的最优路径吧,感觉这种做法似乎需要经过比较严格的考察才成

对,这也是个问题
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

29

帖子

0

威望

2088

eV
积分
2117

Level 5 (御坂)

10#
 楼主 Author| 发表于 Post on 2021-6-24 09:35:47 | 只看该作者 Only view this author
wzkchem5 发表于 2021-6-23 21:48
对,这也是个问题

看了几位前辈的回答和指导,我又迷茫了,
1.想要考察有机小分子分解温度和/或分解位点,如果不用AIMD,还有什么好的办法嘛?
2.如果想要尝试用AIMD,我是不是要尝试设置很高的温度,跑一下试试,看分解位点能不能和实验值一致?
生命不息,折腾不止!

1万

帖子

0

威望

9004

eV
积分
20764

Level 6 (一方通行)

11#
发表于 Post on 2021-6-24 15:20:12 | 只看该作者 Only view this author
胖鱼加加 发表于 2021-6-24 02:35
看了几位前辈的回答和指导,我又迷茫了,
1.想要考察有机小分子分解温度和/或分解位点,如果不用AIMD, ...

可以用高温下的AIMD找到几种可能的分解途径,再用过渡态理论计算这几种分解途径的速率
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

29

帖子

0

威望

2088

eV
积分
2117

Level 5 (御坂)

12#
 楼主 Author| 发表于 Post on 2021-7-12 16:50:54 | 只看该作者 Only view this author
wzkchem5 发表于 2021-6-24 15:20
可以用高温下的AIMD找到几种可能的分解途径,再用过渡态理论计算这几种分解途径的速率

前辈,用高温AIMD找可能的反应位点我会,
但是如何用过渡态理论计算热分解途径的反应速率呢?是要找热分解反应的过渡态吗?麻烦前辈指点。
生命不息,折腾不止!

1万

帖子

0

威望

9004

eV
积分
20764

Level 6 (一方通行)

13#
发表于 Post on 2021-7-12 16:56:08 | 只看该作者 Only view this author
胖鱼加加 发表于 2021-7-12 09:50
前辈,用高温AIMD找可能的反应位点我会,
但是如何用过渡态理论计算热分解途径的反应速率呢?是要找热分 ...

对,找热分解反应的过渡态,然后用Eyring方程换算
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

46

帖子

0

威望

448

eV
积分
494

Level 3 能力者

14#
发表于 Post on 2021-7-13 14:20:06 | 只看该作者 Only view this author
之前用reaxff反应动力学模拟过材料的热分解温度。
首先分子的数目必须要多,最好上百个。不然会有很大的偶然性,你看试验的TGA曲线就知道每个分子的分解温度不可能都一样。
立场只能挑选一些合适的然后一个个去尝试,看哪个和实测更加符合。

29

帖子

0

威望

2088

eV
积分
2117

Level 5 (御坂)

15#
 楼主 Author| 发表于 Post on 2021-7-22 10:02:45 | 只看该作者 Only view this author
wzkchem5 发表于 2021-7-12 16:56
对,找热分解反应的过渡态,然后用Eyring方程换算

前辈,我还有些小问题搞不清楚:
(1)热分解反应过渡态是不是就是有机小分子均裂成两个自由基过程的过渡态?
(2)那我需要通过柔性扫描或者刚性扫描,扫描键长获得解离曲线,对不对?
(3)看sob老师的帖子里 http://sobereva.com/474 ,HF的解离曲线能量也没有最高点;这样的情况是不是就没有过渡态呢?
(4)均裂反应没有过渡态的情况该怎么办?
生命不息,折腾不止!

本版积分规则 Credits rule

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

GMT+8, 2025-8-18 05:25 , Processed in 0.302485 second(s), 23 queries , Gzip On.

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