| 9 | 1/1 | 返回列表 |
| 查看: 2250 | 回復(fù): 8 | |||
[交流]
【求助】為什么要單獨計算態(tài)密度呢?
|
|
我看到資料說Vasp計算的步驟為: 1。生成4個輸入文件 2。優(yōu)化晶格參數(shù) 3。固定晶格參數(shù),一個離子步自洽計算求出能態(tài)密度和費米能級 4。固定電荷密度,做非自洽計算能量。 我想問的是: 優(yōu)化晶格參數(shù)的時候應(yīng)該計算費米能級的吧。否則怎么算電子密度,怎么確定hamiltonian呢,怎么自洽計算呢? 那如果優(yōu)化的時候每一步都有計算費米能級的話,費米能級為什么不用優(yōu)化晶格參數(shù)最后一步的費米能級呢(也就是第2步最后的費米能級)?為什么還要用步驟3來計算費米能級呢? |
» 搶金幣啦!回帖就可以得到:
+5/400
+1/83
+1/71
+1/67
+2/54
+1/38
+1/35
+1/28
+1/27
+1/16
+1/15
+1/7
+1/7
+1/6
+1/5
+1/5
+1/4
+1/2
+1/1
+1/1
榮譽版主 (知名作家)
|
晶格優(yōu)化的時候態(tài)密度和費米能級都不是準確的。只有靜態(tài)自洽的時候才是準確的。 靜態(tài)自洽可以得到電荷密度,同時也可以得到態(tài)密度。 在這個基礎(chǔ)上做非自洽計算再來計算態(tài)密度,有一個條件:增加k點網(wǎng)格密度。從原理上講,計算態(tài)密度是否需要減小SIGMA或者設(shè)置ISMEAR=-5,這個我不清楚。 也可以直接在靜態(tài)自洽的時候計算態(tài)密度。 |
|
首先非常感謝您的回答:) 可是我看見那個優(yōu)化的流程圖是每移動一次離子,就要對電子進行一次自洽計算,自洽后再算作用在離子上的力,如果力足夠小就不再移動離子,否則就繼續(xù)移動離子優(yōu)化晶格啊。 我的理解是:最后一次離子移動后,程序是要進行電子自洽計算的,在此基礎(chǔ)上再計算力,而計算的力足夠小,所以不再移動離子。 照上面所說,最后一步其實是已經(jīng)進行過電子自洽計算的。 上面您講的不知道我是不是可以這么理解: 離子優(yōu)化時,雖然進行了電子自洽計算,但是離子優(yōu)化時為了減小計算量,所以K點網(wǎng)絡(luò)不夠密,從而計算的態(tài)密度和費米能級不夠準確。 在優(yōu)化后的靜態(tài)自洽計算中,需要增加K點網(wǎng)絡(luò),從而增加電子密度、態(tài)密度、費米能級的準確性。 之后才進行非自洽的能帶結(jié)構(gòu)計算? 哦,終于知道答案了。 原來基組增大對幾何優(yōu)化沒有幫助,所以在幾何優(yōu)化時沒有必要選擇很大的基組。 但是單點能計算的時候可以選擇大的基組,這樣可以增加精度,但是能量需要用Freq做零點能矯正。 優(yōu)化后加Freq計算單點能和ZPE,并不會出現(xiàn)虛頻。如果出現(xiàn)了虛頻,則說明優(yōu)化得不好,需要重新優(yōu)化。 優(yōu)化出現(xiàn)虛頻這類問題,原因可能是采用了部分優(yōu)化而不是全局優(yōu)化。部分優(yōu)化的結(jié)果計算ZPE是沒有意義的! [ Last edited by Gina88 on 2011-1-24 at 11:44 ] |
|
只有在平衡條件下計算的物理量才是可信的,所以這里的要求優(yōu)化晶格參數(shù)。當然,如果你知道了平衡時的晶格參數(shù),就可以跳過優(yōu)化離子那一步,而只進行靜態(tài)的自洽計算和非自恰的計算。 其實優(yōu)化的最后一步也計算了自洽的電荷密度。但不是很準確的,因為我們優(yōu)化時一般設(shè)置ISMEAR=0或者1,而要得到準確的自洽的電荷密度,要增加K點網(wǎng)格和改進K點抽樣方法,也就是增加K點和設(shè)置ISMEAR=-5,所以優(yōu)化晶格參數(shù)后要做一次自洽靜態(tài)計算。 |
|
非常感謝! 還想問一個非常弱的問題: 就是優(yōu)化的時候如果我選用第一原理MD優(yōu)化,這個優(yōu)化的精度需要選高一些呢還是選低一些呢?前面看見有人說精度不影響優(yōu)化結(jié)果。 MD優(yōu)化后需不需要再選用其它的優(yōu)化方法再優(yōu)化一次呢? 我的理解是MD優(yōu)化需要加溫度啊,這樣可以找到勢能的全局最小值啊,而其它不加溫度的優(yōu)化只能找到勢能的局部最小值啊。但是MD優(yōu)化據(jù)說很慢,這樣我可以先選MD優(yōu)化,并且精度低一些,找到全局最小值的那個勢阱,然后用其它的優(yōu)化方法并且精度高一些,在這個勢阱中找到勢能最小值? 還有,常?匆娙思疫M行自旋極化計算,對于自選極化計算,那優(yōu)化的時候要不要選上自選極化? 如果優(yōu)化的時候不選自選極化,是不是后面計算的自旋極化性質(zhì)就不對啊? [ Last edited by Gina88 on 2011-1-30 at 16:42 ] |
|
本帖內(nèi)容被屏蔽 |
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +11 | 大米飯! 2026-03-15 | 14/700 |
|
|
[考研] 311求調(diào)劑 +8 | 冬十三 2026-03-15 | 8/400 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +3 | 困于星晨 2026-03-17 | 3/150 |
|
|
[考研] 289求調(diào)劑 +6 | 步川酷紫123 2026-03-11 | 6/300 |
|
|
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +5 | 薛云鵬 2026-03-13 | 5/250 |
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化,一志愿:哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
[考研] 331求調(diào)劑(0703有機化學(xué) +5 | ZY-05 2026-03-13 | 6/300 |
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
[考研] 招收0805(材料)調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|
|
[考研] 279求調(diào)劑 +3 | Dizzy123@ 2026-03-10 | 3/150 |
|
|
[考研] 0703化學(xué)一志愿211 總分320求調(diào)劑 +5 | 瑪卡巴卡啊哈 2026-03-11 | 5/250 |
|
|
[考研] 工科,求調(diào)劑 +3 | 我887 2026-03-11 | 3/150 |
|
|
[考研] 材料與化工085600調(diào)劑求老師收留 +9 | jiaanl 2026-03-11 | 9/450 |
|
|
[考研] 工科278分求調(diào)劑 +5 | 周慢熱啊 2026-03-12 | 7/350 |
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|