| 查看: 1542 | 回复: 1 | |||
[求助]
【求助】偶分布函数g(r)截断半径如何选取？横坐标如何取到10A？
|
|
各位前辈好，我计算金属二元体系液态下的偶分布函数g(r)，构型共100个原子，盒子边长为12.20310 A，想要实现g(r)横坐标最大到 10 A。 具体步骤如下： 1、在Materials Studio中用Amorphous Cell模块构建室温下100个原子的初始结构，根据密度确定晶格常数 2、VASP在熔点以上做恒温动力学（IBRION=0，不改变晶胞大小和形状） 3、分割XDATCAR文件得到多个构型，计算g(r) 用来计算g(r)的构型如下图：  晶格常数为 12.20310 A 如下图：  计算g(r)时，默认截断半径取盒子边长（即晶格常数）的一半，所以g(r)曲线横坐标最大到 6 A 左右，如下图中的红色曲线所示：  而文献中的g(r)横坐标最大到 10 A，我就试着取g(r)截断半径为 10 A，如上图黑色曲线所示。6 A 之前基本重合，但到最后g(r)逐渐减小为0，而不是1，所以肯定是错的。 也试着用一个动力学之后的构型建立了2*2*2的超胞（800个原子）算g(r)，还是不对。 g(r)截断半径是不是只能取盒子边长的一半？盒子边长是不是不能随便乱改？而初始结构盒子边长由密度确定，又该如何使g(r)最大半径到10呢？ 另外，各位前辈有分析g(r)的程序可以推荐吗？ 我理解不到位的地方，也请多多指正。 谢谢！ |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理计算-tg |
» 猜你喜欢
93BebMhtakh前后11位开头都是大写
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有86人回复
-
感觉是下周放榜了
已经有11人回复
-
filecode=后面第一个是大写字母
已经有10人回复
-
时间戳变了，能看出什么问题？
已经有4人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
2楼2016-10-29 19:50:11
