| 5 | 1/1 | 返回列表 |
| 查看: 2268 | 回復(fù): 8 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
[交流]
【求助】為什么要單獨計算態(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/285
+1/236
+2/154
+2/66
+2/46
+2/36
+1/33
+1/31
+1/28
+1/22
+1/16
+1/15
+1/11
+1/10
+1/7
+1/5
+1/5
+1/5
+1/2
+1/2
|
只有在平衡條件下計算的物理量才是可信的,所以這里的要求優(yōu)化晶格參數(shù)。當(dāng)然,如果你知道了平衡時的晶格參數(shù),就可以跳過優(yōu)化離子那一步,而只進行靜態(tài)的自洽計算和非自恰的計算。 其實優(yōu)化的最后一步也計算了自洽的電荷密度。但不是很準(zhǔn)確的,因為我們優(yōu)化時一般設(shè)置ISMEAR=0或者1,而要得到準(zhǔn)確的自洽的電荷密度,要增加K點網(wǎng)格和改進K點抽樣方法,也就是增加K點和設(shè)置ISMEAR=-5,所以優(yōu)化晶格參數(shù)后要做一次自洽靜態(tài)計算。 |
榮譽版主 (知名作家)
|
晶格優(yōu)化的時候態(tài)密度和費米能級都不是準(zhǔn)確的。只有靜態(tài)自洽的時候才是準(zhǔn)確的。 靜態(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)密度和費米能級不夠準(zhǔn)確。 在優(yōu)化后的靜態(tài)自洽計算中,需要增加K點網(wǎng)絡(luò),從而增加電子密度、態(tài)密度、費米能級的準(zhǔn)確性。 之后才進行非自洽的能帶結(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)化的時候如果我選用第一原理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 ] |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +5 | Zhangbod 2026-03-21 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 269專碩求調(diào)劑 +5 | 金恩貝 2026-03-21 | 5/250 |
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +10 | 楓橋ZL 2026-03-18 | 12/600 |
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
[考研] 302求調(diào)劑 +12 | 呼呼呼。。。。 2026-03-17 | 12/600 |
|
|
[考研] 279分求調(diào)劑 一志愿211 +14 | chaojifeixia 2026-03-19 | 15/750 |
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3 | 研究生何瑤明 2026-03-18 | 3/150 |
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|
[考研] 311求調(diào)劑 +5 | 冬十三 2026-03-18 | 5/250 |
|
|
[考研] 290求調(diào)劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
[考研]
|
簡木ChuFront 2026-03-19 | 8/400 |
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 3/150 |
|
|
[考研] 319求調(diào)劑 +3 | 小力氣珂珂 2026-03-20 | 3/150 |
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
[考研] 281求調(diào)劑(0805) +14 | 煙汐憶海 2026-03-16 | 25/1250 |
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|