应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » DFT+U中对U和库伦关联的的理解
51/1返回列表
查看: 4312  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

jugengfans

金虫 (正式写手)

[交流] DFT+U中对U和库伦关联的的理解

在castep和vasp里面，现在对于过渡金属体系都可以采用DFT+U的方法修正带隙。
虽然大家都可以参照文献，选取对于体系合适的U值进行计算，但是对这个U的理论上的理解又是什么呢？

1.电子在格点间迁移的时候，如果某个原子位已经有电子占据，那么，这个迁移的电子跳跃过来，就需要克服这个电子对它强烈的库伦排斥。
这句话本身，应该怎样理解？某个原子位，已经有电子占据，是指包括s,p轨道的占据？还是仅考虑d轨道的占据？我认为某个原子位有d轨道的占据，这样说可能更准确一些，因为s,p轨道的非局域性通常很强，占据不占据不影响3d轨道电子的输运。

2.LDA和GGA的能量泛函形式中，到底包括不包括电子的库伦相互作用项？
我认为包括，因为在KS方程中，哈密跟量里面包括有电子之间的经典库伦相互作用能
{n（r1）n（r2）/[(r1-r2)^2]}，而把其它的关联项放入了Exc中，但是传统的GGA和LDA形式中，没有考虑电子跟自身相互作用也被这一项包括进去引起的误差（对其的修正称之为自相互作用修正），另外也没有考虑多电子体系的关联项f【n（r1，r2）】。
因此，我认为，DFT只是低估了电子之间的库伦相互作用能，并不是完全没考虑电子之间的库伦相互作用。所以，在DFT+U的计算中，需要多指定一个参数J，叫做double-counting part，去处掉已经包括在DFT计算中的能量。

但是，如果说库伦关联能，那我们应该说传统DFT中完全没有包含库伦关联能，还是说低估呢？在请各位虫友解答。


3.在Hubbard模型中，哈密跟量里面的库伦关联U，是加在同一个格点的自旋向上和自旋向下电子之间的库伦关联。
但是在DFT+U计算中，是不是这个U也包含同方向自旋（比如两个自旋向上的3d轨道电子）之间的库伦关联效应？虽然同向自旋之间有交换作用，导致它们之间不可能靠的很近，所以库伦排斥会比自旋反向的电子之间的排斥小的多。所以我的理解是可以不包括。但是为什么在计算过程中，我考虑U以后，3d轨道自旋向上的三个能级的劈裂，也会越来越明显呢？那个U到底是以什么样的形式加在哈密跟量中的？

4. 在KS方程中，理论上可以把所有的关联效应全部放入Exc中，那么通过KS方程得到的粒子数密度n应该与原来多电子的体系的粒子数密度n一模一样，那么单电子近似，岂不是没有不成立的时候了？所有的错误只不过是Exc不准确引起的？


5.另外，HK定理只针对体系的基态有效，而对体系的激发态无法给出正确的结果。那么，我们算那些空带，带隙还有什么意义？如果现在得到一个能带图，费米面上有一些空带，但是一旦受外场激发，电子跃迁到了上面的空带，那这个时候HK定理的条件就不满足了，空带的位置说不定也发生了变化，不是所有电子呆在基态时候的空带位置。那么能带理论给出的带隙还有用么？

6.有人说DFT+U是个dirty的做法，是不是就是指引入U对这个带隙的修正？那除此之外，还有什么更好的办法么？能带理论是建立在周期性边界条件导致block波的基础上的，这是单电子的block波，那么多电子体系如果关联作用很强的情况下，能带理论不是就失效了嘛！！

分享一个DFT+U方法介绍的PPT，不知道是不是以前从小木虫下的，个人认为这个PPT写的比较清楚。[ Last edited by jugengfans on 2012-4-6 at 05:02 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : dft+u.pdf
    • 2012-04-06 03:29:27, 3.41 M

» 收录本帖的淘帖专辑推荐

第一性原理 研究生资料集 理论计算及软件的使用 VASP
网贴收集 田凯 First principles 第一性原理计算
磁性计算

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2012-04-06 02:26:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

quietdistant

至尊木虫 (知名作家)
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
sunyang1988: 金币+1, 谢谢交流 2013-06-06 08:41:50
DFT+U虽然能修正带隙，但同一U值对同一材料的不同体系带来的却是不同的带隙修正，这就更加混乱了
赞一下 回复此楼
9楼2013-06-06 08:12:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

jugengfans

金虫 (正式写手)
现在没金币就没人关注么 = =
用LDA+U的应该不在少数吧
赞一下 回复此楼
2楼2012-04-06 14:02:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)
★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
WDD880227: 金币+2, 感谢您详细的解答~~辛苦了~~ 2012-04-06 14:52:22
jugengfans: 金币+5, 谢谢，只是像知道更多。另外对于3，hubbard模型里面貌似只有自旋向上和自旋向下电子自旋的乘积。是不是因为同向自旋受交换作用，不可能离得很近，导致其库伦作用比反向时弱得多，可以忽略？ 2012-04-06 22:59:17
1.指的是d电子之间的相互作用
2.LDA和GGA都已经考虑了库伦作用，但是基于平均场的考虑，没有考虑很局域的强关联
3.库伦作用是静电作用，应该和自旋无关吧？
4.单电子近似无法计算多体效应，比如kondo效应
5.HK定理第一定理是普适的，并不局限于基态
6.看近似程度来判断计算的精度，一般体系DFT的描述还是挺好的
赞一下(2人) 回复此楼

» 本帖已获得的红花（最新10朵）

jugengfans
3楼2012-04-06 14:19:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)
引用回帖:
2楼: Originally posted by jugengfans at 2012-04-06 14:02:16:
现在没金币就没人关注么 = =
用LDA+U的应该不在少数吧

其实您也可以反过来问一个问题：求助了这么多问题，为什么连几个金币都不舍得给呢？
赞一下 回复此楼
4楼2012-04-06 14:21:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
黄河honour12楼
2014-12-13 23:51   回复  
up [ 发自小木虫客户端 ]
查看全部 12 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭