| 5 | 1/1 | 返回列表 |
| 查看: 2256 | 回復(fù): 8 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
Gina88木蟲 (正式寫手)
|
[交流]
【求助】為什么要單獨(dú)計(jì)算態(tài)密度呢?
|
||
|
我看到資料說Vasp計(jì)算的步驟為: 1。生成4個(gè)輸入文件 2。優(yōu)化晶格參數(shù) 3。固定晶格參數(shù),一個(gè)離子步自洽計(jì)算求出能態(tài)密度和費(fèi)米能級(jí) 4。固定電荷密度,做非自洽計(jì)算能量。 我想問的是: 優(yōu)化晶格參數(shù)的時(shí)候應(yīng)該計(jì)算費(fèi)米能級(jí)的吧。否則怎么算電子密度,怎么確定hamiltonian呢,怎么自洽計(jì)算呢? 那如果優(yōu)化的時(shí)候每一步都有計(jì)算費(fèi)米能級(jí)的話,費(fèi)米能級(jí)為什么不用優(yōu)化晶格參數(shù)最后一步的費(fèi)米能級(jí)呢(也就是第2步最后的費(fèi)米能級(jí))?為什么還要用步驟3來計(jì)算費(fèi)米能級(jí)呢? |
» 搶金幣啦!回帖就可以得到:
+1/468
+1/97
+1/85
+1/80
+1/55
+1/42
+1/36
+1/32
+1/25
+1/16
+1/14
+1/10
+1/8
+1/7
+1/5
+1/4
+1/4
+1/4
+1/1
+1/1
|
你好 “因?yàn)槲覀儍?yōu)化時(shí)一般設(shè)置ISMEAR=0或者1,而要得到準(zhǔn)確的自洽的電荷密度,要增加K點(diǎn)網(wǎng)格和改進(jìn)K點(diǎn)抽樣方法,也就是增加K點(diǎn)和設(shè)置ISMEAR=-5,所以優(yōu)化晶格參數(shù)后要做一次自洽靜態(tài)計(jì)算。" 意思是靜態(tài)自洽的時(shí)候已經(jīng)需要增加K點(diǎn)了嗎, ISMEAR=-5,那態(tài)密度的時(shí)候還需要繼續(xù)增加K嗎 |
榮譽(yù)版主 (知名作家)
|
晶格優(yōu)化的時(shí)候態(tài)密度和費(fèi)米能級(jí)都不是準(zhǔn)確的。只有靜態(tài)自洽的時(shí)候才是準(zhǔn)確的。 靜態(tài)自洽可以得到電荷密度,同時(shí)也可以得到態(tài)密度。 在這個(gè)基礎(chǔ)上做非自洽計(jì)算再來計(jì)算態(tài)密度,有一個(gè)條件:增加k點(diǎn)網(wǎng)格密度。從原理上講,計(jì)算態(tài)密度是否需要減小SIGMA或者設(shè)置ISMEAR=-5,這個(gè)我不清楚。 也可以直接在靜態(tài)自洽的時(shí)候計(jì)算態(tài)密度。 |
木蟲 (正式寫手)
|
首先非常感謝您的回答:) 可是我看見那個(gè)優(yōu)化的流程圖是每移動(dòng)一次離子,就要對(duì)電子進(jìn)行一次自洽計(jì)算,自洽后再算作用在離子上的力,如果力足夠小就不再移動(dòng)離子,否則就繼續(xù)移動(dòng)離子優(yōu)化晶格啊。 我的理解是:最后一次離子移動(dòng)后,程序是要進(jìn)行電子自洽計(jì)算的,在此基礎(chǔ)上再計(jì)算力,而計(jì)算的力足夠小,所以不再移動(dòng)離子。 照上面所說,最后一步其實(shí)是已經(jīng)進(jìn)行過電子自洽計(jì)算的。 上面您講的不知道我是不是可以這么理解: 離子優(yōu)化時(shí),雖然進(jìn)行了電子自洽計(jì)算,但是離子優(yōu)化時(shí)為了減小計(jì)算量,所以K點(diǎn)網(wǎng)絡(luò)不夠密,從而計(jì)算的態(tài)密度和費(fèi)米能級(jí)不夠準(zhǔn)確。 在優(yōu)化后的靜態(tài)自洽計(jì)算中,需要增加K點(diǎn)網(wǎng)絡(luò),從而增加電子密度、態(tài)密度、費(fèi)米能級(jí)的準(zhǔn)確性。 之后才進(jìn)行非自洽的能帶結(jié)構(gòu)計(jì)算? 哦,終于知道答案了。 原來基組增大對(duì)幾何優(yōu)化沒有幫助,所以在幾何優(yōu)化時(shí)沒有必要選擇很大的基組。 但是單點(diǎn)能計(jì)算的時(shí)候可以選擇大的基組,這樣可以增加精度,但是能量需要用Freq做零點(diǎn)能矯正。 優(yōu)化后加Freq計(jì)算單點(diǎn)能和ZPE,并不會(huì)出現(xiàn)虛頻。如果出現(xiàn)了虛頻,則說明優(yōu)化得不好,需要重新優(yōu)化。 優(yōu)化出現(xiàn)虛頻這類問題,原因可能是采用了部分優(yōu)化而不是全局優(yōu)化。部分優(yōu)化的結(jié)果計(jì)算ZPE是沒有意義的! [ Last edited by Gina88 on 2011-1-24 at 11:44 ] |
金蟲 (著名寫手)
|
只有在平衡條件下計(jì)算的物理量才是可信的,所以這里的要求優(yōu)化晶格參數(shù)。當(dāng)然,如果你知道了平衡時(shí)的晶格參數(shù),就可以跳過優(yōu)化離子那一步,而只進(jìn)行靜態(tài)的自洽計(jì)算和非自恰的計(jì)算。 其實(shí)優(yōu)化的最后一步也計(jì)算了自洽的電荷密度。但不是很準(zhǔn)確的,因?yàn)槲覀儍?yōu)化時(shí)一般設(shè)置ISMEAR=0或者1,而要得到準(zhǔn)確的自洽的電荷密度,要增加K點(diǎn)網(wǎng)格和改進(jìn)K點(diǎn)抽樣方法,也就是增加K點(diǎn)和設(shè)置ISMEAR=-5,所以優(yōu)化晶格參數(shù)后要做一次自洽靜態(tài)計(jì)算。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 354求調(diào)劑 +4 | Tyoumou 2026-03-18 | 7/350 |
|
|
[考研] 0703化學(xué)調(diào)劑 +7 | 妮妮ninicgb 2026-03-15 | 11/550 |
|
|
[考研] 0854可跨調(diào)劑,一作一項(xiàng)核心論文五項(xiàng)專利,省、國級(jí)證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | pupcoco 2026-03-17 | 7/350 |
|
|
[考研] 環(huán)境工程調(diào)劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
[考研] 265求調(diào)劑 +3 | 梁梁校校 2026-03-17 | 3/150 |
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
[基金申請(qǐng)] 被我言中:新模板不強(qiáng)調(diào)格式了,假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
[考研] 材料與化工求調(diào)劑 +6 | 為學(xué)666 2026-03-16 | 6/300 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 274求調(diào)劑 +5 | 時(shí)間點(diǎn) 2026-03-13 | 5/250 |
|
|
[考研] 283求調(diào)劑 +10 | 小樓。 2026-03-12 | 14/700 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|
|
[考研] 學(xué)碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生,2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|
|
[考研] 321求調(diào)劑(食品/專碩) +3 | xc321 2026-03-12 | 6/300 |
|