| 查看: 1547 | 回复: 3 | |||
[交流]
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文件 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
申博发邮件
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2012-01-05 16:38:53
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
sunyang1988(金币+3): 谢谢指点 2012-01-06 23:47:48
kmw.8668(金币+5): 多谢了！十分感激！ 2012-01-07 15:10:52
小木虫(金币+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可以设置小一点。 |
3楼2012-01-06 22:30:26
4楼2012-01-07 15:24:43
