24小时热门版块排行榜    

CyRhmU.jpeg
查看: 4294  |  回复: 14

xiaojie7783

金虫 (小有名气)

[交流] 【求助】U值的确定已有3人参与

请教大家,计算氧化物掺杂的时候想要用GGA+U,这个时候的元素的U值和J值应该怎么确定呢?是不是要计算一下掺杂元素的氧化物的能隙值,当此时的能隙值和实验值接近的时候对应的U值就是掺杂计算时应该选用的U值呢,还是别的办法呢?谢谢大家!
        另外还有一个问题,就是J的值要怎么确定呢?再用GGA+U计算的时候,LDAUJ这个参数是不是也要设置上呢?又该怎么确定呢?
         如果我不选用自旋的话,不设置ISPIN的值,那么此时LDAUTYPE能选择为2吗?还是应该为4呢?
回复此楼
成功其实很简单,就是当你坚持不住的时候,再坚持一下!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

锐利的碎片

木虫 (正式写手)

star watcher


lvjian8596(金币+1):感谢交流! 2010-07-16 10:02:19
j一般固定0.8左右就行了。U值可以参考文献中的大小。
2楼2010-07-16 08:42:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

★ ★
lvjian8596(金币+2):感谢交流! 2010-07-16 15:28:22
xiaojie7783(金币+8):谢谢,我还想请问:如果我不打开自旋(用GGA+U算)的话,那么可以不设置ISPIN的值然后将LDAUTYPE能选择为2吗?还是也应该设置为为4呢? 2010-07-16 23:47:10
J的值的意义与LDAUTYPE有关的。当LDAUTYPE为2时,U-J的差值才有有意义,即有效的U参数。

至于有效U参数的值怎么选取,通过线性响应的方式来估测,类似pwscf里面的。要么将U值从小测到大,并计算结果同实验结果(比如电子结构等)进行比较,然后选取一个合适的值。再要不然,参考文献中的值。

[ Last edited by valenhou001 on 2010-7-17 at 09:18 ]
3楼2010-07-16 15:22:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaojie7783

金虫 (小有名气)

引用回帖:
Originally posted by valenhou001 at 2010-07-16 15:22:38:
J的值的意义与LDAUTYPE有关的。当LDAUTYPE为1或4时,U-J的差值才有有意义,即有效的U参数。

至于有效U参数的值怎么选取,通过线性响应的方式来估测,类似pwscf里面的。要么将U值从小测到大,并计算结果同实验结 ...

谢谢,
我看说明书里面说的是:in Dudarev’s approach the parametersU and J do not enter seperately, only the difference (U ¡J) is meaningfull.
并且对于LDAUTYPE,
1 Rotationally invariant LSDA+U according to Liechtenstein et al.
4 Idem 1., but LDA+U instead of LSDA+U (i.e. no LSDA exchange splitting)
2 Dudarev’s approach to LSDA+U (Default)
那是不是说LDAUTYPE=2时只有U-J的值有意义呢?

我还想请问一下:如果我不打开自旋(用GGA+U算)的话,那么可以设置ISPIN的值为1然后将LDAUTYPE能选择为2吗?还是说非自旋计算时只能将LDAUTYPE设置为4呢?

[ Last edited by xiaojie7783 on 2010-7-16 at 23:58 ]
成功其实很简单,就是当你坚持不住的时候,再坚持一下!
4楼2010-07-16 23:46:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaojie7783

金虫 (小有名气)

引用回帖:
Originally posted by 锐利的碎片 at 2010-07-16 08:42:52:
j一般固定0.8左右就行了。U值可以参考文献中的大小。

谢谢,我还想请问:如果我不打开自旋(用GGA+U算)的话,那么可以设置ISPIN的值为1然后将LDAUTYPE能选择为2吗?还是说非自旋计算时只能将LDAUTYPE设置为4呢?
谢谢
成功其实很简单,就是当你坚持不住的时候,再坚持一下!
5楼2010-07-16 23:56:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)


gavinliu7390(金币+1):谢谢交流! 2010-07-17 12:00:43
有关U,J的意义,前一个回帖中说的有点小误,我已经更正过来了。

既然LSDA+U里面U,J参数在目前来说,用的时候,究竟什么样的值才有意义,难以说的清楚,所以尽量减少不明确性,推荐采用默认的方式,即LDAUTYPE=2,要给定的参数就是U-J的值了。LDAUTYPE=2的效果和LDAUTYPE=1或4的都是差不多的。
6楼2010-07-17 09:38:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)


gavinliu7390(金币+1):谢谢交流! 2010-07-17 12:00:48
引用回帖:
Originally posted by xiaojie7783 at 2010-07-16 23:56:15:
谢谢,我还想请问:如果我不打开自旋(用GGA+U算)的话,那么可以设置ISPIN的值为1然后将LDAUTYPE能选择为2吗?还是说非自旋计算时只能将LDAUTYPE设置为4呢?
谢谢

对非自旋选用LDAUTYPE=4,而自旋极化是选用LDAUTYPE=2这样不行。为了比较自旋极化和非自旋极化的稳定性,最好要么都是用LDAUTYPE=2,要么LDAUTYPE=1对自旋极化,LDAUTYPE=4对非自旋极化。
7楼2010-07-17 10:08:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

尹珺加油

金虫 (小有名气)

弱弱的问,在modify里设置的u值为u-j吗?
幸福源自我们的挚爱!
8楼2011-04-05 19:13:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fyh987

金虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
545594楼: Originally posted by valenhou001 at 2010-07-17 10:08:46
对非自旋选用LDAUTYPE=4,而自旋极化是选用LDAUTYPE=2这样不行。为了比较自旋极化和非自旋极化的稳定性,最好要么都是用LDAUTYPE=2,要么LDAUTYPE=1对自旋极化,LDAUTYPE=4对非自旋极化。...

弱弱的问下LZ,加U是对体系里面所有需加U的原子加U吗
努力挑战自己
9楼2012-05-29 20:43:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hanyanli0475

铜虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
544562楼: Originally posted by valenhou001 at 2010-07-17 09:38:09
有关U,J的意义,前一个回帖中说的有点小误,我已经更正过来了。

既然LSDA+U里面U,J参数在目前来说,用的时候,究竟什么样的值才有意义,难以说的清楚,所以尽量减少不明确性,推荐采用默认的方式,即LDAUTYPE=2 ...

侯老师您好,不好意思打扰了,现在有个问题,实在想不明白,想请教您。我最近用http://prl.aps.org/abstract/PRL/v97/i10/e103001
文献上的方法,对U值进行优化(pwscf),我的材料是InN,但是对铟电子加U时,Uout出来的都是0.1多,都特别小,而且拟合曲线的斜率是正值,不知道哪里弄错了呢??求侯老师指点迷津,谢谢哈!这样提问有点唐突,现在有些着急了,请见谅。

我输出的Uout是在这里找的,对吗?
one-electron contribution =     -79.26562683 Ry
     hartree contribution      =      68.78255490 Ry
     xc contribution           =     -94.90442606 Ry
     ewald contribution        =    -207.18697592 Ry
     Hubbard energy            =       0.16900435 Ry
     smearing contrib. (-TS)   =       0.00027725 Ry
什么时候能算完!!!!
10楼2012-06-13 14:31:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaojie7783 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见