| 5 | 1/1 | 返回列表 |
| 查看: 2267 | 回復(fù): 8 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
Gina88木蟲 (正式寫手)
|
[交流]
【求助】為什么要單獨(dú)計(jì)算態(tài)密度呢?
|
||
|
我看到資料說(shuō)Vasp計(jì)算的步驟為: 1。生成4個(gè)輸入文件 2。優(yōu)化晶格參數(shù) 3。固定晶格參數(shù),一個(gè)離子步自洽計(jì)算求出能態(tài)密度和費(fèi)米能級(jí) 4。固定電荷密度,做非自洽計(jì)算能量。 我想問(wèn)的是: 優(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來(lái)計(jì)算費(fèi)米能級(jí)呢? |
» 搶金幣啦!回帖就可以得到:
+5/1775
+1/63
+1/36
+2/36
+1/32
+1/31
+1/23
+1/17
+2/16
+1/15
+1/10
+1/9
+1/6
+1/5
+1/3
+1/3
+1/3
+1/1
+1/1
+1/1
木蟲 (正式寫手)
|
首先非常感謝您的回答:) 可是我看見那個(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)離子。 照上面所說(shuō),最后一步其實(shí)是已經(jīng)進(jìn)行過(guò)電子自洽計(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ì)算? 哦,終于知道答案了。 原來(lái)基組增大對(duì)幾何優(yōu)化沒(méi)有幫助,所以在幾何優(yōu)化時(shí)沒(méi)有必要選擇很大的基組。 但是單點(diǎn)能計(jì)算的時(shí)候可以選擇大的基組,這樣可以增加精度,但是能量需要用Freq做零點(diǎn)能矯正。 優(yōu)化后加Freq計(jì)算單點(diǎn)能和ZPE,并不會(huì)出現(xiàn)虛頻。如果出現(xiàn)了虛頻,則說(shuō)明優(yōu)化得不好,需要重新優(yōu)化。 優(yōu)化出現(xiàn)虛頻這類問(wèn)題,原因可能是采用了部分優(yōu)化而不是全局優(yōu)化。部分優(yōu)化的結(jié)果計(jì)算ZPE是沒(méi)有意義的! [ Last edited by Gina88 on 2011-1-24 at 11:44 ] |
榮譽(yù)版主 (知名作家)
|
晶格優(yōu)化的時(shí)候態(tài)密度和費(fèi)米能級(jí)都不是準(zhǔn)確的。只有靜態(tài)自洽的時(shí)候才是準(zhǔn)確的。 靜態(tài)自洽可以得到電荷密度,同時(shí)也可以得到態(tài)密度。 在這個(gè)基礎(chǔ)上做非自洽計(jì)算再來(lái)計(jì)算態(tài)密度,有一個(gè)條件:增加k點(diǎn)網(wǎng)格密度。從原理上講,計(jì)算態(tài)密度是否需要減小SIGMA或者設(shè)置ISMEAR=-5,這個(gè)我不清楚。 也可以直接在靜態(tài)自洽的時(shí)候計(jì)算態(tài)密度。 |
金蟲 (著名寫手)
|
只有在平衡條件下計(jì)算的物理量才是可信的,所以這里的要求優(yōu)化晶格參數(shù)。當(dāng)然,如果你知道了平衡時(shí)的晶格參數(shù),就可以跳過(guò)優(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ì)算。 |
木蟲 (正式寫手)
|
非常感謝! 還想問(wèn)一個(gè)非常弱的問(wèn)題: 就是優(yōu)化的時(shí)候如果我選用第一原理MD優(yōu)化,這個(gè)優(yōu)化的精度需要選高一些呢還是選低一些呢?前面看見有人說(shuō)精度不影響優(yōu)化結(jié)果。 MD優(yōu)化后需不需要再選用其它的優(yōu)化方法再優(yōu)化一次呢? 我的理解是MD優(yōu)化需要加溫度啊,這樣可以找到勢(shì)能的全局最小值啊,而其它不加溫度的優(yōu)化只能找到勢(shì)能的局部最小值啊。但是MD優(yōu)化據(jù)說(shuō)很慢,這樣我可以先選MD優(yōu)化,并且精度低一些,找到全局最小值的那個(gè)勢(shì)阱,然后用其它的優(yōu)化方法并且精度高一些,在這個(gè)勢(shì)阱中找到勢(shì)能最小值? 還有,常?匆娙思疫M(jìn)行自旋極化計(jì)算,對(duì)于自選極化計(jì)算,那優(yōu)化的時(shí)候要不要選上自選極化? 如果優(yōu)化的時(shí)候不選自選極化,是不是后面計(jì)算的自旋極化性質(zhì)就不對(duì)啊? [ Last edited by Gina88 on 2011-1-30 at 16:42 ] |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +3 | vv迷 2026-03-21 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 考研調(diào)劑 +3 | 呼呼?~+123456 2026-03-21 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +11 | 妮妮ninicgb 2026-03-15 | 15/750 |
|
|
[考研] 326求調(diào)劑 +4 | mlpqaz03 2026-03-15 | 4/200 |
|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
[考研] 265求調(diào)劑 +9 | 梁梁校校 2026-03-17 | 9/450 |
|
|
[考研] 求調(diào)劑 +3 | Ma_xt 2026-03-17 | 3/150 |
|
|
[考研] 一志愿南昌大學(xué),327分,材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂(lè)29 2026-03-18 | 5/250 |
|
|
[考研] 323求調(diào)劑 +3 | 洼小桶 2026-03-18 | 3/150 |
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
[論文投稿]
申請(qǐng)回稿延期一個(gè)月,編輯同意了。但系統(tǒng)上的時(shí)間沒(méi)變,給編輯又寫郵件了,沒(méi)回復(fù)
10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
[考研] 一志愿,福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭(zhēng) 2026-03-15 | 3/150 |
|
|
[考研] 機(jī)械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|