 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2873 | 回復(fù): 12 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
劉仕晨木蟲 (正式寫手)
|
[求助]
關(guān)于pwscf的ibrav=0和不為0時的問題
|
|
|
pwscf我是初學(xué)者，由于之前用的是vasp，所以在寫in文件的時候也沒有多想，一直是用的ibrav=0，然后下面給出晶格常數(shù)，然后用xcrysden，檢查自己建的結(jié)構(gòu)有沒有問題，后來師兄偶然的一次聊天，我突然發(fā)現(xiàn)同樣的體系當(dāng)ibrav不為0時會得到完全不同的結(jié)果。我用example8里的Ni做了測試，因為我當(dāng)初是因為學(xué)習(xí)畫費米面而學(xué)習(xí)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看會得到同樣的結(jié)構(gòu)，應(yīng)該是沒有設(shè)錯的，但是計算的結(jié)果，卻是4倍的關(guān)系，無論是體積還是能量，這個我就不能理解了，我覺得我第二種設(shè)法，沒有設(shè)成超胞啊，不過最離譜的是費米面，算出來的完全不同，這讓我懷疑我之前算的費米面都是錯的，我就很奇怪了，一樣的結(jié)果不是應(yīng)該得到差不多的費米面么，為什么會這樣呢？不知道大家有沒有注意過這個問題，希望知道的能指導(dǎo)一下，謝謝！ |
回復(fù)此樓» 猜你喜歡
一志愿西電085401數(shù)一英一299求調(diào)劑 六級521
已經(jīng)有3人回復(fù)
-
生物考研337分求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑：085600材料與化工，考材科基，總分319
已經(jīng)有14人回復(fù)
-
285求調(diào)劑
已經(jīng)有10人回復(fù)
-
化學(xué)工程085602 305分求調(diào)劑
已經(jīng)有29人回復(fù)
-
理學(xué)07化學(xué) 303求調(diào)劑
已經(jīng)有5人回復(fù)
-
081200-11408-276學(xué)碩求調(diào)劑
已經(jīng)有3人回復(fù)
-
土木304求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿中國科學(xué)院大學(xué)265求調(diào)劑
已經(jīng)有7人回復(fù)
-
286求調(diào)劑
已經(jīng)有8人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 新手求助pwscf的進行磁性運算和LDA+U運算
已經(jīng)有35人回復(fù)
- 編譯pwscf提示找不到F77的庫
已經(jīng)有10人回復(fù)
- Pwscf編譯出錯
已經(jīng)有4人回復(fù)
- pwscf計算integrated charge出錯及不收斂
已經(jīng)有13人回復(fù)
- pwscf并行計算時總是隔一段時間重新開始
已經(jīng)有5人回復(fù)
- PWscf計算homo lumo或者valence conduction band，問題已解決
已經(jīng)有14人回復(fù)
- 自旋體系的聲子計算應(yīng)該用什么軟件呢？
已經(jīng)有5人回復(fù)
- [求助]PWscf安裝
已經(jīng)有11人回復(fù)
- PWscf中CP的溫度單位
已經(jīng)有5人回復(fù)
- 【已解決】【求助】為何PWscf不能得到類似于CASTEP那樣的收斂測試曲線呢？
已經(jīng)有14人回復(fù)
- 【求助】求助pwscf一些參數(shù)的物理意義和設(shè)定經(jīng)驗
已經(jīng)有13人回復(fù)
- 【求助】pwscf結(jié)構(gòu)優(yōu)化結(jié)果
已經(jīng)有4人回復(fù)
- 【原創(chuàng)】PWSCF中的nbnd和degauss
已經(jīng)有11人回復(fù)
- 【求助】關(guān)于PWSCF介電常數(shù)和聲子譜的計算
已經(jīng)有7人回復(fù)
- 【求助】pwscf中的問題
已經(jīng)有4人回復(fù)
- 【求助】PWSCF計算介電極化時為何出現(xiàn)如下提示？
已經(jīng)有19人回復(fù)
- 【分享】pwscf模擬水分子
已經(jīng)有5人回復(fù)
- 【求助】用PWSCF中的LDA+U做反鐵磁絕緣體系的弛豫計算
已經(jīng)有23人回復(fù)
- 【求助】請教PWSCF晶體建模和自洽計算，不收斂
已經(jīng)有21人回復(fù)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
sunyang1988: 金幣+1, 謝謝交流 2013-06-13 00:19:48
劉仕晨: 金幣+10, ★★★★★最佳答案 2013-06-13 16:19:28
感謝參與，應(yīng)助指數(shù) +1
sunyang1988: 金幣+1, 謝謝交流 2013-06-13 00:19:48
劉仕晨: 金幣+10, ★★★★★最佳答案 2013-06-13 16:19:28
|
設(shè)置ibrav=2時，相當(dāng)于POSCAR設(shè)置為： 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 這是建立的結(jié)構(gòu)為FCC的primitive cell。 而LZ給的第二個情況，則是典型的FCC的conventional cell，原子數(shù)量是primitive cell的四倍，自然做出來的總能是四倍的關(guān)系。至于費米面，我想LZ說的差別大不是在于費米面本身的形貌，而是顯示出的原子模型吧。使用ibrav=2做費米面時，如果沒記錯，似乎應(yīng)該是顯示primitive cell的模型的。 |
lsdwrk
3樓2013-06-12 20:45:35
劉仕晨
木蟲 (正式寫手)
- 應(yīng)助: 30 (小學(xué)生)
- 金幣: 3276.8
- 散金: 6
- 紅花: 11
- 帖子: 319
- 在線: 179.2小時
- 蟲號: 2201333
- 注冊: 2012-12-23
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
2樓2013-06-12 19:44:12
樓上解釋很到位啊，果然大神啊 4樓2013-06-12 22:23:27
劉仕晨
木蟲 (正式寫手)
- 應(yīng)助: 30 (小學(xué)生)
- 金幣: 3276.8
- 散金: 6
- 紅花: 11
- 帖子: 319
- 在線: 179.2小時
- 蟲號: 2201333
- 注冊: 2012-12-23
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
|
，謝謝，謝謝，這樣我就明白了，我把費米面的圖進行了給出，第一個是ibrav=2的，第二個是ibrav=0的，也就是說其實這兩張圖都是對的，雖然ibrav=0的時候沒有ibrav=2時好看。不過我糾結(jié)的是，那計算體系的時候到底應(yīng)該是用那種建模的方法呢？比如我算的那些體系的費米面怎么辦，我都是用ibrav=0算的，還是需要全部重算么，而且比如CrO2，我之前是用的ibrav=0的，我發(fā)現(xiàn)ibrav=6時，和ibrav=0時，給出的原子數(shù)是相同的，不知道對不對： &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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +5 | hanamiko 2026-03-30 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 考研生物與醫(yī)藥調(diào)劑 +3 | 鐵憨憨123425 2026-03-31 | 3/150 |
|
|
|
[考研] 環(huán)境工程 085701，267求調(diào)劑 +11 | minht 2026-03-29 | 11/550 |
|
|
|
[考研] 本科211總分289，08工學(xué)真心求調(diào)劑 +3 | utopiaE 2026-03-30 | 3/150 |
|
|
|
[考研] 346求調(diào)劑 一志愿070303有機化學(xué) +7 | 蘿卜燉青菜 2026-03-28 | 8/400 |
|
|
|
[考研] 085600 286分 材料求調(diào)劑 +11 | 麻辣魷魚 2026-03-27 | 12/600 |
|
|
|
[考研] 抱歉 +3 | 田洪有 2026-03-30 | 3/150 |
|
|
|
[有機交流] 考研調(diào)劑 +8 | watb 2026-03-26 | 8/400 |
|
|
|
[考研] 一志愿211，335分，0856，求調(diào)劑院校和導(dǎo)師 +7 | 傾____蕭 2026-03-27 | 8/400 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-26 | 7/350 |
|
|
|
[考研] 349求調(diào)劑 +6 | 李木子啊哈哈 2026-03-25 | 6/300 |
|
|
|
[考研] 本科新能源科學(xué)與工程，一志愿華理能動285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
|
[考研] 0856，材料與化工321分求調(diào)劑 +12 | 大饞小子 2026-03-27 | 13/650 |
|
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
|
[考研] 085405 考的11408求各位老師帶走 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
|
[考研] 315分求調(diào)劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +4 | 零八# 2026-03-27 | 4/200 |
|
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
|
[考博] 26申博 +3 | 加油沖�。� 2026-03-26 | 3/150 |
|
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|