| 查看: 1528 | 回复: 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 |
» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
fileCode有新解读？
已经有16人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
基金中了
已经有14人回复
-
关于filecode
已经有7人回复
-
娱乐
已经有6人回复
2楼2016-10-29 19:50:11
100