 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3950 | 回復(fù): 11 | ||||||||||||
| 【獎勵】 本帖被評價(jià)8次，作者xinqing98增加金幣 6 個 | ||||||||||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||||||||||
[資源]
以石墨烯為例Pwscf計(jì)算態(tài)密度的詳細(xì)步驟和要點(diǎn)
|
||||||||||||
|
使用PWSCF計(jì)算態(tài)密度需要以下幾個步驟： 1. 構(gòu)建好結(jié)構(gòu)以后進(jìn)行scf計(jì)算，以附件中的Graphene.scf.in為例，其中重要的是在calculation='scf' 這是告訴系統(tǒng)進(jìn)行自洽計(jì)算。另外，在K_POINTS的設(shè)置中，要寫成“K_POINTS automatic”，然后下面寫上 12 12 1 1 1 1，前面空一個空格，這是產(chǎn)生K點(diǎn)的網(wǎng)格大小，后面寫1 1 1，而不是像以前那樣寫0 0 0。 &control calculation='scf' prefix='Graphene', pseudo_dir = './', outdir='./' / &system ibrav=4, celldm(1) =4.6478043778273781, celldm(3) =8.1316939322482502, nat=2, ntyp=1, ecutwfc = 27.0, / &electrons conv_thr = 1.0d-6 mixing_beta = 0.7 / ATOMIC_SPECIES C 12.0107 C.pbe-rrkjus.UPF ATOMIC_POSITIONS crystal C 0.6666666666666667 0.3333333333333333 0.5000000000000000 C 0.3333333333333333 0.6666666666666667 0.5000000000000000 K_POINTS automatic 12 12 1 1 1 1 構(gòu)建好以后，執(zhí)行： pw.x < Graphene.scf.in > Graphene.scf.out & 2. Graphene.nscf.in 在這個文件中，calculation='nscf'，這是告訴系統(tǒng)做非自洽計(jì)算。另外還要在 &system中加入nbnd = 這是繪制的能帶的數(shù)目，這一項(xiàng)參數(shù)在計(jì)算態(tài)密度的這個步驟中一定要添加，如果 不添加，輸入的態(tài)密度結(jié)果只有費(fèi)米能級以下的態(tài)密度。 另外，在K_POINTS的設(shè)置中，要寫成“K_POINTS automatic”，K點(diǎn)的網(wǎng)格可以給大一些。 &control calculation='nscf' prefix='Graphene', pseudo_dir = './', outdir='./' / &system ibrav=4, celldm(1) =4.6478043778273781, celldm(3) =8.1316939322482502, nat=2, ntyp=1, ecutwfc = 27.0, nbnd = 8, / &electrons conv_thr = 1.0d-6 mixing_beta = 0.7 / ATOMIC_SPECIES C 12.0107 C.pbe-rrkjus.UPF ATOMIC_POSITIONS crystal C 0.6666666666666667 0.3333333333333333 0.5000000000000000 C 0.3333333333333333 0.6666666666666667 0.5000000000000000 K_POINTS automatic 12 12 1 1 1 1 執(zhí)行：pw.x < Graphene.nscf.in > Graphene.nscf.out & 3. 編輯dos.in文件，注意其中的prefix和outdir的數(shù)值一定要和前兩步完全一樣，否則將無法讀取前兩步 的結(jié)果。 &dos prefix = 'Graphene' outdir = './' fildos = 'Graphene.dos' Emin=-22.0, Emax=15.0, DeltaE=0.1, ngauss=0, degauss=0.025 / 注意：其中的每一個字符都要寫對，否則計(jì)算一開始就出錯，要仔細(xì)檢查。 執(zhí)行：dos.x < dos.in > dos.out & 這步計(jì)算完以后，就會自動生成一個名位“Graphene.dos”的文件，用gnuplot這個文件的圖形即得到 體系的總態(tài)密度。 4. 如果要計(jì)算體系的分波態(tài)密度，則需要重新編輯另外一個文件：pdos.in，該文件的形式和dos.in的形式 大體相同，但要注意第四行的“fildos”要改成“filpdos”，表示計(jì)算的是分波態(tài)密度,這里如果不換過來， 計(jì)算立即報(bào)錯，我用了很長時間才找到這個錯誤，真不容易！還有就是文件的第一行最好改成&projwfc。 &projwfc prefix = 'Graphene' outdir = './' filpdos = 'Graphene.dos' Emin=-22.0, Emax=15.0, DeltaE=0.1, ngauss=0, degauss=0.025 / 執(zhí)行：projwfc.x < pdos.in > pdos.out & 計(jì)算完以后，就生成每個原子各個電子態(tài)的態(tài)密度文件。 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 精華網(wǎng)帖收集 | VASP and MS | 石墨烯模擬計(jì)算 | vasp 第一性原理 |
| 第一性 | 量子化學(xué)計(jì)算與使用經(jīng)驗(yàn) | 科研 | pwscf |
| 方法不錯 | 化工專業(yè)課資料 |
» 猜你喜歡
309求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷
已經(jīng)有8人回復(fù)
-
求調(diào)劑
已經(jīng)有5人回復(fù)
-
求化學(xué)調(diào)劑
已經(jīng)有5人回復(fù)
-
085600，材料與化工321分求調(diào)劑
已經(jīng)有10人回復(fù)
-
0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分
已經(jīng)有5人回復(fù)
-
一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑
已經(jīng)有3人回復(fù)
-
085600，專業(yè)課化工原理，320分求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料與化工272求調(diào)劑
已經(jīng)有16人回復(fù)
-
290求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 石墨烯的總態(tài)密度與分波態(tài)密度
已經(jīng)有25人回復(fù)
- 用wien2k計(jì)算反鐵磁詳細(xì)步驟，舉例說明！
已經(jīng)有12人回復(fù)
4樓2013-11-10 16:31:07
3樓2013-08-23 15:30:44
5樓2013-11-27 10:19:57
6樓2014-01-16 05:57:08
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600，材料與化工321分求調(diào)劑 +10 | 大饞小子 2026-03-28 | 10/500 |
|
|---|---|---|---|---|
|
|
[考研] 求調(diào)劑 +4 | QiMing7 2026-03-25 | 5/250 |
|
|
|
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +8 | 柳溪邊 2026-03-26 | 8/400 |
|
|
[考研] 307求調(diào)劑 +9 | 超級伊昂大王 2026-03-24 | 10/500 |
|
|
|
[考研] 375求調(diào)劑 +4 | 雨夏整夜 2026-03-29 | 4/200 |
|
|
|
[考研] 298求調(diào)劑 +3 | 種圣賜 2026-03-29 | 3/150 |
|
|
|
[考研] 085600，專業(yè)課化工原理，321分求調(diào)劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
|
[考研] 085701求調(diào)劑初試286分 +4 | secret0328 2026-03-28 | 4/200 |
|
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-03-28 | 3/150 |
|
|
|
[考研] 材料277求調(diào)劑 +7 | min3 2026-03-24 | 7/350 |
|
|
|
[考研] 328求調(diào)劑 +7 | 嗯滴的基本都 2026-03-27 | 7/350 |
|
|
|
[考研] 求調(diào)劑推薦 材料 304 +15 | 荷包蛋hyj 2026-03-26 | 15/750 |
|
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
|
[考研] 一志愿吉大071010，316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|
|
|
[考研] 359求調(diào)劑 +4 | 王了個楠 2026-03-25 | 4/200 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
|
[考研] 282求調(diào)劑 +3 | wcq131415 2026-03-24 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +6 | 研研，接電話 2026-03-24 | 7/350 |
|
|
|
[考研]
 085404電子信息284分求調(diào)劑
+4
|
13659058978 2026-03-24 | 4/200 |
|