卡开发发 发表于 2025-12-18 21:28
这问题我经验有限，或许可以看下动能截断是否充足、超胞扩展是否足够，另外我不确定你用的是有限位移还是 ...

好的老师，我大概明白了，这是个异质结构，z方向残余力一直收敛不下去，没办法放弃他了

zhy996zhy 发表于 2025-12-16 00:01
卡开老师，时隔半年冒昧打扰你一下，这是我新计算的一个体相材料，g点处出现下沉包，且我提高优化精度到1 ...

这问题我经验有限，或许可以看下动能截断是否充足、超胞扩展是否足够，另外我不确定你用的是有限位移还是线性响应法，前者的差分步长也会存在影响。

如图所示

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

卡开老师，时隔半年冒昧打扰你一下，这是我新计算的一个体相材料，g点处出现下沉包，且我提高优化精度到10-9然后再次计算声子谱后虚频几乎没有变化，如图所示，之前的材料在高精度优化后虚频都消失了，这个似乎对结构优化没有反应，我尝试了波恩修正，发现修正后结果更差，对此您有什么建议吗

卡开发发 发表于 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
卡开老师，您觉得这两个虚频问题大吗，是否可以忽略，审稿人会不会一直追问

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