| 4 | 1/1 | 返回列表 |
| 查看: 1210 | 回復(fù): 3 | |||
pgi1706木蟲 (小有名氣)
|
[交流]
晶格常數(shù)和能量是這種形式? 已有2人參與
|
|
暈了,算出來(lái)的晶格常數(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é)果,我想得到比較精確的能量。 原來(lái)采用單G點(diǎn)優(yōu)化時(shí),a=10.894時(shí)是最低能量,但是采用4*4*4后,能量卻是現(xiàn)在這個(gè)模樣? 可有人遇到這種情況,能否交流一下 ?謝謝 |
銀蟲 (小有名氣)
木蟲 (正式寫手)
|
說(shuō)兩點(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í),這也是和上面說(shuō)的收斂測(cè)試密切相關(guān)的。一般來(lái)說(shuō),能量收斂到meV的量級(jí),晶格常數(shù)收斂到0.01埃以下,就已經(jīng)是比較不錯(cuò)的精度了。從k點(diǎn)數(shù)目來(lái)看,你使用的參數(shù)并非精度極高的計(jì)算吧,如果a是原子單位的話,你所追究的這些0.001的差別可能已經(jīng)超出了你的計(jì)算參數(shù)所能達(dá)到的精度極限了,能量也是同理。因此你所說(shuō)的這些”能量此起彼伏“可能都是在誤差范圍之內(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 來(lái)計(jì)算的,難道是Ecut,之前1*1*1的晶胞優(yōu)化后的Ecut=400eV,2*2*2的超晶胞,還要優(yōu)化Ecut? |
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +10 | 楓橋ZL 2026-03-18 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +3 | 13341 2026-03-20 | 3/150 |
|
|
[考研] 一志愿深大,0703化學(xué),總分302,求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
[考研] 求助 +5 | 夢(mèng)里的無(wú)言 2026-03-21 | 6/300 |
|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
[考研] 265求調(diào)劑 +9 | 梁梁校校 2026-03-17 | 9/450 |
|
|
[考研] 307求調(diào)劑 +10 | 冷笙123 2026-03-17 | 10/500 |
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
[考研] 材料專業(yè)求調(diào)劑 +6 | hanamiko 2026-03-18 | 6/300 |
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8 | 小材化本科 2026-03-18 | 8/400 |
|
|
[考研] 323求調(diào)劑 +3 | 洼小桶 2026-03-18 | 3/150 |
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽(yáng)花天天開心 2026-03-16 | 4/200 |
|
|
[考研] 機(jī)械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|