应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 【求助】MD模拟中的截断半径？
51/1返回列表
查看: 1468  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

ssttssqq9118

新虫 (初入文坛)

[交流] 【求助】MD模拟中的截断半径？

MD模拟中截断半径选取是模拟原胞中边长的一半还是模拟的超级晶胞边长的一半？选取经验势是否考虑次近邻的相互作用？

[ Last edited by lei0736 on 2009-11-25 at 12:14 ]
回复此楼

» 收录本帖的淘帖专辑推荐

分子动力学

» 猜你喜欢
1楼 2009-05-15 10:46:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lei0736

荣誉版主 (职业作家)

优秀版主
★ ★ ★
gwdavid(金币+3,VIP+0):3q 5-29 15:36
引用回帖:
Originally posted by ssttssqq9118 at 2009-5-15 10:46:
MD模拟中截断半径选取是模拟原胞中边长的一半还是模拟的超级晶胞边长的一半？选取经验势是否考虑次近邻的相互作用？

取势的结尾吧  盒子的一半必须大于势的截尾 当然盒子越大越好
考不考虑次近邻是由你选取的势形式决定的啊
你选的那种有次近邻就有 没有就没有。。。
比如原型EAM 修正EAM 和更新的考虑第二近邻的MEAM
赞一下(10人) 回复此楼
2楼2009-05-15 11:02:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ssttssqq9118

新虫 (初入文坛)
那如果是Buckingham势，应该怎么取？
赞一下 回复此楼
3楼2009-05-29 11:42:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
lei0736(金币+2,VIP+0):谢谢 5-30 11:33
不管是什么势，都要有一个截断半径的。截断半径是势函数的一部分。文献上如果没有截断半径，则要自己选择，但在周期性边界条件的情况下，截断半径不应超过MD胞长度的一半。
赞一下(12人) 回复此楼
4楼2009-05-29 15:43:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)
★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
lei0736(金币+8,VIP+0):谢谢 5-30 15:17
引用回帖:
Originally posted by 老虎大王 at 2009-5-29 15:43:
不管是什么势，都要有一个截断半径的。截断半径是势函数的一部分。文献上如果没有截断半径，则要自己选择，但在周期性边界条件的情况下，截断半径不应超过MD胞长度的一半。

我晕。我说错了，把势函数截断与RDF的截断混在一起说了。当时回这个帖子的时候，我刚上传了一个计算RDF（径向分布函数，或者双体分布函数等）的程序。还在回味之中。呵呵。

应该是，在周期性边界条件下，求算RDF时，截断距离不应超过元胞最短边长的一半。势函数的截断距离则不是这样。呵呵。为什么呢？因为——

RDF在几何上与球壳相联系（归一化因子里面包含球壳体积），如果你截断距离超过了元胞最短边长的一半，就会造成在这个距离以外的原子无法形成一个完整的球壳，归一化就归错掉了。这一点，你画个正方形，再画个内切圆，就知道了。中心一点的原子与内切圆以外的原子（比如靠近顶点的原子）之间的配位情况是不能统计的，因为它们构不成完整球壳。

但势函数截断则不同，势函数讲究的是，要把所有的在截断距离以内的近邻原子都统计在内，只要你在截断距离之内，我就计算你的能量贡献，截断距离大过了MD胞的半边长，只是说明MD胞较小，但不会出现漏计的情况。这是没有问题的。截断距离通常是势函数的一部分，而势函数是先于模拟而确定的。你可以根据势函数的截断选取MD胞的大小，却不应根据MD胞的大小去选择势的截断。所以正确的说法是：最好让Md胞的最短边长不小于文献给出的势函数的截断距离。

要注意的是，对于模拟离子体系的长程库仑势和短程BMH势之类的势函数，有时候文献上没有截断距离，这时候，一般可以把截断距离取为MD胞边长的一半。如果取的截断大于边长的一半会出现什么情况呢？那样的话，当第二次模拟，放大MD胞时，可能计算得到的能量会不相同。但是问题不会太大。因为短程排斥力的作用距离都是很短的，而且随距离衰减得很厉害，只要你的MD胞不取得过分小，是不会有多大影响的。也就是说，对于BMH之类的短程势，你的截断半径取10ang，与取15ang，实际上能量计算结果不会有什么区别。至于长程库仑势以及色散项(r^-6之类)等，通常的可以处理这类长程力的程序，一般都会有长程修正项。所以早期的MD模拟，对离子体系也只用216个原子，Md胞相当小，但结果也可以是比较安全的。

[ Last edited by 老虎大王 on 2009-5-30 at 14:49 ]
赞一下(12人) 回复此楼
5楼2009-05-30 14:29:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ssttssqq9118 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭