 |
 |
 |
| 24小時(shí)熱門版塊排行榜 |
- >論壇更新日志 (4280)
- >導(dǎo)師招生 (586)
- >考研 (499)
- >基金申請(qǐng) (319)
- >文獻(xiàn)求助 (219)
- >考博 (159)
- >蟲(chóng)友互識(shí) (155)
- >休閑灌水 (77)
- >碩博家園 (50)
- >招聘信息布告欄 (49)
- >論文投稿 (49)
- >博后之家 (42)
- >論文道賀祈福 (31)
- >綠色求助(高懸賞) (21)
- >SciFinder/Reaxys (20)
- >找工作 (20)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1619 | 回復(fù): 8 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
Emerald01新蟲(chóng) (初入文壇)
|
[交流]
求助：vasp計(jì)算rutile TiO2的電子DOS時(shí)，為何單胞和supercell的結(jié)果不一樣？ 已有2人參與
|
||
|
我一般都是用單胞來(lái)計(jì)算電子dos。這次因?yàn)橄肟匆粋�(gè)特定的晶格distortion對(duì)rutile的影響，而這個(gè)distortion需要在一個(gè)2*2*2的supercell實(shí)現(xiàn)，所以我就嘗試用一個(gè)supercell來(lái)算dos了。 優(yōu)化的單胞展開(kāi)成2*2*2的supercell（這個(gè)supercell沒(méi)有任何distortion，我是想看一看單胞和8倍的supercell得到的東西是不是一樣），開(kāi)始計(jì)算，想著結(jié)果應(yīng)該和單胞出來(lái)的沒(méi)有什么區(qū)別啊（我覺(jué)得dos的值應(yīng)該大約是8倍，除此之外分布應(yīng)該沒(méi)啥區(qū)別），但是奇怪的是supercell產(chǎn)生的dos很明顯的少了一個(gè)高能級(jí)的conduction band peak，我不知道是咋回事。比較單胞和supercell的區(qū)別，只有kpoint的密度低了一倍，這應(yīng)該是沒(méi)有問(wèn)題的，以下是我對(duì)這個(gè)超胞的計(jì)算文件 （1）先以4*4*4的 K point 優(yōu)化結(jié)構(gòu)，進(jìn)而得到chgcar INCAR: SYSTEM = rutile Relax ISTART = 0 ICHARG = 2 ENCUT = 500 PREC = Accurate NELMIN = 4 # do a minimum of four electronic steps NELM = 40 # maximum of 40 electronic steps EDIFF = 1E-8 EDIFFG = -1E-6 # high accuracy ADDGRID = .TRUE. IBRION = -1 ISMEAR = 0 SIGMA = 0.05 LREAL = .False. LWAVE = .False. LPLANE = .TRUE. NPAR = 8 # twice the number of nodes ------------------------------------------------------------------------------ (2) 拷貝chgcar到新directory，增加KPOINT 到 11*11*11 INCAR SYSTEM = rutile Relax ISTART = 0 ICHARG = 11 ENCUT = 500 PREC = Accurate NELMIN = 4 # do a minimum of four electronic steps NELM = 40 # maximum of 40 electronic steps EDIFF = 1E-8 EDIFFG = -1E-6 # high accuracy ADDGRID = .TRUE. IBRION = -1 ISMEAR = 0 SIGMA = 0.05 LORBIT = 10 #EMIN = -15 #EMAX = 15 #NEDOS = 1000 LREAL = .False. LWAVE = .False. LPLANE = .TRUE. NPAR = 8 # twice the number of nodes ----------------------------------------------------------------------- 單胞的計(jì)算和上面一樣，除了KPOINT 分別是8*8*8 和 21*21*21，我也改過(guò)ISMEAR 到-1，木有啥變化。兩個(gè)系統(tǒng)內(nèi)力都收斂的很小，是優(yōu)化的很好的結(jié)構(gòu)。二者的能量是7.999倍，非常好的8倍關(guān)系。附件中給出了元胞和supercell的比較，很明顯，最上面的一個(gè)peak在超胞計(jì)算中不見(jiàn)了 本人菜鳥(niǎo)，實(shí)在是想不出來(lái)什么地方不對(duì)了，望高人指點(diǎn)，謝謝啦! dos_rutile.png |
回復(fù)此樓» 猜你喜歡
294求調(diào)劑材料與化工專碩
已經(jīng)有10人回復(fù)
-
288求調(diào)劑
已經(jīng)有6人回復(fù)
-
287求調(diào)劑
已經(jīng)有4人回復(fù)
-
招收調(diào)劑碩士
已經(jīng)有7人回復(fù)
-
085601材料工程專碩求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿武漢理工材料工程專碩調(diào)劑
已經(jīng)有3人回復(fù)
-
【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠(chéng)心求調(diào)劑
已經(jīng)有3人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有19人回復(fù)
-
材料考研調(diào)劑
已經(jīng)有3人回復(fù)
-
346求調(diào)劑[0856]
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- VASP計(jì)算TiO2表面能問(wèn)題 求助
已經(jīng)有8人回復(fù)
- VASP計(jì)算TiO2能帶竟然沒(méi)有帶隙
已經(jīng)有3人回復(fù)
- VESTA畫電荷密度圖的疑問(wèn)
已經(jīng)有20人回復(fù)
- TiO2（101）晶面的優(yōu)化
已經(jīng)有13人回復(fù)
- 【轉(zhuǎn)帖】用Material Studio構(gòu)建復(fù)雜模型及轉(zhuǎn)換成poscar的2種方法
已經(jīng)有56人回復(fù)
Emerald01
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 帖子: 4
- 在線: 24分鐘
- 蟲(chóng)號(hào): 3513440
- 注冊(cè): 2014-11-02
|
遇見(jiàn)大神了  果然是專家，一語(yǔ)點(diǎn)破玄機(jī)，我發(fā)現(xiàn)1個(gè)單胞的時(shí)候vasp by default的nbands = 32, 但是8個(gè)超胞的情況下nbands = 216，分明小于8倍的關(guān)系，結(jié)果導(dǎo)致上層空帶沒(méi)有被考慮進(jìn)去，所以少了一個(gè)Peak。。。我把nbands提高到256，也就是正好8倍，結(jié)果就一樣了 這個(gè)問(wèn)題我是不可能想到的，我以為限定了能量區(qū)間，該算的都能算了呢 這個(gè)問(wèn)題vasp online manul說(shuō)的好像不太對(duì)啊，上面說(shuō)如果非spin polarized的計(jì)算，default的nbands = nelect/2+nion/2，我用的ti_pv所以是10個(gè)電子，o是6個(gè)，一共44個(gè)電子/cell，這個(gè)怎么也算不出nbands=32��？ |
6樓2014-11-26 13:01:07
jpchou
專家顧問(wèn) (著名寫手)
-
專家經(jīng)驗(yàn): +20 - 1ST強(qiáng)帖: 5
- 應(yīng)助: 231 (大學(xué)生)
- 金幣: 2342
- 散金: 4367
- 紅花: 103
- 帖子: 1061
- 在線: 757.3小時(shí)
- 蟲(chóng)號(hào): 585915
- 注冊(cè): 2008-08-08
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
- 管轄: 第一性原理
2樓2014-11-25 16:05:27
Emerald01
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 帖子: 4
- 在線: 24分鐘
- 蟲(chóng)號(hào): 3513440
- 注冊(cè): 2014-11-02
3樓2014-11-26 01:20:55
jpchou
專家顧問(wèn) (著名寫手)
-
專家經(jīng)驗(yàn): +20 - 1ST強(qiáng)帖: 5
- 應(yīng)助: 231 (大學(xué)生)
- 金幣: 2342
- 散金: 4367
- 紅花: 103
- 帖子: 1061
- 在線: 757.3小時(shí)
- 蟲(chóng)號(hào): 585915
- 注冊(cè): 2008-08-08
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
- 管轄: 第一性原理
4樓2014-11-26 03:40:40
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 288求調(diào)劑 +6 | 于海海海海 2026-03-19 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
|
[考研] 332求調(diào)劑 +3 | ydfyh 2026-03-17 | 3/150 |
|
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +9 | rare12345 2026-03-18 | 9/450 |
|
|
|
[考研] 354求調(diào)劑 +4 | Tyoumou 2026-03-18 | 7/350 |
|
|
|
[考研] 一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +4 | 1孫悟空 2026-03-17 | 4/200 |
|
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯！ 2026-03-15 | 18/900 |
|
|
|
[考研] 299求調(diào)劑 +5 | △小透明* 2026-03-17 | 5/250 |
|
|
|
[考研] 環(huán)境工程調(diào)劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號(hào) 2026-03-16 | 16/800 |
|
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +8 | 一鴨鴨喲 2026-03-14 | 10/500 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
|
[考研] 085600調(diào)劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
|
[考研] 080500，材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識(shí)可樂(lè) 2026-03-14 | 5/250 |
|
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
|
[考研] 289求調(diào)劑 +3 | 李政瑩 2026-03-12 | 3/150 |
|