24小时热门版块排行榜

返回列表

kmw.8668

金虫 (正式写手)

应助: 9 (幼儿园)
金币: 1731.6
帖子: 668
在线: 235.9小时
虫号: 1008177

[交流] H2O的计算

计算一个H2O或者O2分子的能量的时候，遇到点问题，请大家帮帮忙：
1，K点设为1x1x1，改变ENCUT的值去计算。
2，ENCUT的值不变，改变K点去计算。
应该用以上哪种算法呢？如果改变ENCUT，k点取1x1x1计算一个H2O分子的时候则没问题。看一篇文献上说K点是描述原子间的相互作用，计算单个原子的时候可以不变k点，我想对于分子需要改变K点吧？可是如果改变k点，计算H2O分子的结果总是不收敛，而且计算过程中占的内存好像还比较大。

SYSTEM=I
NSW=200
IBRION=2
ISIF=3
EDIFF=1E-8
EDIFFG=-0.005
PREC=Accurate
ENCUT=600
!ADDGRID = .True.
NPAR=1
NSIM=8
LPLANE= .TRUE.
LSCALU= .FALSE.这是我 INCAR文件

回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2012-01-05 11:16:21

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

llh2010

至尊木虫 (著名写手)

应助: 34 (小学生)
金币: 10007.8
帖子: 1614
在线: 350.1小时
虫号: 1471268

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
liliangfang(金币+1): 多谢指教 2012-01-05 16:39:56
kmw.8668(金币+5): 谢谢，把你所说的那份VASP手册给我传一份吧。。。。。邮箱 physics201111@yahoo.cn 2012-01-07 15:14:19

对于原子与分子的计算，VASP手册上明确说了只需要一个K点即可，原因是没有使用布洛赫定理。

赞一下(2人)

回复此楼

2楼2012-01-05 16:38:53

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

zhoudayan2

木虫 (小有名气)

应助: 1 (幼儿园)
金币: 1550.3
帖子: 109
在线: 67小时
虫号: 1098641

★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
sunyang1988(金币+3): 谢谢指点 2012-01-06 23:47:48
kmw.8668(金币+5): 多谢了！十分感激！ 2012-01-07 15:10:52

1: 改变K点算，3 3 3 ，5 5 5，算出来的结果应该跟1 1 1，的不会相差很大，用 1 1 1 够了，不要改变。
2：改变ENCUT算，基本上不会使计算结果改变，ENCUT 的设置看手册，你设置600，足够大了，如果改变ENCUT算有变化，只能说你设置的不够大，是不准确的。
3，你出现问题的原因，应该不在你说的两个问题中任何一个，算O2是会有问题，我算过很多了，跟计算方法，与你所用赝势是很相关的，但是想与实验值吻合是有点困难的，这可能是VASP的缺陷，但是算H2O不应该存在问题，你的EDIFF=1E-8,EDIFFG=-0.005可以设置小一点。

赞一下(2人)

回复此楼

3楼2012-01-06 22:30:26

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