24小时热门版块排行榜    

查看: 5324  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 没有任何消息-是不是就凉了 +7 图啦图啦 2026-08-24 8/400 2026-08-24 22:00 by maomao_da
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-24 21:59 by BZKMTicpDhFj
[基金申请] 人气不行了 +10 fansofjerry 2026-08-21 10/500 2026-08-24 21:03 by zhanghaozhu
[基金申请] 我面上完蛋了 +9 且听虎啸 2026-08-20 10/500 2026-08-24 20:08 by 苦寒来香
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 19/950 2026-08-24 16:52 by iaeyuan
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天基金会出结果吗?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
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
请填处理意见