24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知(长期有效)
查看: 1382  |  回复: 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的回帖

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的回帖

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的回帖

huangyc

木虫 (正式写手)

谢谢各位的讨论。不过我的问题好象仍未解决...
6楼2009-03-09 20:12:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

y1ding

铁杆木虫 (著名写手)


zdhlover(金币+1,VIP+0):谢谢 3-13 19:49
引用回帖:
Originally posted by nkleof at 2009-3-5 19:35:
有个问题请教一下用过L(S)DA的各位,vasp中是要加U就必须对所有原子都加吗?
我看到user guide里面有一句话
NB: LDAUL, LDAUU and LDAUJ must be specified for all atomic species!
没看懂,请指点一下。而且 ...

不是对于原子加U,VASP是对于同一类原子加相同的U。
7楼2009-03-09 20:47:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)

这个+U,加得人心烦意乱。比如我研究金属A的结合能,对A要用LDA+U,那么我要计算晶体A的能量和孤立A的能量,这样的话,是不是要对晶体A和孤立A加同样的U啊?
然后,我要计算晶体AB的结合能,这个时候,对晶体AB中的A,照道理应该是加不同的U了,可是这个时候同样要用到孤立A的能量,这时候对孤立A的U又是怎么个加法啊?

同样的原子,不同的环境,竟然要加不同的U,真是太烦了。
8楼2009-03-09 21:09:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zeoliters

木虫 (正式写手)

Vasp&MedeA存档
9楼2009-06-13 21:46:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 huangyc 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见