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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP 计算CeO2用DFT+U方法的输入如何设置
北京石油化工学院2026年研究生招生接收调剂公告
101/1返回列表
查看: 3669  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0855求调剂材料 +10 红桃灼灼 2026-04-04 10/500 2026-04-05 18:08 by 猪会飞
[考研] 调剂 +4 不逢春 2026-04-05 4/200 2026-04-05 18:07 by ms629
[考研] 085500机械专硕初试288求调剂 +3 GZJguo666- 2026-04-05 3/150 2026-04-05 18:06 by jkddd
[考研] 312求调剂 +3 Say Never 2026-04-04 3/150 2026-04-05 11:19 by guoweigw
[考研] 359求调剂22408 +3 123456789qw 2026-03-31 3/150 2026-04-05 10:09 by zhq0425
[考研] 材料化工306分找合适调剂 +14 沧海轻舟e 2026-04-04 14/700 2026-04-05 09:53 by 朱云虎202
[考研] 求 +5 化工专硕323分 2026-04-04 5/250 2026-04-05 08:02 by 544594351
[考博] 申博 +7 IQwQl 2026-04-04 7/350 2026-04-04 23:32 by mumin1990
[考研] 一志愿郑州大学材料与化工085600，求调剂 +24 吃的不少 2026-04-02 24/1200 2026-04-04 23:20 by 永字号
[考研] 283分求调剂 +7 小聂爱学习 2026-04-03 7/350 2026-04-04 21:51 by hemengdong
[考研] 321求调剂 +6 认真求上学 2026-04-03 6/300 2026-04-04 19:51 by dongzh2009
[考研] 359求调剂 +7 hhhhaaaa$ 2026-04-04 7/350 2026-04-04 18:49 by imissbao
[考研] 294求调剂 +6 Grey_Ey 2026-04-03 6/300 2026-04-03 20:46 by 欣喜777
[考研] 321求调剂 +17 y-yh 2026-04-01 20/1000 2026-04-03 12:57 by y-yh
[考研] 273求调剂 +20 李芷新1 2026-03-31 20/1000 2026-04-03 09:58 by linyelide
[考研] 一志愿武汉理工0856，初试334 +3 26考研材料 2026-04-02 3/150 2026-04-02 21:22 by dongzh2009
[考研] 085602化学工程268分蹲调剂 +8 月照花林。 2026-04-01 8/400 2026-04-01 22:08 by 无际的草原
[考研] 0817化工学硕调剂 +11 努力上岸中！ 2026-03-31 11/550 2026-04-01 20:30 by 赖春艳
[考研] 358求调剂 +3 王向阳花 2026-03-31 3/150 2026-04-01 09:56 by zzchen2000
[考研] 考研生物与医药调剂 +7 铁憨憨123425 2026-03-31 7/350 2026-04-01 08:45 by JourneyLucky
信息提示
关闭
请填处理意见
关闭