24小时热门版块排行榜    

查看: 1484  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

kmw.8668

金虫 (正式写手)


[交流] 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文件
回复此楼

» 猜你喜欢

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

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

查看全部散金贴

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

kmw.8668

金虫 (正式写手)


引用回帖:
: Originally posted by llh2010 at 2012-01-05 16:38:53:
对于原子与分子的计算,VASP手册上明确说了只需要一个K点即可,原因是没有使用布洛赫定理。

我只看到计算单原子的,分子的没见过,希望你能传给我一份,多谢了
4楼2012-01-07 15:24:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

llh2010

至尊木虫 (著名写手)


★ ★
小木虫(金币+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楼2012-01-05 16:38:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhoudayan2

木虫 (小有名气)


★ ★ ★ ★
小木虫(金币+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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026国自然放榜时间 +3 布布和一二 2026-08-08 3/150 2026-08-08 18:15 by Lightingo
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 好奇怪的filecode +3 布布和一二 2026-08-08 4/200 2026-08-08 17:46 by zhanghaozhu
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:27 by oEVWOejN9taj
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 17:02 by oEVWOejN9taj
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 16:47 by oEVWOejN9taj
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 16:47 by oEVWOejN9taj
[基金申请] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 16:27 by oEVWOejN9taj
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 15:07 by oEVWOejN9taj
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:47 by oEVWOejN9taj
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:27 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
信息提示
请填处理意见