| 查看: 1572 | 回复: 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 |
» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有13人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有6人回复
-
申请科研项目纯靠实力能成功吗？
已经有11人回复
-
极限拉扯，如何抉择
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
2楼2016-10-29 19:50:11
200