 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2869 | 回復(fù): 12 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
劉仕晨木蟲 (正式寫手)
|
[求助]
關(guān)于pwscf的ibrav=0和不為0時(shí)的問題
|
|
|
pwscf我是初學(xué)者，由于之前用的是vasp，所以在寫in文件的時(shí)候也沒有多想，一直是用的ibrav=0，然后下面給出晶格常數(shù)，然后用xcrysden，檢查自己建的結(jié)構(gòu)有沒有問題，后來師兄偶然的一次聊天，我突然發(fā)現(xiàn)同樣的體系當(dāng)ibrav不為0時(shí)會得到完全不同的結(jié)果。我用example8里的Ni做了測試，因?yàn)槲耶?dāng)初是因?yàn)閷W(xué)習(xí)畫費(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 然后是我寫的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看會得到同樣的結(jié)構(gòu)，應(yīng)該是沒有設(shè)錯(cuò)的，但是計(jì)算的結(jié)果，卻是4倍的關(guān)系，無論是體積還是能量，這個(gè)我就不能理解了，我覺得我第二種設(shè)法，沒有設(shè)成超胞啊，不過最離譜的是費(fèi)米面，算出來的完全不同，這讓我懷疑我之前算的費(fèi)米面都是錯(cuò)的，我就很奇怪了，一樣的結(jié)果不是應(yīng)該得到差不多的費(fèi)米面么，為什么會這樣呢？不知道大家有沒有注意過這個(gè)問題，希望知道的能指導(dǎo)一下，謝謝！ |
回復(fù)此樓» 猜你喜歡
277跪求調(diào)劑
已經(jīng)有11人回復(fù)
-
一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑
已經(jīng)有18人回復(fù)
-
求調(diào)劑
已經(jīng)有6人回復(fù)
-
0703化學(xué)
已經(jīng)有14人回復(fù)
-
面上5B能上會嗎？
已經(jīng)有7人回復(fù)
-
0703本科鄭州大學(xué)求調(diào)劑
已經(jīng)有7人回復(fù)
-
286分調(diào)劑
已經(jīng)有3人回復(fù)
-
0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可
已經(jīng)有15人回復(fù)
-
317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑
已經(jīng)有11人回復(fù)
-
332求調(diào)劑
已經(jīng)有12人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 新手求助pwscf的進(jìn)行磁性運(yùn)算和LDA+U運(yùn)算
已經(jīng)有35人回復(fù)
- 編譯pwscf提示找不到F77的庫
已經(jīng)有10人回復(fù)
- Pwscf編譯出錯(cuò)
已經(jīng)有4人回復(fù)
- pwscf計(jì)算integrated charge出錯(cuò)及不收斂
已經(jīng)有13人回復(fù)
- pwscf并行計(jì)算時(shí)總是隔一段時(shí)間重新開始
已經(jīng)有5人回復(fù)
- PWscf計(jì)算homo lumo或者valence conduction band，問題已解決
已經(jīng)有14人回復(fù)
- 自旋體系的聲子計(jì)算應(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)驗(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中的問題
已經(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ù)
- 【求助】請教PWSCF晶體建模和自洽計(jì)算，不收斂
已經(jīng)有21人回復(fù)
樓上解釋很到位啊，果然大神啊 
4樓2013-06-12 22:23:27
劉仕晨
木蟲 (正式寫手)
- 應(yīng)助: 30 (小學(xué)生)
- 金幣: 3276.8
- 散金: 6
- 紅花: 11
- 帖子: 319
- 在線: 179.2小時(shí)
- 蟲號: 2201333
- 注冊: 2012-12-23
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
2樓2013-06-12 19:44:12
【答案】應(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的四倍，自然做出來的總能是四倍的關(guān)系。至于費(fèi)米面，我想LZ說的差別大不是在于費(fèi)米面本身的形貌，而是顯示出的原子模型吧。使用ibrav=2做費(fèi)米面時(shí)，如果沒記錯(cuò)，似乎應(yīng)該是顯示primitive cell的模型的。 |
» 本帖已獲得的紅花（最新10朵）
lsdwrk3樓2013-06-12 20:45:35
劉仕晨
木蟲 (正式寫手)
- 應(yīng)助: 30 (小學(xué)生)
- 金幣: 3276.8
- 散金: 6
- 紅花: 11
- 帖子: 319
- 在線: 179.2小時(shí)
- 蟲號: 2201333
- 注冊: 2012-12-23
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
|
，謝謝，謝謝，這樣我就明白了，我把費(fèi)米面的圖進(jìn)行了給出，第一個(gè)是ibrav=2的，第二個(gè)是ibrav=0的，也就是說其實(shí)這兩張圖都是對的，雖然ibrav=0的時(shí)候沒有ibrav=2時(shí)好看。不過我糾結(jié)的是，那計(jì)算體系的時(shí)候到底應(yīng)該是用那種建模的方法呢？比如我算的那些體系的費(fèi)米面怎么辦，我都是用ibrav=0算的，還是需要全部重算么，而且比如CrO2，我之前是用的ibrav=0的，我發(fā)現(xiàn)ibrav=6時(shí)，和ibrav=0時(shí)，給出的原子數(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ā)表 | |
|---|---|---|---|---|
|
[考研] 297求調(diào)劑 +13 | 田洪有 2026-03-26 | 13/650 |
|
|---|---|---|---|---|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-26 | 7/350 |
|
|
|
[考研] 307求調(diào)劑 +9 | 超級伊昂大王 2026-03-24 | 10/500 |
|
|
[考研] 求調(diào)劑，一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +7 | @taotao 2026-03-29 | 7/350 |
|
|
|
[考研] 一志愿211 初試270分 求調(diào)劑 +8 | 谷雨上岸 2026-03-23 | 9/450 |
|
|
|
[考研] 320分，材料與化工專業(yè)，求調(diào)劑 +9 | 一定上岸aaa 2026-03-27 | 13/650 |
|
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +7 | hanamiko 2026-03-27 | 7/350 |
|
|
|
[考研] 081200-314 +3 | LILIQQ 2026-03-27 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +4 | 零八# 2026-03-27 | 4/200 |
|
|
|
[考研] 08開頭275求調(diào)劑 +4 | 拉誰不重要 2026-03-26 | 4/200 |
|
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
|
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
|
[考研] 359求調(diào)劑 +4 | 王了個(gè)楠 2026-03-25 | 4/200 |
|
|
|
[考研]
材料調(diào)劑
 5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
|
[考研] 309求調(diào)劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
|
[考研] 0854AI CV方向招收調(diào)劑 +4 | 章小魚567 2026-03-23 | 4/200 |
|
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|