| 查看: 4385 | 回复: 19 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
c_chi木虫 (正式写手)
|
[求助]
pwscf收敛速度，参数调节与赝势已有4人参与
|
|
|
大家好，我现在使用pwscf一段时间了，也碰到各种问题，有些能够解决，但还有些棘手的问题。 我需要计算MnO2（3个原子），发现scf很难收敛，感觉比siesta和vasp慢，而且有时不容易收敛，会出现很多“ c_bands: 2 eigenvalues not converged”，试过调节'ecutwfc', 'ecutrho', 'nbnd', 'mixing_beta' 等参数，也换过赝势，但是还是不行。所以我有以下一些问题： 1。在pwscf中有哪些参数可以加速scf收敛，大家是如何调节的，调节后对精度有没有影响？对于离子步如何加速收敛？ 2. “ c_bands: 2 eigenvalues not converged“是对结果没有太大影响，但是过多也会影响速度，特别是加电荷后这种warnning特别多，是结构造成的还是有什么参数不恰当？ 3.如何去选择赝势？大家一般使用官方的PSL还是其他版本，或者是自己生成？对于这种过渡金属氧化物是选择pbe还是pw91合适 4.还有个问题，就是我如果想监控算AIMD过程中的电荷密度，大家有什么便捷的方法吗？ 在此谢谢各位了！！ 我的参数如下： &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 |
回复此楼» 猜你喜欢
导师想让我从独立一作变成了共一第一
已经有9人回复
-
博士读完未来一定会好吗
已经有23人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有11人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
-
申请2026年博士
已经有6人回复
12楼2014-06-03 16:12:43
c_chi
木虫 (正式写手)
- 应助: 17 (小学生)
- 金币: 2448.7
- 散金: 228
- 红花: 2
- 帖子: 551
- 在线: 493.6小时
- 虫号: 1253512
- 注册: 2011-04-02
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
2楼2014-06-03 09:58:24
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
c_chi: 金币+10, ★★★很有帮助 2014-06-03 14:56:42
c_chi: 金币+10, ★★★很有帮助 2014-06-03 14:56:42
|
1.MnO2应该是有磁性的，你的input中nspin=1明显不合理，应设置为nspin=2； 2.我自己的经验，过渡金属氧化物选择pbe赝势较好； 3.etot_conv_thr和forc_conv_thr用default值即可吧，没必要加高精度； 4.我不太清楚你这算的是slab还是bulk，如果算的是slab，那么设置mixing_mode = 'local-TF'比较好，如果算的是bulk，那么k-points一般用k×k×k； 5. 在relax的前几个scf可以使electron_maxstep大一些，我一般设为400； 6.可以尝试给Mn加上U，使其3d电子局域化； 7.你的最后一个问题，我真的不懂。 以上是我的一些个人看法，仅供参考，你有什么想法大家继续讨论。我也是用pwscf的，也算过一些Mn掺杂催化剂，应该有一些可以讨论的东西。 |
3楼2014-06-03 11:37:35
gongjijun
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 4161.8
- 散金: 48
- 红花: 5
- 帖子: 507
- 在线: 846小时
- 虫号: 935467
- 注册: 2009-12-30
- 专业: 凝聚态物性 II ：电子结构
4楼2014-06-03 12:59:09