 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 4635 | 回復(fù): 19 | |||
c_chi主管區(qū)長
|
[求助]
pwscf收斂速度，參數(shù)調(diào)節(jié)與贗勢 已有4人參與
|
||
|
大家好，我現(xiàn)在使用pwscf一段時(shí)間了，也碰到各種問題，有些能夠解決，但還有些棘手的問題。 我需要計(jì)算MnO2（3個原子），發(fā)現(xiàn)scf很難收斂，感覺比siesta和vasp慢，而且有時(shí)不容易收斂，會出現(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ù)此樓» 收錄本帖的淘帖專輯推薦
 第一性 |
» 猜你喜歡
271求調(diào)劑
已經(jīng)有12人回復(fù)
-
322求調(diào)劑
已經(jīng)有3人回復(fù)
-
環(huán)境工程 085701，267求調(diào)劑
已經(jīng)有15人回復(fù)
-
材料0856 英一數(shù)二 323 求調(diào)劑
已經(jīng)有3人回復(fù)
-
301求調(diào)劑
已經(jīng)有3人回復(fù)
-
348求調(diào)劑
已經(jīng)有3人回復(fù)
-
350求調(diào)劑
已經(jīng)有3人回復(fù)
-
309分085801求調(diào)劑
已經(jīng)有7人回復(fù)
-
生物學(xué)308分求調(diào)劑（一志愿華東師大）
已經(jīng)有3人回復(fù)
-
求調(diào)劑：一志愿：南京大學(xué) 專業(yè)：0705 總分320 ，本科985，四六級已過
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- pwscf計(jì)算收斂太慢
已經(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)行磁性運(yùn)算和LDA+U運(yùn)算
已經(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)驗(yàn)
已經(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ì)算，不收斂
已經(jīng)有21人回復(fù)
- 【討論】什么是贗勢？為什么用贗勢？
已經(jīng)有18人回復(fù)
c_chi
版主
- 應(yīng)助: 17 (小學(xué)生)
- 金幣: 2449.7
- 散金: 228
- 紅花: 2
- 帖子: 551
- 在線: 493.6小時(shí)
- 蟲號: 1253512
- 注冊: 2011-04-02
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
10樓2014-06-03 16:01:33
c_chi
實(shí)習(xí)版主
- 應(yīng)助: 17 (小學(xué)生)
- 金幣: 2449.7
- 散金: 228
- 紅花: 2
- 帖子: 551
- 在線: 493.6小時(shí)
- 蟲號: 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小時(shí)
- 蟲號: 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)驗(yàn)，過渡金屬氧化物選擇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
實(shí)習(xí)版主
- 應(yīng)助: 11 (小學(xué)生)
- 金幣: 4161.8
- 散金: 48
- 紅花: 5
- 帖子: 507
- 在線: 846小時(shí)
- 蟲號: 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小時(shí)
- 蟲號: 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小時(shí)
- 蟲號: 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小時(shí)
- 蟲號: 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小時(shí)
- 蟲號: 1253512
- 注冊: 2011-04-02
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
|
謝謝你的指點(diǎn)哈，我這里還有幾點(diǎn)疑問 1. 我覺得在做結(jié)構(gòu)優(yōu)化的時(shí)候可以不加磁性，然后在算能量和其他性質(zhì)的時(shí)候再加磁性參數(shù)，這樣計(jì)算量會小一些 2. 我是用psl生成的贗勢，有兩種kjpaw和rrkjus，應(yīng)該是分別對應(yīng)PAW和超軟贗勢，用的時(shí)候PAW贗勢報(bào)錯用不了，超軟贗勢有比較慢，而且出現(xiàn)很多c_band的警告，很難對角化，不知道你是如何選擇的。另外，我生成的Mn的pbe贗勢只有PAW贗勢，沒有超軟的，不知道是不是生成出錯了。關(guān)于贗勢最后一點(diǎn)，可不可以通過混合使用贗勢在保證一定精度的情況下來提高速度，如果可以，一般有什么規(guī)律？ 3. 確實(shí)我的精度高了點(diǎn)，用默認(rèn)的對角化不了的情況少了很多。不過感覺計(jì)算速度比siesta和VASP慢一些，可能是siesta是LCAO的方法，而VASP有比較成熟的PAW勢 4. 我算的是slab，那么mixing_mode的選擇是不是不會影響到最后的精度，只是和計(jì)算收斂速度關(guān)系比較大？ 5. 在relax里面可以分段設(shè)置scf最大步數(shù)嗎？還是直接將electron_maxstep設(shè)置大一點(diǎn)。我看到有的帖子說，scf會在100步內(nèi)收斂，如果收斂不了一般就是參數(shù)有問題 6. 加U應(yīng)該對結(jié)構(gòu)和能量的影響比較小，在算帯隙的時(shí)候用的比較多，我是不是一般可以不加考慮 7.那你有用過pwscf做第一原理分子動力學(xué)嗎？相對VASP和siesta，速度和準(zhǔn)確度你覺得怎么樣。不知道是不是參數(shù)的問題，我試過幾個原子的小體系還是蠻快的，一旦大了起來就受不了了。也可能因?yàn)锳IMD本來就一個好計(jì)算量大的工作，看來還是得優(yōu)化一下參數(shù) 8. mixing_beta和degauss還有沒有調(diào)節(jié)的可能去提高速度，我試了下感覺不明顯，但是我在siesta里面減小類似的參數(shù)感覺收斂速度明顯提高 9.除了我上面的一些參數(shù)，你一般還會調(diào)節(jié)其他哪些參數(shù) 問題多了些哈。。。對于pwscf我覺得功能是蠻強(qiáng)大，就是贗勢的問題和沒有向其他軟件一樣有比較系統(tǒng)的說明文檔，只有參數(shù)說明，這有點(diǎn)頭疼。 |
8樓2014-06-03 14:56:30
c_chi
版主
- 應(yīng)助: 17 (小學(xué)生)
- 金幣: 2449.7
- 散金: 228
- 紅花: 2
- 帖子: 551
- 在線: 493.6小時(shí)
- 蟲號: 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ā)表 | |
|---|---|---|---|---|
|
[考研] 環(huán)境工程 085701，267求調(diào)劑 +14 | minht 2026-03-29 | 15/750 |
|
|---|---|---|---|---|
|
|
[考研] 350求調(diào)劑 +3 | 阿佳～ 2026-03-31 | 3/150 |
|
|
|
[考研] 考研調(diào)劑 +6 | Amber00 2026-03-31 | 6/300 |
|
|
[考研] 08工科275分求調(diào)劑 +6 | AaAa7420 2026-03-31 | 6/300 |
|
|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +10 | 李李不服輸 2026-03-25 | 10/500 |
|
|
|
[考研] 349求調(diào)劑 +4 | zwjjjjjj 2026-03-31 | 4/200 |
|
|
|
[考研] 考研材料工程351分調(diào)劑 +4 | 整個好的 2026-03-31 | 4/200 |
|
|
|
[考研] 085600，321分求調(diào)劑 +7 | 大饞小子 2026-03-31 | 7/350 |
|
|
|
[考研] 一志愿華東師范大學(xué)有機(jī)化學(xué)專業(yè)，初試351分，復(fù)試被刷求調(diào)劑! +9 | 真名有冰 2026-03-29 | 10/500 |
|
|
|
[考研] 290分調(diào)劑求助 +10 | 吉祥止止陳 2026-03-25 | 10/500 |
|
|
|
[考研] 0710生物學(xué)求調(diào)劑！ +5 | 敘述文 2026-03-31 | 5/250 |
|
|
|
[考研] 329求調(diào)劑，一志愿西北工業(yè)大學(xué)，材料工程（085601） +6 | 小小機(jī)靈蟲 2026-03-29 | 12/600 |
|
|
|
[考研] 本科211安全工程，初試290分，求調(diào)劑 +3 | 2719846834 2026-03-28 | 3/150 |
|
|
|
[考研] 南京大學(xué)化學(xué)調(diào)劑 +11 | 景隨風(fēng) 2026-03-29 | 16/800 |
|
|
|
[有機(jī)交流]
 甲基亞磺磺酸鈉和甲基磺酸酯反應(yīng)機(jī)理
 10+3
|
kaobao456 2026-03-29 | 4/200 |
|
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +16 | 李李不服輸 2026-03-26 | 16/800 |
|
|
|
[考研] 291求調(diào)劑 +5 | Y-cap 2026-03-29 | 6/300 |
|
|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +7 | 閆旭東 2026-03-28 | 9/450 |
|
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|