| 查看: 4342 | 回复: 8 | ||||
cos60木虫 (著名写手)
|
[求助]
CASTEP LDA+U 的U值能设置成负的吗？ 已有2人参与
|
|
CASTEP LDA+U 的U值能设置成负的吗？ 谢谢！ |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理 | 刺客 |
» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有73人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
-
招聘 金属平磨液，抛光液研发工程师
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- CASTEP对晶胞几何优化可以设置不同的温度吗？还有压力
已经有8人回复
- LDA+U或GGA+U，U值选取问题
已经有14人回复
- 在castep，如何设置超软赝势，降低截断能？请大家看看内容中提到的如何设置。
已经有12人回复
- 【求助】castep 优化铁时 initial spin该如何设置
已经有26人回复
- 【求助】CASTEP中可以实现GGA+U算法吗
已经有6人回复
- 【求助】CASTEP中只能有LDA+U没有GGA+U嘛？
已经有6人回复
- 【求助】请教MS CASTEP模块中的LDA+U 计算
已经有6人回复
sunboyccc
铁杆木虫 (正式写手)
- 应助: 79 (初中生)
- 金币: 5418.4
- 红花: 10
- 帖子: 352
- 在线: 182.6小时
- 虫号: 674686
- 注册: 2008-12-14
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2012-03-28 07:24:01
sunboyccc
铁杆木虫 (正式写手)
- 应助: 79 (初中生)
- 金币: 5418.4
- 红花: 10
- 帖子: 352
- 在线: 182.6小时
- 虫号: 674686
- 注册: 2008-12-14
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2012-03-28 07:28:44
shmilylxc
木虫 (正式写手)
- 应助: 139 (高中生)
- 金币: 3765.6
- 散金: 111
- 红花: 5
- 帖子: 533
- 在线: 387.1小时
- 虫号: 1566205
- 注册: 2012-01-06
- 专业: 半导体材料
4楼2012-03-28 12:34:38
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
贺仪: 金币+5, 谢谢指教！ 2012-03-28 14:57:28
感谢参与，应助指数 +1
贺仪: 金币+5, 谢谢指教！ 2012-03-28 14:57:28
|
还是老话，先弄清楚U的物理含义，U是用于弥补分数占据与整数占据之间的误差的，并且是物理自洽的，不是用来fit to band gap的，尽管合理的U可以很大程度上使band gap贴近实验值。但是仅仅用U调整到完全一样，似乎意义也不大，比如简单的+U没有考虑GW中那样的声子矫正作用之类的。而且，U作为一个有物理意义的数值，也是通过自洽来实现的。比较逗的是，在各个DFT软件中，都推荐Cococcioni的自洽方法确定U值，唯独CASTEP，不但不提供求此值必要的参数，并且还在help中写明，由于直接fit band gap得到的U比自洽得到的U给的带隙之类的结果更好，所以不提供此功能…… 如果仅仅是想调整带隙，干嘛还搞+U啊？直接用剪刀算符就好了～跟自己直接在Origin上修改个数据一样简单有效，只不过物理意义几乎没有～  |
5楼2012-03-28 12:44:36
sunboyccc
铁杆木虫 (正式写手)
- 应助: 79 (初中生)
- 金币: 5418.4
- 红花: 10
- 帖子: 352
- 在线: 182.6小时
- 虫号: 674686
- 注册: 2008-12-14
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★
WDD880227: 金币+2, 感谢您的详细分析~~~ 2012-03-29 09:28:11
WDD880227: 金币+2, 感谢您的详细分析~~~ 2012-03-29 09:28:11
|
1.从唯象角度，U的物理意义可以这样理解，（一家之言不一定对，呵呵）当电子想去它不该去的地方时，U是一个势垒，这样，这样就可以克服DFT方法中特有的缺陷，不该有相互作用的不同壳层就不会产生作用。 2.此外，对于掺杂的情况，不加U时可能没有带隙，这时剪刀算符并不起作用，而加U才产生带隙。 3.事实上CASTEP与VASP两种软件加U并无本质区别，新版VASP中有一种简单的只设U值的加U方法，就是CASTEP的矩阵对角线方法，当然VASP中能够设U与J的方法更好。 4.利用U值调整带隙，并且使掺杂能级位于带隙中，对于过渡金属掺杂的磁性体系中，电子的正确的跃迁非常重要。例如可以使Co原子自旋向上、向下能级劈裂，向下的能级实际应该位于带隙中，如果不加U或加U后带隙不正确，向下的能级可能位于导带中，那么本来应该跃迁到自旋向下能级的电子就会跃迁到导带中，结果可想而知。 一点愚见，呵呵 |
6楼2012-03-29 08:23:25
7楼2015-03-17 18:33:14
qnclqq
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 3862.6
- 红花: 4
- 帖子: 467
- 在线: 132.8小时
- 虫号: 1464347
- 注册: 2011-10-28
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
我使用QE、abinit等都采用Cococcioni教授的方法计算过TiO2的Ti-3d的U值，并不能对带隙有什么好的修正。而且查阅当时他们FeO计算QE的输入文件，对总自旋进行了严格的限定才能得到正确的带隙。 就近两年出版的书籍中，有几种观点： 1. KS态本来就不能用于描述激发态，所以修正的物理意义不是那么明显。 Interacting Electrons: Theory and Computational Approaches 2. SIE并不是造成DFT描述带隙不准确的根本原因，这与1类似，KS态的物理意义。 Density Functional Theory in Quantum Chemistry 3. 有些论文中认为+U对于任何原子的电子都是适用的，反正是对当作库仑作用的修正。 所以，个人的看法，MS没有这种方法，或者说用带隙来拟合U的参数都是没有问题的。关键是在拟合带隙的同时，还要对结构参数进行拟合。因为如果考虑晶格中不同格点的同种或者异种原子的不同电子都+U，那么符合的带隙参数应该一个多维矩阵。从中选取一组各方面参数都合适的应该也是一种方法。 个人愚见，请指教。 |
8楼2017-04-16 20:36:09
9楼2017-12-26 20:24:27
25