24小时热门版块排行榜

返回列表

当前主题已经存档。

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

素婕求学

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 9
帖子: 65
在线: 1.8小时
虫号: 535493
注册: 2008-03-29
性别: MM
专业: 分子动力学模拟/纳米流

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

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

以一个系统包含了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 ]

回复此楼

» 猜你喜欢

1楼 2008-08-15 04:03:56

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

素婕求学

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 9
帖子: 65
在线: 1.8小时
虫号: 535493
注册: 2008-03-29
性别: MM
专业: 分子动力学模拟/纳米流

引用回帖:

Originally posted by acridine at 2008-8-15 19:48:
>但是这里同样出现了一个问题，那就是，g原子变成了g'原子后，它的位置变到了模拟盒左边，这时，它离a,b,c 这些原子很近，按照距离来确定成键的原则，这时a,b,c 可能与g'原子有了相互作用，甚至形成了键，但这 ...

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

前面的回复不太清楚，贴个图来说明一下吧！
如下图所示，原子12离开了模拟盒，其镜像原子12' 从模拟盒左边又回到了模拟盒。原子12最初时和10，11，2这三个原子都成键，当它离开模拟盒后，其镜像原子的位置处于一个石墨原胞之内，12' 可能与2，1，4分别成键，那这样的话，不是就不满足最初的模型了吗？而且，当我用visualization的软件处理时，也同样会出错，因为原子12的坐标换成了12'，显示原子空间图时，原子12的位置就变成了如图示12’的位置，这样的结构图看起来就很不对了。
哎，问题好多呀，希望我能说明白了，看明白了的高手如果知道问题所在的话，麻烦指点一下，感激不尽。。。。。。。。。。。。。。。。