| 2 | 1/1 | 返回列表 |
| 查看: 1440 | 回復: 1 | |||
[求助]
雜化泛函計算出能帶需要的費米能級應該取哪一個?
|
|
在雜化泛函計算能帶結(jié)構(gòu)的時候,通常有以下順序 針對雜化泛函特殊K點的PBE計算——〉 針對雜化泛函特殊K點的HSEscf計算,利用前一步的波函數(shù)(ISTART=1)——〉 針對雜化泛函特殊K點的HSEnoscf計算(ISTART=1,ICHARG=11) 經(jīng)過測試發(fā)現(xiàn),利用第二步或者第三步出的能帶結(jié)構(gòu)均可。。沒啥太大改變。。 但是。出現(xiàn)與PBE計算不同的是,比如純半導體的體系。在PBE計算出能帶結(jié)構(gòu)后,取PBE靜態(tài)自洽的或者DOS的費米能級,那么價帶頂正好是在能量零點處。。。 而用HSE方法計算能帶發(fā)現(xiàn)無論取HSEscf還是HSEnoscf步的費米能級都達不到這樣。。 后來卻發(fā)現(xiàn)一個很奇怪的現(xiàn)象,就是利用 普通K點(自動生成)PBE的電荷和波函數(shù),直接進行HSE的非自洽計算 (即ISTART=1,ICHARG=11),在前N步電子步都會出現(xiàn)不收斂的情況 entering main loop N E dE d eps ncg rms ort N E dE rms(c) ROT: 1 -0.507082286949E+02 -0.50708E+02 ROT: 2 -0.506938482460E+02 0.14380E-01 。。。。。。。 ROT: ** -0.680892037469E+02 0.88443E+01 ROT: ** -0.769335096047E+02 -0.88443E+01 ROT: ** -0.680892037470E+02 0.88443E+01 ROT: ** -0.769335096047E+02 -0.88443E+01 ROT: ** -0.680892037469E+02 0.88443E+01 ROT: ** -0.769335096047E+02 -0.88443E+01 ROT: ** -0.680892037469E+02 0.88443E+01 ROT: ** -0.769335096047E+02 -0.88443E+01 ROT: ** -0.680892037470E+02 0.88443E+01 ROT: ** -0.769335096047E+02 -0.88443E+01 ROT: ** -0.680892037470E+02 0.88443E+01 ROT: ** -0.769335096047E+02 -0.88443E+01 ROT: ** -0.680892037470E+02 0.88443E+01 ROT: ** -0.769335096047E+02 -0.88443E+01 ROT: ** -0.680892037470E+02 0.88443E+01 ROT: ** -0.769335096046E+02 -0.88443E+01 ROT: ** -0.680892037469E+02 0.88443E+01 ROT: ** -0.769335096047E+02 -0.88443E+01 ROT: ** -0.680892037469E+02 0.88443E+01 ROT: ** -0.769335096047E+02 -0.88443E+01 ROT: ** -0.680892037470E+02 0.88443E+01 gam= 0.000 g(H,U,f)= 0.949E+00-0.182E+01 0.157E-14 ort(H,U,f) = 0.000E+00 0.000E+00 0.000E+00 SDA: 1 -0.237059317507E+02 -0.23706E+02 0.34776E+00 864 -0.869E+00 0.000E+00 。。。。。。。 ROT: 1 -0.107404641718E+03 0.47400E-01 ROT: 2 -0.107404640662E+03 0.10556E-05 ROT: 3 -0.107404641256E+03 -0.59404E-06 gam= 0.382 g(H,U,f)= 0.108E-05 0.878E-05 0.111E-16 ort(H,U,f) = 0.141E-05 0.598E-05 0.827E-16 DMP: 19 -0.243777224012E+02 -0.40814E-05 -0.50766E-05 864 0.987E-05 0.739E-05 1 F= -.24377722E+02 E0= -.24377722E+02 d E =-.153305E-14 mag= 0.0000 writing wavefunctions 直到最后收斂,這樣得到的費米能級卻是合理的。。。但是覺得這樣的計算不合理啊 HSE應該先用波函數(shù)進行自洽的啊。。不知道大家有什么看法? |
雜化泛函方面的研究 |
木蟲 (正式寫手)
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 南京大學化學調(diào)劑 +9 | 景隨風 2026-03-29 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 一志愿北京工業(yè)大學,324分求調(diào)劑 +6 | 零八# 2026-03-28 | 6/300 |
|
|
[考研] 材料專碩調(diào)劑 +5 | 椰椰。 2026-03-29 | 5/250 |
|
|
[考研] 275求調(diào)劑 +15 | Micky11223 2026-03-25 | 20/1000 |
|
|
[考研] 本科雙非材料,跨考一志愿華電085801電氣,283求調(diào)劑,任何專業(yè)都可以 +6 | 芝士雪baoo 2026-03-28 | 8/400 |
|
|
[考研] 266分,求材料冶金能源化工等調(diào)劑 +7 | 哇呼哼呼哼 2026-03-27 | 9/450 |
|
|
[考研] 070300化學求調(diào)劑 +4 | 起個名咋這么難 2026-03-27 | 4/200 |
|
|
[考研] 化學調(diào)劑 +4 | 愛吃番茄的旭 2026-03-24 | 5/250 |
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
[碩博家園] 北京林業(yè)大學碩導招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
[考研] 325求調(diào)劑 +5 | 李嘉圖·S·路 2026-03-23 | 5/250 |
|
|
[考研] 325求調(diào)劑 +3 | Aoyijiang 2026-03-23 | 3/150 |
|
|
[考研]
材料調(diào)劑
5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|
|
[考研] 中國科學院深圳先進技術(shù)研究院-光纖傳感課題組招生-中國科學院大學、深圳理工大學聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
[考研] 一志愿哈工大,085400,320,求調(diào)劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
[考研] 293求調(diào)劑 +7 | 加一一九 2026-03-24 | 7/350 |
|
|
[考研] 化工專碩求調(diào)劑 +3 | question挽風 2026-03-24 | 3/150 |
|
|
[考研] 接收2026碩士調(diào)劑(學碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|