 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3949 | 回復: 11 | ||||||||||||
| 【獎勵】 本帖被評價8次，作者xinqing98增加金幣 6 個 | ||||||||||||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||||||||||||
[資源]
以石墨烯為例Pwscf計算態(tài)密度的詳細步驟和要點
|
||||||||||||
|
使用PWSCF計算態(tài)密度需要以下幾個步驟： 1. 構建好結構以后進行scf計算，以附件中的Graphene.scf.in為例，其中重要的是在calculation='scf' 這是告訴系統(tǒng)進行自洽計算。另外，在K_POINTS的設置中，要寫成“K_POINTS automatic”，然后下面寫上 12 12 1 1 1 1，前面空一個空格，這是產生K點的網格大小，后面寫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 構建好以后，執(zhí)行： pw.x < Graphene.scf.in > Graphene.scf.out & 2. Graphene.nscf.in 在這個文件中，calculation='nscf'，這是告訴系統(tǒng)做非自洽計算。另外還要在 &system中加入nbnd = 這是繪制的能帶的數(shù)目，這一項參數(shù)在計算態(tài)密度的這個步驟中一定要添加，如果 不添加，輸入的態(tài)密度結果只有費米能級以下的態(tài)密度。 另外，在K_POINTS的設置中，要寫成“K_POINTS automatic”，K點的網格可以給大一些。 &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ù)值一定要和前兩步完全一樣，否則將無法讀取前兩步 的結果。 &dos prefix = 'Graphene' outdir = './' fildos = 'Graphene.dos' Emin=-22.0, Emax=15.0, DeltaE=0.1, ngauss=0, degauss=0.025 / 注意：其中的每一個字符都要寫對，否則計算一開始就出錯，要仔細檢查。 執(zhí)行：dos.x < dos.in > dos.out & 這步計算完以后，就會自動生成一個名位“Graphene.dos”的文件，用gnuplot這個文件的圖形即得到 體系的總態(tài)密度。 4. 如果要計算體系的分波態(tài)密度，則需要重新編輯另外一個文件：pdos.in，該文件的形式和dos.in的形式 大體相同，但要注意第四行的“fildos”要改成“filpdos”，表示計算的是分波態(tài)密度,這里如果不換過來， 計算立即報錯，我用了很長時間才找到這個錯誤，真不容易！還有就是文件的第一行最好改成&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 & 計算完以后，就生成每個原子各個電子態(tài)的態(tài)密度文件。 |
回復此樓» 收錄本帖的淘帖專輯推薦
 精華網帖收集 | VASP and MS | 石墨烯模擬計算 | vasp 第一性原理 |
| 第一性 | 量子化學計算與使用經驗 | 科研 | pwscf |
| 方法不錯 | 化工專業(yè)課資料 |
» 猜你喜歡
291求調劑
已經有7人回復
-
材料與化工272求調劑
已經有13人回復
-
283求調劑（080500）
已經有5人回復
-
材料科學與工程求調劑
已經有5人回復
-
一志愿北京化工大學材料與化工（085600）296求調劑
已經有17人回復
-
考研調劑
已經有7人回復
-
343求調劑085601
已經有3人回復
-
08工科求調劑286
已經有3人回復
-
337求調劑
已經有3人回復
-
329求調劑，一志愿西北工業(yè)大學，材料工程（085601）
已經有10人回復
» 本主題相關商家推薦: (我也要在這里推廣)
» 本主題相關價值貼推薦，對您同樣有幫助:
- 石墨烯的總態(tài)密度與分波態(tài)密度
已經有25人回復
- 用wien2k計算反鐵磁詳細步驟，舉例說明！
已經有12人回復
3樓2013-08-23 15:30:44
4樓2013-11-10 16:31:07
5樓2013-11-27 10:19:57
6樓2014-01-16 05:57:08
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 調劑求院校招收 +6 | 鶴鯨鴿 2026-03-28 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 083000學碩274求調劑 +8 | Li李魚 2026-03-26 | 8/400 |
|
|
[考研] 一志愿中南大學化學0703總分337求調劑 +5 | niko- 2026-03-27 | 5/250 |
|
|
|
[考研] 0703一志愿9，初試成績：338，四六級已過，有科研經歷，求調劑！ +4 | Zuhui0306 2026-03-25 | 4/200 |
|
|
|
[考研] 調劑 +3 | 好好讀書。 2026-03-28 | 3/150 |
|
|
|
[考研] 292求調劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
|
[考研] 0703化學/290求調劑/本科經歷豐富/工科也可 +9 | 丹青奶蓋 2026-03-26 | 10/500 |
|
|
|
[考研] 086502化學工程342求調劑 +6 | 阿姨復古不過 2026-03-27 | 6/300 |
|
|
|
[考研] 材料求調劑 一志愿哈工大總分298分，前三科223分 +5 | dongfang59 2026-03-27 | 5/250 |
|
|
|
[考研] 275求調劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|
|
|
[考研] 265求調劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
|
[考研] 266求調劑 +11 | 陽陽哇塞 2026-03-27 | 12/600 |
|
|
|
[考研] 化學調劑 +4 | 愛吃番茄的旭 2026-03-24 | 5/250 |
|
|
|
[考研] 考研調劑 +10 | 呼呼？~+123456 2026-03-24 | 10/500 |
|
|
|
[考研] 333求調劑 +3 | question挽風 2026-03-23 | 3/150 |
|
|
|
[考研] 286求調劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
|
[考研] 321求調劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
|
[考研] 調劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
|
[考研] 總分322求生物學/生化與分子/生物信息學相關調劑 +5 | 星沉uu 2026-03-26 | 6/300 |
|
|
|
[考研] 中國科學院深圳先進技術研究院-光纖傳感課題組招生-中國科學院大學、深圳理工大學聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|