24小时热门版块排行榜

返回列表

lqtsdu

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 243.5
帖子: 13
在线: 43小时
虫号: 4063557
注册: 2015-09-10
性别: MM
专业: 金属材料的微观结构

[求助] 【求助】偶分布函数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)的程序可以推荐吗？

我理解不到位的地方，也请多多指正。

谢谢！

回复此楼