VASP+phono3py计算晶格热导率时关于k点以及小虚频对热导率计算结果影响求助

xuewuqiongxiao

请教各位老师：
计算描述：使用vasp6.4.3结合phono3py计算UO3的晶格热导率：
过程描述：
       1、计算过程中使用原胞（4.2*4.2*4.2A）进行结构优化（PBE；k点为15 15 15；EDIFF = 1E-010；EDIFFG = -1E-05）
       2、使用结构优化的结构进行频率分析（PBE；k点为15 15 15；EDIFF = 1E-010；EDIFFG = -1E-05）
       3、使用phono3py -d --dim="2 2 2" -c POSCAR进行扩胞并产生随机位移（共计754个位置文件POSCAR）


         4、使用脚本进行批量任务提交，提取二阶、三阶力常数进行晶格热导率计算
问题描述：
        1、在进行频率分析时，出现虚频，为小虚频，尝试包括以下方法进行虚频消除：（1）将原胞扩胞为2*2*2后进行结构优化、频率分析（出现较大虚频，后使用vaspkit以及ASE进行虚频的消除后出现3个小虚频（2）使用更加严格的收敛限，增加k点数以及使用脚本对原子位置加微小扰动，破坏对称性。仍然出现3个小虚频，提问：如图所示的小虚频对晶格热导率的影响是否可以忽略？
       2、在进行热导率计算时，使用phono3py扩胞为原胞的2*2*2，共得到754个结构，提问：请问是否可以减小K点数（因计算资源有限），可以减小到多少不会对热导率产生显著影响能够算出较为合理的值？