| 9 | 1/1 | 返回列表 |
| 查看: 2248 | 回復(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/920
+1/86
+1/81
+1/78
+1/69
+1/49
+1/36
+1/31
+1/21
+5/20
+1/18
+1/15
+1/8
+1/7
+1/7
+1/6
+1/2
+1/2
+1/2
+1/1
榮譽版主 (知名作家)
|
晶格優(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)化的流程圖是每移動一次離子,就要對電子進(jìn)行一次自洽計算,自洽后再算作用在離子上的力,如果力足夠小就不再移動離子,否則就繼續(xù)移動離子優(yōu)化晶格啊。 我的理解是:最后一次離子移動后,程序是要進(jìn)行電子自洽計算的,在此基礎(chǔ)上再計算力,而計算的力足夠小,所以不再移動離子。 照上面所說,最后一步其實是已經(jīng)進(jìn)行過電子自洽計算的。 上面您講的不知道我是不是可以這么理解: 離子優(yōu)化時,雖然進(jìn)行了電子自洽計算,但是離子優(yōu)化時為了減小計算量,所以K點網(wǎng)絡(luò)不夠密,從而計算的態(tài)密度和費米能級不夠準(zhǔn)確。 在優(yōu)化后的靜態(tài)自洽計算中,需要增加K點網(wǎng)絡(luò),從而增加電子密度、態(tài)密度、費米能級的準(zhǔn)確性。 之后才進(jì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)化晶格參數(shù)。當(dāng)然,如果你知道了平衡時的晶格參數(shù),就可以跳過優(yōu)化離子那一步,而只進(jìn)行靜態(tài)的自洽計算和非自恰的計算。 其實優(yōu)化的最后一步也計算了自洽的電荷密度。但不是很準(zhǔn)確的,因為我們優(yōu)化時一般設(shè)置ISMEAR=0或者1,而要得到準(zhǔn)確的自洽的電荷密度,要增加K點網(wǎng)格和改進(jìn)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(jìn)行自旋極化計算,對于自選極化計算,那優(yōu)化的時候要不要選上自選極化? 如果優(yōu)化的時候不選自選極化,是不是后面計算的自旋極化性質(zhì)就不對啊? [ Last edited by Gina88 on 2011-1-30 at 16:42 ] |
|
本帖內(nèi)容被屏蔽 |
金蟲 (著名寫手)
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 344求調(diào)劑 +4 | knight344 2026-03-16 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[基金申請] 被我言中:新模板不強(qiáng)調(diào)格式了,假專家開始管格式了 +3 | beefly 2026-03-14 | 3/150 |
|
|
[考研] 一志愿蘇州大學(xué)材料工程(085601)專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[考研] 085601專碩,總分342求調(diào)劑,地區(qū)不限 +3 | share_joy 2026-03-16 | 3/150 |
|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +6 | SUICHILD 2026-03-12 | 6/300 |
|
|
[基金申請] 國自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
[考研] 085600調(diào)劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
[考研] 一志愿華中師范071000,325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化,一志愿:哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
[考研] 【0703化學(xué)調(diào)劑】-一志愿華中師范大學(xué)-六級475 +5 | Becho359 2026-03-11 | 5/250 |
|
|
[考研] 材料371求調(diào)劑 +9 | 鱷魚? 2026-03-11 | 11/550 |
|
|
[考研] 329求調(diào)劑 +3 | miaodesi 2026-03-12 | 4/200 |
|
|
[考研] 求調(diào)劑 +7 | 18880831720 2026-03-11 | 7/350 |
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|
|
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生,2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-12 | 4/200 |
|