vasp+phonopy计算声子谱虚频求助
正在计算Ta2C（六方结构）的热膨胀系数，首先要计算11组体积不同（改变POSCAR基矢缩放系数）的声子谱
然后在缩放比较小的时候就在G点附近产生了虚频
缩放0.95：
0.95.jpg
缩放1.1：
1.1.jpg
缩放1.05没有虚频
扩胞3 3 3 ，INCAR参照phonopy手册， K点取3 3 3
计算声子谱的INCAR：
PREC = Accurate
IBRION = -1
ENCUT = 500
EDIFF = 1.0e-08
ISMEAR = 0; SIGMA = 0.01
IALGO = 38
LREAL = .FALSE.
ADDGRID = .TRUE.
LWAVE = .FALSE.
LCHARG = .FALSE.
之前优化的INCAR：
SYSTEM=TaC
ENCUT=550
ISMEAR = 1
SIGMA = 0.01
NSW=41
IBRION = 2
ISIF = 3
EDIFF=1.0E-4
EDIFFG=-0.001
PREC =Accurate
师兄说只要有虚频结果就不可靠，请问是哪个参数还是优化的问题呢？小白只能仰仗各位了
返回小木虫查看更多
今日热帖
京公网安备 11010802022153号
结构优化的时候精度可以再高一些吧。十的负七或者负八
请问EDIFFG需要改吗？
这周我增加K点到5 5 5算了两个点，还是有这样的虚频。
重新优化了一下原胞，这是原先的优化：
Ta C
1.00000000000000
3.1244339623514183 -0.0000000000000000 0.0000000000000000
-1.5622169811757092 2.7058391838478966 -0.0000000000000000
-0.0000000000000000 -0.0000000000000000 4.9634578427587872
2 1
Direct
0.3333333429999996 0.6666666870000029 0.2526072086980557
0.6666666269999979 0.3333333129999971 0.7473927913019447
-0.0000000000000000 -0.0000000000000000 -0.0000000000000000
0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
EDIFF=10E-8，EDIFFG=-10E-6得到的原胞：
Ta C
1.00000000000000
3.1237671313434077 -0.0000000000000000 -0.0000000000000000
-1.5618835656717038 2.7052616912549250 -0.0000000000000000
0.0000000000000000 0.0000000000000000 4.9634839229430820
Ta C
2 1
Direct
0.3333333429999996 0.6666666870000029 0.2526028572060771
0.6666666269999979 0.3333333129999971 0.7473971427939231
0.0000000000000000 -0.0000000000000000 0.0000000000000000
0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
因为计算很费时，想先问您一下这个差别算大吗？别的参数是否需要修改呢？
——————————————————————————————————————————
另外IBRION=8/6 的话，就是参照phonopy手册里的VASP-DFPT & phonopy calculation 部分？这是手册给的INCAR，把IBRION改成6就可以了吗？
PREC = Accurate
ENCUT = 500
61
phonopy manual, Release 1.7.4
IBRION = 8
EDIFF = 1.0e-08
IALGO = 38
ISMEAR = 0; SIGMA = 0.1
LREAL = .FALSE.
ADDGRID = .TRUE.
LWAVE = .FALSE.
LCHARG = .FALSE
，
一个办法是缩放系数变小点把。比如0.99
请问缩放系数在哪里设置啊？
我用phonopy除了改变K点路径，都用的是默认设置
请问问题解决了吗？这11个体积用什么优化的呢？是3吗？