博致远 发表于 2016-10-10 21:01
感谢老师的指点，在细问一下细节；
1、对于价态的影响，我又进行了以下测试，同样fine精度下对于+1、+2 ...

1、MS显然还没有支持约束密度泛函理论的方案，按说“价态”的设置对第一性原理计算无影响。castep的初猜是随机数，据说确实有人发生过前后两次计算结果有些许差异的情况，保险起见我很建议你前后测试的时候采用的random seeds采用相同的编号看看情况如何。

2、带电的周期性系统的能量计算和中心计算能量要进行相比的话需要进行修正，有人对此作了一些工作，可以查文献找找，这是掺杂计算的常规情形之一。

3、VCA≡virtual crystal approximation。

卡开发发 发表于 2016-10-10 12:03
MS的properties中的价态对第一性原理的模块做计算都没有影响，这样后续问题自然迎刃而解。如果存在变价问题 ...

感谢老师的指点，在细问一下细节；
1、对于价态的影响，我又进行了以下测试，同样fine精度下对于+1、+2价铜的CZTS模型进行了一下结果优化和布居分析，结果二者晶格常数分别为5.47795和5.47473，而且布局分析发现Cu的原子态分布也分别是-0.050和-0.040，不知道是否算是合理的误差范围；

2、对于整体胞带电手动计算的问题，这个正确性或者说可信度有多少？因为我本来直接想用这种掺杂相跟其他可能第二相做个能量对比，但您这么一提带电体系对于能量影响大，所以想问一下这个有可比性么.

PS:那个VCA的中文是？
麻烦老师

zyj19831206 发表于 2016-10-10 10:51
这个是substitution.pl那个脚本吧？

对的

MS的properties中的价态对第一性原理的模块做计算都没有影响，这样后续问题自然迎刃而解。如果存在变价问题，需要改变的是晶格整体（或者说system）的charge，针对周期性引进电荷引起发散的问题程序会根据jellium近似来进行抵消，但电荷引起能量的变化可能需要手工计算（公式就不写了，可以直接参考那些价态不同元素掺杂的文章），至于原子带多少电量可以通过布居分析。价态变化的情况我不确定VCA能不能给出好结果，最好找找文献确认一下。

这个是substitution.pl那个脚本吧？