 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 9 | 1/1 | 返回列表 |
| 查看: 1618 | 回復: 8 | |||
[交流]
求助：vasp計算rutile TiO2的電子DOS時，為何單胞和supercell的結果不一樣？ 已有2人參與
|
|
我一般都是用單胞來計算電子dos。這次因為想看一個特定的晶格distortion對rutile的影響，而這個distortion需要在一個2*2*2的supercell實現(xiàn)，所以我就嘗試用一個supercell來算dos了。 優(yōu)化的單胞展開成2*2*2的supercell（這個supercell沒有任何distortion，我是想看一看單胞和8倍的supercell得到的東西是不是一樣），開始計算，想著結果應該和單胞出來的沒有什么區(qū)別�。ㄎ矣X得dos的值應該大約是8倍，除此之外分布應該沒啥區(qū)別），但是奇怪的是supercell產(chǎn)生的dos很明顯的少了一個高能級的conduction band peak，我不知道是咋回事。比較單胞和supercell的區(qū)別，只有kpoint的密度低了一倍，這應該是沒有問題的，以下是我對這個超胞的計算文件 （1）先以4*4*4的 K point 優(yōu)化結構，進而得到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 ----------------------------------------------------------------------- 單胞的計算和上面一樣，除了KPOINT 分別是8*8*8 和 21*21*21，我也改過ISMEAR 到-1，木有啥變化。兩個系統(tǒng)內(nèi)力都收斂的很小，是優(yōu)化的很好的結構。二者的能量是7.999倍，非常好的8倍關系。附件中給出了元胞和supercell的比較，很明顯，最上面的一個peak在超胞計算中不見了 本人菜鳥，實在是想不出來什么地方不對了，望高人指點，謝謝啦! dos_rutile.png |
回復此樓» 猜你喜歡
一志愿985，本科211，0817化學工程與技術319求調(diào)劑
已經(jīng)有10人回復
-
一志愿吉林大學材料學碩321求調(diào)劑
已經(jīng)有8人回復
-
一志愿中海洋材料工程專碩330分求調(diào)劑
已經(jīng)有6人回復
-
一志愿西安交通大學材料工程專業(yè) 282分求調(diào)劑
已經(jīng)有3人回復
-
332求調(diào)劑
已經(jīng)有3人回復
-
294求調(diào)劑材料與化工專碩
已經(jīng)有7人回復
-
281求調(diào)劑（0805）
已經(jīng)有17人回復
-
267一志愿南京工業(yè)大學0817化工求調(diào)劑
已經(jīng)有10人回復
-
0817調(diào)劑
已經(jīng)有3人回復
-
焦慮
已經(jīng)有13人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- VASP計算TiO2表面能問題 求助
已經(jīng)有8人回復
- VASP計算TiO2能帶竟然沒有帶隙
已經(jīng)有3人回復
- VESTA畫電荷密度圖的疑問
已經(jīng)有20人回復
- TiO2（101）晶面的優(yōu)化
已經(jīng)有13人回復
- 【轉(zhuǎn)帖】用Material Studio構建復雜模型及轉(zhuǎn)換成poscar的2種方法
已經(jīng)有56人回復
jpchou
專家顧問 (著名寫手)
-
專家經(jīng)驗: +20 - 1ST強帖: 5
- 應助: 231 (大學生)
- 金幣: 2342
- 散金: 4367
- 紅花: 103
- 帖子: 1061
- 在線: 757.3小時
- 蟲號: 585915
- 注冊: 2008-08-08
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
- 管轄: 第一性原理
2樓2014-11-25 16:05:27
3樓2014-11-26 01:20:55
jpchou
專家顧問 (著名寫手)
-
專家經(jīng)驗: +20 - 1ST強帖: 5
- 應助: 231 (大學生)
- 金幣: 2342
- 散金: 4367
- 紅花: 103
- 帖子: 1061
- 在線: 757.3小時
- 蟲號: 585915
- 注冊: 2008-08-08
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
- 管轄: 第一性原理
4樓2014-11-26 03:40:40
失敗的路人
木蟲 (正式寫手)
- 應助: 3 (幼兒園)
- 金幣: 2824.7
- 散金: 20
- 紅花: 1
- 帖子: 362
- 在線: 292.9小時
- 蟲號: 2506020
- 注冊: 2013-06-13
- 性別: GG
- 專業(yè): 理論和計算化學
5樓2014-11-26 08:19:41
|
遇見大神了  果然是專家，一語點破玄機，我發(fā)現(xiàn)1個單胞的時候vasp by default的nbands = 32, 但是8個超胞的情況下nbands = 216，分明小于8倍的關系，結果導致上層空帶沒有被考慮進去，所以少了一個Peak。。。我把nbands提高到256，也就是正好8倍，結果就一樣了 這個問題我是不可能想到的，我以為限定了能量區(qū)間，該算的都能算了呢 這個問題vasp online manul說的好像不太對啊，上面說如果非spin polarized的計算，default的nbands = nelect/2+nion/2，我用的ti_pv所以是10個電子，o是6個，一共44個電子/cell，這個怎么也算不出nbands=32�。� |
6樓2014-11-26 13:01:07
jpchou
專家顧問 (著名寫手)
-
專家經(jīng)驗: +20 - 1ST強帖: 5
- 應助: 231 (大學生)
- 金幣: 2342
- 散金: 4367
- 紅花: 103
- 帖子: 1061
- 在線: 757.3小時
- 蟲號: 585915
- 注冊: 2008-08-08
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
- 管轄: 第一性原理
7樓2014-11-26 16:07:57
jpchou
專家顧問 (著名寫手)
-
專家經(jīng)驗: +20 - 1ST強帖: 5
- 應助: 231 (大學生)
- 金幣: 2342
- 散金: 4367
- 紅花: 103
- 帖子: 1061
- 在線: 757.3小時
- 蟲號: 585915
- 注冊: 2008-08-08
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
- 管轄: 第一性原理
8樓2014-11-26 16:09:18
感謝感謝 9樓2014-11-27 04:09:15
相關版塊跳轉(zhuǎn)
我要訂閱樓主 Emerald01 的主題更新
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[教師之家] 焦慮 +9 | 水冰月月野兔 2026-03-13 | 13/650 |
|
|---|---|---|---|---|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-16 | 10/500 |
|
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
|
[考研] 265求調(diào)劑 +3 | 梁梁校校 2026-03-17 | 3/150 |
|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠意在機?/a> 2026-03-16 | 3/150 |
|
|
|
[考研] 290求調(diào)劑 +6 | 孔志浩 2026-03-12 | 11/550 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽風就是雨； 2026-03-16 | 3/150 |
|
|
|
[考研] 274求調(diào)劑 +5 | 時間點 2026-03-13 | 5/250 |
|
|
|
[基金申請]
今年的國基金是打分制嗎？
 50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy！ 2026-03-13 | 6/300 |
|
|
|
[考研] 考研調(diào)劑 +4 | 芬達46 2026-03-12 | 4/200 |
|
|
|
[考研] 274求調(diào)劑 +3 | S.H1 2026-03-12 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
|
[考研] 081200-11408-276學碩求調(diào)劑 +3 | 崔wj 2026-03-12 | 4/200 |
|