| 5 | 1/1 | 返回列表 |
| 查看: 3354 | 回復(fù): 9 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
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í),那是心理安慰而已。 建議樓主先好好看一下幫助文檔~ |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +4 | 脫穎而出 2026-03-16 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 324分 085600材料化工求調(diào)劑 +3 | llllkkkhh 2026-03-18 | 3/150 |
|
|
[考研] 287求調(diào)劑 +3 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
[考研] 一志愿天大材料與化工(085600)總分338 +5 | 蔡大美女 2026-03-13 | 5/250 |
|
|
[考研] 材料080500調(diào)劑求收留 +4 | 一顆meteor 2026-03-13 | 4/200 |
|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +10 | SUICHILD 2026-03-12 | 10/500 |
|
|
[考研] 304求調(diào)劑 +6 | 司空. 2026-03-18 | 6/300 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 297求調(diào)劑 +8 | 戲精丹丹丹 2026-03-17 | 8/400 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林,給自己招師弟中~ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15 | 7/350 |
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
[考研] 26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy! 2026-03-13 | 6/300 |
|