 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2867 | 回復(fù): 12 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
劉仕晨木蟲 (正式寫手)
|
[求助]
關(guān)于pwscf的ibrav=0和不為0時(shí)的問(wèn)題
|
|
|
pwscf我是初學(xué)者，由于之前用的是vasp，所以在寫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í)畫費(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看會(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ù)此樓» 猜你喜歡
317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑
已經(jīng)有11人回復(fù)
-
332求調(diào)劑
已經(jīng)有12人回復(fù)
-
材料專碩 085600求調(diào)劑
已經(jīng)有5人回復(fù)
-
0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可
已經(jīng)有14人回復(fù)
-
283求調(diào)劑（080500）
已經(jīng)有14人回復(fù)
-
材料工程專碩求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料專碩調(diào)劑
已經(jīng)有9人回復(fù)
-
求調(diào)劑 有機(jī)化學(xué)考研356分
已經(jīng)有11人回復(fù)
-
環(huán)境科學(xué)與工程334分求調(diào)劑
已經(jīng)有8人回復(fù)
-
考研調(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í)間重新開始
已經(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)助: 30 (小學(xué)生)
- 金幣: 3276.8
- 散金: 6
- 紅花: 11
- 帖子: 319
- 在線: 179.2小時(shí)
- 蟲號(hào): 2201333
- 注冊(cè): 2012-12-23
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
9樓2013-06-13 18:40:20
劉仕晨
木蟲 (正式寫手)
- 應(yīng)助: 30 (小學(xué)生)
- 金幣: 3276.8
- 散金: 6
- 紅花: 11
- 帖子: 319
- 在線: 179.2小時(shí)
- 蟲號(hào): 2201333
- 注冊(cè): 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的四倍，自然做出來(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
樓上解釋很到位啊，果然大神啊 4樓2013-06-12 22:23:27
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 296求調(diào)劑 +10 | 彼岸t 2026-03-29 | 10/500 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿中南大學(xué)化學(xué)0703總分337求調(diào)劑 +6 | niko- 2026-03-27 | 6/300 |
|
|
|
[考研] 327求調(diào)劑 +4 | 小卡不卡. 2026-03-29 | 4/200 |
|
|
[考研] 085600 286分 材料求調(diào)劑 +10 | 麻辣魷魚 2026-03-27 | 11/550 |
|
|
|
[考研] 307求調(diào)劑 +9 | 超級(jí)伊昂大王 2026-03-24 | 10/500 |
|
|
|
[考研] 08工科求調(diào)劑286 +3 | tgs_001 2026-03-28 | 3/150 |
|
|
|
[考研] 085701求調(diào)劑初試286分 +4 | secret0328 2026-03-28 | 4/200 |
|
|
|
[考研] 295材料工程專碩求調(diào)劑 +7 | 1428151015 2026-03-27 | 7/350 |
|
|
|
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +4 | 小張zxy張 2026-03-26 | 8/400 |
|
|
|
[考研] 一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑 +7 | Micky11223 2026-03-25 | 7/350 |
|
|
|
[考研] 085602 化工專碩 338分 求調(diào)劑 +12 | 路癡小琪 2026-03-27 | 12/600 |
|
|
|
[考研] 394求調(diào)劑 +3 | 好事多磨靜候佳?/a> 2026-03-26 | 5/250 |
|
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
|
[考研] 085600，材料與化工321分調(diào)劑 +4 | 大饞小子 2026-03-27 | 6/300 |
|
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
|
[考研] 343求調(diào)劑 +4 | 贈(zèng)我一本書 2026-03-23 | 4/200 |
|
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|
|
|
[考研] 一志愿吉大化學(xué)322求調(diào)劑 +4 | 17501029541 2026-03-23 | 6/300 |
|