| 查看: 3592 | 回复: 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. 请大家帮忙看看这个输入对不对啊？谢谢啦 |
回复此楼» 猜你喜欢
0703化学调剂，求各位老师收留
已经有8人回复
-
070303一志愿西北大学学硕310找调剂
已经有8人回复
-
085600材料与化工 求调剂
已经有14人回复
-
本人考085602 化学工程 专硕
已经有13人回复
-
材料与化工求调剂
已经有3人回复
-
283求调剂
已经有13人回复
-
312求调剂
已经有4人回复
-
285求调剂
已经有6人回复
-
材料与化工专硕调剂
已经有3人回复
-
0703化学调剂 ，六级已过，有科研经历
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 【原创】关于结构弛豫设置及收敛判据选择的一点体会
已经有140人回复
emilyoyang
木虫 (正式写手)
- 应助: 232 (大学生)
- 金币: 4340.4
- 红花: 42
- 帖子: 667
- 在线: 189.9小时
- 虫号: 1630291
- 注册: 2012-02-20
- 性别: MM
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
cuiwen0314: 金币+10, ★★★很有帮助 2013-04-13 15:31:50
cuiwen0314: 金币+10, ★★★很有帮助 2013-04-13 15:39:20
感谢参与，应助指数 +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
cuiwen0314
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 400.7
- 红花: 1
- 帖子: 247
- 在线: 71.2小时
- 虫号: 1102818
- 注册: 2010-09-19
- 性别: MM
- 专业: 无机合成和制备化学
|
感谢您的回复！有一点我不明白的是如果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
wuli8
荣誉版主 (知名作家)
…………
- 1ST强帖: 2
- 应助: 35 (小学生)
- 贵宾: 12.924
- 金币: 20189.4
- 散金: 15888
- 红花: 88
- 沙发: 4
- 帖子: 7840
- 在线: 1114.6小时
- 虫号: 465889
- 注册: 2007-11-23
- 专业: 物理学I
- 管辖: 计算模拟
4楼2013-04-13 21:54:35
cuiwen0314
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 400.7
- 红花: 1
- 帖子: 247
- 在线: 71.2小时
- 虫号: 1102818
- 注册: 2010-09-19
- 性别: MM
- 专业: 无机合成和制备化学
|
谢谢您的解释！ 那您能给我解答下面这个问题吗？ 感谢您的回复！有一点我不明白的是如果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
emilyoyang
木虫 (正式写手)
- 应助: 232 (大学生)
- 金币: 4340.4
- 红花: 42
- 帖子: 667
- 在线: 189.9小时
- 虫号: 1630291
- 注册: 2012-02-20
- 性别: MM
- 专业: 理论和计算化学
6楼2013-04-13 22:56:07
emilyoyang
木虫 (正式写手)
- 应助: 232 (大学生)
- 金币: 4340.4
- 红花: 42
- 帖子: 667
- 在线: 189.9小时
- 虫号: 1630291
- 注册: 2012-02-20
- 性别: MM
- 专业: 理论和计算化学
|
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
cuiwen0314
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 400.7
- 红花: 1
- 帖子: 247
- 在线: 71.2小时
- 虫号: 1102818
- 注册: 2010-09-19
- 性别: MM
- 专业: 无机合成和制备化学
8楼2013-04-15 09:01:18
bjmaowei25
铜虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 362.1
- 红花: 1
- 帖子: 121
- 在线: 293.5小时
- 虫号: 1531554
- 注册: 2011-12-10
- 性别: GG
- 专业: 结构陶瓷
|
请问楼主，你是选取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
10楼2013-11-28 13:54:18
20