VASP5.0以上版本有VDW校正,使用PBE赝势,在INCAR文件中加入LVDW=.T.
我简单设置后计算得到的层间距跟实验或其他文献结果差挺多的,为3.2A左右,实验为3.354A,其他文献结果与此接近,为3.34…。您说的“最要注意以的是与交换关联相关的那个关键词的设置”是指什么?烦请老师具体说一下是什么关键词,如何设置。谢谢!
我的参数都用的手册上默认的:
LVDW=.TRUE.
VDW_RADIUS=30.0
VDW_SCALING=0.75
VDW_D=20.0
VDWC6=1.75
VDW_R0=1.452
其他的参数是否做了检查,ENCUT, K-mesh, 是不是用的GGA=PE?
采用ms导出来的模型,包含4个C原子
ENCUT=520eV
Kpoints采用21*21*7
GGA=PE
不知道我的问题出在哪?您那边有做过考虑vdw的计算吗?大概是什么情况,方便给介绍一下吗?谢谢!
虽然是按c/a的比率按scale的,但k点网格沿着c方向的分割数显得有点过小。
引致新浪"有志者"的一片博文:在VASP5.2.11以上版本中可以看到有计算范德华力的INCAR关键词设置项,具体说明可以在http://cms.mpi.univie.ac.at/vasp/vasp/vasp.html上面查到。
在阅读文献时,也经常注意到因为采用广义梯度近似的密度泛函理论计算没有考虑到范德华作用,如果体系中范德华力成分显著并影响较大的话,忽略范德华力的计算的精准度是需要质疑的。
因此,最近关注了一下范德华力的计算,也注意到版内已经有些帖子在讨论这个问题了,比如http://emuch.net/html/201103/2885194.html和http://emuch.net/html/201111/3766493.html。
为了方便虫友一起学习提高,将范德华力的理论计算方法和应用的三篇文档分享给大家,见附件,前两篇的作者都是出自一个研究组,主要是方法基础;第三篇来自另外一个研究组,计算了金属表面石墨烯的性质,属于方法应用。
再根据文献提供一个列表,都是属于范德华力泛函的应用计算文章。
CODE:[Copyto clipboard]本人还没做过考虑范德华力的体系计算,主要纯属抛砖引玉,引起大家对范德华力计算的重视,互相交流这方面的学习心得。
A. Puzder, M. Dion, D.C. Langreth
J. Chem. Phys., 124 (2006), p. 164105
T. Thonhauser, A. Puzder, D.C. Langreth
J. Chem. Phys., 124 (2006), p. 164106
S.D. Chakarova, E. Schroder
J. Chem. Phys., 122 (2005), p. 054102
E. Ziambaras, J. Kleis, E. Schröuml, P.Hyldgaard
Phys. Rev. B, 76 (2007), p. 155425
J. Kleis, B.I. Lundqvist, D.C. Langreth, E. Schroder
Phys. Rev. B, 76 (2007), p. 100201
V.R. Cooper, T. Thonhauser, A. Puzder, E. Schroder, B.I. Lundqvist,D.C. Langreth
J. Am. Chem. Soc., 130 (2008), p. 1304
P.G. Moses, J.J. Mortensen, B.I. Lundqvist, J.K. Norskov
J. Chem. Phys., 130 (2009), p. 104709
L. Romaner, D. Nabok, P. Puschnig, E. Zojer, C.Ambrosch-Draxl
New J. Phys., 11 (2009), p. 053010
S.D. Chakarova-Kack, E. Schroder, B.I. Lundqvist, D.C.Langreth
Phys. Rev. Lett., 96 (2006), p. 146107
D.C. Langreth, B.I. Lundqvist, S.D. Chakarova-Kack, V.R. Cooper, M.Dion, P. Hyldgaard, A. Kelkkanen, J. Kleis, L.Z. Kong, S. Li, P.G.Moses, E. Murray, A. Puzder, H. Rydberg, E. Schroder, T.Thonhauser
J. Phys.: Condens. Matter, 21 (2009), p. 084203
欢迎讨论。
各位朋友,有用过vasp5.2算范德华力的吗?我看了下说明书,给了几个参数的设置,并且要用PBE势,那么,是不是在优化时,就要加入这些参数呢?还是怎样?并且,说明书里,还提到了在优化时,ISIF>=3的都不可以用的。大伙都是怎么处理的啊?踊跃讨论下啊。:
既然manual都说得这么清楚了,那就按照manual说的做吧。这些参数应该在优化的时候就加入。
不过如果涉及晶格常数在优化的时候有变化啥的,那估计还得用ISIF=3。
wang-lf09(站内联系TA)1楼: Originally posted by wanglei200 at 2011-10-312126:
各位朋友,有用过vasp5.2算范德华力的吗?我看了下说明书,给了几个参数的设置,并且要用PBE势,那么,是不是在优化时,就要加入这些参数呢?还是怎样?并且,说明书里,还提到了在优化时,ISIF>=3的都不可以用...你做什么体系?要用VDW-DF泛函。vasp里面还有DFT-D的方法考虑范德华力,不知道你们有没有考虑过?
wanglei200(站内联系TA)3楼: Originally posted by wang-lf09 at 2011-11-011133:
你做什么体系?要用VDW-DF泛函。vasp里面还有DFT-D的方法考虑范德华力,不知道你们有没有考虑过?DFT-D的方法?可否告知一二。感谢这位朋友。:)
buct2010(站内联系TA)1楼: Originally posted by wanglei200 at 2011-10-312126:
各位朋友,有用过vasp5.2算范德华力的吗?我看了下说明书,给了几个参数的设置,并且要用PBE势,那么,是不是在优化时,就要加入这些参数呢?还是怎样?并且,说明书里,还提到了在优化时,ISIF>=3的都不可以用...这个范德华力是在吸附的时候用的到,其他情况,比如晶格常数优化是没用的。因为有些气体吸附在体系表面时(弱相互作用),算到的吸附能与实验值偏差大,原因是没有考虑长程的范德华作用力。最近在看这方面 的文献,希望对你有帮助:THEJOURNAL OF CHEMICAL PHYSICS 134, 104707 (2011)
wang-lf09(站内联系TA)4楼: Originally posted by wanglei200 at 2011-11-011344:
DFT-D的方法?可否告知一二。感谢这位朋友。:)vasp的手册说明上不就有吗?应该是叫DFT-D2。你看了VDW-DF,应该也看到这个了呀,是挨着的。