 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 2866 | 回復: 12 | ||
劉仕晨木蟲 (正式寫手)
|
[求助]
關于pwscf的ibrav=0和不為0時的問題
|
|
|
pwscf我是初學者，由于之前用的是vasp，所以在寫in文件的時候也沒有多想，一直是用的ibrav=0，然后下面給出晶格常數，然后用xcrysden，檢查自己建的結構有沒有問題，后來師兄偶然的一次聊天，我突然發(fā)現同樣的體系當ibrav不為0時會得到完全不同的結果。我用example8里的Ni做了測試，因為我當初是因為學習畫費米面而學習pwscf的。 首先是官方給出的輸入文件，ibrav=2，顯示如下（省略了60個K點）： &control calculation='scf' restart_mode='from_scratch', prefix='ni', pseudo_dir = '/home/liushichen/tools/pwscf/espresso-4.3.2/pseudo/', outdir='/home/liushichen/work/pwscf/Ni/tmp/' / &system ibrav=2, celldm(1) =6.48, nat=1, ntyp=1, nspin = 2, starting_magnetization(1)=0.7, ecutwfc = 24.0, ecutrho = 288.0, occupations='smearing', smearing='methfessel-paxton', degauss=0.02 / &electrons conv_thr = 1.0e-10 mixing_beta = 0.7 / ATOMIC_SPECIES Ni 58.69 Ni.pz-nd-rrkjus.UPF ATOMIC_POSITIONS Ni 0.0 0.0 0.0 K_POINTS 60 然后是我寫的ibrav=0時的in文件，顯示如下： &control calculation='scf' restart_mode='from_scratch', prefix='ni', pseudo_dir = '/home/liushichen/tools/pwscf/espresso-4.3.2/pseudo/', outdir='/home/liushichen/work/pwscf/Ni/tmp0/' / &system ibrav=0, nat=4, ntyp=1, nspin = 2, starting_magnetization(1)=0.7, ecutwfc = 24.0, ecutrho = 288.0, occupations='smearing', smearing='methfessel-paxton', degauss=0.02 / &electrons conv_thr = 1.0e-10 mixing_beta = 0.7 / CELL_PARAMETERS cubic 6.480000000 0.000000000 0.000000000 0.000000000 6.480000000 0.000000000 0.000000000 0.000000000 6.480000000 ATOMIC_SPECIES Ni 58.69 Ni.pz-nd-rrkjus.UPF ATOMIC_POSITIONS Ni 0.0 0.0 0.0 Ni 0.5 0.5 0.0 Ni 0.0 0.5 0.5 Ni 0.5 0.0 0.5 K_POINTS 60 這兩個輸入文件用xcrysden看會得到同樣的結構，應該是沒有設錯的，但是計算的結果，卻是4倍的關系，無論是體積還是能量，這個我就不能理解了，我覺得我第二種設法，沒有設成超胞啊，不過最離譜的是費米面，算出來的完全不同，這讓我懷疑我之前算的費米面都是錯的，我就很奇怪了，一樣的結果不是應該得到差不多的費米面么，為什么會這樣呢？不知道大家有沒有注意過這個問題，希望知道的能指導一下，謝謝！ |
回復此樓» 猜你喜歡
309求調劑
已經有9人回復
-
環(huán)境科學與工程334分求調劑
已經有4人回復
-
085404求調劑，總分309，本科經歷較為豐富
已經有6人回復
-
285求調劑
已經有5人回復
-
求化學調劑
已經有7人回復
-
材料工程專碩求調劑
已經有4人回復
-
317求調劑
已經有10人回復
-
318求調劑
已經有4人回復
-
化學308分調劑
已經有3人回復
-
一志愿河北工業(yè)大學0817化工278分求調劑
已經有18人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 新手求助pwscf的進行磁性運算和LDA+U運算
已經有35人回復
- 編譯pwscf提示找不到F77的庫
已經有10人回復
- Pwscf編譯出錯
已經有4人回復
- pwscf計算integrated charge出錯及不收斂
已經有13人回復
- pwscf并行計算時總是隔一段時間重新開始
已經有5人回復
- PWscf計算homo lumo或者valence conduction band，問題已解決
已經有14人回復
- 自旋體系的聲子計算應該用什么軟件呢？
已經有5人回復
- [求助]PWscf安裝
已經有11人回復
- PWscf中CP的溫度單位
已經有5人回復
- 【已解決】【求助】為何PWscf不能得到類似于CASTEP那樣的收斂測試曲線呢？
已經有14人回復
- 【求助】求助pwscf一些參數的物理意義和設定經驗
已經有13人回復
- 【求助】pwscf結構優(yōu)化結果
已經有4人回復
- 【原創(chuàng)】PWSCF中的nbnd和degauss
已經有11人回復
- 【求助】關于PWSCF介電常數和聲子譜的計算
已經有7人回復
- 【求助】pwscf中的問題
已經有4人回復
- 【求助】PWSCF計算介電極化時為何出現如下提示？
已經有19人回復
- 【分享】pwscf模擬水分子
已經有5人回復
- 【求助】用PWSCF中的LDA+U做反鐵磁絕緣體系的弛豫計算
已經有23人回復
- 【求助】請教PWSCF晶體建模和自洽計算，不收斂
已經有21人回復
【答案】應助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應助指數 +1
sunyang1988: 金幣+1, 謝謝交流 2013-06-13 00:19:48
劉仕晨: 金幣+10, ★★★★★最佳答案 2013-06-13 16:19:28
感謝參與，應助指數 +1
sunyang1988: 金幣+1, 謝謝交流 2013-06-13 00:19:48
劉仕晨: 金幣+10, ★★★★★最佳答案 2013-06-13 16:19:28
|
設置ibrav=2時，相當于POSCAR設置為： lattice_param 0.0 0.5 0.5 0.5 0.0 0.5 0.5 0.5 0.0 1 0.0 0.0 0.0 這是建立的結構為FCC的primitive cell。 而LZ給的第二個情況，則是典型的FCC的conventional cell，原子數量是primitive cell的四倍，自然做出來的總能是四倍的關系。至于費米面，我想LZ說的差別大不是在于費米面本身的形貌，而是顯示出的原子模型吧。使用ibrav=2做費米面時，如果沒記錯，似乎應該是顯示primitive cell的模型的。 |
» 本帖已獲得的紅花（最新10朵）
lsdwrk
3樓2013-06-12 20:45:35
劉仕晨
木蟲 (正式寫手)
- 應助: 30 (小學生)
- 金幣: 3276.8
- 散金: 6
- 紅花: 11
- 帖子: 319
- 在線: 179.2小時
- 蟲號: 2201333
- 注冊: 2012-12-23
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
2樓2013-06-12 19:44:12
樓上解釋很到位啊，果然大神啊 4樓2013-06-12 22:23:27
劉仕晨
木蟲 (正式寫手)
- 應助: 30 (小學生)
- 金幣: 3276.8
- 散金: 6
- 紅花: 11
- 帖子: 319
- 在線: 179.2小時
- 蟲號: 2201333
- 注冊: 2012-12-23
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
|
，謝謝，謝謝，這樣我就明白了，我把費米面的圖進行了給出，第一個是ibrav=2的，第二個是ibrav=0的，也就是說其實這兩張圖都是對的，雖然ibrav=0的時候沒有ibrav=2時好看。不過我糾結的是，那計算體系的時候到底應該是用那種建模的方法呢？比如我算的那些體系的費米面怎么辦，我都是用ibrav=0算的，還是需要全部重算么，而且比如CrO2，我之前是用的ibrav=0的，我發(fā)現ibrav=6時，和ibrav=0時，給出的原子數是相同的，不知道對不對： &SYSTEM ibrav = 6, celldm(1) = 8.3503, celldm(3) = 0.65897, nat = 6, ntyp = 2, ecutwfc = 300 , ecutrho = 1200 , occupations = 'smearing' , degauss = 0.02 , smearing = 'gaussian' , nspin = 2 , starting_magnetization(1) = 0.5, lda_plus_u = .true. , Hubbard_U(1) = 3, / &ELECTRONS conv_thr = 1.0e-10 , mixing_beta = 0.7 , / ATOMIC_SPECIES Cr 52.00000 Cr.pz-hgh.UPF O 16.00000 O.pz-hgh.UPF ATOMIC_POSITIONS crystal Cr 0.000000000 0.000000000 0.000000000 Cr 0.500000000 0.500000000 0.500000000 O 0.303000000 0.303000000 0.000000000 O 0.697000000 0.697000000 0.000000000 O 0.803000000 0.197000000 0.500000000 O 0.197000000 0.803000000 0.500000000 ibrav=0： &SYSTEM ibrav = 0, nat = 6, ntyp = 2, ecutwfc = 300 , ecutrho = 1200 , occupations = 'smearing' , degauss = 0.02 , smearing = 'gaussian' , nspin = 2 , starting_magnetization(1) = 0.5, lda_plus_u = .true. , Hubbard_U(1) = 3, / &ELECTRONS conv_thr = 1.0e-10 , mixing_beta = 0.7 , / CELL_PARAMETERS cubic 8.350300000 0.000000000 0.000000000 0.000000000 8.350300000 0.000000000 0.000000000 0.000000000 5.502600000 ATOMIC_SPECIES Cr 52.00000 Cr.pz-hgh.UPF O 16.00000 O.pz-hgh.UPF ATOMIC_POSITIONS crystal Cr 0.000000000 0.000000000 0.000000000 Cr 0.500000000 0.500000000 0.500000000 O 0.303000000 0.303000000 0.000000000 O 0.697000000 0.697000000 0.000000000 O 0.803000000 0.197000000 0.500000000 O 0.197000000 0.803000000 0.500000000 Bands_FS_down.png Bands_FS_up.png |
5樓2013-06-13 13:17:14
6樓2013-06-13 14:44:29
劉仕晨
木蟲 (正式寫手)
- 應助: 30 (小學生)
- 金幣: 3276.8
- 散金: 6
- 紅花: 11
- 帖子: 319
- 在線: 179.2小時
- 蟲號: 2201333
- 注冊: 2012-12-23
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
7樓2013-06-13 16:19:21
一個下午愣是沒能擠進去的怨念飄過…… 8樓2013-06-13 16:31:12
劉仕晨
木蟲 (正式寫手)
- 應助: 30 (小學生)
- 金幣: 3276.8
- 散金: 6
- 紅花: 11
- 帖子: 319
- 在線: 179.2小時
- 蟲號: 2201333
- 注冊: 2012-12-23
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
9樓2013-06-13 18:40:20
weiyongkai
木蟲 (職業(yè)作家)
- 應助: 13 (小學生)
- 金幣: 3192
- 散金: 209
- 紅花: 3
- 帖子: 3186
- 在線: 228.1小時
- 蟲號: 1104094
- 注冊: 2010-09-20
- 專業(yè): 凝聚態(tài)物性I:結構、力學和
10樓2013-07-06 11:26:33
相關版塊跳轉
我要訂閱樓主 劉仕晨 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經歷 +8 | cccchenso 2026-03-29 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 一志愿北京工業(yè)大學，324分求調劑 +6 | 零八# 2026-03-28 | 6/300 |
|
|
|
[考研] 289求調劑 +5 | BrightLL 2026-03-29 | 5/250 |
|
|
|
[考研] 070305高分子化學與物理 304分求調劑 +12 | c297914 2026-03-28 | 12/600 |
|
|
|
[考研] 一志愿鄭州大學，080500學碩，總分317分求調劑 +8 | 舉個栗子oi 2026-03-24 | 9/450 |
|
|
|
[考研] 0856材料化工調劑 總分330 +14 | zhubinhao 2026-03-27 | 14/700 |
|
|
|
[考研] 壓國家一區(qū)線，求導師收留，有恩必謝！ +7 | 迷人的哈哈 2026-03-28 | 7/350 |
|
|
|
[考研] 復試調劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
|
[考研] 070300求調劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 266分，求材料冶金能源化工等調劑 +7 | 哇呼哼呼哼 2026-03-27 | 9/450 |
|
|
|
[考研] 286求調劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
|
[考研] 085602 307分 求調劑 +7 | 不知道叫什么！ 2026-03-26 | 7/350 |
|
|
|
[考研] 一志愿南昌大學324求調劑 +7 | hanamiko 2026-03-27 | 7/350 |
|
|
|
[考研] 291求調劑 +6 | HanBeiNingZC 2026-03-24 | 6/300 |
|
|
|
[考研] 一志愿 南京郵電大學 288分 材料考研 求調劑 +3 | jl0720 2026-03-26 | 3/150 |
|
|
|
[考研] 309求調劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
|
[考研] 材料與化工304求B區(qū)調劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
|
[考研] 考研一志愿蘇州大學初始315（英一）求調劑 +3 | sbdksD 2026-03-24 | 4/200 |
|
|
|
[考研]
 材料考研調劑生
+3
|
黃粱一夢千年 2026-03-24 | 3/150 |
|
|
|
[考研] 材料專碩331求調劑 +4 | 鮮當牛 2026-03-24 | 4/200 |
|