 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 12 | 1/1 | 返回列表 |
| 查看: 3948 | 回復(fù): 11 | ||||||||||||
| 【獎(jiǎng)勵(lì)】 本帖被評價(jià)8次，作者xinqing98增加金幣 6 個(gè) | ||||||||||||
[資源]
以石墨烯為例Pwscf計(jì)算態(tài)密度的詳細(xì)步驟和要點(diǎn)
|
||||||||||||
|
使用PWSCF計(jì)算態(tài)密度需要以下幾個(gè)步驟： 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，前面空一個(gè)空格，這是產(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 在這個(gè)文件中，calculation='nscf'，這是告訴系統(tǒng)做非自洽計(jì)算。另外還要在 &system中加入nbnd = 這是繪制的能帶的數(shù)目，這一項(xiàng)參數(shù)在計(jì)算態(tài)密度的這個(gè)步驟中一定要添加，如果 不添加，輸入的態(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 / 注意：其中的每一個(gè)字符都要寫對，否則計(jì)算一開始就出錯(cuò)，要仔細(xì)檢查。 執(zhí)行：dos.x < dos.in > dos.out & 這步計(jì)算完以后，就會自動(dòng)生成一個(gè)名位“Graphene.dos”的文件，用gnuplot這個(gè)文件的圖形即得到 體系的總態(tài)密度。 4. 如果要計(jì)算體系的分波態(tài)密度，則需要重新編輯另外一個(gè)文件：pdos.in，該文件的形式和dos.in的形式 大體相同，但要注意第四行的“fildos”要改成“filpdos”，表示計(jì)算的是分波態(tài)密度,這里如果不換過來， 計(jì)算立即報(bào)錯(cuò)，我用了很長時(shí)間才找到這個(gè)錯(cuò)誤，真不容易！還有就是文件的第一行最好改成&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ì)算完以后，就生成每個(gè)原子各個(gè)電子態(tài)的態(tài)密度文件。 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 精華網(wǎng)帖收集 | VASP and MS | 石墨烯模擬計(jì)算 | vasp 第一性原理 |
| 第一性 | 量子化學(xué)計(jì)算與使用經(jīng)驗(yàn) | 科研 | pwscf |
| 方法不錯(cuò) | 化工專業(yè)課資料 |
» 猜你喜歡
329求調(diào)劑
已經(jīng)有8人回復(fù)
-
070305高分子化學(xué)與物理 304分求調(diào)劑
已經(jīng)有8人回復(fù)
-
349求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑
已經(jīng)有7人回復(fù)
-
291求調(diào)劑
已經(jīng)有4人回復(fù)
-
0856材料化工調(diào)劑 總分330
已經(jīng)有14人回復(fù)
-
352分-085602-一志愿985
已經(jīng)有5人回復(fù)
-
各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩
已經(jīng)有8人回復(fù)
-
求調(diào)劑，一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分
已經(jīng)有4人回復(fù)
-
279求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 石墨烯的總態(tài)密度與分波態(tài)密度
已經(jīng)有25人回復(fù)
- 用wien2k計(jì)算反鐵磁詳細(xì)步驟，舉例說明！
已經(jīng)有12人回復(fù)
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
6樓2014-01-16 05:57:08
8樓2015-04-10 20:10:34
11樓2016-06-17 13:24:38
12樓2016-07-10 14:58:05
簡單回復(fù)
卡崔娜2樓
2013-08-21 08:02
回復(fù)
五星好評 
吳聲豪9樓
2015-09-21 09:30
回復(fù)
五星好評 thx
吳聲豪10樓
2015-10-11 10:17
回復(fù)
thx
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 xinqing98 的主題更新
| 12 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 291求調(diào)劑 +4 | Y-cap 2026-03-29 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 298求調(diào)劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
[考研] 壓國家一區(qū)線，求導(dǎo)師收留，有恩必謝！ +7 | 迷人的哈哈 2026-03-28 | 7/350 |
|
|
|
[考研]
 一志愿廈門大學(xué)化學(xué)學(xué)碩307求調(diào)劑
+10
|
y7czhao 2026-03-26 | 10/500 |
|
|
|
[考研] 352分 化工與材料 +5 | 海納百川Ly 2026-03-27 | 5/250 |
|
|
|
[考研] 279 分 求調(diào)劑 +4 | 睡個(gè)好覺_16 2026-03-24 | 4/200 |
|
|
|
[考研] 312求調(diào)劑 +9 | 上岸吧ZJY 2026-03-22 | 13/650 |
|
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
|
[考研] 081700 調(diào)劑 267分 +11 | 迷人的哈哈 2026-03-23 | 11/550 |
|
|
|
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4 | 小魚愛有機(jī) 2026-03-25 | 4/200 |
|
|
|
[考研] 材料與化工328分調(diào)劑 +6 | 。，。，。，。i 2026-03-23 | 6/300 |
|
|
|
[考研] 材料專碩 335 分求調(diào)劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|
|
|
[考研] 0854電子信息求調(diào)劑 +7 | α____ 2026-03-22 | 9/450 |
|
|
|
[考研] 293求調(diào)劑 +7 | 加一一九 2026-03-24 | 7/350 |
|
|
|
[考研] 材料調(diào)劑 +3 | iwinso 2026-03-23 | 3/150 |
|
|
|
[考研] 277分求調(diào)劑，跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|
|
|
[考博] 26申博自薦 +3 | whh869393 2026-03-24 | 3/150 |
|
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|