频率分析计算时虚频过大怎么办
13 f/i= 0.295923 THz 1.859342 2PiTHz 9.870942 cm-1 1.223841 meV
X Y Z dx dy dz
0.000000 0.000000 15.000000 -0.008433 0.022563 -0.307981
0.935885 1.621000 12.403390 -0.015800 0.032033 -0.401723
1.871770 0.000000 17.596610 -0.013088 0.028841 -0.438809
0.935885 1.621000 16.290960 -0.016631 0.043391 -0.534236
1.871770 0.000000 13.709040 -0.018835 0.046475 -0.508204
14 f/i= 3.550214 THz 22.306655 2PiTHz 118.422402 cm-1 14.682512 meV
X Y Z dx dy dz
0.000000 0.000000 15.000000 -0.549241 0.821441 0.025603
0.935885 1.621000 12.403390 -0.072545 0.065059 0.009396
1.871770 0.000000 17.596610 -0.070262 0.064003 0.008797
0.935885 1.621000 16.290960 -0.039250 0.009443 0.015130
1.871770 0.000000 13.709040 -0.044662 0.010918 0.015682
15 f/i= 4.741184 THz 29.789735 2PiTHz 158.148854 cm-1 19.607966 meV
X Y Z dx dy dz
0.000000 0.000000 15.000000 0.830121 0.552150 0.013299
0.935885 1.621000 12.403390 0.034219 0.033258 0.000083
1.871770 0.000000 17.596610 0.029014 0.034399 -0.002635
0.935885 1.621000 16.290960 -0.031509 -0.004526 0.002676
1.871770 0.000000 13.709040 -0.021512 -0.005668 0.005238
该结构虚频过大,已经把第一个原子位置移动过重新进行优化,频率计算,结果虚频还是很大,该怎么办? 返回小木虫查看更多
可能本身就是不稳定的吧
不擅长声子计算,但算过几次,解决方法
0. 再次优化
1. 扩胞
2. 好像有个什么冷冻原子法一类的可以消虚频,但我不会
你这问题说得不清楚,到底是虚频多了还是不该有虚频?
如果你是优化结构的,首先要保证模型是合理的,初始结构不合理优化未必收敛到你要的结构,其次,如果初始结构合理,可以提高收敛精度继续优化试一下,或者按照相应虚频给出的dx,dy,dz,取其一半,添加到POSCAR相应原子坐标上去,再优化。注意坐标形式,如果POSCAR里是分数坐标,则取dx.dy.dz一半后须转换为分数坐标单位再相加。可以从最大虚频尝试一下,
如果是过渡态消除多余虚频,也可以仿照微调坐标来做
你的体系里不会只有5个原子吧?如果是,你需要固定其中一个
提高精度再算