| 10 | 1/1 | 返回列表 |
| 查看: 3351 | 回復(fù): 9 | ||
gynan2008金蟲 (正式寫手)
|
[求助]
鈦酸鋇晶格參數(shù)哪個正確?
|
|
文獻里面有看到立方相晶格參數(shù)是3.95*****,前幾天用abinit計算時,別人給的常數(shù)是7.901******,我用不同的acell計算總能,也是七點幾左右。 兩個問題,1哪個正確? 2、我用不同acel去計算總能,是不是就是在優(yōu)化晶格參數(shù)? 3.nband的設(shè)置我還是不明白,不是說價電子數(shù)/2么,可是設(shè)置的要比物質(zhì)的價電子數(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 #定義原胞中原子的總個數(shù),本例子中有2個原子 typat 1 2 3 3 3 #定義每類原子的個數(shù),本例子中第一類原子有2個 #下面定義了原胞中原子的坐標 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 #定義了要計算的能帶的數(shù)目,最好按這樣來設(shè)置: # nband > 原胞中總的價電子數(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 #定義了平面波的切斷動能 #Definition of the k-point grid ngkpt 8 8 8 #下面定義了k點網(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 #電子自洽收斂的標準 enunit 2 rprim: 1 0 0 #下面三行定義了原胞的基矢,本例子中Si是fcc結(jié)構(gòu) 0 1 0 0 0 1 我算了一下,結(jié)果是下面的樣子,參數(shù)設(shè)置有問題么? 沒辦法傳圖片了,是一條收斂的曲線 |

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

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

金蟲 (正式寫手)

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

銀蟲 (初入文壇)
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 281求調(diào)劑(0805) +4 | 煙汐憶海 2026-03-16 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 331求調(diào)劑(0703有機化學(xué) +7 | ZY-05 2026-03-13 | 8/400 |
|
|
[考研] 299求調(diào)劑 +5 | △小透明* 2026-03-17 | 5/250 |
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +5 | 慕寒mio 2026-03-16 | 5/250 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林,給自己招師弟中~ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
[考研] 08工科 320總分 求調(diào)劑 +4 | 梨花珞晚風(fēng) 2026-03-17 | 4/200 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研]
|
笨笨兔子 2026-03-12 | 3/150 |
|
|
[考研] 材料080500調(diào)劑求收留 +3 | 一顆meteor 2026-03-13 | 3/150 |
|
|
[考研] [0860]321分求調(diào)劑,ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
[考研] 333求調(diào)劑 +3 | 球球古力 2026-03-11 | 3/150 |
|
|
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|