 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 2870 | 回復(fù): 12 | ||
劉仕晨木蟲(chóng) (正式寫(xiě)手)
|
[求助]
關(guān)于pwscf的ibrav=0和不為0時(shí)的問(wèn)題
|
|
|
pwscf我是初學(xué)者，由于之前用的是vasp，所以在寫(xiě)in文件的時(shí)候也沒(méi)有多想，一直是用的ibrav=0，然后下面給出晶格常數(shù)，然后用xcrysden，檢查自己建的結(jié)構(gòu)有沒(méi)有問(wèn)題，后來(lái)師兄偶然的一次聊天，我突然發(fā)現(xiàn)同樣的體系當(dāng)ibrav不為0時(shí)會(huì)得到完全不同的結(jié)果。我用example8里的Ni做了測(cè)試，因?yàn)槲耶?dāng)初是因?yàn)閷W(xué)習(xí)畫(huà)費(fèi)米面而學(xué)習(xí)pwscf的。 首先是官方給出的輸入文件，ibrav=2，顯示如下（省略了60個(gè)K點(diǎn)）： &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 然后是我寫(xiě)的ibrav=0時(shí)的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 這兩個(gè)輸入文件用xcrysden看會(huì)得到同樣的結(jié)構(gòu)，應(yīng)該是沒(méi)有設(shè)錯(cuò)的，但是計(jì)算的結(jié)果，卻是4倍的關(guān)系，無(wú)論是體積還是能量，這個(gè)我就不能理解了，我覺(jué)得我第二種設(shè)法，沒(méi)有設(shè)成超胞啊，不過(guò)最離譜的是費(fèi)米面，算出來(lái)的完全不同，這讓我懷疑我之前算的費(fèi)米面都是錯(cuò)的，我就很奇怪了，一樣的結(jié)果不是應(yīng)該得到差不多的費(fèi)米面么，為什么會(huì)這樣呢？不知道大家有沒(méi)有注意過(guò)這個(gè)問(wèn)題，希望知道的能指導(dǎo)一下，謝謝！ |
回復(fù)此樓» 猜你喜歡
一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑
已經(jīng)有22人回復(fù)
-
322求調(diào)劑：一志愿湖南大學(xué) 材料與化工（085600），已過(guò)六級(jí)。
已經(jīng)有7人回復(fù)
-
332求調(diào)劑
已經(jīng)有5人回復(fù)
-
310求調(diào)劑
已經(jīng)有10人回復(fù)
-
26考研-291分-廈門(mén)大學(xué)（085601）-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑
已經(jīng)有5人回復(fù)
-
085600 材料與化工 329分求調(diào)劑
已經(jīng)有15人回復(fù)
-
化學(xué)0703 調(diào)劑 306分 一志愿211
已經(jīng)有8人回復(fù)
-
342求調(diào)劑
已經(jīng)有4人回復(fù)
-
334分 一志愿武理 材料求調(diào)劑
已經(jīng)有14人回復(fù)
-
材料科學(xué)與工程 317求調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 新手求助pwscf的進(jìn)行磁性運(yùn)算和LDA+U運(yùn)算
已經(jīng)有35人回復(fù)
- 編譯pwscf提示找不到F77的庫(kù)
已經(jīng)有10人回復(fù)
- Pwscf編譯出錯(cuò)
已經(jīng)有4人回復(fù)
- pwscf計(jì)算integrated charge出錯(cuò)及不收斂
已經(jīng)有13人回復(fù)
- pwscf并行計(jì)算時(shí)總是隔一段時(shí)間重新開(kāi)始
已經(jīng)有5人回復(fù)
- PWscf計(jì)算homo lumo或者valence conduction band，問(wèn)題已解決
已經(jīng)有14人回復(fù)
- 自旋體系的聲子計(jì)算應(yīng)該用什么軟件呢？
已經(jīng)有5人回復(fù)
- [求助]PWscf安裝
已經(jīng)有11人回復(fù)
- PWscf中CP的溫度單位
已經(jīng)有5人回復(fù)
- 【已解決】【求助】為何PWscf不能得到類似于CASTEP那樣的收斂測(cè)試曲線呢？
已經(jīng)有14人回復(fù)
- 【求助】求助pwscf一些參數(shù)的物理意義和設(shè)定經(jīng)驗(yàn)
已經(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ì)算
已經(jīng)有7人回復(fù)
- 【求助】pwscf中的問(wèn)題
已經(jīng)有4人回復(fù)
- 【求助】PWSCF計(jì)算介電極化時(shí)為何出現(xiàn)如下提示？
已經(jīng)有19人回復(fù)
- 【分享】pwscf模擬水分子
已經(jīng)有5人回復(fù)
- 【求助】用PWSCF中的LDA+U做反鐵磁絕緣體系的弛豫計(jì)算
已經(jīng)有23人回復(fù)
- 【求助】請(qǐng)教PWSCF晶體建模和自洽計(jì)算，不收斂
已經(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時(shí)，相當(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給的第二個(gè)情況，則是典型的FCC的conventional cell，原子數(shù)量是primitive cell的四倍，自然做出來(lái)的總能是四倍的關(guān)系。至于費(fèi)米面，我想LZ說(shuō)的差別大不是在于費(fèi)米面本身的形貌，而是顯示出的原子模型吧。使用ibrav=2做費(fèi)米面時(shí)，如果沒(méi)記錯(cuò)，似乎應(yīng)該是顯示primitive cell的模型的。 |
» 本帖已獲得的紅花（最新10朵）
lsdwrk
3樓2013-06-12 20:45:35
劉仕晨
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 30 (小學(xué)生)
- 金幣: 3276.8
- 散金: 6
- 紅花: 11
- 帖子: 319
- 在線: 179.2小時(shí)
- 蟲(chóng)號(hào): 2201333
- 注冊(cè): 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
劉仕晨
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 30 (小學(xué)生)
- 金幣: 3276.8
- 散金: 6
- 紅花: 11
- 帖子: 319
- 在線: 179.2小時(shí)
- 蟲(chóng)號(hào): 2201333
- 注冊(cè): 2012-12-23
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
|
，謝謝，謝謝，這樣我就明白了，我把費(fèi)米面的圖進(jìn)行了給出，第一個(gè)是ibrav=2的，第二個(gè)是ibrav=0的，也就是說(shuō)其實(shí)這兩張圖都是對(duì)的，雖然ibrav=0的時(shí)候沒(méi)有ibrav=2時(shí)好看。不過(guò)我糾結(jié)的是，那計(jì)算體系的時(shí)候到底應(yīng)該是用那種建模的方法呢？比如我算的那些體系的費(fèi)米面怎么辦，我都是用ibrav=0算的，還是需要全部重算么，而且比如CrO2，我之前是用的ibrav=0的，我發(fā)現(xiàn)ibrav=6時(shí)，和ibrav=0時(shí)，給出的原子數(shù)是相同的，不知道對(duì)不對(duì)： &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
劉仕晨
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 30 (小學(xué)生)
- 金幣: 3276.8
- 散金: 6
- 紅花: 11
- 帖子: 319
- 在線: 179.2小時(shí)
- 蟲(chóng)號(hào): 2201333
- 注冊(cè): 2012-12-23
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
7樓2013-06-13 16:19:21
一個(gè)下午愣是沒(méi)能擠進(jìn)去的怨念飄過(guò)…… 8樓2013-06-13 16:31:12
劉仕晨
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 30 (小學(xué)生)
- 金幣: 3276.8
- 散金: 6
- 紅花: 11
- 帖子: 319
- 在線: 179.2小時(shí)
- 蟲(chóng)號(hào): 2201333
- 注冊(cè): 2012-12-23
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
9樓2013-06-13 18:40:20
weiyongkai
木蟲(chóng) (職業(yè)作家)
- 應(yīng)助: 13 (小學(xué)生)
- 金幣: 3192
- 散金: 209
- 紅花: 3
- 帖子: 3186
- 在線: 228.1小時(shí)
- 蟲(chóng)號(hào): 1104094
- 注冊(cè): 2010-09-20
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
10樓2013-07-06 11:26:33
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 劉仕晨 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分 +7 | Shern—- 2026-03-28 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 332求調(diào)劑 +13 | 032500 2026-03-25 | 13/650 |
|
|
|
[考研] 348求調(diào)劑 +6 | 小懶蟲(chóng)不懶了 2026-03-28 | 6/300 |
|
|
|
[考研] 311求調(diào)劑 +10 | lin0039 2026-03-26 | 10/500 |
|
|
|
[考研] 337求調(diào)劑 +6 | 《樹(shù)》 2026-03-29 | 6/300 |
|
|
|
[考研] 085600，專業(yè)課化工原理，320分求調(diào)劑 +4 | 大饞小子 2026-03-29 | 4/200 |
|
|
|
[考研] 281求調(diào)劑 +4 | 亞克西good 2026-03-26 | 6/300 |
|
|
|
[考研] 275求調(diào)劑 +15 | Micky11223 2026-03-25 | 20/1000 |
|
|
|
[考研] 071000生物學(xué)求調(diào)劑，初試成績(jī)343 +7 | 小小甜面團(tuán) 2026-03-25 | 7/350 |
|
|
|
[考研] 295材料工程專碩求調(diào)劑 +7 | 1428151015 2026-03-27 | 7/350 |
|
|
|
[考研] 調(diào)劑 +3 | 好好讀書(shū)。 2026-03-28 | 3/150 |
|
|
|
[考研] 藥學(xué)105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
|
[考研] 292求調(diào)劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +9 | 多久上課 2026-03-27 | 9/450 |
|
|
|
[考研] 322求調(diào)劑 +4 | 我真的很想學(xué)習(xí) 2026-03-23 | 4/200 |
|
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|
|
|
[考研] 調(diào)劑求收留 +7 | 果然有我 2026-03-26 | 7/350 |
|
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
|
[考研] 302求調(diào)劑 +4 | 錦衣衛(wèi)藤椒 2026-03-25 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|