24小时热门版块排行榜    

查看: 4202  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

ft1071

新虫 (初入文坛)

[求助] vasp+phonopy计算声子谱虚频求助 已有2人参与

正在计算Ta2C(六方结构)的热膨胀系数,首先要计算11组体积不同(改变POSCAR基矢缩放系数)的声子谱
然后在缩放比较小的时候就在G点附近产生了虚频
缩放0.95:
vasp+phonopy计算声子谱虚频求助
0.95.jpg
缩放1.1:
vasp+phonopy计算声子谱虚频求助-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

师兄说只要有虚频结果就不可靠,请问是哪个参数还是优化的问题呢?小白只能仰仗各位了
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ft1071

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 人人余yan at 2016-02-02 22:04:23
结构优化的时候精度可以再高一些吧。十的负七或者负八

请问EDIFFG需要改吗?
3楼2016-02-03 02:19:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

人人余yan

金虫 (著名写手)

结构优化的时候精度可以再高一些吧。十的负七或者负八

发自小木虫Android客户端
2楼2016-02-02 22:04:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

人人余yan

金虫 (著名写手)

引用回帖:
3楼: Originally posted by ft1071 at 2016-02-03 02:19:02
请问EDIFFG需要改吗?...

精度高些更准确,结构太大的话算的慢,也可以是一下ibrion=8/6,为了消除虚频的话

发自小木虫IOS客户端
4楼2016-02-03 05:27:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ft1071

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by 人人余yan at 2016-02-03 05:27:52
精度高些更准确,结构太大的话算的慢,也可以是一下ibrion=8/6,为了消除虚频的话
...

这周我增加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.
5楼2016-02-03 17:12:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见