计算化学公社

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

[Gaussian/gview] 关于ADMP计算中Random Number Generator Seed对结果的影响

[复制链接 Copy URL]

20

帖子

0

威望

246

eV
积分
266

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 高章磊 于 2023-10-2 10:46 编辑

最近一段时间一直在用ADMP算一个分子的震动,分子模型如图所示:(氢原子这里省略了,没有显示)
主要是研究上下两个五元环在进行ADMP计算中的运动情况。
之前只关注了NKE对于计算结果的影响,前段时间仔细看了一下out文件,发现Random Number Generator Seed也会对计算结果产生很大的影响,之前的计算这个参数都是默认的398465,通过IOp(1/44=x)的方式可以在输入文件中修改这个参数为x。这个参数我查了一下能设置的范围很广,很多回答中也都说这个参数是完全随机产生新的初始能量分配,我想问一下真的是完全随机的吗?因为在其他条件相同的情况下,不同的Random Number Generator Seed会产生完全不同的计算结果。而我们和实验组合作,始终没办法得到合适的计算数据,所以很困扰。请问有人指点一下吗?感谢万分。

20

帖子

0

威望

246

eV
积分
266

Level 3 能力者

2#
 楼主 Author| 发表于 Post on 2023-9-27 12:32:45 | 只看该作者 Only view this author
本帖最后由 高章磊 于 2023-10-2 10:47 编辑

前面图片好像上传有问题,重新上传一下。

1万

帖子

0

威望

7402

eV
积分
18167

Level 6 (一方通行)

3#
发表于 Post on 2023-9-27 17:47:31 | 只看该作者 Only view this author
仔细看http://sobereva.com/88
正确的做法是随便取很多个不同的Random Number Generator Seed,然后统计有多少条轨迹得到想要的结果,多少条没有得到。一般而言不需要所有轨迹都跑出想要的结果,对于某些研究目的(例如研究实验上唯一一个能够观测到的反应途径,但实验灵敏度不高,不排除有少许副反应),只要绝大部分轨迹的结果满足预期就行了;有的情况下(例如研究主反应途径)只要过半满足预期就行了;有的情况下(例如探究一个反应是否能发生)只要有一条满足预期就行了。只取一个Random Number Generator Seed做计算在绝大多数情况下是不行的,哪怕这个Random Number Generator Seed是默认的、你之前用过的或前人用过的,也不行
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

5万

帖子

99

威望

5万

eV
积分
112469

管理员

公社社长

4#
发表于 Post on 2023-9-28 05:55:09 | 只看该作者 Only view this author
高章磊 发表于 2023-9-27 12:32
前面图片好像上传有问题,重新上传一下。

有别人回复之前若需要对帖子进行修改、补充,应直接编辑原帖,不要通过回帖进行补充,这点在置顶的新社员必读贴里明确说了。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

20

帖子

0

威望

246

eV
积分
266

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2023-9-28 12:57:53 | 只看该作者 Only view this author
sobereva 发表于 2023-9-28 05:55
有别人回复之前若需要对帖子进行修改、补充,应直接编辑原帖,不要通过回帖进行补充,这点在置顶的新社员 ...

了解了,不好意思啊,没注意,刚改过原帖了。

本版积分规则 Credits rule

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

GMT+8, 2024-11-26 14:53 , Processed in 1.274431 second(s), 20 queries , Gzip On.

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