应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » vasp阶段能测试求助，谢谢啊！急求。
161/2返回列表12下一页
查看: 2658  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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,}
教程中第二列显示的是能量值，请问这是什么问题啊？谢谢啊！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-11-14 20:16:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

甲斐之虎

银虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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左面那个键），而不是'(单引号)，你试试看
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

cuicui2397
别人笑我太疯颠，我笑他人看不穿。不见五陵豪杰墓，无花无酒锄作田。
2楼2013-11-14 20:30:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

甲斐之虎

银虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
sunyang1988: 金币+1, 谢谢交流 2013-11-15 18:40:17
cuicui2397: 金币+12, ★★★很有帮助 2013-11-23 10:11:24
赝势是不是不一样，你有没有选正确的赝势？侯老师的那个教程，有没有说选什么赝势？
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

cuicui2397
别人笑我太疯颠，我笑他人看不穿。不见五陵豪杰墓，无花无酒锄作田。
4楼2013-11-14 23:06:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cuicui2397

新虫 (小有名气)
引用回帖:
7楼: Originally posted by 甲斐之虎 at 2013-11-15 12:48:07
谢谢你的花，但是请别送花了，浪费金币，留着以后提问题吧，有的用呢！
我说的不是这个意思，而是你采用了何种方法？PAW，UUSP？何种泛函91，PBE，LDA？这些不一样，得出结果也不一样...

嗯，呵呵，我自己选用的是paw_pbe赝势，而侯的教程选用的是uspp、LDA。之前没看清楚，给过来之后，结果和教程中一致了。谢谢你啊！不过我还想请教一下，针对具体体系，选择赝势是paw或是uspp，交换关联函数选择是LDA或是GGA，有什么根据吗？谢谢！
赞一下(1人) 回复此楼
8楼2013-11-15 16:03:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

微6凉

铁虫 (小有名气)
请问，我按照您这样的做法，总会出现这样一些结果Try `rm --help' for more information.
'un_ecut: line 4: syntax error near unexpected token `do
'un_ecut: line 4: `do  这是怎么回事
赞一下(1人) 回复此楼
14楼2015-08-03 11:50:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

cuicui2397

新虫 (小有名气)
送红花一朵
引用回帖:
2楼: Originally posted by 甲斐之虎 at 2013-11-14 20:30:12
那个，有可能是符号有问题，那个E='grep "TOTEN" OUTCAR |tail-1|awk '{printf"%12.6f \n",$5}''
是`（ESC下面那个键，数字1左面那个键），而不是'(单引号)，你试试看

你好，首先谢谢你，我把符号改过之后，第二列出现了能量值。不过和教程上差别很大，你觉得是什么原因造成的啊？我的结果是这样的：
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cuicui2397

新虫 (小有名气)
送红花一朵
引用回帖:
4楼: Originally posted by 甲斐之虎 at 2013-11-14 23:06:04
赝势是不是不一样，你有没有选正确的赝势？侯老师的那个教程，有没有说选什么赝势？

针对这个例子来说，侯的教程里没有说具体选什么赝势，而Si有两个赝势，分别是Si  245和Si_h 380。之前侯的教程里说一般优先选择Si 245，所以我选择了245那个，不知道对不？那我再试试380那个赝势。
赞一下 回复此楼
5楼2013-11-15 09:54:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cuicui2397

新虫 (小有名气)
引用回帖:
4楼: Originally posted by 甲斐之虎 at 2013-11-14 23:06:04
赝势是不是不一样，你有没有选正确的赝势？侯老师的那个教程，有没有说选什么赝势？

改成Si 380的赝势，貌似更不对了，能量值是这样的：
150  -11.407850
200  -10.938620
250  -10.825592
300  -10.813574
350  -10.830185
400  -10.849576
请问你一般是根据什么选择赝势啊？谢谢！
赞一下 回复此楼
6楼2013-11-15 10:29:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

甲斐之虎

银虫 (正式写手)
★
sunyang1988: 金币+1, 谢谢交流 2013-11-15 18:40:33
引用回帖:
6楼: Originally posted by cuicui2397 at 2013-11-15 10:29:02
改成Si 380的赝势，貌似更不对了，能量值是这样的：
150  -11.407850
200  -10.938620
250  -10.825592
300  -10.813574
350  -10.830185
400  -10.849576
请问你一般是根据什么选择赝势啊？谢谢！...

谢谢你的花，但是请别送花了，浪费金币，留着以后提问题吧，有的用呢！
我说的不是这个意思，而是你采用了何种方法？PAW，UUSP？何种泛函91，PBE，LDA？这些不一样，得出结果也不一样
赞一下 回复此楼
别人笑我太疯颠，我笑他人看不穿。不见五陵豪杰墓，无花无酒锄作田。
7楼2013-11-15 12:48:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

忠义二郎

铜虫 (初入文坛)
引用回帖:
2楼: Originally posted by 甲斐之虎 at 2013-11-14 20:30:12
那个，有可能是符号有问题，那个E='grep "TOTEN" OUTCAR |tail-1|awk '{printf"%12.6f \n",$5}''
是`（ESC下面那个键，数字1左面那个键），而不是'(单引号)，你试试看

请教一个问题，同样是侯博士的这个测试Si的一个例子，我没有UUSP的赝势文件，所有就选了LDA、PBE的几个Si的赝势文件试一下，可是不论用哪个POTCAR，在comment中总是
150   0.0000
200   0.0000

请问是个什么情况？
赞一下 回复此楼
http://sdustzk.lofter.com
9楼2013-12-02 10:30:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

忠义二郎

铜虫 (初入文坛)
请问这个例子中你的INCAR是怎么写的？
赞一下 回复此楼
http://sdustzk.lofter.com
10楼2013-12-02 10:34:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cuicui2397 的主题更新
161/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 焦虑 +9 水冰月月野兔 2026-03-13 13/650 2026-03-19 09:50 by otani
[考研] 271材料工程求调剂 +5 .6lL 2026-03-18 5/250 2026-03-19 03:07 by 无懈可击111
[考研] 一志愿中国海洋大学，生物学，301分，求调剂 +4 1孙悟空 2026-03-17 4/200 2026-03-18 17:59 by fivewind
[考研] 0854可跨调剂，一作一项核心论文五项专利，省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 生物学071000 329分求调剂 +3 我爱生物生物爱� 2026-03-17 3/150 2026-03-18 10:12 by macy2011
[考研] 301求调剂 +9 yy要上岸呀 2026-03-17 9/450 2026-03-18 08:58 by 无际的草原
[考研] 085600材料与化工求调剂 +5 绪幸与子 2026-03-17 5/250 2026-03-17 16:40 by laoshidan
[考研] 26考研求调剂 +6 丶宏Sir 2026-03-13 6/300 2026-03-17 16:13 by 醉在风里
[考研] 一志愿南京大学，080500材料科学与工程，调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 283求调剂 +3 听风就是雨； 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 326求调剂 +4 诺贝尔化学奖觊� 2026-03-15 7/350 2026-03-16 17:11 by 诺贝尔化学奖觊�
[基金申请] 今年的国基金是打分制吗？ 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 0703一志愿211 285分求调剂 +5 ly3471z 2026-03-13 5/250 2026-03-16 16:16 by 哦哦123
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 材料与化工 323 英一+数二+物化，一志愿：哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
信息提示
关闭
请填处理意见
关闭