24小时热门版块排行榜    

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

cuiwen0314

银虫 (小有名气)

[求助] VASP 计算CeO2用DFT+U方法的输入如何设置

我在重复一篇文献,是用vasp计算有关CeO2吸附小分子的研究,总的电子和电子密度的计算是使用DFT+U方法,其中Ce的U=4.5eV,方法是GGA=PE,截断能为415eV,我在计算晶面时的输入为
SYSTEM = GeO2_cell
ISTART=0
ICHARG=2
#ISPIN=2
ENCUT=415
PREC=Accurate
LREAL=Auto;
ISMEAR=0;
SIGMA=0.05;
NELMIN=8;
EDIFF=1E-4;
EDIFFG=-0.01;
MAXMIX=80;
GGA=PE;
ISIF=2;      
NSW=200;     
IBRION=2;
LPLANE = .TRUE.
NCORE=4
LSCALU= .FALSE.
NSIM=4
NPAR=8
LDAU=.TRUE.
LDAL=2 2
LDAUU=0 2
LDAUJ=0 0
LMAXMIX=4
LWAVE=.FALSE.
请大家帮忙看看这个输入对不对啊?谢谢啦
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

emilyoyang

木虫 (正式写手)

引用回帖:
4楼: Originally posted by wuli8 at 2013-04-13 21:54:35
LDAU        = .FALSE.不对

Setting LDAU=.TRUE. in the INCAR file switches on the L(S)DA+U....

Yes, U are right. I mad a mistak, because i just copy the setting suggested in the manual without changing the defaut valuse
6楼2013-04-13 22:56:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

emilyoyang

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
cuiwen0314: 金币+10, ★★★很有帮助 2013-04-13 15:31:50
cuiwen0314: 金币+10, ★★★很有帮助 2013-04-13 15:39:20
there are some problems on the DFT+U settings, U may use these setting
LDAU        = .FALSE.
LDAUTYPE        = 2
LDAUL= 3 -1
LDAUU= 4.5 0
LDAUJ= 0 0
LASPH = .TRUE.
LMAXMIX = 6
Moreover, U should double check your NCORE and NPAR  tags in VASP maunal, these two tags have connections with each other, U shall make sure the specific meaning of ur settings in ur INCAR
2楼2013-04-13 12:17:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cuiwen0314

银虫 (小有名气)

引用回帖:
2楼: Originally posted by emilyoyang at 2013-04-13 12:17:17
there are some problems on the DFT+U settings, U may use these setting
LDAU        = .FALSE.
LDAUTYPE        = 2
LDAUL= 3 -1
LDAUU= 4.5 0
LDAUJ= 0 0
LASPH = .TRUE.
LMAXMIX = 6
Moreover, U should double che ...

感谢您的回复!有一点我不明白的是如果U不给定的话,那LDAUL=
LDAUU=
LDAUJ=
LMAXMIX =
这几个量是如何设定的呢,我看说明书
LDAUU=U1 U2   ... specifies the effective on-site Coulomb interaction parameters.

LDAUJ=J1 J2   ... specifies the effective on-site Exchange interaction parameters.

那么U1,U2,J1,J2如何得到?
LMAXMIX to 4 (d elements) and 6 (f elements) 在CeO2体系中规定的是将Ce(4f,5s,5p,5d,6s)作为价电子处理的,那是不是就是看最内层的价电子是d电子还是f电子,而加以选择的?那如果是S电子呢?谢谢啦
3楼2013-04-13 15:30:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuli8

荣誉版主 (知名作家)

…………

优秀版主优秀版主

引用回帖:
2楼: Originally posted by emilyoyang at 2013-04-13 12:17:17
there are some problems on the DFT+U settings, U may use these setting
LDAU        = .FALSE.
LDAUTYPE        = 2
LDAUL= 3 -1
LDAUU= 4.5 0
LDAUJ= 0 0
LASPH = .TRUE.
LMAXMIX = 6
Moreover, U should double che ...

LDAU        = .FALSE.不对

Setting LDAU=.TRUE. in the INCAR file switches on the L(S)DA+U.
…………
4楼2013-04-13 21:54:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见