vasp 使用VCA方法时，如何优化radial cutoffs？

Shana

各位老师好，小弟最近在使用VCA计算固溶体体系，参考vasp的VCA使用教程设置了POSCAR和POTCAR，在INCAR中加入了VCA的参数。具体情形如下图所示： .
但是在计算固溶体系晶格常数的时候出现了问题.
当我设置VCA=1.0 1.0 0.0 1.0 1.0的时候,能够得到LiScGeO的晶格常数, VCA=1.0 0.0 1.0 1.0 1.0的时候,得到LiInGeO的晶格常数,其中LiInGeO的晶格常数大于LiScGeO的晶格常数.
当我设置In和Sc混合占据的时候, 哪怕是VCA=1.0 0.9 0.1 1.0 1.0, 得到的晶格常数也是大于LiInGeO和LiScGeO的.晶格常数并没有随着In的比例增加逐渐增加到LiInGeO.而是直接大于两个边界的晶格常数.最后得到的结果如下图所示: ,实验结果是近似于线性变化的,基本满足Vegard's law ,而VCA计算的结果,则呈现非线性变化,不仅和实验不符,也和基本的预期变化趋势不符.
请问各位老师,这个应该怎么调整,才能让晶格常数至少按照线性趋势变化?
vaspwiki里提到需要优化radial cutoffs, 这个应该设置哪些参数才能达到?