应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP 计算CeO2用DFT+U方法的输入如何设置
101/1返回列表
查看: 3891  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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.
请大家帮忙看看这个输入对不对啊？谢谢啦
回复此楼

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-04-13 11:45:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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

cuiwen0314

银虫 (小有名气)
引用回帖:
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....

谢谢您的解释！
那您能给我解答下面这个问题吗？
感谢您的回复！有一点我不明白的是如果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电子呢？谢谢啦
赞一下 回复此楼
5楼2013-04-13 22:08:33
已阅   回复此楼   关注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的回帖

emilyoyang

木虫 (正式写手)
引用回帖:
3楼: Originally posted by cuiwen0314 at 2013-04-13 15:30:57
感谢您的回复！有一点我不明白的是如果U不给定的话，那LDAUL=
LDAUU=
LDAUJ=
LMAXMIX =
这几个量是如何设定的呢，我看说明书
LDAUU=U1 U2   ... specifies the effective on-site Coulomb interaction p ...

Sorry, please change LDAU        = .FALSE. to .TRUE. as suggeted by wuli8 . i mad a mistake.
If not DFT+U, U　do not need these settings in INCAR, please leave it, these settings enjoy its default valuse in INCAR.
there are many ways to get the U values, well, it's another topic we can discuss, which u could reserch it emuch.
the vasp manual give a thorough descriping on the value of LMAXMIX, u can find it by ur self. the default values of LMAXMIX is 2 for s and p
赞一下 回复此楼
7楼2013-04-13 23:09:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cuiwen0314

银虫 (小有名气)
引用回帖:
7楼: Originally posted by emilyoyang at 2013-04-13 23:09:35
Sorry, please change LDAU        = .FALSE. to .TRUE. as suggeted by wuli8 . i mad a mistake.
If not DFT+U, U　do not need these settings in INCAR, please leave it, these settings enjoy its default  ...

Thank you for you explainning. I will research it much.
赞一下 回复此楼
8楼2013-04-15 09:01:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bjmaowei25

铜虫 (小有名气)
请问楼主，你是选取Ce 赝势包四个中那一个呢？

Ce   VRHFIN =Ce : [core= Kr 4d10]       很慢得收敛，如果不调参数几乎不收敛
  PAW_PBE Ce 23Dec2003                  
   12.0000000000000     
parameters from PSCTR are:
   VRHFIN =Ce : [core= Kr 4d10]
   LEXCH  = PE
   EATOM  =  1063.0861 eV,   78.1346 Ry

Ce_3  VRHFIN =Ce : [core=Xe4] s2d1f1  收敛很快，但这里含 4f 电子？
  PAW_PBE Ce_3 06Sep2000                 
   11.0000000000000     
parameters from PSCTR are:
   VRHFIN =Ce : [core=Xe4] s2d1f1
   LEXCH  = PE
   EATOM  =   890.3669 eV,   65.4401 Ry


Ce_GW                                        慢
  PAW Ce_GW 26Mar2009                    
   12.0000000000000     
parameters from PSCTR are:
   VRHFIN =Ce : [core=Kr4d]
   LEXCH  = PE
   EATOM  =  1063.0861 eV,   78.1346 Ry


Ce_h                                                    快
  PAW_PBE Ce_h 03Mar2005                 
   12.0000000000000     
parameters from PSCTR are:
   VRHFIN =Ce : [core=Xe4] s2d1f1
   LEXCH  = PE
   EATOM  =  1063.0861 eV,   78.1346 Ry
赞一下 回复此楼
追求真理
9楼2013-06-05 22:56:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cjcorange

金虫 (初入文坛)
请问下楼主现在做得怎么样了， 能不能请教一下
赞一下 回复此楼
10楼2013-11-28 13:54:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cuiwen0314 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 化学专业申博 +3 赵子羊 2026-05-23 4/200 2026-05-24 18:10 by 工大学长
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 1rx34o113h 2026-05-23 3/150 2026-05-24 17:41 by 0i3mu4vkjz
[基金申请] 评审有感 +16 popular289 2026-05-18 27/1350 2026-05-24 17:34 by hhs666
[教师之家] 论文撤稿了 +4 bjvtcliu 2026-05-24 7/350 2026-05-24 17:29 by bjvtcliu
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 hvkbtfonbv 2026-05-23 4/200 2026-05-24 17:21 by 75ui6h7z2t
[博后之家] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 hvkbtfonbv 2026-05-23 3/150 2026-05-24 17:10 by 75ui6h7z2t
[考博] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 hvkbtfonbv 2026-05-23 3/150 2026-05-24 17:01 by 75ui6h7z2t
[论文投稿] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 a2tycdlnq1 2026-05-23 4/200 2026-05-24 16:16 by hhx1yx9evi
[基金申请] 河北省自然科学基金 +6 Peterchao 2026-05-18 9/450 2026-05-24 16:02 by 130067131
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 pmo95bazuy 2026-05-23 8/400 2026-05-24 15:56 by 1uy1ht2y9r
[基金申请] 西安交大新媒学院副院长用撤稿论文结题 +3 bjvtcliu 2026-05-24 5/250 2026-05-24 10:16 by kudofaye
[教师之家] 某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路！ +4 zju2000 2026-05-21 5/250 2026-05-24 09:35 by songwz
[考博] 26/27申博自荐 10+4 ZXW0202 2026-05-22 9/450 2026-05-24 08:47 by bjvtcliu
[考博] 博士申请 +3 焦晓明 2026-05-21 3/150 2026-05-23 11:26 by mlc840311
[基金申请] 国自然评分 +4 无名者登山 2026-05-20 5/250 2026-05-21 16:35 by swuq
[基金申请] 国自然上会要求 +7 无名者登山 2026-05-18 11/550 2026-05-21 15:50 by draco1987
[基金申请] 提交了我也来说说感想 +9 fummck 2026-05-20 10/500 2026-05-21 14:17 by draco1987
[有机交流] 反应很差，大量原料没有反应 5+3 Mr.Zot 2026-05-19 8/400 2026-05-20 22:19 by Equinoxhua
[考博] 如果工作了想读博，可以边工作边读全日制嘛？ 30+3 铁达火车 2026-05-18 5/250 2026-05-20 09:33 by tfang
[考博] 博士申请 +5 星…… 2026-05-18 6/300 2026-05-18 23:49 by 糊糊涂涂好
信息提示
关闭
请填处理意见
关闭