| 查看: 1551 | 回复: 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 |
» 猜你喜欢
梦想
已经有8人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
2楼2016-10-29 19:50:11
30