LDA+U体系用dynaphopy计算声子谱求助

卡开发发

zhy996zhy 发表于 2025-3-28 11:12
谢谢卡老师，我ALGO用的是all，我觉得可能和这个有关，也就是说他没有启动微扰吗，昨晚我重新跑了一个， ...

我不排除vasp内部可能ALGO等于特定值会强制NSW=0，我有怀疑过这件事，但搜了下源码中NSW和IBRION等被赋值0的情况没有找到。

当然按照道理说OUTCAR开头的INCAR信息应该和你输入是一致的，否则OUTCAT开头INCAR信息的IBRION=8按照正常逻辑也会被打印为IBRION=-1，但实际上这里并没有，而是在后续参数check的部分显示IBRION=-1，这里应该是在检查过NSW=0后被改为了IBRION=-1。

所以或许你应该再检查一下那个任务的INCAR，或者重新提交试试，具体机制我不明确。

zhy996zhy

卡开发发 发表于 2025-3-28 11:54
我不排除vasp内部可能ALGO等于特定值会强制NSW=0，我有怀疑过这件事，但搜了下源码中NSW和IBRION等被赋值 ...

我明白了老师，我再检查一下，重新跑那个OUTCAR 中IBRION是8，NSW是1，也有以下输出：Linear response reoptimize wavefunctions to high precision
DAV:   1    -0.731358150771E+03   -0.22685E-06   -0.44903E-09  4320   0.844E-05
DAV:   2    -0.731358150771E+03    0.36380E-09   -0.32204E-09  4320   0.502E-05
DAV:   3    -0.731358150770E+03    0.21100E-09   -0.18410E-09  4320   0.268E-05
Linear response DOF=          72
Linear response progress:
  Degree of freedom:   1/ 72
       N       E                     dE             d eps       ncg     rms          rms(c)
RMM:   1    -0.248471924827E-01   -0.24847E-01    0.19780E-02 50396   0.157E-01
RMM:   2    -0.230871519418E-01    0.17600E-02   -0.33716E-03 36088   0.451E-02    0.706E-02
RMM:   3    -0.233513401164E-01   -0.26419E-03   -0.79871E-04 41575   0.155E-02    0.385E-02
RMM:   4    -0.236379966169E-01   -0.28666E-03   -0.15534E-04 42126   0.717E-03    0.200E-02
RMM:   5    -0.238251825231E-01   -0.18719E-03   -0.64102E-05 40427   0.493E-03    0.501E-03
RMM:   6    -0.238101303555E-01    0.15052E-04   -0.35627E-05 53008   0.147E-03    0.306E-03
RMM:   7    -0.237959157941E-01    0.14215E-04   -0.21905E-05 52190   0.939E-04    0.621E-04
RMM:   8    -0.237979507496E-01   -0.20350E-05   -0.21008E-05 59352   0.562E-04    0.262E-04
RMM:   9    -0.237976832849E-01    0.26746E-06   -0.20606E-05 63412   0.507E-04    0.865E-05
RMM:  10    -0.237974415404E-01    0.24174E-06   -0.18336E-05 66204   0.492E-04    0.358E-05
RMM:  11    -0.237974015809E-01    0.39959E-07   -0.18550E-05 67839   0.485E-04    0.146E-05
RMM:  12    -0.237973589879E-01    0.42593E-07   -0.16976E-05 68205   0.486E-04    0.125E-05
first order change of Fermi-energy -6.076101350771523E-003

卡开发发

zhy996zhy 发表于 2025-3-28 12:16
我明白了老师，我再检查一下，重新跑那个OUTCAR 中IBRION是8，NSW是1，也有以下输出：Linear response re ...

是啊，要开始声子计算NSW至少是大于0的值，否则IBRION会被自动改成-1。

zhy996zhy

卡开发发 发表于 2025-3-28 12:59
是啊，要开始声子计算NSW至少是大于0的值，否则IBRION会被自动改成-1。

非常感谢卡开老师，希望这次计算顺利吧

zhy996zhy

卡开发发 发表于 2025-3-28 12:59
是啊，要开始声子计算NSW至少是大于0的值，否则IBRION会被自动改成-1。

卡开老师，抱歉想在问个问题，这是我的outcar文件声子谱计算结果有三个虚频，以及绘制的声子谱，请问这种大小的虚频有可能消除吗？如果要消除，可以往那个方向努力呢？

zhy996zhy

卡开老师，这是高对称点：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

zhy996zhy

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

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

卡开发发

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

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

zhy996zhy

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

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

卡开发发

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

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

zhy996zhy

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

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

卡开发发

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

有多少虚频呢？

zhy996zhy

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

zhy996zhy

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

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

卡开发发

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

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