| 查看: 769 | 回复: 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 ] |
回复此楼» 猜你喜欢
夏威夷大学生物胶体与界面实验室诚招全奖博士生和博士后
已经有1人回复
-
200年来人们一直无法破解的难题： 气液临界点状态方程
已经有7人回复
-
物理化学论文润色/翻译怎么收费?
已经有153人回复
-
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答！
已经有10人回复
-
求助立方晶体In(OH)3的CIF卡片
已经有2人回复
-
请问四氢呋喃溶解的聚合物用甲醇沉淀时，如何使沉淀过程加速？
已经有2人回复
-
[电子教材]天津大学物理化学（第七版）上、下册
已经有1人回复
-
阴离子交换膜电解二氧化碳还原
已经有0人回复
-
Chemical Bonding at Surfaces and Interfaces,最经典的一本表面上化学相互作用教材
已经有0人回复
acridine
木虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 4901.6
- 红花: 13
- 帖子: 706
- 在线: 389.5小时
- 虫号: 305478
- 注册: 2006-12-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★ ★
lei0736(金币+3,VIP+0):谢谢
lei0736(金币+3,VIP+0):谢谢
|
>但是这里同样出现了一个问题，那就是，g原子变成了g'原子后，它的位置变到了模拟盒左边，这时，它离a,b,c 这些原子很近，按照距离来确定成键的原则，这时a,b,c 可能与g'原子有了相互作用，甚至形成了键，但这是不符合最初系统的构型的... 没太弄清楼主的疑问在哪。在你的初始构型下，把a放在了左侧，g在右侧，那么正如楼主所说，根据最小镜像原则，它们初始的距离就是很近，之后在运动中成键也是完全可能的，可以说与初始构型并没有矛盾。如果楼主所希望的是a和g远离，那么是你的初始结构设置有问题。 总之关键就在于：在楼主所设置的初始形貌下，a和g是很近的 [ Last edited by acridine on 2008-8-15 at 19:50 ] |
5楼2008-08-15 19:48:18
2楼2008-08-15 04:15:00
3楼2008-08-15 09:33:38
4楼2008-08-15 11:55:33