| 查看: 9865 | 回复: 20 | |||||||||||||||||||||
[求助]
LDA+U,为什么要加U，作用是什么？
|
|||||||||||||||||||||
|
关于LDA+U，查过说处理d或f电子时的方法，为什么要加U呢，还是没弄明白，U起的作用？ 谢谢 |
回复此楼» 收录本帖的淘帖专辑推荐
 VASP and MS | 第一性原理 | 第一原理资料汇编 | 仿真建模与计算 |
| VASP | first principle | ATK和电子输运 | VASP |
| 计算资料 | Photochemistry | 模拟 | 刺客 |
| 第一性原理计算 | 第一性原理计算 方法、理论、参数设置 | 第一性原理计算 | 第一性原理 |
| 稀土 赝势 加U | 其他 | 第一性原理杂记 |
» 猜你喜欢
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有3人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有21人回复
-
能否退出参与的面上项目解除限项
已经有17人回复
-
估计是周四
已经有4人回复
-
我面上完蛋了
已经有9人回复
-
2026国自然函评费到账
已经有18人回复
-
建议基金发布提前给出明确的时间点
已经有15人回复
-
范进中举一文的中心思想
已经有7人回复
-
朋友圈看到的
已经有9人回复
-
让我中一个面上吧！
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教，为何采用LDA+U计算的时候Co变成了低自旋？
已经有7人回复
- 关于磁性问题请教
已经有7人回复
- LDA+U到底能够对能带的计算有多大的改善呢
已经有22人回复
- 【求助】LDA+U优化问题
已经有10人回复
- 【求助】请教PWSCF晶体建模和自洽计算，不收敛
已经有21人回复
韩D小希: 回帖置顶 2012-02-23 11:29:59
uuv2010(金币+2): 多谢分享！ 2012-02-24 10:23:20
uuv2010(金币+2): 多谢分享！ 2012-02-24 10:23:20
|
您好，谢谢您的回答，我在网上找个些资料，大致理解如下，您帮我看下有哪些地方有出入的，谢谢 对于一般的体系，电子间的作用不像d、f电子那样强烈，所以在涉及其交换关联能时，只需要考虑一项，即电子在不同原子间跳跃的动量项，这也就是传统的LDA或GGA近似；但是对于过渡金属之类的电子作用十分强的体系，利用传统的算法得到的结果不够精准，于是Hubbard model 出现了，它除了考虑前面说过的动量项之外顾忌了两个电子在一个轨道或不同轨道间的相互作用，这项就是U，也就是电子间的库仑排斥项；对于U值得确定，我查到有几种方法： 1.计算，即分别计算所需离子 得到和失去一个电子能量，然后利用下面的式子得到U； 2.凭查找文献（个人觉得这种方法不容易，万一是难找到的体系怎么办） 3.凭测试，多试几个U ，但是评判的标准是什么呢，是算出的带隙和实验值一样么，还有人说是是要看加U后d电子的能级与实验值d电子的能级是否对应 ？哪一种后者怎么做呢 以上是个人的一点理解，希望您能指导一下，十分感谢。  |
7楼2012-02-23 11:29:21
wbt11129
木虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 2378.5
- 散金: 1444
- 帖子: 349
- 在线: 122.9小时
- 虫号: 834541
- 注册: 2009-08-26
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
韩D小希: 回帖置顶 2012-02-25 08:10:51
uuv2010(金币+3): 多谢详细解答 2012-02-25 08:42:44
uuv2010(金币+3): 多谢详细解答 2012-02-25 08:42:44
|
1.U 的物理本质是什么？ 2.加U 最早是为处理何种体系而引入的？ 3.为何要加U ? 4. 对哪些体系需要加U ? 5. 如何确定加U 值的大小？ 6加U 带来的典型后果是什么？ 7.加U 的方法会有哪些局限性 我试着回答一下，1，U值描述屏蔽库伦排斥能的大小，U主要是为了准确描述电子间强的库伦排斥引进的，传统的DFT没有考虑。2，没有细看方法类文献，不清楚，但是知道典型的研究材料是NiO和UO2。 3，参1。4，一般在研究镧系和锕系化合物时要对4f或者5f点子加U，对含d电子的过渡金属化合物中和O,N等成键的金属元素d电子加U。5，一般采用和实验电子能谱对比的方法确定U值，参（PRB81,045119），俺的文章，欢迎拍砖  。6，加U以后一般会使得电子能谱更局域化，峰值变平缓（此为加U计算的缺陷），DFT得到金属性，加U会变成Mott绝缘体，更好的计算能隙，原子和电子结构，磁性稳定性等。7，加U方法依赖参数和实验，各种材料需要测试U值，不能提前确定U值，所以从根本上讲，此方法有很大的局限性，比较机械，但是它近几年也算是风光了。最近比较热门的方法是LDA+DMFT。正在挑战加U方法的领军地位。先说到这，慢慢交流。  |
11楼2012-02-25 00:29:11
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
★ ★
fzx2008: 金币+2, 专家考核, 谢谢指教 2012-03-17 21:25:49
fzx2008: 金币+2, 专家考核, 谢谢指教 2012-03-17 21:25:49
|
其实很早的一篇经典的RPB论文就已经讨论了LDA+U了，其中涉及到了U的物理含义，确定方法。17楼的答案已经很好了，我就顺便再给出文章来吧，这篇文章是从作者的PhD毕业论文中汲取精华写成的，篇幅很长，论述很详细，后来作者好像去了MIT，再后来就不清楚了，咱不是狗仔队的…… 本文基于PWscf计算结果，可广泛应用于Abinit和VASP，未见CASTEP的报道，估计跟本人不怎么关心CASTEP有关，不太清楚是否也可用于FPLAPW方法，这个就欢迎大家补充了。不管如何，该文章是使用线性响应确定U值的鼻祖级的了，不是说以前没有，而是没有这么完善的：Phy Rev B, 71, 035105 (2005) |
20楼2012-03-17 20:28:14
dujunhui
金虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 1769.3
- 帖子: 403
- 在线: 136.3小时
- 虫号: 739582
- 注册: 2009-04-03
- 性别: GG
- 专业: 无机材料化学
6楼2012-02-23 09:57:07
fzx2008
荣誉版主 (著名写手)
- 1ST强帖: 9
- 应助: 272 (大学生)
- 贵宾: 1.045
- 金币: 21814.1
- 散金: 710
- 红花: 33
- 沙发: 1
- 帖子: 1808
- 在线: 1601小时
- 虫号: 621341
- 注册: 2008-10-09
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 第一性原理
【答案】应助回帖
感谢参与，应助指数 +1
uuv2010(金币+1): 多谢提示 2012-02-22 22:15:50
韩D小希(金币+3): ★有帮助 谢谢 2012-02-24 16:09:30
uuv2010(金币+1): 多谢提示 2012-02-22 22:15:50
韩D小希(金币+3): ★有帮助 谢谢 2012-02-24 16:09:30
2楼2012-02-22 20:20:47
sunboyccc
铁杆木虫 (正式写手)
- 应助: 79 (初中生)
- 金币: 5418.4
- 红花: 10
- 帖子: 352
- 在线: 182.6小时
- 虫号: 674686
- 注册: 2008-12-14
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2012-02-22 21:48:21
|
您好，谢谢您的回答，我在网上找个些资料，大致理解如下，您帮我看下有哪些地方有出入的，谢谢 对于一般的体系，电子间的作用不像d、f电子那样强烈，所以在涉及其交换关联能时，只需要考虑一项，即电子在不同原子间跳跃的动量项，这也就是传统的LDA或GGA近似；但是对于过渡金属之类的电子作用十分强的体系，利用传统的算法得到的结果不够精准，于是Hubbard model 出现了，它除了考虑前面说过的动量项之外顾忌了两个电子在一个轨道或不同轨道间的相互作用，这项就是U，也就是电子间的库仑排斥项；对于U值得确定，我查到有几种方法： 1.计算，即分别计算所需离子 得到和失去一个电子能量，然后利用下面的式子得到U； 2.凭查找文献（个人觉得这种方法不容易，万一是难找到的体系怎么办） 3.凭测试，多试几个U ，但是评判的标准是什么呢，是算出的带隙和实验值一样么，还有人说是是要看加U后d电子的能级与实验值d电子的能级是否对应 ？哪一种后者怎么做呢 以上是个人的一点理解，希望您能指导一下，十分感谢。  |
9楼2012-02-23 11:32:05
来去无风6988
木虫 (著名写手)
- 应助: 7 (幼儿园)
- 金币: 3228.6
- 帖子: 1361
- 在线: 1351.5小时
- 虫号: 1035874
- 注册: 2010-06-04
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
感谢参与，应助指数 +1
WDD880227(金币+1): 感谢交流 2012-02-23 20:10:40
韩D小希(金币+2): ★有帮助 谢谢 2012-02-24 16:10:24
WDD880227(金币+1): 感谢交流 2012-02-23 20:10:40
韩D小希(金币+2): ★有帮助 谢谢 2012-02-24 16:10:24
|
我来尝试解答一下（详细请参照Hubbard 模型）： U的物理本质：两个电子占据统一轨道需要克服库伦排斥能，导致能量增加，称为Hubbard能U。 加U是为了处理中强关联体系，例如在Mott绝缘体中，通常的能带理论不考虑电子电子相互作用，而在中强关联体系通必须考虑。但在考虑如分数量子Hall效应，以及超导现象等强关联体系时，加U也失去作用，因此DFT不能计算超导材料。 3、4问在以上能找到答案，就不废话了。 我的感觉是加U的大小是根据经验决定的，或者对于体系的认识来估计，然后通过迭代算法来优化。 加U的典型后果一般确实是使能隙变大。如考虑CoO，如果按照通常的能带理论计算，是金属，而实际上它是一种Mott绝缘体，只有通过加U才能准确算得起能带结构。 |
10楼2012-02-23 18:52:49
17楼2012-02-27 20:21:02
3楼2012-02-22 20:33:43
5楼2012-02-23 09:11:50
8楼2012-02-23 11:30:57
20