| 4 | 1/1 | 返回列表 |
| 查看: 1209 | 回復(fù): 3 | |||
pgi1706木蟲 (小有名氣)
|
[交流]
晶格常數(shù)和能量是這種形式? 已有2人參與
|
|
暈了,算出來的晶格常數(shù)和能量是這種形式? a= 10.891 1 F= -.22064566E+03 E0= -.22064566E+03 d E =-.220646E+03 a= 10.892 1 F= -.22064519E+03 E0= -.22064519E+03 d E =-.220645E+03 a= 10.893 1 F= -.22064626E+03 E0= -.22064626E+03 d E =-.220646E+03 a= 10.894 1 F= -.22064570E+03 E0= -.22064570E+03 d E =-.220646E+03 a= 10.895 1 F= -.22064509E+03 E0= -.22064509E+03 d E =-.220645E+03 a= 10.896 1 F= -.22064517E+03 E0= -.22064517E+03 d E =-.220645E+03 a= 10.897 1 F= -.22064448E+03 E0= -.22064448E+03 d E =-.220644E+03 a= 10.898 1 F= -.22064530E+03 E0= -.22064530E+03 d E =-.220645E+03 a= 10.899 1 F= -.22064452E+03 E0= -.22064452E+03 d E =-.220645E+03 a= 10.900 1 F= -.22064591E+03 E0= -.22064591E+03 d E =-.220646E+03 所謂能量此起彼伏,何為最優(yōu)結(jié)構(gòu)? 這是采用4*4*4K點(diǎn)優(yōu)化的結(jié)果,我想得到比較精確的能量。 原來采用單G點(diǎn)優(yōu)化時(shí),a=10.894時(shí)是最低能量,但是采用4*4*4后,能量卻是現(xiàn)在這個(gè)模樣? 可有人遇到這種情況,能否交流一下 ?謝謝 |
銀蟲 (小有名氣)
木蟲 (正式寫手)
|
說兩點(diǎn): 1,k點(diǎn)可能不夠,要做收斂測(cè)試,也就是不斷的增加k點(diǎn)數(shù)目進(jìn)行晶格常數(shù)的計(jì)算,得到一系列的晶格常數(shù)值,當(dāng)k點(diǎn)足夠多時(shí),直到得到的晶格常數(shù)趨近穩(wěn)定,這個(gè)時(shí)候的k數(shù)目才是足夠的。同樣,也應(yīng)該對(duì)Ecut做類似的測(cè)試。 2,也是更重要的,要對(duì)自己計(jì)算精度有一定的認(rèn)識(shí),這也是和上面說的收斂測(cè)試密切相關(guān)的。一般來說,能量收斂到meV的量級(jí),晶格常數(shù)收斂到0.01埃以下,就已經(jīng)是比較不錯(cuò)的精度了。從k點(diǎn)數(shù)目來看,你使用的參數(shù)并非精度極高的計(jì)算吧,如果a是原子單位的話,你所追究的這些0.001的差別可能已經(jīng)超出了你的計(jì)算參數(shù)所能達(dá)到的精度極限了,能量也是同理。因此你所說的這些”能量此起彼伏“可能都是在誤差范圍之內(nèi)的。 |
木蟲 (小有名氣)
|
K點(diǎn)我已經(jīng)優(yōu)化了,如下 1 Found 1 irreducible k-points: free energy TOTEN = -219.066683 eV 2 Found 1 irreducible k-points: free energy TOTEN = -220.632965 eV 3 Found 4 irreducible k-points: free energy TOTEN = -220.644897 eV 4 Found 4 irreducible k-points: free energy TOTEN = -220.645697 eV 5 Found 10 irreducible k-points: free energy TOTEN = -220.645609 eV 所以我采用4*4*4的MP-Kpoints 來計(jì)算的,難道是Ecut,之前1*1*1的晶胞優(yōu)化后的Ecut=400eV,2*2*2的超晶胞,還要優(yōu)化Ecut? |
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 278求調(diào)劑 +7 | 煙火先于春 2026-03-17 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +6 | 墨墨漠 2026-03-20 | 6/300 |
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
[考研] 346求調(diào)劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
[考研] 307求調(diào)劑 +10 | 冷笙123 2026-03-17 | 10/500 |
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8 | 小材化本科 2026-03-18 | 8/400 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 6/300 |
|
|
[考研] 323求調(diào)劑 +3 | 洼小桶 2026-03-18 | 3/150 |
|
|
[考研] 085600材料與化工 +8 | 安全上岸! 2026-03-16 | 8/400 |
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生(接收計(jì)算機(jī)專業(yè)調(diào)劑) +3 | 123阿標(biāo) 2026-03-17 | 3/150 |
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15 | 7/350 |
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷,論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|