| 查看: 4202 | 回复: 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 师兄说只要有虚频结果就不可靠，请问是哪个参数还是优化的问题呢？小白只能仰仗各位了  |
回复此楼» 猜你喜欢
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有6人回复
-
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）
已经有4人回复
-
求个博导看看
已经有18人回复
-
自荐读博
已经有6人回复
-
青基代表作，AAAI之类的A会的special track在国内认可度高吗？还是归为workshop之流？
已经有3人回复
-
上海工程技术大学【激光智能制造】课题组招收硕士
已经有6人回复
-
上海工程技术大学张培磊教授团队招收博士生
已经有4人回复
-
临港实验室与上科大联培博士招生1名
已经有9人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用VASP-DFPT+phonopy计算石墨烯声子谱求助
已经有37人回复
- vasp+phonopy计算一个新结构的声子谱有虚频，怎么消掉！！急！！！
已经有4人回复
3楼2016-02-03 02:19:02
人人余yan
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 2590.4
- 散金: 1321
- 红花: 14
- 帖子: 1137
- 在线: 382.3小时
- 虫号: 2877165
- 注册: 2013-12-17
- 专业: 凝聚态物性 II ：电子结构
2楼2016-02-02 22:04:23
人人余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