| 查看: 1693 | 回复: 5 | ||||
swuhaizhilan金虫 (小有名气)
|
[求助]
sigma测试问题 已有2人参与
|
|
请大家帮忙看一下sigma测试问题啊，测试过程参照侯博士的教程，smear取0，得到的结果一直单调变化，应该选哪一个合理啊，标准是啥？侯老师教程里只说了一句“选择entropy T*S EENTRO值中最小的所对应的SIGMA” 0.01 -0.006668 0.02 -0.019231 0.03 -0.033323 0.04 -0.048926 0.05 -0.065022 0.06 -0.081251 0.07 -0.097711 0.08 -0.114861 0.10 -0.152795 0.12 -0.195451 0.14 -0.241477 0.16 -0.289859 0.18 -0.340186 0.20 -0.392413 0.22 -0.446791 0.24 -0.503891 0.26 -0.564442 0.28 -0.629312 0.30 -0.699572 下面是脚本 rm WAVECAR for i in 0.10 0.12 0.14 0.16 0.18 0.20 0.22 0.24 0.26 0.28 0.30 do cat > INCAR <<! SYSTEM = bcc Fe ENCUT = 500 ISTART = 0 ICHARG = 2 ISMEAR = 0 SIGMA = $i VOSKOWN=1 NPAR=4 SPIN=2 LREAL=Auto EDIFFG=-0.02 ! echo " SIGMA = $i eV ";mpirun -n 12 vasp.5.2.12 TS=$(grep "EENTRO" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5}') echo $i $TS >> comment done |
回复此楼» 收录本帖的淘帖专辑推荐
 科研 |
» 猜你喜欢
情人节自我反思：在爱情中有过遗憾吗？
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 琼脂稀释MIC测定时，抗生素的配制问题？
已经有4人回复
- 参数测试能量怎么都不稳定，求救！
已经有10人回复
- 求教，高效液相跑出来的图谱是这样该怎么解决
已经有11人回复
- 国标法测定植酸
已经有6人回复
- 结构优化 siesta 和 vasp
已经有7人回复
- k点和截断能的测试问题
已经有14人回复
- vasp测试SIGMA
已经有12人回复
- 目前测原油蜡含量的主要方法
已经有6人回复
- VASP的ISMEAR=0，测试SIGMA的时候，TS为0是怎么回事？
已经有5人回复
- MTT法检测hegf生长因子活性系列问题求助
已经有22人回复
- ENCUT测试脚本
已经有17人回复
- 考马斯亮蓝法测蛋白质含量
已经有12人回复
- 怎样使用购买的福林酚试剂测定多酚的含量
已经有11人回复
- VASP计算时的几个细节设置问题，希望有经验的人过来指点一下！
已经有23人回复
- 【求助】关于测试SIGMA的问题
已经有9人回复
- 【求助】思考：关于截断能
已经有11人回复
- 【求助】测试sigma取值时，T*S项的值是负数合理吗？
已经有6人回复
- 【求助】求教肝微粒体体外代谢实验问题
已经有20人回复
嗨王and娜娜
木虫 (正式写手)
博士
- 应助: 110 (高中生)
- 金币: 5888.3
- 红花: 11
- 帖子: 786
- 在线: 185.7小时
- 虫号: 2742524
- 注册: 2013-10-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2014-04-03 11:28:38
mazuju028
木虫 (正式写手)
- 1ST强帖: 2
- 应助: 30 (小学生)
- 贵宾: 0.141
- 金币: 2794.4
- 散金: 528
- 红花: 10
- 帖子: 942
- 在线: 419.4小时
- 虫号: 668835
- 注册: 2008-12-05
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
swuhaizhilan: 金币+15, ★★★很有帮助 2014-04-03 16:43:17
感谢参与，应助指数 +1
swuhaizhilan: 金币+15, ★★★很有帮助 2014-04-03 16:43:17
|
根据ＶＡＳＰ官方说明： SIGMA should be as large as possible keeping the difference between the free energy and the total energy (i.e. the term 'entropy T*S') in the OUTCAR file negligible (1 meV/atom). SIGMA的取值规则为：使'entropy T*S数值小于1 meV/atom的最大值。 从你的结果看出：取0.01差不多。 |
3楼2014-04-03 12:19:17
swuhaizhilan
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 64.2
- 散金: 200
- 红花: 1
- 帖子: 140
- 在线: 271.6小时
- 虫号: 2432678
- 注册: 2013-04-23
- 专业: 凝聚态物性 II ：电子结构
4楼2014-04-03 16:38:55
swuhaizhilan
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 64.2
- 散金: 200
- 红花: 1
- 帖子: 140
- 在线: 271.6小时
- 虫号: 2432678
- 注册: 2013-04-23
- 专业: 凝聚态物性 II ：电子结构
5楼2014-04-03 16:42:23
yaoyuan0711
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 207.7
- 散金: 380
- 红花: 2
- 帖子: 197
- 在线: 37.8小时
- 虫号: 4210071
- 注册: 2015-11-10
- 专业: 材料物理化学
6楼2019-01-04 20:23:21
10