 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 3612 | 回復(fù): 43 | ||||||
jugengfans金蟲 (正式寫手)
|
[求助]
DFT+U考慮同自旋電子之間的相互作用嗎？比如dxy，和dxz軌道之間的占位庫(kù)倫排斥？ 已有2人參與
|
|||||
|
VASP里面怎么處理+U的？ 比如一個(gè)電子占據(jù)xy軌道，一個(gè)占據(jù)xz軌道，它們之間的相互作用在DFT+U里面考慮到?jīng)]有？ 如果按傳統(tǒng)Hubbard模型，H=H0+U*(Nup*Ndn)，是不考慮不同軌道Nup和Ndn的�。。。。�！ 按照傳統(tǒng)的理解，U應(yīng)該考慮的是同一個(gè)軌道已經(jīng)占據(jù)一個(gè)電子，再來一個(gè)電子所要克服的庫(kù)倫能啊！ 但是我覺得一個(gè)軌道占據(jù)電子，別的軌道上電子跑到這個(gè)原子附近應(yīng)該也需要克服這個(gè)庫(kù)倫能啊。。。。 而且我做了實(shí)驗(yàn)，也發(fā)現(xiàn)，加U以后，同自旋方向上軌道的能級(jí)也會(huì)有更大的劈裂？ |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 vasp | 第一性原理和電化學(xué) | VASP算例 | VASP第一性原理計(jì)算 |
» 猜你喜歡
329求調(diào)劑
已經(jīng)有8人回復(fù)
-
070305高分子化學(xué)與物理 304分求調(diào)劑
已經(jīng)有8人回復(fù)
-
349求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑
已經(jīng)有7人回復(fù)
-
291求調(diào)劑
已經(jīng)有4人回復(fù)
-
0856材料化工調(diào)劑 總分330
已經(jīng)有14人回復(fù)
-
352分-085602-一志愿985
已經(jīng)有5人回復(fù)
-
各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩
已經(jīng)有8人回復(fù)
-
求調(diào)劑，一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分
已經(jīng)有4人回復(fù)
-
279求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- DFT+U中對(duì)U和庫(kù)倫關(guān)聯(lián)的的理解
已經(jīng)有11人回復(fù)
KalaShayminS
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 432 (碩士)
- 金幣: 6174.2
- 紅花: 56
- 帖子: 1088
- 在線: 379小時(shí)
- 蟲號(hào): 2408111
- 注冊(cè): 2013-04-08
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
jugengfans: 金幣+10, ★★★★★最佳答案, 謝謝您的耐心回答！學(xué)習(xí)了。也就是說，如果spin up某幾個(gè)軌道全占，spindown全空，這個(gè)時(shí)候即使加多大的U，spin up也不會(huì)發(fā)生任何的改變。對(duì)應(yīng)的spin down也不會(huì)發(fā)生變化。而只有在非整數(shù)占據(jù)的時(shí)候，U的引入先影響spin up能級(jí)和占據(jù)數(shù)，電子云密度變了有效勢(shì)也就變了，這個(gè)時(shí)候spin up 和spin dn都要重新自洽了。那就引出一個(gè)問題。是不是U加到一定程度，引起d軌道100%占據(jù)的時(shí)候，再增大U，對(duì)體系的splitting就沒有影響了呢。比如加到10eV已經(jīng)100%占據(jù)，加到100eV是不是gap還是這么大。 2014-09-03 23:57:57
jugengfans: 金幣+10, ★★★★★最佳答案, 謝謝您的耐心回答！學(xué)習(xí)了。也就是說，如果spin up某幾個(gè)軌道全占，spindown全空，這個(gè)時(shí)候即使加多大的U，spin up也不會(huì)發(fā)生任何的改變。對(duì)應(yīng)的spin down也不會(huì)發(fā)生變化。而只有在非整數(shù)占據(jù)的時(shí)候，U的引入先影響spin up能級(jí)和占據(jù)數(shù)，電子云密度變了有效勢(shì)也就變了，這個(gè)時(shí)候spin up 和spin dn都要重新自洽了。那就引出一個(gè)問題。是不是U加到一定程度，引起d軌道100%占據(jù)的時(shí)候，再增大U，對(duì)體系的splitting就沒有影響了呢。比如加到10eV已經(jīng)100%占據(jù)，加到100eV是不是gap還是這么大。 2014-09-03 23:57:57
|
我覺得LZ有個(gè)概念理解有些不清。哈密頓量只有一個(gè)，是多電子波函數(shù)的方程。所以不存在什么up的哈密頓量和down的哈密頓量。 KS的單電子方程是變分法引出的方程，不是哈密頓量的方程。但其中有一項(xiàng)是哈密頓量變分出來的。注意，Exc和U的變分是自洽求解的，它并未退化成真的單電子方程，從哈密頓量的形式說某條件下，此哈密頓量和U無關(guān)，不表明其變分無關(guān)。 然后看一下U有什么效果。 首先up和down在這里是可以任選的，我們假設(shè)占據(jù)數(shù)比較大的那個(gè)自旋叫up。 用Ueff=U-J這個(gè)模型的時(shí)候別忘了假設(shè)：DFT在整數(shù)占據(jù)下是完全正確的，我們才能用這個(gè)式子。 所以，在這個(gè)框架下，對(duì)于你假設(shè)的體系，如果是整數(shù)占據(jù)——DFT給出來的基態(tài)就是spinup全占，spindown全空——此時(shí)Ueff就不產(chǎn)生任何效果——因?yàn)檫@個(gè)框架已經(jīng)假定它是精確的了�。ㄔ赨，J框架下可以有效果） 但如果不是全占或全空，即使down全空，U也有作用。 因?yàn)椴皇侨既�，意味著up軌道（能帶）能量穿過Fermi面。這時(shí)，這個(gè)U可以改動(dòng)up的能級(jí)。up能級(jí)改動(dòng)——up占據(jù)數(shù)改變——電荷密度改變——KS有效勢(shì)改變——單電子KS方程改變。這時(shí)候就算down一直保持全空，其對(duì)應(yīng)解出來的能級(jí)也會(huì)變。 |
31樓2014-09-03 19:24:07
卡開發(fā)發(fā)
專家顧問 (著名寫手)
Ab Initio Amateur
-
專家經(jīng)驗(yàn): +224 - 1ST強(qiáng)帖: 2
- 應(yīng)助: 47 (小學(xué)生)
- 金幣: 3514.8
- 散金: 14
- 紅花: 223
- 帖子: 2477
- 在線: 1254小時(shí)
- 蟲號(hào): 1369423
- 注冊(cè): 2011-08-16
- 性別: GG
- 專業(yè): 金屬材料的磨損與磨蝕
- 管轄: 第一性原理
2樓2014-08-30 18:39:33
卡開發(fā)發(fā)
專家顧問 (著名寫手)
Ab Initio Amateur
-
專家經(jīng)驗(yàn): +224 - 1ST強(qiáng)帖: 2
- 應(yīng)助: 47 (小學(xué)生)
- 金幣: 3514.8
- 散金: 14
- 紅花: 223
- 帖子: 2477
- 在線: 1254小時(shí)
- 蟲號(hào): 1369423
- 注冊(cè): 2011-08-16
- 性別: GG
- 專業(yè): 金屬材料的磨損與磨蝕
- 管轄: 第一性原理
|
這個(gè)文章前半部分提到校正的構(gòu)造似乎建立在“整數(shù)占據(jù)”這個(gè)問題上，從公式推導(dǎo)的過程也很容易看出。 我認(rèn)為上述問題不沖突，你所謂的“傳統(tǒng)”的U是建立在忽略多中心積分的近似，如果完全些考慮，U實(shí)際上是Hartree-Fock-like的，先不管J，U的部分這樣也包含雙中心的部分，這樣鄰近的能級(jí)是有作用的，多中心積分可忽略是有條件的。如果是這樣，從鄰近能級(jí)拿個(gè)電子過來到已經(jīng)有電子的能級(jí)就肯定也就會(huì)有相應(yīng)的能量變化。具體問題印象中李正中老師的書上有些討論。 上面我也是自己的理解，有些東西學(xué)完也差不多忘了，也希望大牛們來指導(dǎo)討論。 |
» 本帖已獲得的紅花（最新10朵）
jugengfans7樓2014-08-31 00:16:43
卡開發(fā)發(fā)
專家顧問 (著名寫手)
Ab Initio Amateur
-
專家經(jīng)驗(yàn): +224 - 1ST強(qiáng)帖: 2
- 應(yīng)助: 47 (小學(xué)生)
- 金幣: 3514.8
- 散金: 14
- 紅花: 223
- 帖子: 2477
- 在線: 1254小時(shí)
- 蟲號(hào): 1369423
- 注冊(cè): 2011-08-16
- 性別: GG
- 專業(yè): 金屬材料的磨損與磨蝕
- 管轄: 第一性原理
|
這個(gè)理解有點(diǎn)問題，如果你推導(dǎo)過HF你就會(huì)看到，交換部分只體現(xiàn)在自旋平行的電子上，而關(guān)聯(lián)部分則是都存在的。前者主要是受到Pauli不相容原理限制，波函數(shù)反對(duì)稱化造成的，這個(gè)很顯然只存在于自旋平行的電子間，自旋相反的項(xiàng)因正交性抵消；后者則是采用平均場(chǎng)產(chǎn)生的缺陷。 對(duì)于自旋平行的電子，因?yàn)橐呀?jīng)存在Pauli排斥，所以上述相關(guān)的作用就沒那么強(qiáng)，但是自旋相反的電子這種效應(yīng)就不能被忽略。 這里從公式的表面看應(yīng)該是修正了平行電子之間的作用，趨勢(shì)上應(yīng)該和交換作用的方向一致，我沒辦法說這部分就是交換或是關(guān)聯(lián)造成的，但是導(dǎo)致的結(jié)果是電子整體的關(guān)聯(lián)都被修正。 |
17樓2014-09-02 15:47:26
KalaShayminS
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 432 (碩士)
- 金幣: 6174.2
- 紅花: 56
- 帖子: 1088
- 在線: 379小時(shí)
- 蟲號(hào): 2408111
- 注冊(cè): 2013-04-08
【答案】應(yīng)助回帖
|
這里理解有問題，U雖然看起來是作用于同自旋的，但其實(shí)U本身就能使兩個(gè)自旋之間的能差加大，這一步還不用考慮DFT。原因是里面n*n的項(xiàng)在電子占滿同自旋軌道的時(shí)候能量最小——如果電子正好占滿一個(gè)自旋，那么另一個(gè)自旋只能在費(fèi)米面以上，這就拉大了自旋之間的能量差。 舉例：如果是球?qū)ΨQ的原子，那么公式只有對(duì)角元保留， 修正值為(U-J)/2 *Σ_σ (Σ_m nσ - Σ_m nσ*nσ ) 假設(shè)除了一條軌道以外其他軌道都是空的，那么就是(U-J)/2 *Σ_σ (nσ - nσ*nσ )。 如果是兩個(gè)自旋完全極化，能量就是0，其中 Σ_σ nσ*nσ = 0*0+1*1=1 如果電子被平均分在兩個(gè)自旋上，意味著每個(gè)自旋上是0.5個(gè)電子，那么Σ_σ nσ*nσ = 0.5*0.5+0.5*0.5 = 0.5，此時(shí)能量變成(U-J)/4，是個(gè)正數(shù)。 |
18樓2014-09-03 00:12:00
卡開發(fā)發(fā)
專家顧問 (著名寫手)
Ab Initio Amateur
-
專家經(jīng)驗(yàn): +224 - 1ST強(qiáng)帖: 2
- 應(yīng)助: 47 (小學(xué)生)
- 金幣: 3514.8
- 散金: 14
- 紅花: 223
- 帖子: 2477
- 在線: 1254小時(shí)
- 蟲號(hào): 1369423
- 注冊(cè): 2011-08-16
- 性別: GG
- 專業(yè): 金屬材料的磨損與磨蝕
- 管轄: 第一性原理
20樓2014-09-03 06:44:35
卡開發(fā)發(fā)
專家顧問 (著名寫手)
Ab Initio Amateur
-
專家經(jīng)驗(yàn): +224 - 1ST強(qiáng)帖: 2
- 應(yīng)助: 47 (小學(xué)生)
- 金幣: 3514.8
- 散金: 14
- 紅花: 223
- 帖子: 2477
- 在線: 1254小時(shí)
- 蟲號(hào): 1369423
- 注冊(cè): 2011-08-16
- 性別: GG
- 專業(yè): 金屬材料的磨損與磨蝕
- 管轄: 第一性原理
|
總的Hamilton肯定有變化啊，算是Exc比較合適吧，畢竟這也算是采用平均場(chǎng)過后產(chǎn)生差異的部分，并且也是動(dòng)能、平均排斥和核對(duì)電子吸引之外的作用。這個(gè)問題類似于這樣，如果把PBE的交換部分的25%換成HF，計(jì)算得到的能級(jí)和原來的PBE會(huì)一樣嗎？你再好好想想。 |
» 本帖已獲得的紅花（最新10朵）
30樓2014-09-03 15:29:07
provels
金蟲 (正式寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 1534.9
- 散金: 81
- 紅花: 1
- 帖子: 448
- 在線: 216.7小時(shí)
- 蟲號(hào): 639531
- 注冊(cè): 2008-10-28
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
【答案】應(yīng)助回帖
★ ★
jugengfans: 金幣+2, 據(jù)說castep里的+U就是采用你說的這個(gè)簡(jiǎn)化模型的，而沒有其它選項(xiàng)，是不是說明castep的+U是非常粗略的？并不可信？本來一個(gè)計(jì)算任務(wù)打算用castep算的。。 2014-09-14 21:44:18
jugengfans: 金幣+2, 據(jù)說castep里的+U就是采用你說的這個(gè)簡(jiǎn)化模型的，而沒有其它選項(xiàng)，是不是說明castep的+U是非常粗略的？并不可信？本來一個(gè)計(jì)算任務(wù)打算用castep算的。。 2014-09-14 21:44:18
|
anistropic inteaction 是J 控制的 在HATREE， EXCHANGE 項(xiàng)都是 那么VASP 里面幾個(gè)模型有簡(jiǎn)化的 U/EFF之類就應(yīng)該沒有這個(gè)是近MEAN FIELD APPROACH, TYPE=1的 這里 HUND\S RULE 就只考慮了第一條 所以要考慮的話 選好模型 ！ 不能選1 |
34樓2014-09-04 04:46:06
卡開發(fā)發(fā)
專家顧問 (著名寫手)
Ab Initio Amateur
-
專家經(jīng)驗(yàn): +224 - 1ST強(qiáng)帖: 2
- 應(yīng)助: 47 (小學(xué)生)
- 金幣: 3514.8
- 散金: 14
- 紅花: 223
- 帖子: 2477
- 在線: 1254小時(shí)
- 蟲號(hào): 1369423
- 注冊(cè): 2011-08-16
- 性別: GG
- 專業(yè): 金屬材料的磨損與磨蝕
- 管轄: 第一性原理
37樓2014-09-09 13:49:32
jugengfans
金蟲 (正式寫手)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 773.9
- 散金: 151
- 紅花: 3
- 帖子: 537
- 在線: 205.6小時(shí)
- 蟲號(hào): 291609
- 注冊(cè): 2006-10-29
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2014-08-30 19:31:41
卡開發(fā)發(fā)
專家顧問 (著名寫手)
Ab Initio Amateur
-
專家經(jīng)驗(yàn): +224 - 1ST強(qiáng)帖: 2
- 應(yīng)助: 47 (小學(xué)生)
- 金幣: 3514.8
- 散金: 14
- 紅花: 223
- 帖子: 2477
- 在線: 1254小時(shí)
- 蟲號(hào): 1369423
- 注冊(cè): 2011-08-16
- 性別: GG
- 專業(yè): 金屬材料的磨損與磨蝕
- 管轄: 第一性原理
4樓2014-08-30 19:41:54
jugengfans
金蟲 (正式寫手)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 773.9
- 散金: 151
- 紅花: 3
- 帖子: 537
- 在線: 205.6小時(shí)
- 蟲號(hào): 291609
- 注冊(cè): 2006-10-29
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
5樓2014-08-30 23:04:28
jugengfans
金蟲 (正式寫手)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 773.9
- 散金: 151
- 紅花: 3
- 帖子: 537
- 在線: 205.6小時(shí)
- 蟲號(hào): 291609
- 注冊(cè): 2006-10-29
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
|
U的準(zhǔn)確理解，國(guó)內(nèi)的文獻(xiàn)中總是說什么同軌道占據(jù)一個(gè)自旋相反的電子， 但是我看Phys. Rev. B 57, 1505 (1998)里面似乎說， Since the strength of the effective on-site Coulomb interaction between d electrons ~characterized by the Hubbard U) is comparable with the valence bandwidth, the processes associated with electron transfer between two metal ions or resulting from addition or removal of d electrons give rise to large fluctuations of the energy of the system, leading to the localization of carriers and to the formation of band gaps.17 他的意思似乎是說，在兩個(gè)金屬離子之間電子的傳輸?shù)碾娮拥倪^程，或者換句話說，從一個(gè)離子旁邊拿走或加上一個(gè)d電子，會(huì)對(duì)體系的能量造成很大的影響。 他沒有說在同一個(gè)d軌道上，是不是就是說d軌道占據(jù)一個(gè)電子，別的電子占據(jù)其余四個(gè)簡(jiǎn)并軌道之一，也會(huì)有庫(kù)倫排斥？而這個(gè)排斥已經(jīng)考慮在它的模型中了。貌似是這個(gè)意思。 也不知道我的理解是不是有誤。高量一些東西忘光光了。 |
6樓2014-08-30 23:14:07
jugengfans
金蟲 (正式寫手)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 773.9
- 散金: 151
- 紅花: 3
- 帖子: 537
- 在線: 205.6小時(shí)
- 蟲號(hào): 291609
- 注冊(cè): 2006-10-29
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
8樓2014-08-31 00:29:00
KalaShayminS
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 432 (碩士)
- 金幣: 6174.2
- 紅花: 56
- 帖子: 1088
- 在線: 379小時(shí)
- 蟲號(hào): 2408111
- 注冊(cè): 2013-04-08
【答案】應(yīng)助回帖
★ ★
感謝參與，應(yīng)助指數(shù) +1
jugengfans: 金幣+2 2014-08-31 12:40:33
liliangfang: 謝謝交流 2014-08-31 18:33:48
感謝參與，應(yīng)助指數(shù) +1
jugengfans: 金幣+2 2014-08-31 12:40:33
liliangfang: 謝謝交流 2014-08-31 18:33:48
|
single band Hubbard model沒有不同軌道之間的作用（只有一條軌道），如果是multi band的話會(huì)引入各軌道之間的作用。LDA+U屬于后者，例如d軌道5個(gè)分量之間都有作用，U是個(gè)矩陣。 另外，Hubbard model里面U關(guān)聯(lián)的是只有兩中心的項(xiàng)，其他項(xiàng)忽略了，因此在球形近似下U和J就足夠表示整個(gè)矩陣了。 進(jìn)一步近似DFT在整數(shù)占據(jù)下是正確的，那么只有整數(shù)-非整數(shù)的差需要考慮，這時(shí)就只需要一個(gè)數(shù)U-J就夠了。 |
9樓2014-08-31 06:14:14
jugengfans
金蟲 (正式寫手)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 773.9
- 散金: 151
- 紅花: 3
- 帖子: 537
- 在線: 205.6小時(shí)
- 蟲號(hào): 291609
- 注冊(cè): 2006-10-29
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
10樓2014-08-31 12:44:33
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 jugengfans 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 352分-085602-一志愿985 +5 | 海納百川Ly 2026-03-29 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 本科雙非材料，跨考一志愿華電085801電氣，283求調(diào)劑，任何專業(yè)都可以 +6 | 芝士雪baoo 2026-03-28 | 8/400 |
|
|
|
[考研] 2026年華南師范大學(xué)歡迎化學(xué)，化工，生物，生醫(yī)工等專業(yè)優(yōu)秀學(xué)子加入！ +3 | llss0711 2026-03-28 | 5/250 |
|
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
|
[考研] 071000生物學(xué)求調(diào)劑，初試成績(jī)343 +7 | 小小甜面團(tuán) 2026-03-25 | 7/350 |
|
|
|
[考研] 一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑 +7 | Micky11223 2026-03-25 | 7/350 |
|
|
|
[考研]
 一志愿廈門大學(xué)化學(xué)學(xué)碩307求調(diào)劑
+10
|
y7czhao 2026-03-26 | 10/500 |
|
|
|
[考研] 286求調(diào)劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
|
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分，英語(yǔ)已過四六級(jí)，六級(jí)521，求調(diào)劑 +4 | 美味可樂雞翅 2026-03-26 | 4/200 |
|
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
|
[考研] 292求調(diào)劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
|
[考研] 【雙一流院校新能源、環(huán)境材料，材料加工與模擬招收大量調(diào)劑】 +4 | Higraduate 2026-03-22 | 8/400 |
|
|
|
[考研]
 085602化學(xué)工程求調(diào)劑。
+4
|
平樂樂樂 2026-03-26 | 4/200 |
|
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
|
[考研] 281求調(diào)劑 +6 | Koxui 2026-03-24 | 7/350 |
|
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|
|
|
[考研] 0854AI CV方向招收調(diào)劑 +4 | 章小魚567 2026-03-23 | 4/200 |
|
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
|
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4 | 為學(xué)666 2026-03-24 | 4/200 |
|
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|