24小时热门版块排行榜    

查看: 1760  |  回复: 8
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

huangyc

木虫 (正式写手)

[交流] 【求助】关于L(S)DA+U的疑问

1、对AB合金体系,怎么样确定U和J的值。
比如研究CO在AB合金上的吸附。若AB是1:1时,可以选择不同的U值优化晶胞,取晶胞参数(或费米能)与实验数据接近的U值建模并作下一步计算(AB都是过渡金属恐怕很费时间和精力);但AB计量比不是1:1且没有实验数据时,如何确定U值?
2、L(S)DA+U方法适用于强关联体系,比如常见的CO与过渡金属的吸附。但过渡金属中有些原子价轨道是充满的,LDAUL又怎么设置?
举个例子,如Pd 4d105s0,与CO形成电子授受键并不是靠d轨道,LDAUL=1?又如Cd 4d105s2,LDAUL=-1?
关于CO adsorption puzzle的问题,讨论最多的是CO/Pt体系,对于价层充满的过渡金属来说,是否意味着不是强的关联体系,不需要+U处理。
望各位指点迷津,不甚感激!

[ Last edited by mingdong on 2009-9-15 at 09:24 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yzcluster

金虫 (著名写手)

小木虫扫盲人

★ ★ ★ ★ ★ ★ ★ ★
wuli8(金币+3,VIP+0):(*^__^*)谢谢,辛苦了。 3-5 22:46
huangyc(金币+5,VIP+0):谢谢 3-9 20:11
我来说一下。
LDA+U或是GGA+U中的U指的是HUBBARD模型中的自旋相反电子的强关联排斥能。通常的DFT计算是单电子近似的,或者是平均场近似,在这种近似下,的确可以解释很多问题。但是,真正的凝聚态体系,分子或是固体,本质上是多体理论。很多问题,比如MOTT绝缘体,在单电子近似下的结论就是错的。因此必须超越单电子近似,考虑体系的量子多体效应。U就是在这种情况下出现的。
简单的说,或者在HUBBARD模型一级近似下,U就是考虑了同一个原子上自旋相反的局域电子之间的库仑排斥,从而导致了能带的“重正化”。比如所,在普通能带理论中,如果原胞中的电子数是奇数个,则体系一定是金属。但是在考虑过在位格点(ON-SITE)的排斥能后,就出现了绝缘体。
因此,考虑U的情况一般是体系中存在局域化的d,f电子(比如Ce)。这些电子的存在是”标准“的DFT理论无法合理处理的,必须要考虑U的影响。
而通常U是一个试验值,不同的化学环境,同一个原子的U就可能不同,必须通过试验来确定。
呵呵,先说这么多,累死了。
5楼2009-03-05 20:16:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

sss3012

银虫 (正式写手)

★ ★ ★
wuli8(金币+1,VIP+0):(*^__^*)辛苦了 3-4 09:27
huangyc(金币+2,VIP+0):谢谢。我这个是合金,只有晶胞参数等数据,不同比例的什么都没有了,怎么找实验值! 3-5 19:30
LDA+U 中的U最好参考实验值,如果找不到实验值就要自己去试,通常认为带隙或晶格常数不再发生变化时的U值作为近似参考值
2楼2009-03-03 22:46:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chineseli

★ ★ ★ ★
spur(金币+2,VIP+0): 3Q!!谢谢参与! 3-5 18:39
huangyc(金币+2,VIP+0):谢谢。“貌似”我还是不太放心啊。下面说的加拿大的什么方法恐怕难以实现 3-5 19:33
貌似一般算Pd不用加U的。U用的多的是半导体和有f电子的,比如Ce。

看过Emily Carter一个学生做ab initio DFT+U,现在在加拿大的Queen's, 有兴趣的话不妨去了解一下,拿那个方法来试试看~
3楼2009-03-05 14:16:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nkleof

银虫 (正式写手)


huangyc(金币+1,VIP+0):谢谢讨论 3-9 20:10
有个问题请教一下用过L(S)DA的各位,vasp中是要加U就必须对所有原子都加吗?
我看到user guide里面有一句话
NB: LDAUL, LDAUU and LDAUJ must be specified for all atomic species!
没看懂,请指点一下。而且在user guide 里面没有提到如何对某些原子单独设置U值,请问该怎么实现?
先谢谢了!
4楼2009-03-05 19:35:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 什么时候开奖? +9 CrisMessi 2026-08-18 10/500 2026-08-23 12:03 by 丶昵称占用
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 9/450 2026-08-23 10:44 by IXZuIJ2Q7OVy
[基金申请] 让我中一个面上吧! +12 大萍1987 2026-08-20 14/700 2026-08-23 10:30 by wrm
[基金申请] 今天放榜吗? +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一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:43 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 11/550 2026-08-23 02:55 by OEbVnUOu01ol
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 今天基金会出结果吗?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
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 重要消息,中午系统在维护 +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
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见