24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

c_chi

木虫 (正式写手)

应助: 17 (小学生)
金币: 2448.7
散金: 228
红花: 2
帖子: 551
在线: 493.6小时
虫号: 1253512
注册: 2011-04-02
性别: GG
专业: 凝聚态物性I:结构、力学和

[求助] 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

回复此楼

1楼2014-05-29 15:14:35

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

1244464635

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 1105.3
帖子: 314
在线: 90.2小时
虫号: 1239162
注册: 2011-03-20
专业: 分离过程

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
c_chi: 金币+10, ★★★很有帮助 2014-06-03 14:56:42

引用回帖:

2楼: Originally posted by c_chi at 2014-06-03 09:58:24
这些问题还是没有解决，不知道有没没有了解的，谢谢各位赐教一下！！

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掺杂催化剂，应该有一些可以讨论的东西。