| 10 | 1/1 | 返回列表 |
| 查看: 3350 | 回復(fù): 9 | ||
gynan2008金蟲 (正式寫手)
|
[求助]
鈦酸鋇晶格參數(shù)哪個(gè)正確?
|
|
文獻(xiàn)里面有看到立方相晶格參數(shù)是3.95*****,前幾天用abinit計(jì)算時(shí),別人給的常數(shù)是7.901******,我用不同的acell計(jì)算總能,也是七點(diǎn)幾左右。 兩個(gè)問題,1哪個(gè)正確? 2、我用不同acel去計(jì)算總能,是不是就是在優(yōu)化晶格參數(shù)? 3.nband的設(shè)置我還是不明白,不是說價(jià)電子數(shù)/2么,可是設(shè)置的要比物質(zhì)的價(jià)電子數(shù)多很多啊。 |

金蟲 (正式寫手)

鐵桿木蟲 (著名寫手)
金蟲 (正式寫手)
|
#Definition of the unit cell ndtset 51 acell: 3.5 3.5 3.5 acell+ 0.1 0.1 0.1 #Definition of the atom types ntypat 3 #定義原胞中原子的類別的數(shù)目,本例子中只有1類原子 znucl 56 22 8 #定義原胞中原子的核電荷數(shù) natom 5 #定義原胞中原子的總個(gè)數(shù),本例子中有2個(gè)原子 typat 1 2 3 3 3 #定義每類原子的個(gè)數(shù),本例子中第一類原子有2個(gè) #下面定義了原胞中原子的坐標(biāo) xred 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00 5.0000000000E-01 5.0000000000E-01 5.0000000000E-01 5.0000000000E-01 5.0000000000E-01 0.0000000000E+00 5.0000000000E-01 0.0000000000E+00 5.0000000000E-01 0.0000000000E+00 5.0000000000E-01 5.0000000000E-01 #Gives the number of band, explicitely (do not take the default) nband 35 #定義了要計(jì)算的能帶的數(shù)目,最好按這樣來設(shè)置: # nband > 原胞中總的價(jià)電子數(shù)目/2 + 10 #Exchange-correlation functional ixc 1 #定義交換關(guān)聯(lián)函數(shù),本例子中,采用的是Teter Pade 參數(shù)化的LDA形式 #Definition of the planewave basis set ecut 10.0 #定義了平面波的切斷動(dòng)能 #Definition of the k-point grid ngkpt 8 8 8 #下面定義了k點(diǎn)網(wǎng)格取樣的大小 nshiftk 1 shiftk 0.5 0.5 0.5 #Definition of the SCF procedure nstep 30 #電子自洽迭代的最大步數(shù) diemac 12.0 #介電常數(shù)設(shè)置 toldfe 1.0d-10 #電子自洽收斂的標(biāo)準(zhǔn) enunit 2 rprim: 1 0 0 #下面三行定義了原胞的基矢,本例子中Si是fcc結(jié)構(gòu) 0 1 0 0 0 1 我算了一下,結(jié)果是下面的樣子,參數(shù)設(shè)置有問題么? 沒辦法傳圖片了,是一條收斂的曲線 |

鐵桿木蟲 (正式寫手)
西部牛仔
|
3.9*和7.9*都對(duì)的,只是單位不同而已。一個(gè)是埃,一個(gè)是波爾。 用不同的acel計(jì)算總能,是晶胞優(yōu)化的一種方法,不過效率太低而且太復(fù)雜了。但是,你這個(gè)代碼的功能是體積膨脹,不是結(jié)構(gòu)優(yōu)化! nband的設(shè)置按章規(guī)則來就行了,不懂樓主“設(shè)置的要比物質(zhì)的價(jià)電子數(shù)多很多”啥意思?對(duì)于絕緣體,完全沒必要加10個(gè)導(dǎo)帶能級(jí),那是心理安慰而已。 建議樓主先好好看一下幫助文檔~ |

木蟲 (小有名氣)
金蟲 (正式寫手)

金蟲 (正式寫手)

鐵桿木蟲 (正式寫手)
西部牛仔

銀蟲 (初入文壇)
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 344求調(diào)劑 +5 | knight344 2026-03-16 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +9 | yangfz 2026-03-17 | 9/450 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考研] 299求調(diào)劑 +5 | △小透明* 2026-03-17 | 5/250 |
|
|
[考博] 環(huán)境領(lǐng)域全國重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 268求調(diào)劑 +6 | 簡單點(diǎn)0 2026-03-17 | 6/300 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|
[考研] 一志愿蘇州大學(xué)材料工程(085601)專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨; 2026-03-16 | 3/150 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[基金申請(qǐng)] 國自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
[基金申請(qǐng)]
今年的國基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考研] 311求調(diào)劑 +5 | 26研0 2026-03-15 | 5/250 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
[考研] 331求調(diào)劑(0703有機(jī)化學(xué) +5 | ZY-05 2026-03-13 | 6/300 |
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|