| 查看: 4498 | 回复: 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 师兄说只要有虚频结果就不可靠，请问是哪个参数还是优化的问题呢？小白只能仰仗各位了  |
回复此楼» 猜你喜欢
云南大学材料与能源学院解琳课题组钙钛矿博士招生
已经有7人回复
-
西安交大新媒学院副院长用撤稿论文结题
已经有5人回复
-
论文撤稿了
已经有5人回复
-
某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路！
已经有5人回复
-
26/27申博自荐
已经有9人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有7人回复
-
河北省自然科学基金
已经有7人回复
-
揭秘青基评审内幕：几个A才能顺利中标
已经有4人回复
-
青B发送上会通知了吗
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用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
10