请教一下用VASP+Phonopy-QHA,计算热膨胀系数的嘛?我算完声子后，明明mesh.conf,文件里面温度设置的最高温度是2000，但是最后在用phonopy-qha的时候热膨胀系数只到了1000，哪里出错了呢，...

最近被phonopy-qha的这个准简谐近似搞得迷迷糊糊，有没有哪位大神亲自指导一下！下面是我的理解：1.使用phonopy做完声子谱之后，调节POSCAR（原胞，未扩胞）获得不同体积的POSCAR，然后固定...

File"/opt/software/phonopy-1.8.4.1/bin/phonopy-qha",line277,in)File"/opt/software/phonopy-1.8.4.1/bin/phonopy-qha",line113,inread_data_from_thermal_propertiesentropy=np.array...

1.dim="222"pa="01/21/21/201/21/21/20"mp="414141"nunit=32计算QHA的时候，nunit的数值怎样取？选默认值32？我发现取不同的值拟合出来的数值精度上会有差别？2.T=0K时，拟合主来的吉布斯...

取了10个体积点1.为什么虚线很不平滑？这是什么原因造成的？2.为什么热膨胀系数下降？本来应该是缓慢上升的埃eos.png

在Phonopy的准谐振近似QHA计算中，输出一个与熵有关的文件entropy-volume.dat内容如下,(列了一个温度下数据).#temperature298.000000000000000(温度K)#0.000004028697269-0....

我先用phonopy–t-tmax=1000–tstep=100得到thermal_properties.yaml文件，然后再命令phonopy-qhae-v.datthermal_properties-{1.16}.yaml-s画图，最后怎么图中QHA只给出一条来埃急求

想用phonopy-qha产生多个体积进行计算。其默认值为0.95~1.05,间隔为0.01，共十个体积。想扩大范围且产生十个以上的体积，不知道如何设置-cu=后面的内容。

最近想用phonopy计算热力学性质，直接计算出来的是赫姆霍兹自由能，发现想算吉布斯自由能的话需要用到QHA（quasiharmonicapproximation），看手册上最后的命令是phonopy-qhae-v.datthermal_...