24小时热门版块排行榜    

查看: 1566  |  回复: 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的回帖

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的回帖

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的回帖
相关版块跳转 我要订阅楼主 kmw.8668 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 02:25 by YXM1FDCSE8Bn
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 02:25 by YXM1FDCSE8Bn
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 02:02 by YXM1FDCSE8Bn
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 00:54 by YXM1FDCSE8Bn
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 6/300 2026-11-01 00:49 by YXM1FDCSE8Bn
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 00:48 by YXM1FDCSE8Bn
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 00:44 by YXM1FDCSE8Bn
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 4/200 2026-11-01 00:41 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 X5bHxOULxFqb 2026-10-30 15/750 2026-11-01 00:22 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 xvJS1BYDxvBS 2026-10-31 7/350 2026-10-31 23:18 by YXM1FDCSE8Bn
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-10-31 23:10 by YXM1FDCSE8Bn
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 X5bHxOULxFqb 2026-10-30 19/950 2026-10-31 22:44 by YXM1FDCSE8Bn
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 X5bHxOULxFqb 2026-10-30 24/1200 2026-10-31 22:39 by YXM1FDCSE8Bn
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 4/200 2026-10-31 21:54 by YXM1FDCSE8Bn
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 9/450 2026-10-31 21:34 by YXM1FDCSE8Bn
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +8 zuesli 2026-10-27 14/700 2026-10-30 22:44 by 好运连连公主
[基金申请] 论文的挂名怎么评判贡献?基金申请中,排名第一的挂名文章 +3 biily 2026-10-25 3/150 2026-10-30 18:33 by 旁观2020
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
[文学芳草园] 闲聊 +3 myrtle 2026-10-25 6/300 2026-10-29 15:28 by myrtle
[考博] 申博求助 20+3 Katrinae 2026-10-26 4/200 2026-10-27 09:56 by 北京莱茵编辑
信息提示
请填处理意见