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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP 计算CeO2用DFT+U方法的输入如何设置
101/1返回列表
查看: 4187  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 9/450 2026-08-23 10:59 by Tide man
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 10/500 2026-08-23 07:17 by OEbVnUOu01ol
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 00:55 by OEbVnUOu01ol
[教师之家] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 今天基金会出结果吗？20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
关闭
请填处理意见
关闭