phono3py --dim="1 1 1" --sym-fc -c POSCAR 第五步:计算热导率。 phono3py --fc3 --fc2 --dim="1 1 1" --mesh="19 19 19" -c POSCAR --br | tee out.txt 打开out.txt下拉到底就可以看到计算的热热导率了,直接复制到originlab里作图就可以了。 这是笔者选用不同赝势利用该方法计算的...
当我们查找不到时就需要自己计算出Relaxation time。 晶格热导率的计算 计算晶格热导率我们需要用到的软件包括Phonopy,Thirdorder,ShengBTE。 其中Phonopy用于计算声子谱及二阶力常数,Thirdorder用于计算三阶力常数,ShengBTE用于结合前面两者的结果计算晶格热导率。 1、Phonopy计算声子谱及二阶力常数 计算声子谱及二阶力...
将每个计算目录job-*中的vasprun.xml文件复制到上级目录，改名并编号。接着，执行三条命令以获取二阶力常数矩阵FORCE_CONSTANTS和声子谱数据551.dat。声子谱数据551.dat可以导入OriginLab软件绘制，FORCE_CONSTANTS文件则可用于进一步计算晶格热导率。
看到的都是算晶体的，非晶体既然是无序那怎么可以得到一个原胞呢，没有原胞就没法进行计算呀 ...
计算热导率所需的3个文件，1是CONTROL文件，2是2阶力常数，3是3阶力常数 2阶力常数是通过phonopy对...
其中Phonopy用于计算声子谱及二阶力常数,Thirdorder用于计算三阶力常数,ShengBTE用于结合前面两者的结果计算晶格热导率。 1、Phonopy计算声子谱及二阶力常数 计算声子谱及二阶力常数的具体流程如下: (1)对初始结构进行高精度的结构优化 这一步中INCAR的主要参数是EDIFFG,一般情况下应达到EDIFFG=-1E-8的标准。考虑优...
最后我们通过phonopy里面的QHA模块 去得到此材料的非偕性质 phonopy-qha -p -s e-v.dat *00/thermal_properties.yaml 就可以得到很多热力学性质比如热膨胀系数,格林爱森参数,自由能,等压热容等等 (4)我们想要进一步得到材料的非偕性质热导率等要计算三阶力常数,和步骤(2)里面扩胞一样,我们通过thirdorder扩胞 ...
其中Phonopy用于计算声子谱及二阶力常数,Thirdorder用于计算三阶力常数,ShengBTE用于结合前面两者的结果计算晶格热导率。 1、Phonopy计算声子谱及二阶力常数 计算声子谱及二阶力常数的具体流程如下: (1)对初始结构进行高精度的结构优化 这一步中INCAR的主要参数是EDIFFG,一般情况下应达到EDIFFG=-1E-8的标准。考虑优...
最后我们通过phonopy里面的QHA模块 去得到此材料的非偕性质 phonopy-qha -p -s e-v.dat *00/thermal_properties.yaml 就可以得到很多热力学性质比如热膨胀系数,格林爱森参数,自由能,等压热容等等 (4)我们想要进一步得到材料的非偕性质热导率等要计算三阶力常数,和步骤(2)里面扩胞一样,我们通过thirdorder扩胞 ...
其中Phonopy用于计算声子谱及二阶力常数,Thirdorder用于计算三阶力常数,ShengBTE用于结合前面两者的结果计算晶格热导率。 1、Phonopy计算声子谱及二阶力常数 计算声子谱及二阶力常数的具体流程如下: (1)对初始结构进行高精度的结构优化 这一步中INCAR的主要参数是EDIFFG,一般情况下应达到EDIFFG=-1E-8的标准。考虑优...