| 查看: 3077 | 回复: 15 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
cuicui2397新虫 (小有名气)
|
[求助]
vasp阶段能测试求助，谢谢啊！急求。
|
|
|
看的侯博士的教程，里面有测试Si的截断能，但我自己去测试老是出现问题，求解啊！谢谢。 POSCAR文件： Si-Diamond 5.430 0.0 0.5 0.5 0.5 0.0 0.5 0.5 0.5 0.0 2 Direct 0.0 0.0 0.0 0.25 0.25 0.25 KPOINTS文件： Automatic generation 0 Monhkorst-Pack 9 9 9 0.0 0.0 0.0 用来确定ENCUT的脚本程序为run_ecut，其内容是： #!/bin/sh rm -rf WAVECAR for i in 150 200 250 300 350 400 do cat > INCAR <<! SYSTEM = Si-Diamond ENCUT = $i ISTART = 0; ICHARG = 2 ISMEAR = -5 PREC = Accurate ! echo "ENCUT = $i eV"; mpirun -np 2 vasp E='grep "TOTEN" OUTCAR |tail-1|awk '{printf"%12.6f \n",$5}'' echo $i $E >>comment done 同时选择了赝势文件POTCAR，把以上四个文件POSCAR、KPOINTS、run_ecut和POTCAR放入我的计算目录下，在终端输入sh.run_ecut命令，在生成的comment里是这样的： 150 grep "TOTEN" OUTCAR |tail-1|awk {printf%12.6f\n,} 200 grep "TOTEN" OUTCAR |tail-1|awk {printf%12.6f\n,} 250 grep "TOTEN" OUTCAR |tail-1|awk {printf%12.6f\n,} 300 grep "TOTEN" OUTCAR |tail-1|awk {printf%12.6f\n,} 350 grep "TOTEN" OUTCAR |tail-1|awk {printf%12.6f\n,} 400 grep "TOTEN" OUTCAR |tail-1|awk {printf%12.6f\n,} 教程中第二列显示的是能量值，请问这是什么问题啊？谢谢啊！ |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有4人回复
-
跳槽后在研项目怎么办？
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
什么时候开奖？
已经有10人回复
cuicui2397
新虫 (小有名气)
- 应助: 5 (幼儿园)
- 金币: 1611.8
- 红花: 2
- 帖子: 211
- 在线: 130.4小时
- 虫号: 2198831
- 注册: 2012-12-21
- 专业: 凝聚态物性 II ：电子结构
11楼2013-12-02 11:16:46
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
cuicui2397(fzx2008代发): 金币+3 2013-11-14 21:25:06
sunyang1988: 金币+1, 谢谢交流 2013-11-15 18:40:10
cuicui2397: 金币+5, ★★★★★最佳答案 2013-11-23 10:10:40
感谢参与，应助指数 +1
cuicui2397(fzx2008代发): 金币+3 2013-11-14 21:25:06
sunyang1988: 金币+1, 谢谢交流 2013-11-15 18:40:10
cuicui2397: 金币+5, ★★★★★最佳答案 2013-11-23 10:10:40
|
那个，有可能是符号有问题，那个E='grep "TOTEN" OUTCAR |tail-1|awk '{printf"%12.6f \n",$5}'' 是`（ESC下面那个键，数字1左面那个键），而不是'(单引号)，你试试看 |
» 本帖已获得的红花（最新10朵）
2楼2013-11-14 20:30:12
cuicui2397
新虫 (小有名气)
- 应助: 5 (幼儿园)
- 金币: 1611.8
- 红花: 2
- 帖子: 211
- 在线: 130.4小时
- 虫号: 2198831
- 注册: 2012-12-21
- 专业: 凝聚态物性 II ：电子结构
送红花一朵 |
你好，首先谢谢你，我把符号改过之后，第二列出现了能量值。不过和教程上差别很大，你觉得是什么原因造成的啊？我的结果是这样的： 150 -10.744031 200 -10.795431 250 -10.822115 300 -10.836639 350 -10.840896 400 -10.842717 而教程中是这样写的： 150 -11.900655 200 -11.938864 250 -11.944599 300 -11.945248 350 -11.945503 400 -11.945622 总能变化在0.001eV左右就足够了，在这个例子中，因此，选择ENCUT=250eV。 我这样算出的没什么参考价值啊，求解，谢谢！ |
3楼2013-11-14 21:18:24
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
sunyang1988: 金币+1, 谢谢交流 2013-11-15 18:40:17
cuicui2397: 金币+12, ★★★很有帮助 2013-11-23 10:11:24
sunyang1988: 金币+1, 谢谢交流 2013-11-15 18:40:17
cuicui2397: 金币+12, ★★★很有帮助 2013-11-23 10:11:24
| 赝势是不是不一样，你有没有选正确的赝势？侯老师的那个教程，有没有说选什么赝势？ |
» 本帖已获得的红花（最新10朵）
4楼2013-11-14 23:06:04
