| 查看: 4776 | 回复: 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文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有6人回复
-
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有6人回复
-
现代”学阀”该如何界定
已经有14人回复
-
求合成方法
已经有8人回复
-
师弟论文见刊大半年才想起来申请专利，还能抢救一下吗？
已经有4人回复
-
有机合成以后会不会被AI改变？做科研的虫友怎么看
已经有8人回复
-
我的奶奶
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用VASP-DFPT+phonopy计算石墨烯声子谱求助
已经有37人回复
- vasp+phonopy计算一个新结构的声子谱有虚频，怎么消掉！！急！！！
已经有4人回复
|
这周我增加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
人人余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