| 查看: 1947 | 回复: 11 | ||
Gina88木虫 (正式写手)
|
[求助]
ISMEAR=1的测试问题
|
|
|
一个体系有70多个原子， 我计算的时候ISMEAR=1，然后测试sigma，当sigma=0.2的时候，得到TS为-0.036318 但是几何优化的时候sigma=0.2时，得到TS=-0.089262，然后采用sigma=0.15，得到TS=-0.063869 那我sigma应该选0.2呢还是应该选0.15呢？ 大家有遇到静态测试盒几何优化时，ISMEAR和sigma取相同的值，却得到不同数值TS的情况吗？ 还有大家测试ENCUT和Kpoints的时候，最终精度是不是都达到了delta-total energy小于1meV呢？我觉得我自己测试的时候好像很难到达1meV的精度啊，即使到达1meV的精度，那ENCUT和KPOINTs取值也比较大，因而计算量很大。但是手册上及侯老师的确这么说的。 然后查看文献的时候，发现有的人用的delta-total energy的精度是10meV。 所以对着这个参数测试的精度犯糊涂了， 不知道诸位在具体计算的时候都是怎么测试的，还望不吝赐教！ 先谢谢各位大侠了！ |
回复此楼» 猜你喜欢
276求调剂。有半年电池和半年高分子实习经历
已经有4人回复
-
284求调剂
已经有6人回复
-
306求调剂
已经有7人回复
-
291求调剂
已经有5人回复
-
070300化学求调剂
已经有4人回复
-
070300，一志愿北航320求调剂
已经有5人回复
-
085600材料与化工
已经有9人回复
-
0703 调剂
已经有6人回复
-
招08考数学
已经有11人回复
-
276求调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用VASP-DFPT+phonopy计算石墨烯声子谱求助
已经有37人回复
- HSE计算的一些问题，希望大家帮我解决谢谢
已经有13人回复
- 关于Si算例结果分析
已经有5人回复
- 超晶格结构原子弛豫优化出问题
已经有5人回复
- 求助一下如何超算中心提交vasp任务
已经有8人回复
- 求大神帮忙解释截断能测试中的问题！
已经有7人回复
- sigma测试问题
已经有5人回复
- NEB过渡态问题，搜索的结果为什么和插入结果差距这么大？
已经有6人回复
- 关于VASP中对ENCUT收敛性测试的疑问
已经有5人回复
- vasp阶段能测试求助，谢谢啊！急求。
已经有15人回复
- NEB的计算出现负值，测试的其他能量都比初始的00和终止的10低
已经有17人回复
- VASP计算ZnO的性质，为什么ENCUT的测试曲线是这个样子的？
已经有10人回复
- 【重新开帖】KPOINTS的测试脚本问题
已经有2人回复
- ENCUT测试
已经有12人回复
- 求助：带电缺陷HSE06总能计算结果匪夷所思，请赐教！
已经有11人回复
- [求助成功]VASP 计算 graphene 的错误
已经有15人回复
- 紧急求助vasp5.2 计算氧化铈表面吸附遇到的若干问题
已经有4人回复
- 利用GGA+U计算ZnO带隙的问题
已经有3人回复
- vasp4.6问题
已经有2人回复
- 为什么我进行ENCUT测试时是这个样子的？问题出在哪里？
已经有8人回复
- 【求助】关于测试SIGMA的问题
已经有9人回复
wzbhit
至尊木虫 (著名写手)
- 应助: 84 (初中生)
- 金币: 11959.7
- 散金: 26
- 红花: 36
- 帖子: 1162
- 在线: 256.7小时
- 虫号: 1128043
- 注册: 2010-10-21
- 专业: 理论和计算化学
3楼2012-11-04 11:54:45
colinchen2
金虫 (小有名气)
- 应助: 15 (小学生)
- 金币: 1652
- 帖子: 59
- 在线: 33.7小时
- 虫号: 1536228
- 注册: 2011-12-13
- 性别: GG
- 专业: 半导体物理
2楼2012-11-04 10:21:38
Gina88
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 9 (幼儿园)
- 金币: 3328.8
- 散金: 124
- 红花: 8
- 帖子: 395
- 在线: 113小时
- 虫号: 1047684
- 注册: 2010-06-26
- 专业: 理论和计算化学
4楼2012-11-04 23:33:41
Gina88
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 9 (幼儿园)
- 金币: 3328.8
- 散金: 124
- 红花: 8
- 帖子: 395
- 在线: 113小时
- 虫号: 1047684
- 注册: 2010-06-26
- 专业: 理论和计算化学
5楼2012-11-04 23:43:04
colinchen2
金虫 (小有名气)
- 应助: 15 (小学生)
- 金币: 1652
- 帖子: 59
- 在线: 33.7小时
- 虫号: 1536228
- 注册: 2011-12-13
- 性别: GG
- 专业: 半导体物理
6楼2012-11-05 07:26:43
wzbhit
至尊木虫 (著名写手)
- 应助: 84 (初中生)
- 金币: 11959.7
- 散金: 26
- 红花: 36
- 帖子: 1162
- 在线: 256.7小时
- 虫号: 1128043
- 注册: 2010-10-21
- 专业: 理论和计算化学
7楼2012-11-05 09:20:31
Gina88
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 9 (幼儿园)
- 金币: 3328.8
- 散金: 124
- 红花: 8
- 帖子: 395
- 在线: 113小时
- 虫号: 1047684
- 注册: 2010-06-26
- 专业: 理论和计算化学
8楼2012-11-05 09:48:31
wzbhit
至尊木虫 (著名写手)
- 应助: 84 (初中生)
- 金币: 11959.7
- 散金: 26
- 红花: 36
- 帖子: 1162
- 在线: 256.7小时
- 虫号: 1128043
- 注册: 2010-10-21
- 专业: 理论和计算化学
9楼2012-11-05 11:44:51
Gina88
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 9 (幼儿园)
- 金币: 3328.8
- 散金: 124
- 红花: 8
- 帖子: 395
- 在线: 113小时
- 虫号: 1047684
- 注册: 2010-06-26
- 专业: 理论和计算化学
|
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=Si111_relax ISTART=0;ICHARG=2 PREC=Accurate ENCUT=300 GGA=91;VOSKOWN=1 ISMEAR=1;SIGMA=$i EDIFF=1E-6 LREAL=A ALGO=Veryfast;IALGO=48 NGX=170 NGY=170 NGZ=80 NBANDS=200 !ISIF=4 !ISYM=2 !NSW=200 !IBRION=2 !POTIM=0.5 LWAVE=.FALSE. LCHARG=.FALSE. ! |
10楼2012-11-05 21:34:45
40