24小时热门版块排行榜    

查看: 4200  |  回复: 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的回帖

人人余yan

金虫 (著名写手)

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

发自小木虫Android客户端
2楼2016-02-02 22:04:23
已阅   回复此楼   关注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的回帖

人人余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的回帖

小虫虫LG

木虫 (正式写手)

【答案】应助回帖

★ ★ ★
ft1071: 金币+3 2016-02-14 19:28:46
一个办法是缩放系数变小点把。比如0.99
6楼2016-02-11 16:28:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

navyzhang

木虫 (正式写手)

引用回帖:
6楼: Originally posted by 小虫虫LG at 2016-02-11 16:28:02
一个办法是缩放系数变小点把。比如0.99

请问缩放系数在哪里设置啊?
我用phonopy除了改变K点路径,都用的是默认设置
7楼2016-06-24 15:54:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Katherinexhl

新虫 (初入文坛)

蛋蛋小童鞋: 应助指数-1 2016-08-04 10:38:10
引用回帖:
7楼: Originally posted by navyzhang at 2016-06-24 15:54:52
请问缩放系数在哪里设置啊?
我用phonopy除了改变K点路径,都用的是默认设置...

请问问题解决了吗?这11个体积用什么优化的呢?是3吗?
8楼2016-08-01 17:07:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Katherinexhl

新虫 (初入文坛)

【答案】应助回帖

请问您是用ISIF=?优化的?
9楼2016-08-05 20:01:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ft1071 的主题更新
信息提示
请填处理意见