24小时热门版块排行榜

>论坛更新日志 (1544)
>论文投稿 (35)
>硕博家园 (31)
>考研 (31)
>论文道贺祈福 (28)
>找工作 (27)
>导师招生 (24)
>虫友互识 (20)
>公派出国 (20)
>基金申请 (19)
>考博 (19)
>教师之家 (16)
>博后之家 (13)
>文献求助 (13)
>休闲灌水 (6)
>海外博后 (5)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

Gina88

木虫 (正式写手)

1ST强帖: 1
应助: 9 (幼儿园)
金币: 3328.8
散金: 124
红花: 8
帖子: 395
在线: 113小时
虫号: 1047684
注册: 2010-06-26
专业: 理论和计算化学

[求助] 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。
所以对着这个参数测试的精度犯糊涂了，
不知道诸位在具体计算的时候都是怎么测试的，还望不吝赐教！
先谢谢各位大侠了！

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

用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人回复

1楼 2012-11-04 10:18:29

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

Gina88

木虫 (正式写手)

1ST强帖: 1
应助: 9 (幼儿园)
金币: 3328.8
散金: 124
红花: 8
帖子: 395
在线: 113小时
虫号: 1047684
注册: 2010-06-26
专业: 理论和计算化学

引用回帖:

9楼: Originally posted by wzbhit at 2012-11-05 11:44:51
变化趋势比较诡异，看不出来原因，把INCAR贴出来看看...

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.
!