24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

sheath

至尊木虫 (知名作家)

应助: 0 (幼儿园)
金币: 25367.5
红花: 2
帖子: 5848
在线: 416.3小时
虫号: 445940
注册: 2007-10-28
性别: GG
专业: 凝聚态物性I:结构、力学和

[交流] 【讨论】自编程序中正则系综的体系温度为什么不能保持恒定？

在自编的分子动力学程序中，计算原胞的原子体系为正则系综。原子位置在三维方向应用周期性边界条件，恒温采用Nose-Hoover算法。
   但是，程序执行后，出现了体系数目减少，体系温度不断升高，甚至到达上万度的错误结果。

   经过仔细检查程序结构和每个子程序，初步分析认为：在每个时间步，周期性边界条件使跑出计算原胞范围的原子重新回到原胞内的相应位置。而这可能使部分原子重叠，尽管这种情况的几率很小。随后进行“更新原子存在与对应数据”操作，这将可能造成部分重叠原子数据的丢失，当然，原子数目也可能减小。
   也就是说，周期性边界条件造成原子重叠，而原子重叠是导致原子数目减少的关键原因。如果这个分析对的话，如何解决？若不是，那么原因又在哪里？

敬请各位同行不吝赐教！！！

[ Last edited by cadick on 2009-12-14 at 02:30 ]

回复此楼

» 猜你喜欢

求国际会议网站已经有1人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有254人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有19人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
26申博推荐：南京航空航天大学国际前沿科学研究院光学方向招收博士生！已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复

» 本主题相关商家推荐: (我也要在这里推广)

1楼2009-11-08 16:46:45

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

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

引用回帖:

Originally posted by sheath at 2009-11-9 17:20:

我的程序在计算原子位置和相互作用力都使用了周期性边界条件。

计算能量和相互作用力都要用周期性边界条件
我是说在使用周期性边界条件时，只调用而不改变粒子的原坐标

引用回帖:

Originally posted by sheath at 2009-11-9 17:20:

"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."

第一句没有问题，第二句值得商榷
１、在周期性边界条件下，同在原胞下的两个粒子是不能通过直接计算它们之间距离来计算能量和受力作用的。把原胞外的粒子放回原胞的操作虚耗了计算量，还可能出错
２、对于高分子体系，这样的操作会破坏共价键

Allen 的＂Computer simulation of liquid＂２９－３１页讲了关于周期性边界条件的问题，建议看下

[ Last edited by dellus on 2009-11-10 at 13:20 ]