 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 4633 | 回復(fù): 19 | |||
c_chi木蟲 (正式寫手)
|
[求助]
pwscf收斂速度，參數(shù)調(diào)節(jié)與贗勢 已有4人參與
|
||
|
大家好，我現(xiàn)在使用pwscf一段時間了，也碰到各種問題，有些能夠解決，但還有些棘手的問題。 我需要計算MnO2（3個原子），發(fā)現(xiàn)scf很難收斂，感覺比siesta和vasp慢，而且有時不容易收斂，會出現(xiàn)很多“ c_bands: 2 eigenvalues not converged”，試過調(diào)節(jié)'ecutwfc', 'ecutrho', 'nbnd', 'mixing_beta' 等參數(shù)，也換過贗勢，但是還是不行。所以我有以下一些問題： 1。在pwscf中有哪些參數(shù)可以加速scf收斂，大家是如何調(diào)節(jié)的，調(diào)節(jié)后對精度有沒有影響？對于離子步如何加速收斂？ 2. “ c_bands: 2 eigenvalues not converged“是對結(jié)果沒有太大影響，但是過多也會影響速度，特別是加電荷后這種warnning特別多，是結(jié)構(gòu)造成的還是有什么參數(shù)不恰當(dāng)？ 3.如何去選擇贗勢？大家一般使用官方的PSL還是其他版本，或者是自己生成？對于這種過渡金屬氧化物是選擇pbe還是pw91合適 4.還有個問題，就是我如果想監(jiān)控算AIMD過程中的電荷密度，大家有什么便捷的方法嗎？ 在此謝謝各位了�。� 我的參數(shù)如下： &control calculation='relax', restart_mode='from_scratch', pseudo_dir = './', outdir= './', prefix='MnO2', tprnfor = .true. tstress = .true. etot_conv_thr = 1.0D-6, !Default: 1.0D-4 (a.u) forc_conv_thr = 1.0D-4, !Default: 1.0D-3 (a.u) nstep = 100, !number of ionic + electronic steps verbosity = 'high' / &system ibrav = 4, a=2.75, c=18, nat = 3, ntyp=2, ecutwfc = 30, ecutwfc = 300, occupations='smearing', smearing='mp', degauss=0.03, nspin=1, / &electrons electron_maxstep=100, diagonalization = 'david' mixing_mode = 'plain' mixing_beta = 0.7, conv_thr = 1.D-6, / &ions ion_dynamics = 'bfgs' pot_extrapolation = 'second_order' wfc_extrapolation = 'second_order' / ATOMIC_SPECIES Mn 54.938 Mn.pbe-sp-van_mit.UPF O 15.999 O.pbe-n-rrkjus_psl.1.0.0.UPF ATOMIC_POSITIONS {crystal} Mn 0.00000 0.00000 0.57292 O 0.66679 0.33353 0.51863 O 0.33321 0.66647 0.62720 K_POINTS {automatic} 8 8 1 0 0 0 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 第一性 |
» 猜你喜歡
各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩
已經(jīng)有12人回復(fù)
-
359求調(diào)劑
已經(jīng)有7人回復(fù)
-
070300化學(xué)專業(yè)279調(diào)劑
已經(jīng)有10人回復(fù)
-
環(huán)境工程297分求調(diào)劑一志愿杭高院
已經(jīng)有9人回復(fù)
-
一志愿武理材料工程302調(diào)劑環(huán)化或化工
已經(jīng)有9人回復(fù)
-
309分085801求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩
已經(jīng)有5人回復(fù)
-
材料調(diào)劑
已經(jīng)有4人回復(fù)
-
材料調(diào)劑
已經(jīng)有3人回復(fù)
-
311（085601）求調(diào)劑
已經(jīng)有12人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- pwscf計算收斂太慢
已經(jīng)有5人回復(fù)
- pwscf的AIMD參數(shù)
已經(jīng)有8人回復(fù)
- 關(guān)于VASP如何選擇贗勢的問題
已經(jīng)有4人回復(fù)
- 說說vasp中的參數(shù)NPAR的使用
已經(jīng)有3人回復(fù)
- 新手求助pwscf的進(jìn)行磁性運算和LDA+U運算
已經(jīng)有35人回復(fù)
- 尋找一個帖子：參照VASP贗勢信息生成pwscf PAW贗勢的帖子！
已經(jīng)有10人回復(fù)
- 關(guān)于pwscf測試各參數(shù)收斂的交流
已經(jīng)有5人回復(fù)
- pwscf中的動力學(xué)矩陣的單位
已經(jīng)有9人回復(fù)
- PWSCF贗勢混用問題
已經(jīng)有9人回復(fù)
- 關(guān)于PWscf中的位置坐標(biāo)設(shè)置問題
已經(jīng)有15人回復(fù)
- vasp 中的贗勢文件（）
已經(jīng)有3人回復(fù)
- 求PWSCF使用說明書
已經(jīng)有3人回復(fù)
- 【求助】pwscf 贗勢選擇
已經(jīng)有6人回復(fù)
- 【已解決】【求助】為何PWscf不能得到類似于CASTEP那樣的收斂測試曲線呢？
已經(jīng)有14人回復(fù)
- 【求助】求助pwscf一些參數(shù)的物理意義和設(shè)定經(jīng)驗
已經(jīng)有13人回復(fù)
- 【求助】求VASP的贗勢庫。。
已經(jīng)有11人回復(fù)
- 【求助】請問vasp有B3LYP贗勢嗎？
已經(jīng)有7人回復(fù)
- 【討論】討論使用vasp超軟贗勢參數(shù)生成pwscf超軟贗勢文件的可行性
已經(jīng)有6人回復(fù)
- 【求助】PWSCF中將abinit贗勢轉(zhuǎn)化成pwscf贗勢能的程序怎么用？
已經(jīng)有7人回復(fù)
- 【求助】PWSCF中Ti的贗勢是不是有問題呢？
已經(jīng)有17人回復(fù)
- 【求助】請教PWSCF晶體建模和自洽計算，不收斂
已經(jīng)有21人回復(fù)
- 【討論】什么是贗勢？為什么用贗勢？
已經(jīng)有18人回復(fù)
c_chi
木蟲 (正式寫手)
- 應(yīng)助: 17 (小學(xué)生)
- 金幣: 2449.7
- 散金: 228
- 紅花: 2
- 帖子: 551
- 在線: 493.6小時
- 蟲號: 1253512
- 注冊: 2011-04-02
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
10樓2014-06-03 16:01:33
c_chi
木蟲 (正式寫手)
- 應(yīng)助: 17 (小學(xué)生)
- 金幣: 2449.7
- 散金: 228
- 紅花: 2
- 帖子: 551
- 在線: 493.6小時
- 蟲號: 1253512
- 注冊: 2011-04-02
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
2樓2014-06-03 09:58:24
1244464635
金蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1069.1
- 散金: 100
- 帖子: 333
- 在線: 92.5小時
- 蟲號: 1239162
- 注冊: 2011-03-20
- 專業(yè): 分離過程
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
c_chi: 金幣+10, ★★★很有幫助 2014-06-03 14:56:42
c_chi: 金幣+10, ★★★很有幫助 2014-06-03 14:56:42
|
1.MnO2應(yīng)該是有磁性的，你的input中nspin=1明顯不合理，應(yīng)設(shè)置為nspin=2； 2.我自己的經(jīng)驗，過渡金屬氧化物選擇pbe贗勢較好； 3.etot_conv_thr和forc_conv_thr用default值即可吧，沒必要加高精度； 4.我不太清楚你這算的是slab還是bulk，如果算的是slab，那么設(shè)置mixing_mode = 'local-TF'比較好，如果算的是bulk，那么k-points一般用k×k×k； 5. 在relax的前幾個scf可以使electron_maxstep大一些，我一般設(shè)為400； 6.可以嘗試給Mn加上U，使其3d電子局域化； 7.你的最后一個問題，我真的不懂。 以上是我的一些個人看法，僅供參考，你有什么想法大家繼續(xù)討論。我也是用pwscf的，也算過一些Mn摻雜催化劑，應(yīng)該有一些可以討論的東西。 |
3樓2014-06-03 11:37:35
gongjijun
木蟲 (正式寫手)
- 應(yīng)助: 11 (小學(xué)生)
- 金幣: 4161.8
- 散金: 48
- 紅花: 5
- 帖子: 507
- 在線: 846小時
- 蟲號: 935467
- 注冊: 2009-12-30
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
4樓2014-06-03 12:59:09
c_chi
木蟲 (正式寫手)
- 應(yīng)助: 17 (小學(xué)生)
- 金幣: 2449.7
- 散金: 228
- 紅花: 2
- 帖子: 551
- 在線: 493.6小時
- 蟲號: 1253512
- 注冊: 2011-04-02
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
5樓2014-06-03 13:27:58
gongjijun
木蟲 (正式寫手)
- 應(yīng)助: 11 (小學(xué)生)
- 金幣: 4161.8
- 散金: 48
- 紅花: 5
- 帖子: 507
- 在線: 846小時
- 蟲號: 935467
- 注冊: 2009-12-30
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
6樓2014-06-03 13:30:32
c_chi
木蟲 (正式寫手)
- 應(yīng)助: 17 (小學(xué)生)
- 金幣: 2449.7
- 散金: 228
- 紅花: 2
- 帖子: 551
- 在線: 493.6小時
- 蟲號: 1253512
- 注冊: 2011-04-02
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
7樓2014-06-03 14:14:44
c_chi
木蟲 (正式寫手)
- 應(yīng)助: 17 (小學(xué)生)
- 金幣: 2449.7
- 散金: 228
- 紅花: 2
- 帖子: 551
- 在線: 493.6小時
- 蟲號: 1253512
- 注冊: 2011-04-02
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
|
謝謝你的指點哈，我這里還有幾點疑問 1. 我覺得在做結(jié)構(gòu)優(yōu)化的時候可以不加磁性，然后在算能量和其他性質(zhì)的時候再加磁性參數(shù)，這樣計算量會小一些 2. 我是用psl生成的贗勢，有兩種kjpaw和rrkjus，應(yīng)該是分別對應(yīng)PAW和超軟贗勢，用的時候PAW贗勢報錯用不了，超軟贗勢有比較慢，而且出現(xiàn)很多c_band的警告，很難對角化，不知道你是如何選擇的。另外，我生成的Mn的pbe贗勢只有PAW贗勢，沒有超軟的，不知道是不是生成出錯了。關(guān)于贗勢最后一點，可不可以通過混合使用贗勢在保證一定精度的情況下來提高速度，如果可以，一般有什么規(guī)律？ 3. 確實我的精度高了點，用默認(rèn)的對角化不了的情況少了很多。不過感覺計算速度比siesta和VASP慢一些，可能是siesta是LCAO的方法，而VASP有比較成熟的PAW勢 4. 我算的是slab，那么mixing_mode的選擇是不是不會影響到最后的精度，只是和計算收斂速度關(guān)系比較大？ 5. 在relax里面可以分段設(shè)置scf最大步數(shù)嗎？還是直接將electron_maxstep設(shè)置大一點。我看到有的帖子說，scf會在100步內(nèi)收斂，如果收斂不了一般就是參數(shù)有問題 6. 加U應(yīng)該對結(jié)構(gòu)和能量的影響比較小，在算帯隙的時候用的比較多，我是不是一般可以不加考慮 7.那你有用過pwscf做第一原理分子動力學(xué)嗎？相對VASP和siesta，速度和準(zhǔn)確度你覺得怎么樣。不知道是不是參數(shù)的問題，我試過幾個原子的小體系還是蠻快的，一旦大了起來就受不了了。也可能因為AIMD本來就一個好計算量大的工作，看來還是得優(yōu)化一下參數(shù) 8. mixing_beta和degauss還有沒有調(diào)節(jié)的可能去提高速度，我試了下感覺不明顯，但是我在siesta里面減小類似的參數(shù)感覺收斂速度明顯提高 9.除了我上面的一些參數(shù)，你一般還會調(diào)節(jié)其他哪些參數(shù) 問題多了些哈。。。對于pwscf我覺得功能是蠻強大，就是贗勢的問題和沒有向其他軟件一樣有比較系統(tǒng)的說明文檔，只有參數(shù)說明，這有點頭疼。 |
8樓2014-06-03 14:56:30
c_chi
木蟲 (正式寫手)
- 應(yīng)助: 17 (小學(xué)生)
- 金幣: 2449.7
- 散金: 228
- 紅花: 2
- 帖子: 551
- 在線: 493.6小時
- 蟲號: 1253512
- 注冊: 2011-04-02
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
9樓2014-06-03 15:05:31
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 c_chi 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑 +16 | Micky11223 2026-03-25 | 17/850 |
|
|---|---|---|---|---|
|
[考研] 375求調(diào)劑 +7 | 雨夏整夜 2026-03-29 | 7/350 |
|
|
|
[考研] 張芳銘-中國農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +9 | 手機用戶 2026-03-26 | 9/450 |
|
|
|
[考研] 336材料求調(diào)劑 +10 | 陳瀅瑩 2026-03-26 | 12/600 |
|
|
|
[考研] 322求調(diào)劑 +3 | 熹僖XX 2026-03-31 | 3/150 |
|
|
|
[考研] 263求調(diào)劑 +3 | DDDDuu 2026-03-27 | 3/150 |
|
|
|
[考研] 本科211安全工程，初試290分，求調(diào)劑 +3 | 2719846834 2026-03-28 | 3/150 |
|
|
|
[考研] 085600 295分求調(diào)劑 +13 | W55j 2026-03-30 | 15/750 |
|
|
|
[考研] 22408 359分調(diào)劑 +4 | Qshers 2026-03-27 | 8/400 |
|
|
|
[考研] 085601一志愿中山大學(xué)深圳材料工程330求調(diào)劑 +5 | pipiver 2026-03-30 | 5/250 |
|
|
|
[考研] 11408總分309，一志愿東南大學(xué)求調(diào)劑，不挑專業(yè) +5 | 天賦帶到THU 2026-03-29 | 6/300 |
|
|
|
[考研] 生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績338 +7 | YYYYYNNNNN 2026-03-26 | 9/450 |
|
|
|
[考研] 295材料工程專碩求調(diào)劑 +10 | 1428151015 2026-03-27 | 10/500 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +6 | 丹青奶蓋 2026-03-26 | 8/400 |
|
|
|
[考研] 324求調(diào)劑 +9 | hanamiko 2026-03-26 | 11/550 |
|
|
|
[考研] 2026年華南師范大學(xué)歡迎化學(xué)，化工，生物，生醫(yī)工等專業(yè)優(yōu)秀學(xué)子加入！ +3 | llss0711 2026-03-28 | 6/300 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +6 | 崔wj 2026-03-26 | 6/300 |
|
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
|
[考研] 085600，材料與化工321分調(diào)劑 +4 | 大饞小子 2026-03-27 | 6/300 |
|
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|