卡开发发 发表于 2025-4-14 08:10
原则上说，这么大的虚频肯定是可忽略的，剩下的看文章怎么写和图怎么出。但先不要急，如果这么看，一个是 ...

好的谢谢卡开老师，我计划先再尝试下波恩修正，这个是重新优化过结构算的声子谱，力收敛到10-7之后一直波动，感觉很难下降了，您说的方式我也尝试一下

zhy996zhy 发表于 2025-4-13 13:12
卡开老师，我重新优化结构，力收敛到-7，计算声子谱putcar中结果如下，
550 f/i=    0.000188 THz     0 ...

原则上说，这么大的虚频肯定是可忽略的，剩下的看文章怎么写和图怎么出。但先不要急，如果这么看，一个是之前的收敛精度可能并不高，而且你之前可能没有重新优化过结构，建议多优化几次来消除一下Pulay Stress，再去算声子比较好，我没办法了解到你全部的计算流程。

卡开发发 发表于 2025-4-6 18:47
确实有可能，但也不排除扩更大胞能够缓解。

卡开老师，我重新优化结构，力收敛到-7，计算声子谱putcar中结果如下，
550 f/i=    0.000188 THz     0.001183 2PiTHz    0.006280 cm-1     0.000779 meV
551 f/i=    0.000510 THz     0.003202 2PiTHz    0.016998 cm-1     0.002107 meV
552 f/i=    0.000652 THz     0.004095 2PiTHz    0.021741 cm-1     0.002696 meV，
想询问下您觉得这个虚频够小可以忽略吗？

卡开发发 发表于 2025-4-6 18:47
确实有可能，但也不排除扩更大胞能够缓解。

好的谢谢老师，唉主要放图上不好看，只看outcar理论上够小了

zhy996zhy 发表于 2025-4-6 15:13
就是利用phonopy得出的band.dat绘图后有最大-0.21THZ的虚频，outcar中只有那三个，我不太清楚是否有可能 ...

确实有可能，但也不排除扩更大胞能够缓解。

卡开发发 发表于 2025-4-5 04:25
有多少虚频呢？

就是利用phonopy得出的band.dat绘图后有最大-0.21THZ的虚频，outcar中只有那三个，我不太清楚是否有可能是phonopy本身插值的问题

数量吗，还是挺多的大概有 50 多条，就在图片中那两个下沉区域

zhy996zhy 发表于 2025-4-4 11:11
谢谢卡开老师，我重新看了下，outcar 没给全，最大虚频在-6cm-1，如果是这个量级呢

有多少虚频呢？

卡开发发 发表于 2025-4-4 09:21
按道理不太会，其实按道理这个虚频足够小了。

谢谢卡开老师，我重新看了下，outcar 没给全，最大虚频在-6cm-1，如果是这个量级呢

zhy996zhy 发表于 2025-4-4 08:58
卡开老师，您觉得这两个虚频问题大吗，是否可以忽略，审稿人会不会一直追问

按道理不太会，其实按道理这个虚频足够小了。

卡开发发 发表于 2025-4-3 22:44
那你可能得考虑对精度和超胞尺寸测试一下了。

卡开老师，您觉得这两个虚频问题大吗，是否可以忽略，审稿人会不会一直追问

zhy996zhy 发表于 2025-4-3 21:35
卡开老师，我重新确定了两个虚频最低点对应坐标，一个在G点和Z点之间，另一个在高对称点D上

那你可能得考虑对精度和超胞尺寸测试一下了。

zhy996zhy 发表于 2025-3-28 14:41
非常感谢卡开老师，希望这次计算顺利吧

卡开老师，我重新确定了两个虚频最低点对应坐标，一个在G点和Z点之间，另一个在高对称点D上

卡开老师，这是高对称点：High-symmetry points (in fractional coordinates).
   0.0000000000   0.0000000000   0.0000000000     GAMMA
   0.0000000000   0.5000000000   0.0000000000     Z
   0.0000000000   0.0000000000   0.5000000000     B
   0.0000000000   0.0000000000  -0.5000000000     B_2
   0.5000000000   0.0000000000   0.0000000000     Y
  -0.5000000000   0.0000000000   0.0000000000     Y_2
   0.5000000000   0.5000000000   0.0000000000     C
  -0.5000000000   0.5000000000   0.0000000000     C_2
   0.0000000000   0.5000000000   0.5000000000     D
   0.0000000000   0.5000000000  -0.5000000000     D_2
  -0.5000000000   0.0000000000   0.5000000000     A
  -0.5000000000   0.5000000000   0.5000000000     E
  -0.5000000000   0.0000000000   0.5000000000     H
  -0.5000000000   0.0000000000   0.5000000000     H_2
  -0.5000000000   0.0000000000  -0.5000000000     H_4
  -0.5000000000   0.5000000000   0.5000000000     M
  -0.5000000000   0.5000000000   0.5000000000     M_2
  -0.5000000000   0.5000000000  -0.5000000000     M_4