24小时热门版块排行榜    

查看: 5355  |  回复: 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-30 11/550 2026-08-31 23:16 by hittle2008
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 EZ67DG6CvcKL 2026-08-31 3/150 2026-08-31 22:23 by o3EMxHUAQDtS
[考博] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 EZ67DG6CvcKL 2026-08-31 3/150 2026-08-31 22:14 by o3EMxHUAQDtS
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[文学芳草园] 梦想 +5 myrtle 2026-08-26 5/250 2026-08-31 14:40 by hahaboy
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见