| 查看: 1492 | 回复: 14 | |||
didadida321金虫 (正式写手)
|
[求助]
encut测试时，为什么得到的comment中没有对应的能量？
|
| 在对encut测试时，为什么得到的comment中没有对应的能量？刚学习vasp，望各位帮帮忙，谢了。 |
回复此楼» 猜你喜欢
» 本主题相关价值贴推荐，对您同样有帮助:
未知胶乳的测试表征问题
已经有14人回复
- DSC测试（差示扫描量热法）的来看看
已经有37人回复
- 多普勒过程有没有能量的传递，做功的过程？
已经有18人回复
- 关于收敛性测试问题
已经有6人回复
- 保水剂扫描电镜和透射电镜的测试
已经有7人回复
- 求问为什么使用DMOL3模块优化石墨烯结构的时候没有能量曲线？
已经有13人回复
- 求问APPLIED ENVIRONMENT MICROBIOLOGY的comment letter投稿周期？
已经有4人回复
- 物质和能量(原创)
已经有14人回复
- 三重串联四级杆的碰撞能量优化
已经有12人回复
- 单原子在金红石结构中的能量
已经有4人回复
- 老调重谈，CASTEP计算能量上升全攻略
已经有3人回复
- 高压下氢气的总能量该怎么算
已经有5人回复
liliangfang
荣誉版主 (著名写手)
- 1ST强帖: 7
- 应助: 138 (高中生)
- 贵宾: 0.952
- 金币: 10648.7
- 散金: 4575
- 红花: 42
- 沙发: 3
- 帖子: 1856
- 在线: 825.2小时
- 虫号: 1275010
- 注册: 2011-04-23
- 性别: GG
- 专业: 微/纳机械系统
- 管辖: 第一性原理
2楼2012-05-04 09:00:39
liliangfang
荣誉版主 (著名写手)
- 1ST强帖: 7
- 应助: 138 (高中生)
- 贵宾: 0.952
- 金币: 10648.7
- 散金: 4575
- 红花: 42
- 沙发: 3
- 帖子: 1856
- 在线: 825.2小时
- 虫号: 1275010
- 注册: 2011-04-23
- 性别: GG
- 专业: 微/纳机械系统
- 管辖: 第一性原理
【答案】应助回帖
★
sunyang1988: 金币+1, 谢谢帮助 2012-05-04 12:08:10
sunyang1988: 金币+1, 谢谢帮助 2012-05-04 12:08:10
|
#!/bin/sh rm WAVECAR for i in 250 300 350 400 450 500 550 600 do cat > INCAR < System = Si ENCUT = $i ISTART = 0 ICHARG = 2 ISMEAR = 0 SIGMA = 0.1 PREC = Accurate ! echo "ENCUT = $i eV" ; mpdboot ; mpirun -np 4 vasp E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f\n",$5}'` KP=`grep "irreducible" OUTCAR | tail -1 | awk '{printf "%5i \n", $2 }'` echo $i $KP $E >> comment done 这个你最好不是用复制的，自己打进去，这个mpdboot ; mpirun -np 4 vasp 根据你需要修改 |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : encut
2012-05-04 09:05:15, 391 bytes
3楼2012-05-04 09:05:16
didadida321
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 944.7
- 散金: 30
- 红花: 1
- 帖子: 867
- 在线: 166小时
- 虫号: 1702537
- 注册: 2012-03-19
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2012-05-04 09:48:55
techfan
铁虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 31.4
- 帖子: 165
- 在线: 61.5小时
- 虫号: 779793
- 注册: 2009-05-26
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2012-05-04 10:15:02
liliangfang
荣誉版主 (著名写手)
- 1ST强帖: 7
- 应助: 138 (高中生)
- 贵宾: 0.952
- 金币: 10648.7
- 散金: 4575
- 红花: 42
- 沙发: 3
- 帖子: 1856
- 在线: 825.2小时
- 虫号: 1275010
- 注册: 2011-04-23
- 性别: GG
- 专业: 微/纳机械系统
- 管辖: 第一性原理
【答案】应助回帖
★ ★
sunyang1988: 金币+2, 谢谢帮助 2012-05-04 12:08:00
sunyang1988: 金币+2, 谢谢帮助 2012-05-04 12:08:00
|
这是SIGMA，KPOINTS，晶格常数的脚本，你自己修改一下就可以了 |
6楼2012-05-04 10:37:11
techfan
铁虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 31.4
- 帖子: 165
- 在线: 61.5小时
- 虫号: 779793
- 注册: 2009-05-26
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
7楼2012-05-04 12:42:02
techfan
铁虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 31.4
- 帖子: 165
- 在线: 61.5小时
- 虫号: 779793
- 注册: 2009-05-26
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
8楼2012-05-04 13:49:06
liliangfang
荣誉版主 (著名写手)
- 1ST强帖: 7
- 应助: 138 (高中生)
- 贵宾: 0.952
- 金币: 10648.7
- 散金: 4575
- 红花: 42
- 沙发: 3
- 帖子: 1856
- 在线: 825.2小时
- 虫号: 1275010
- 注册: 2011-04-23
- 性别: GG
- 专业: 微/纳机械系统
- 管辖: 第一性原理
【答案】应助回帖
★
sunyang1988: 金币+1, 谢谢交流 2012-05-04 18:31:34
sunyang1988: 金币+1, 谢谢交流 2012-05-04 18:31:34
|
#!/bin/sh rm WAVECAR for i in 0.01 0.02 0.04 0.05 0.06 0.08 0.10 0.12 0.14 0.16 0.18 0.20 do cat > INCAR < System = Si ENCUT = 350 ISTART = 0 ICHARG = 2 ISMEAR = 0 SIGMA = $i PREC = Accurate! echo "SIGMA = $i eV" ; vasp TS=`grep "EENTRO" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'` KP=`grep "irreducible" OUTCAR | tail -1 | awk '{printf "%5i \n", $2 }'` echo $i $TS >>comment done 红色的地方需要你自己修改，如果不行那就不是脚本的问题啦 |
9楼2012-05-04 15:42:25
techfan
铁虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 31.4
- 帖子: 165
- 在线: 61.5小时
- 虫号: 779793
- 注册: 2009-05-26
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
10楼2012-05-04 15:51:44
