linuxprobe 发表于 2024-2-2 17:24
运行命令错了吧？phonopy --dim="2 2 2" -c POSCAR-unitcell band.conf

https://zhuanlan.zhihu.com/p/481680637

运行命令错了吧？phonopy --dim="2 2 2" -c POSCAR-unitcell band.conf

wmg166 发表于 2024-2-1 19:18
计算过程应该没错，Log文件提示：
Degree of freedom:  18/ 18
       N       E                      ...

如果没额外做固定，有可能两个程序对对称性的识别（尤其是声子计算等价的原子）判断不一致，但具体我很难说，实在不行给phonopy发邮件看看。

卡开发发 发表于 2024-2-1 13:09
问题应该在这：

224原子做出来力常数的矩阵是(28,28)，程序认为维度不匹配，但原因不明，到底是前后两次 ...

计算过程应该没错，Log文件提示：
Degree of freedom:  18/ 18
       N       E                     dE             d eps       ncg     rms          rms(c)
RMM:   1    -0.353956298061E-01   -0.35396E-01   -0.21148E-02 76507   0.391E-01
RMM:   2    -0.341378790009E-01    0.12578E-02   -0.12596E-03 46744   0.166E-01    0.117E-01
RMM:   3    -0.349997938971E-01   -0.86191E-03   -0.36472E-03 66756   0.131E-01    0.268E-01
RMM:   4    -0.338741490041E-01    0.11256E-02   -0.14831E-03 55473   0.101E-01    0.539E-02
RMM:   5    -0.339155098934E-01   -0.41361E-04   -0.54067E-04 63896   0.523E-02    0.729E-02
RMM:   6    -0.338575129263E-01    0.57997E-04   -0.13108E-04 55897   0.328E-02    0.336E-02
RMM:   7    -0.338414852637E-01    0.16028E-04   -0.27214E-05 51368   0.174E-02    0.148E-02
RMM:   8    -0.338323703054E-01    0.91150E-05   -0.32378E-05 64122   0.125E-02    0.732E-03
RMM:   9    -0.338306850702E-01    0.16852E-05   -0.36758E-06 56699   0.523E-03    0.442E-03
RMM:  10    -0.338283559383E-01    0.23291E-05   -0.15780E-06 56855   0.397E-03    0.265E-03
RMM:  11    -0.338275219068E-01    0.83403E-06   -0.62293E-07 62076   0.202E-03    0.135E-03
RMM:  12    -0.338270200208E-01    0.50189E-06   -0.21477E-07 54079   0.134E-03    0.517E-04
RMM:  13    -0.338267865625E-01    0.23346E-06   -0.70327E-08 68040   0.593E-04    0.505E-04
RMM:  14    -0.338264659903E-01    0.32057E-06   -0.21780E-07 62124   0.314E-04    0.311E-04
RMM:  15    -0.338263374331E-01    0.12856E-06   -0.36478E-07 65957   0.242E-04    0.172E-04
RMM:  16    -0.338262736454E-01    0.63788E-07   -0.48133E-07 69957   0.182E-04    0.559E-05
force on displaced ion      25 direction 3 :  -0.700  -0.056  -4.384
Linear response finished

问题应该在这：

Number of atoms in supercell (224) is not consistent with the matrix shape of
force constants (28, 28) read from FORCE_CONSTANTS.

224原子做出来力常数的矩阵是(28,28)，程序认为维度不匹配，但原因不明，到底是前后两次对称性识别的不一样还是计算流程出现问题暂时不好判断。