 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3960 | 回復: 11 | ||||||||||||
| 【獎勵】 本帖被評價8次，作者xinqing98增加金幣 6 個 | ||||||||||||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||||||||||||
[資源]
以石墨烯為例Pwscf計算態(tài)密度的詳細步驟和要點
|
||||||||||||
|
使用PWSCF計算態(tài)密度需要以下幾個步驟： 1. 構(gòu)建好結(jié)構(gòu)以后進行scf計算，以附件中的Graphene.scf.in為例，其中重要的是在calculation='scf' 這是告訴系統(tǒng)進行自洽計算。另外，在K_POINTS的設置中，要寫成“K_POINTS automatic”，然后下面寫上 12 12 1 1 1 1，前面空一個空格，這是產(chǎn)生K點的網(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)做非自洽計算。另外還要在 &system中加入nbnd = 這是繪制的能帶的數(shù)目，這一項參數(shù)在計算態(tài)密度的這個步驟中一定要添加，如果 不添加，輸入的態(tài)密度結(jié)果只有費米能級以下的態(tài)密度。 另外，在K_POINTS的設置中，要寫成“K_POINTS automatic”，K點的網(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 / 注意：其中的每一個字符都要寫對，否則計算一開始就出錯，要仔細檢查。 執(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)密度文件。 |
回復此樓» 收錄本帖的淘帖專輯推薦
 精華網(wǎng)帖收集 | VASP and MS | 石墨烯模擬計算 | vasp 第一性原理 |
| 第一性 | 量子化學計算與使用經(jīng)驗 | 科研 | pwscf |
| 方法不錯 | 化工專業(yè)課資料 |
» 猜你喜歡
309分085801求調(diào)劑
已經(jīng)有7人回復
-
生物學308分求調(diào)劑（一志愿華東師大）
已經(jīng)有3人回復
-
求調(diào)劑：一志愿：南京大學 專業(yè)：0705 總分320 ，本科985，四六級已過
已經(jīng)有3人回復
-
環(huán)境工程 085701，267求調(diào)劑
已經(jīng)有14人回復
-
308求調(diào)劑
已經(jīng)有12人回復
-
求調(diào)劑：085600材料與化工，考材科基，總分319
已經(jīng)有21人回復
-
311求調(diào)劑
已經(jīng)有8人回復
-
294分080500材料科學與工程求調(diào)劑
已經(jīng)有13人回復
-
288求調(diào)劑 一志愿哈工大 材料與化工
已經(jīng)有19人回復
-
288資源與環(huán)境專碩求調(diào)劑，不限專業(yè)，有學上就行
已經(jīng)有23人回復
» 本主題相關商家推薦: (我也要在這里推廣)
» 本主題相關價值貼推薦，對您同樣有幫助:
- 石墨烯的總態(tài)密度與分波態(tài)密度
已經(jīng)有25人回復
- 用wien2k計算反鐵磁詳細步驟，舉例說明！
已經(jīng)有12人回復
7樓2014-05-30 11:03:18
3樓2013-08-23 15:30:44
4樓2013-11-10 16:31:07
5樓2013-11-27 10:19:57
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 070300化學專業(yè)279調(diào)劑 +10 | 哈哈哈^_^ 2026-03-31 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +10 | 李李不服輸 2026-03-25 | 10/500 |
|
|
|
[考研] 一志愿西交大080500材料學碩349 +6 | jqx1258 2026-03-31 | 7/350 |
|
|
|
[考研] 080200學碩，機械工程專業(yè)277分，求帶走！ +4 | 瓶子PZ 2026-03-31 | 4/200 |
|
|
|
[考研] 求調(diào)劑推薦 材料 304 +18 | 荷包蛋hyj 2026-03-26 | 18/900 |
|
|
|
[考研] 362求調(diào)劑 +9 | 西南交材料專碩3 2026-03-31 | 9/450 |
|
|
|
[考研] 287求調(diào)劑 +17 | land xuxu 2026-03-26 | 17/850 |
|
|
|
[考研] 一志愿中海洋320化學工程與技術(shù)學碩求調(diào)劑 +8 | 披星河 2026-03-30 | 8/400 |
|
|
|
[考研] 085601一志愿中山大學深圳材料工程330求調(diào)劑 +5 | pipiver 2026-03-30 | 5/250 |
|
|
|
[考研] 一志愿食品科學與工程083200求調(diào)劑 +4 | XQTJZ 2026-03-30 | 4/200 |
|
|
|
[考研] 英一數(shù)一408，總分284，二戰(zhàn)真誠求調(diào)劑 +3 | 12.27 2026-03-30 | 5/250 |
|
|
|
[考研]
 304材料求調(diào)劑
+5
|
鐘llll 2026-03-26 | 5/250 |
|
|
|
[考研] 一志愿南開大學0710生物學359求調(diào)劑 +5 | 兔兔兔111223314 2026-03-29 | 7/350 |
|
|
|
[考研] 332求調(diào)劑 +6 | Lyy930824@ 2026-03-29 | 6/300 |
|
|
|
[考研] 求調(diào)劑，一志愿 南京航空航天大學大學 ，080500材料科學與工程學碩 +6 | @taotao 2026-03-26 | 7/350 |
|
|
|
[考研] 一志愿中南大學化學0703總分337求調(diào)劑 +6 | niko- 2026-03-27 | 6/300 |
|
|
|
[考研] 337求調(diào)劑 +6 | 《樹》 2026-03-29 | 6/300 |
|
|
|
[考研] 321求調(diào)劑 +7 | 璞玉~~ 2026-03-25 | 8/400 |
|
|
|
[考研] 本科新能源科學與工程，一志愿華理能動285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
|
[考研] 各位老師您好：本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|