| 查看: 4811 | 回复: 8 | ||
[求助]
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 师兄说只要有虚频结果就不可靠，请问是哪个参数还是优化的问题呢？小白只能仰仗各位了  |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有14人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有16人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有15人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有17人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有21人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有21人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用VASP-DFPT+phonopy计算石墨烯声子谱求助
已经有37人回复
- vasp+phonopy计算一个新结构的声子谱有虚频，怎么消掉！！急！！！
已经有4人回复
人人余yan
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 2590.4
- 散金: 1321
- 红花: 14
- 帖子: 1137
- 在线: 382.3小时
- 虫号: 2877165
- 注册: 2013-12-17
- 专业: 凝聚态物性 II ：电子结构
2楼2016-02-02 22:04:23
3楼2016-02-03 02:19:02
人人余yan
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 2590.4
- 散金: 1321
- 红花: 14
- 帖子: 1137
- 在线: 382.3小时
- 虫号: 2877165
- 注册: 2013-12-17
- 专业: 凝聚态物性 II ：电子结构
4楼2016-02-03 05:27:52
|
这周我增加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
小虫虫LG
木虫 (正式写手)
- 应助: 114 (高中生)
- 金币: 1775.9
- 散金: 1871
- 红花: 25
- 帖子: 975
- 在线: 678.5小时
- 虫号: 2128430
- 注册: 2012-11-15
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2016-02-11 16:28:02
navyzhang
木虫 (正式写手)
- 应助: 20 (小学生)
- 金币: 1411.8
- 散金: 1197
- 红花: 5
- 帖子: 640
- 在线: 187.1小时
- 虫号: 883307
- 注册: 2009-10-25
- 性别: GG
- 专业: 理论和计算化学
7楼2016-06-24 15:54:52
Katherinexhl
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 7
- 帖子: 27
- 在线: 39.2小时
- 虫号: 4037098
- 注册: 2015-08-25
- 专业: 凝聚态物性 II ：电子结构
8楼2016-08-01 17:07:47
Katherinexhl
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 7
- 帖子: 27
- 在线: 39.2小时
- 虫号: 4037098
- 注册: 2015-08-25
- 专业: 凝聚态物性 II ：电子结构
9楼2016-08-05 20:01:17
50