24小时热门版块排行榜    

查看: 770  |  回复: 9
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

素婕求学

铁虫 (小有名气)

[交流] 【求助】继续请教周期性边界问题

在我发过的帖子里边我有请教过这个问题,得到了各位好心的大侠的悉心指教,在此再次叩谢!!!!!!!!!!!!!!!
只是,小女子太笨拙,在计算中仍然有迷惑的地方,所以还请各位大侠能再次赐教~~~~~~~~~~~~~~~~~~~~~

以一个系统包含了a,b,c,d,e,f,g个原子为例(a,b,c三个原子在模拟盒左边,d原子在模拟盒中间,e,f,g三个原子在模拟盒右边),这7个原子之间的关系是:a-b,a-c,b-c,c-d 成键,而d-e,e-f,f-g成键,原子之间成键是按照作用距离来判断的。当系统运动一段时间后,g原子从模拟盒右边出来了,按照周期性边界条件,g原子的镜像原子g'从左边进来,g原子的坐标被g'原子的坐标替换,从而对g原子的关注转移到g'上了。在计算力时,则按照最小镜像原则,即如果f-g'的距离大于半个模拟盒边长的长度,则按照f原子与g'原子的镜像原子的距离来计算,这样可以继续保持f-g 成键的事实。但是这里同样出现了一个问题,那就是,g原子变成了g'原子后,它的位置变到了模拟盒左边,这时,它离a,b,c 这些原子很近,按照距离来确定成键的原则,这时a,b,c 可能与g'原子有了相互作用,甚至形成了键,但这是不符合最初系统的构型的,而且我想象不出来g原子因为离开了模拟盒,而能通过周期性边界条件和最小镜像的使用能使得它回到原来的位置去呢,那么如果不能了,这样的构型实质性的发生了变化呢!这样的结果怎么可能正确呢?实在迷惑不解中。。。。。。。。。

[ Last edited by 素婕求学 on 2008-8-15 at 04:16 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

acridine

木虫 (正式写手)

★ ★ ★
lei0736(金币+3,VIP+0):谢谢一直关注
有图就清晰多了,从图上看问题也有可能不是在边界的处理上
12原子运动到了12‘的位置,这个过程其实就相当于6号原子运动到了右侧环中间。而这个过程不是楼主所希望的出现。出现这个现象的可能也有很多种:势函数部分,分子动力学部分等等。当然,也可能是构型自身就是不稳定的。
既然是处理碳的,强烈建议楼主先跑个c60的动力学,不用周期条件或将边界取很大,这样就可以排除处理边界函数的影响。如果计算后c60的结构没有破坏,那么再跑个金刚石,来确定边界是不是出问题了。如果都没问题,就是构型自身的事了
8楼2008-08-16 08:58:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

素婕求学

铁虫 (小有名气)

自己先顶起来~~~~~~~
2楼2008-08-15 04:15:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mengsk

金虫 (著名写手)

我也有相似的疑惑。
奋斗中的小木虫
3楼2008-08-15 09:33:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

素婕求学

铁虫 (小有名气)

引用回帖:
Originally posted by mengsk at 2008-8-15 09:33:
我也有相似的疑惑。

---------------------------------------------------------------------------------------------------------

握手!!!!

希望知道的虫友能慷慨赐教呀~~~~~~~~~~~~
4楼2008-08-15 11:55:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见