24小时热门版块排行榜    

查看: 4437  |  回复: 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
回复此楼

» 收录本帖的淘帖专辑推荐

第一性

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

gongjijun

木虫 (正式写手)

【答案】应助回帖

★ ★
c_chi: 金币+2 2014-06-03 15:06:27
a=2.75?
你确定吗?这个可是不是以A为单位的!!
4楼2014-06-03 12:59:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gongjijun

木虫 (正式写手)

【答案】应助回帖

ibrav = 4,   a=2.75,  c=18,

三个原子的mno,这个参数你确定对头吗?
6楼2014-06-03 13:30:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 c_chi 的主题更新
信息提示
请填处理意见