| 查看: 1455 | 回复: 5 | ||
jk射手座新虫 (著名写手)
|
[求助]
超晶格结构原子弛豫优化出问题 已有1人参与
|
|
今天突然碰上这么个问题，优化2层、4层结构的钛酸铅和铌酸钾的超晶格没问题，但是3:3 4:4 5:5 的结构就出现优化后的能量 的值成千上万、有正有负。。显然是出错误的，已经找过很多原因，; 现在还是在不停的改参数测试中，之余来这里求教一下有没有遇到过类似情况的同僚；nohup中有提示，虽然就是那个老生常谈的 Segmentation Fault ，但是我不理解为什么我n：n （n=1~5）的钛酸铅和铝酸镧结构都没任何问题，怎么铌酸钾的这么尴尬呢？ 附带说一下，用的都是PAW赝势 INCAR: SYSTEM=PTOKNO ENCUT=450 PREC=Normal ISMEAR=0 SIGMA=0.05 IBRION=7 ISTART=0 ICHARG=2 EDIFF=1E-3 ISIF=2 NSW=200 KPOINTS: Auto-generate 0 Monkhorst-Pack 6 6 2 0 0 0 脚本： #!/bin/bash #for ps in 0.5 10 20 30 50 100 150 200 250 300 350 400 450 500 550 600; do for c in 23.5026 23.5224 23.74 23.76; do cat > POSCAR <<! data_PTOKNO 1.0 3.89 0.00 0.00 0.00 3.89 0.00 0.00 0.00 $c 3 3 3 3 18 Selective dynamics Direct 0.000 0.000 0.000 F F T 0.000 0.000 0.167 F F T 0.000 0.000 0.333 F F T 0.500 0.500 0.083 F F T 0.500 0.500 0.250 F F T 0.500 0.500 0.417 F F T 0.000 0.000 0.500 F F T 0.000 0.000 0.667 F F T 0.000 0.000 0.833 F F T 0.500 0.500 0.583 F F T 0.500 0.500 0.750 F F T 0.500 0.500 0.917 F F T 0.500 0.000 0.083 F F T 0.000 0.500 0.083 F F T 0.500 0.500 0.000 F F T 0.500 0.000 0.250 F F T 0.000 0.500 0.250 F F T 0.500 0.500 0.167 F F T 0.500 0.000 0.417 F F T 0.000 0.500 0.417 F F T 0.500 0.500 0.333 F F T 0.500 0.000 0.583 F F T 0.000 0.500 0.583 F F T 0.500 0.500 0.500 F F T 0.500 0.000 0.750 F F T 0.000 0.500 0.750 F F T 0.500 0.500 0.667 F F T 0.500 0.000 0.917 F F T 0.000 0.500 0.917 F F T 0.500 0.500 0.833 F F T ! /public/jzh/vasp5 cp CONTCAR CONTCAR$c cp OUTCAR OUTCAR$c grep "TOTEN" OUTCAR | tail -1 >>E.dat done 麻烦懂的帮看一下，我自己现在也在找解决方法，共同讨论  |
回复此楼» 猜你喜欢
一志愿211，化学学硕，310分，本科重点双非，求调剂
已经有13人回复
-
285求调剂
已经有8人回复
-
22408 266求调剂
已经有6人回复
-
264求调剂
已经有8人回复
-
材料工程专业日语生求调剂
已经有9人回复
-
297求调剂
已经有15人回复
-
277、学硕，求调剂 数一104，
已经有12人回复
-
机械专硕273请求调剂
已经有3人回复
-
080100力学316求调剂
已经有3人回复
-
288环境专硕,求调材料方向
已经有39人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 掺杂表面的分子吸附模型建立
已经有8人回复
- 大家知道什么是亚晶格结构嘛？求交流
已经有3人回复
- 如何做vasp掺杂
已经有15人回复
jk射手座
新虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 9882.6
- 散金: 1955
- 红花: 2
- 帖子: 2281
- 在线: 175.2小时
- 虫号: 2777009
- 注册: 2013-11-04
- 性别: GG
- 专业: 零件成形制造
2楼2014-09-02 16:12:20
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
fzx2008: 金币+3, 谢谢指导 2014-09-02 19:16:22
jk射手座: 金币+5, ★有帮助 2014-09-04 09:33:52
感谢参与，应助指数 +1
fzx2008: 金币+3, 谢谢指导 2014-09-02 19:16:22
jk射手座: 金币+5, ★有帮助 2014-09-04 09:33:52
|
為什麼你的 IBRION = 7 ?? 根據手冊上指出 IBRION=7 and IBRION=8 is only supported starting from VASP.5.1. It determines the Hessian matrix (matrix of second derivatives) using density functional perturbation theory. 你要算 Hessian matrix ? 可是看你的腳本 你似乎是在確認 z 方向的長度，不是嬤? |
3楼2014-09-02 17:38:52
jk射手座
新虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 9882.6
- 散金: 1955
- 红花: 2
- 帖子: 2281
- 在线: 175.2小时
- 虫号: 2777009
- 注册: 2013-11-04
- 性别: GG
- 专业: 零件成形制造
4楼2014-09-03 08:40:47
5楼2014-09-03 13:41:52
jk射手座
新虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 9882.6
- 散金: 1955
- 红花: 2
- 帖子: 2281
- 在线: 175.2小时
- 虫号: 2777009
- 注册: 2013-11-04
- 性别: GG
- 专业: 零件成形制造
6楼2014-09-03 14:59:22
