| 查看: 1192 | 回复: 17 | |||
| 当前主题已经存档。 | |||
sheath至尊木虫 (知名作家)
|
[交流]
【讨论】自编程序中正则系综的体系温度为什么不能保持恒定？
|
||
|
在自编的分子动力学程序中，计算原胞的原子体系为正则系综。原子位置在三维方向应用周期性边界条件，恒温采用Nose-Hoover算法。 但是，程序执行后，出现了体系数目减少，体系温度不断升高，甚至到达上万度的错误结果。 经过仔细检查程序结构和每个子程序，初步分析认为：在每个时间步，周期性边界条件使跑出计算原胞范围的原子重新回到原胞内的相应位置。而这可能使部分原子重叠，尽管这种情况的几率很小。随后进行“更新原子存在与对应数据”操作，这将可能造成部分重叠原子数据的丢失，当然，原子数目也可能减小。 也就是说，周期性边界条件造成原子重叠，而原子重叠是导致原子数目减少的关键原因。如果这个分析对的话，如何解决？若不是，那么原因又在哪里？ 敬请各位同行不吝赐教！！！ [ Last edited by cadick on 2009-12-14 at 02:30 ] |
回复此楼» 猜你喜欢
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
UJN物理学专业调剂
已经有12人回复
-
物理学I论文润色/翻译怎么收费?
已经有274人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
-
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组招收2026硕士研究生1名
已经有0人回复
-
山西大同大学物理学硕士研究生招收调剂生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
［调剂信息］211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
一志愿物理所凝聚态物理352求调剂
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
zeoliters
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 贵宾: 0.902
- 金币: 2088.4
- 散金: 16
- 帖子: 445
- 在线: 60.6小时
- 虫号: 772181
- 注册: 2009-05-16
- 专业: 计算化学，量子化学
2楼2009-11-08 18:17:44
sheath
至尊木虫 (知名作家)
- 应助: 0 (幼儿园)
- 金币: 26023.5
- 红花: 2
- 帖子: 5954
- 在线: 416.3小时
- 虫号: 445940
- 注册: 2007-10-28
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
3楼2009-11-08 18:42:41
sheath
至尊木虫 (知名作家)
- 应助: 0 (幼儿园)
- 金币: 26023.5
- 红花: 2
- 帖子: 5954
- 在线: 416.3小时
- 虫号: 445940
- 注册: 2007-10-28
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
4楼2009-11-09 10:14:31
dellus
木虫 (职业作家)
- 应助: 0 (幼儿园)
- 贵宾: 0.08
- 金币: 5305.6
- 散金: 35
- 红花: 1
- 沙发: 2
- 帖子: 3469
- 在线: 285.9小时
- 虫号: 92581
- 注册: 2005-09-10
- 性别: GG
5楼2009-11-09 11:20:51
6楼2009-11-09 15:56:58
sheath
至尊木虫 (知名作家)
- 应助: 0 (幼儿园)
- 金币: 26023.5
- 红花: 2
- 帖子: 5954
- 在线: 416.3小时
- 虫号: 445940
- 注册: 2007-10-28
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
★ ★ ★ ★
dongdong3881(金币+2,VIP+0):多谢交流 11-10 14:09
zeoliters(金币+2,VIP+0):谢谢回帖交流！ 11-10 22:08
dongdong3881(金币+2,VIP+0):多谢交流 11-10 14:09
zeoliters(金币+2,VIP+0):谢谢回帖交流！ 11-10 22:08
|
这个我以前的确没有考虑过， 不过，试过后还是出现原来的问题。 我的程序在计算原子位置和相互作用力都使用了周期性边界条件。 依据是被有些同行称为“分子动力学圣经”中的一段话： "The wraparound effect of the periodic boundaries must be taken into account in both the integration of the equations of motion and the interaction computations. After each integration step the coordinates must be examined, and if an atom is found to have moved outside the region its coordinates must be adjusted to bring it back inside." 来源：D.C. Rapaport, The Art of Molecular Dynamics Simulation (2ed), Cambridge University Press,2004, P16-17 [ Last edited by sheath on 2009-11-9 at 17:28 ] |
7楼2009-11-09 17:20:47
sheath
至尊木虫 (知名作家)
- 应助: 0 (幼儿园)
- 金币: 26023.5
- 红花: 2
- 帖子: 5954
- 在线: 416.3小时
- 虫号: 445940
- 注册: 2007-10-28
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
8楼2009-11-09 17:39:44
dellus
木虫 (职业作家)
- 应助: 0 (幼儿园)
- 贵宾: 0.08
- 金币: 5305.6
- 散金: 35
- 红花: 1
- 沙发: 2
- 帖子: 3469
- 在线: 285.9小时
- 虫号: 92581
- 注册: 2005-09-10
- 性别: GG
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
dongdong3881(金币+2,VIP+0):多谢回帖交流 11-10 14:09
zeoliters(金币+3,VIP+0):多谢回帖交流！ 11-10 22:09
小木虫(金币+0.5):给个红包，谢谢回帖交流
dongdong3881(金币+2,VIP+0):多谢回帖交流 11-10 14:09
zeoliters(金币+3,VIP+0):多谢回帖交流！ 11-10 22:09
|
计算能量和相互作用力都要用周期性边界条件 我是说在使用周期性边界条件时，只调用而不改变粒子的原坐标 第一句没有问题，第二句值得商榷 １、在周期性边界条件下，同在原胞下的两个粒子是不能通过直接计算它们之间距离来计算能量和受力作用的。把原胞外的粒子放回原胞的操作虚耗了计算量，还可能出错 ２、对于高分子体系，这样的操作会破坏共价键 Allen 的 ＂Computer simulation of liquid＂ ２９－３１页讲了关于周期性边界条件的问题，建议看下 [ Last edited by dellus on 2009-11-10 at 13:20 ] |
9楼2009-11-10 13:13:13
zeoliters
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 贵宾: 0.902
- 金币: 2088.4
- 散金: 16
- 帖子: 445
- 在线: 60.6小时
- 虫号: 772181
- 注册: 2009-05-16
- 专业: 计算化学，量子化学
10楼2009-11-11 05:11:09
