| 查看: 1147 | 回复: 17 | |||
| 当前主题已经存档。 | |||
sheath至尊木虫 (知名作家)
|
[交流]
【讨论】自编程序中正则系综的体系温度为什么不能保持恒定？
|
||
|
在自编的分子动力学程序中，计算原胞的原子体系为正则系综。原子位置在三维方向应用周期性边界条件，恒温采用Nose-Hoover算法。 但是，程序执行后，出现了体系数目减少，体系温度不断升高，甚至到达上万度的错误结果。 经过仔细检查程序结构和每个子程序，初步分析认为：在每个时间步，周期性边界条件使跑出计算原胞范围的原子重新回到原胞内的相应位置。而这可能使部分原子重叠，尽管这种情况的几率很小。随后进行“更新原子存在与对应数据”操作，这将可能造成部分重叠原子数据的丢失，当然，原子数目也可能减小。 也就是说，周期性边界条件造成原子重叠，而原子重叠是导致原子数目减少的关键原因。如果这个分析对的话，如何解决？若不是，那么原因又在哪里？ 敬请各位同行不吝赐教！！！ [ Last edited by cadick on 2009-12-14 at 02:30 ] |
回复此楼» 猜你喜欢
投稿chemical physical letters不送审？
已经有2人回复
-
南开大学物理学院张书辉副教授招收凝聚态物理理论方向博士生、硕士生
已经有22人回复
-
物理学I论文润色/翻译怎么收费?
已经有177人回复
-
杂志要求提供official language-editing certificate
已经有4人回复
-
深圳大学2026年秋博士招生-物理学-活性胶体方向-高永祥课题组
已经有16人回复
-
论物质与能量的统一模型及物理现象解释
已经有0人回复
-
基于基元I统一理论的数学相关应用推导
已经有0人回复
-
基元I统一理论：宇宙本质、层级演化与修炼文明的本源规律
已经有1人回复
-
基元I理论下三大核心空间现象精准推导与细节解析
已经有0人回复
-
基于基元 I 统一理论的反重力理论推导
已经有0人回复
-
基于基元I统一理论的量子力学本源推导
已经有1人回复
» 本主题相关商家推荐: (我也要在这里推广)
zeoliters
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 贵宾: 0.902
- 金币: 2088.4
- 散金: 16
- 帖子: 445
- 在线: 60.6小时
- 虫号: 772181
- 注册: 2009-05-16
- 专业: 计算化学，量子化学
2楼2009-11-08 18:17:44
sheath
至尊木虫 (知名作家)
- 应助: 0 (幼儿园)
- 金币: 25661.5
- 红花: 2
- 帖子: 5896
- 在线: 416.3小时
- 虫号: 445940
- 注册: 2007-10-28
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
3楼2009-11-08 18:42:41
sheath
至尊木虫 (知名作家)
- 应助: 0 (幼儿园)
- 金币: 25661.5
- 红花: 2
- 帖子: 5896
- 在线: 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 (幼儿园)
- 金币: 25661.5
- 红花: 2
- 帖子: 5896
- 在线: 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 (幼儿园)
- 金币: 25661.5
- 红花: 2
- 帖子: 5896
- 在线: 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