| 查看: 4568 | 回复: 15 | |||
[交流]
Gaussian计算的随机性 已有8人参与
|
|
同一个版本Gaussian软件，在三台不同的服务器上进行优化计算，同一个输入文件，前27步的能量完全一样，从28步后能量就开始不一样了，难道Gaussian计算也有随机性，请教各位高手是否也遇到类似的问题。 20 -2679.23298795 -2679.23298795 -2679.23298795 21 -2679.23319064 -2679.23319065 -2679.23319065 22 -2679.23339199 -2679.23339199 -2679.23339199 23 -2679.23355888 -2679.23355888 -2679.23355889 24 -2679.23367177 -2679.23367177 -2679.23367177 25 -2679.23374824 -2679.23374824 -2679.23374825 26 -2679.23381226 -2679.23381226 -2679.23381226 27 -2679.23384277 -2679.23384277 -2679.23384277 28 -2679.23388423 -2679.23388425 -2679.23388425 29 -2679.23390839 -2679.23390840 -2679.23390840 30 -2679.23399369 -2679.23399347 -2679.23399351 31 -2679.23397296 -2679.23417440 -2679.23413908 32 -2679.23400632 -2679.23419135 -2679.23414835 33 -2679.23413415 -2679.23422701 -2679.23419652 |
回复此楼» 猜你喜欢
中南大学化学化工学院招收电池储能、理论计算方向博士生
已经有16人回复
-
请问哪里可以代烧带压力的炉子？
已经有6人回复
-
物理化学论文润色/翻译怎么收费?
已经有259人回复
-
关于合成MIL-125的求助
已经有0人回复
-
双非本硕光催化方向想读博，请各位老师看看我！
已经有2人回复
-
瑞典林雪平大学范立洲课题组现招聘博士后、岗位制博士生
已经有0人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
寻求物理化学学报最新的投稿模板！
已经有4人回复
-
投稿咨询
已经有7人回复
-
自科面上 散金币 祈福
已经有61人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- gaussian计算
已经有4人回复
- 该物质是对称结构，计算结果却是HOMO在左，LUMO在右。
已经有9人回复
- 使用Multiwfn做拓扑分析以及计算孤对电子角度
已经有18人回复
- 【sobereva个人文集】制作动画分析电子结构特征
已经有29人回复
- 【求助】50金币求助，金属间界面和金属氧化物间界面模型怎么建立？
已经有30人回复
lihb734
铁杆木虫 (职业作家)
站在计算化学入门的门槛上
- QC强帖: 2
- 应助: 132 (高中生)
- 金币: 7495.4
- 散金: 9378
- 红花: 74
- 帖子: 3107
- 在线: 1177.3小时
- 虫号: 448403
- 注册: 2007-11-01
- 性别: GG
- 专业: 理论和计算化学
2楼2011-05-05 17:47:35
hairan
木虫 (著名写手)
瞌睡木虫
- QC强帖: 3
- 应助: 61 (初中生)
- 贵宾: 1.06
- 金币: 4974.3
- 散金: 8072
- 红花: 44
- 帖子: 2901
- 在线: 733.4小时
- 虫号: 432069
- 注册: 2007-08-12
- 性别: GG
- 专业: 理论和计算化学
3楼2011-05-05 18:09:50
lihb734
铁杆木虫 (职业作家)
站在计算化学入门的门槛上
- QC强帖: 2
- 应助: 132 (高中生)
- 金币: 7495.4
- 散金: 9378
- 红花: 74
- 帖子: 3107
- 在线: 1177.3小时
- 虫号: 448403
- 注册: 2007-11-01
- 性别: GG
- 专业: 理论和计算化学
4楼2011-05-05 18:13:48
hairan
木虫 (著名写手)
瞌睡木虫
- QC强帖: 3
- 应助: 61 (初中生)
- 贵宾: 1.06
- 金币: 4974.3
- 散金: 8072
- 红花: 44
- 帖子: 2901
- 在线: 733.4小时
- 虫号: 432069
- 注册: 2007-08-12
- 性别: GG
- 专业: 理论和计算化学
★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
gmy1990(金币+1): 非常感谢你对论坛的贡献，奖励一下！ 2011-05-05 18:49:10
375642546(金币+3): 给力！ 2011-05-06 09:00:36
小木虫(金币+0.5):给个红包，谢谢回帖
gmy1990(金币+1): 非常感谢你对论坛的贡献，奖励一下！ 2011-05-05 18:49:10
375642546(金币+3): 给力！ 2011-05-06 09:00:36
|
我的感觉是类似蝴蝶效应，就是在某一步计算出现了一个看似微不足道的差别，但这种差别不断累积下来最后导致了能量的不尽一致。 不过楼主如果能继续优化完，可能最终的优化结构和总能量也不会有太大差别。 我加的那个关键词是让计算积分时的误差尽量小，看看能不能通过减小一点这方面的误差来保持能量的相同。 另外，总能量的差别还不到0.0001 a.u.，不算很大的差别吧？ [ Last edited by hairan on 2011-5-5 at 18:26 ] |
5楼2011-05-05 18:25:23
coolrainbow
木虫 (著名写手)
未来国家冻凉
- QC强帖: 11
- 应助: 1 (幼儿园)
- 贵宾: 0.1
- 金币: 3979.5
- 散金: 10
- 红花: 50
- 帖子: 1137
- 在线: 272.3小时
- 虫号: 239560
- 注册: 2006-04-07
- 性别: GG
- 专业: 理论和计算化学
6楼2011-05-05 19:44:29
7楼2011-05-05 20:12:44
sobereva
至尊木虫 (著名写手)
本人已永久离开小木虫
- QC强帖: 17
- 应助: 226 (大学生)
- 贵宾: 5.4
- 金币: 2180.5
- 红花: 217
- 帖子: 1241
- 在线: 366.8小时
- 虫号: 281271
- 注册: 2006-09-23
- 专业: 宇宙学
★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
375642546(金币+3): 2011-05-06 16:57:38
375642546(金币+3): 2011-05-06 16:58:04
小木虫(金币+0.5):给个红包，谢谢回帖
375642546(金币+3): 2011-05-06 16:57:38
375642546(金币+3): 2011-05-06 16:58:04
|
可能的原因非常多。例如根据内存大小的不同自动使用不同积分储存方式（在写入/读取时引入数值差异），并行计算时数据累积运算顺序存在的一定随机性等等。实际上前27步也并非相同，只是在有限的数值输出精度中没显示出来罢了，随着累积运算的进行差异被逐渐放大到可感知的地步。 有兴趣的话建议你看看J Comput Aided Mol Des (2008) 22:39–51。虽然是讨论数值误差对分子力学优化影响的，但对理解这个问题也有益。 |
8楼2011-05-05 23:20:31
9楼2011-05-06 12:21:35
ghcacj
荣誉版主 (著名写手)
阿超
- QC强帖: 2
- 应助: 26 (小学生)
- 贵宾: 2.764
- 金币: 2543.5
- 散金: 8899
- 红花: 35
- 沙发: 3
- 帖子: 1376
- 在线: 1199小时
- 虫号: 513281
- 注册: 2008-02-27
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
10楼2011-05-06 13:32:21
40