版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

ft1071

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 3.5
帖子: 8
在线: 24.2小时
虫号: 2408439
注册: 2013-04-09
专业: 金属材料表面科学与工程

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

正在计算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

师兄说只要有虚频结果就不可靠，请问是哪个参数还是优化的问题呢？小白只能仰仗各位了

回复此楼

» 猜你喜欢

论文终于录用啦！满足毕业条件了已经有17人回复
不自信的我已经有5人回复
磺酰氟产物，毕不了业了！已经有4人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复

» 本主题相关价值贴推荐，对您同样有帮助:

用VASP-DFPT+phonopy计算石墨烯声子谱求助已经有37人回复
vasp+phonopy计算一个新结构的声子谱有虚频，怎么消掉！！急！！！已经有4人回复

1楼 2016-02-02 12:50:42

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

人人余yan

金虫 (著名写手)

应助: 1 (幼儿园)
金币: 2590.4
散金: 1321
红花: 14
帖子: 1137
在线: 382.3小时
虫号: 2877165
注册: 2013-12-17
专业: 凝聚态物性 II ：电子结构

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

发自小木虫Android客户端

赞一下(2人)

回复此楼

2楼2016-02-02 22:04:23

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

ft1071

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 3.5
帖子: 8
在线: 24.2小时
虫号: 2408439
注册: 2013-04-09
专业: 金属材料表面科学与工程

引用回帖:

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

金虫 (著名写手)

应助: 1 (幼儿园)
金币: 2590.4
散金: 1321
红花: 14
帖子: 1137
在线: 382.3小时
虫号: 2877165
注册: 2013-12-17
专业: 凝聚态物性 II ：电子结构

引用回帖:

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

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 3.5
帖子: 8
在线: 24.2小时
虫号: 2408439
注册: 2013-04-09
专业: 金属材料表面科学与工程

引用回帖:

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