24小时热门版块排行榜    

查看: 3499  |  回复: 43
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

jugengfans

金虫 (正式写手)

[求助] DFT+U考虑同自旋电子之间的相互作用吗?比如dxy,和dxz轨道之间的占位库伦排斥? 已有2人参与

VASP里面怎么处理+U的?
比如一个电子占据xy轨道,一个占据xz轨道,它们之间的相互作用在DFT+U里面考虑到没有?
如果按传统Hubbard模型,H=H0+U*(Nup*Ndn),是不考虑不同轨道Nup和Ndn的啊!!!!!


按照传统的理解,U应该考虑的是同一个轨道已经占据一个电子,再来一个电子所要克服的库伦能啊!
但是我觉得一个轨道占据电子,别的轨道上电子跑到这个原子附近应该也需要克服这个库伦能啊。。。。
而且我做了实验,也发现,加U以后,同自旋方向上轨道的能级也会有更大的劈裂?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

KalaShayminS

铁杆木虫 (著名写手)

【答案】应助回帖

引用回帖:
32楼: Originally posted by jugengfans at 2014-09-03 23:59:25
谢谢,看来是我理解有问题,我以为SPIN UP和SPIN DOWN的哈密根是分开的呢。...

继续回复一下评论的内容。抱歉这里需要修改错误,那一帖讨论的是对哈密顿量,也就是总能量的影响,同时间接影响了down,这个只是个二级的效应。
U对能级的直接影响:
U/2*Σ_iσ(n_i-n_in_i)对波函数变分可得
U(1/2-n_iσ)
就是说现在KS单电子方程从
H_KS ψ_iσ = ε_iσ ψ_iσ 变成
[H_KS + U(1/2-n_iσ)] ψ_iσ = ε_iσ ψ_iσ
假设波函数不变,那么对于占据态能级变成
[H_KS - U/2] ψ_iσ = ε_iσ ψ_iσ
对于未占据态态变成
[H_KS + U/2] ψ_iσ = ε_iσ ψ_iσ
就是说KS方程解出来的单电子能级,占据态下降U/2,未占据态上升U/2。
33楼2014-09-04 01:10:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 44 个回答

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

看vasp手册的公式似乎默认的LDAUTYPE=2的情形n上的自旋指标都是相同的,似乎这里Ueff确实影响的是相同自旋方向的。
不一定挂在论坛,计算问题问题欢迎留言。
2楼2014-08-30 18:39:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jugengfans

金虫 (正式写手)

引用回帖:
2楼: Originally posted by 卡开发发 at 2014-08-30 18:39:33
看vasp手册的公式似乎默认的LDAUTYPE=2的情形n上的自旋指标都是相同的,似乎这里Ueff确实影响的是相同自旋方向的。

也就是说只要有d电子,别的电子无论自旋正负,它们之间的排斥都会考虑到?这样的处理倒是比较合理吧。只是这还算不算传统意义上的U了呢。
3楼2014-08-30 19:31:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
3楼: Originally posted by jugengfans at 2014-08-30 19:31:41
也就是说只要有d电子,别的电子无论自旋正负,它们之间的排斥都会考虑到?这样的处理倒是比较合理吧。只是这还算不算传统意义上的U了呢。...

vasp手册公式(6.57)下面的公式里面的括号似乎没有体现自旋相反的项吧?

我的理解是,U还是原来的U,不过可能还有其他的约束条件,Phys. Rev. B 57, 1505 (1998)
不一定挂在论坛,计算问题问题欢迎留言。
4楼2014-08-30 19:41:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见