| 查看: 5256 | 回复: 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 |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性 |
» 猜你喜欢
课题组会帮你通过考核么？
已经有6人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有9人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- pwscf计算收敛太慢
已经有5人回复
- pwscf的AIMD参数
已经有8人回复
- 关于VASP如何选择赝势的问题
已经有4人回复
- 说说vasp中的参数NPAR的使用
已经有3人回复
- 新手求助pwscf的进行磁性运算和LDA+U运算
已经有35人回复
- 寻找一个帖子：参照VASP赝势信息生成pwscf PAW赝势的帖子！
已经有10人回复
- 关于pwscf测试各参数收敛的交流
已经有5人回复
- pwscf中的动力学矩阵的单位
已经有9人回复
- PWSCF赝势混用问题
已经有9人回复
- 关于PWscf中的位置坐标设置问题
已经有15人回复
- vasp 中的赝势文件（）
已经有3人回复
- 求PWSCF使用说明书
已经有3人回复
- 【求助】pwscf 赝势选择
已经有6人回复
- 【已解决】【求助】为何PWscf不能得到类似于CASTEP那样的收敛测试曲线呢？
已经有14人回复
- 【求助】求助pwscf一些参数的物理意义和设定经验
已经有13人回复
- 【求助】求VASP的赝势库。。
已经有11人回复
- 【求助】请问vasp有B3LYP赝势吗？
已经有7人回复
- 【讨论】讨论使用vasp超软赝势参数生成pwscf超软赝势文件的可行性
已经有6人回复
- 【求助】PWSCF中将abinit赝势转化成pwscf赝势能的程序怎么用？
已经有7人回复
- 【求助】PWSCF中Ti的赝势是不是有问题呢？
已经有17人回复
- 【求助】请教PWSCF晶体建模和自洽计算，不收敛
已经有21人回复
- 【讨论】什么是赝势？为什么用赝势？
已经有18人回复
1244464635
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1084.1
- 散金: 100
- 帖子: 336
- 在线: 92.6小时
- 虫号: 1239162
- 注册: 2011-03-20
- 专业: 分离过程
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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
c_chi
木虫 (正式写手)
- 应助: 17 (小学生)
- 金币: 2449.7
- 散金: 228
- 红花: 2
- 帖子: 551
- 在线: 493.6小时
- 虫号: 1253512
- 注册: 2011-04-02
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
2楼2014-06-03 09:58:24
gongjijun
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 4161.8
- 散金: 48
- 红花: 5
- 帖子: 507
- 在线: 846小时
- 虫号: 935467
- 注册: 2009-12-30
- 专业: 凝聚态物性 II ：电子结构
4楼2014-06-03 12:59:09
c_chi
木虫 (正式写手)
- 应助: 17 (小学生)
- 金币: 2449.7
- 散金: 228
- 红花: 2
- 帖子: 551
- 在线: 493.6小时
- 虫号: 1253512
- 注册: 2011-04-02
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
5楼2014-06-03 13:27:58
