求助计算WSe2加SOC后，自洽计算总能变得特别大

lcq845402558

在下刚开始接触QE。正在重复文献。
求助计算加SOC的单层WSe2后，自洽计算总能变得特别大。
输入文件：
&CONTROL
                 calculation = 'scf',
                restart_mode = 'from_scratch' ,
                      outdir = './' ,
                  pseudo_dir = './ ,
               etot_conv_thr = 1.0d-4 ,
               forc_conv_thr = 2.5d-4 ,
/
&SYSTEM
                       ibrav = 0,
                   celldm(1) = 6.28724957,
                         nat = 3,
                        ntyp = 2,
                     ecutwfc = 50 ,
                     ecutrho = 400 ,
                       nosym = .false. ,
                    noncolin = .true.
                    lspinorb = .true.
/
&ELECTRONS
            electron_maxstep = 100,
                    conv_thr = 1.D-10 ,
/
&IONS
                ion_dynamics = 'bfgs' ,
/
&CELL
               cell_dynamics = 'bfgs' ,
/
ATOMIC_SPECIES
    Se   78.96000  Se.rel-pbe-rrkjus.UPF
    W   183.85000  W.rel-pbe-spn-rrkjus.UPF
CELL_PARAMETERS
   0.999217236   0.000000000  -0.000000000
  -0.499608618   0.865347511  -0.000000000
   0.000000000   0.000000000   7.027222139
ATOMIC_POSITIONS (crystal)
W        0.666666667   0.333333333   0.156845540
Se       0.333333333   0.666666667   0.084951844
Se       0.333333333   0.666666667   0.228739171
K_POINTS automatic
  6 6 1   0 0 0


输出文件中的能量：
scf.pwo:     total energy              =    -213.83800240 Ry
scf.pwo:     total energy              =    -213.85547507 Ry
scf.pwo:     total energy              =   61765.46529166 Ry
scf.pwo:     total energy              =   65917.91273216 Ry
scf.pwo:     total energy              =   38057.62754787 Ry
scf.pwo:     total energy              =   62515.64751687 Ry
scf.pwo:     total energy              =   80283.01487620 Ry
scf.pwo:     total energy              =  103947.29660675 Ry
scf.pwo:     total energy              =   79334.20426000 Ry
scf.pwo:     total energy              =  102232.48951395 Ry
scf.pwo:     total energy              =  -96965.60956842 Ry


后面的kill了，不知道哪里出了问题。
不加SOC的WSe2；
和加了SOC的WS2。
都job done， 能带与文献一致。

截断能和K点均测试过，赝势为超软，真空层设的为20埃。
希望各位老师给予指导！！

谢谢！

卡开发发

可以把赝势换成PSLib0.23的版本看看。实在不行在考虑试试降低一些混合系数mixing_beta，看看能否减小振荡，实在不行把mixing_ndim也降下来。

lcq845402558

卡开发发 发表于 2017-2-16 21:02
可以把赝势换成PSLib0.23的版本看看。实在不行在考虑试试降低一些混合系数mixing_beta，看看能否减小振荡， ...

感谢您的回答，W的加rel的超软赝势确实不太好找，所以我将mixing_beta调到了0.1。
scf成功job done。但是算bands的时候出错了。

出错信息：
Error in routine cdiaghg (105):
     S matrix not positive definite

不知道是不是赝势的问题，小木虫有人说提高截断能，但是我的超软赝势截断能已经达到50Ry。

希望给予指导。

卡开发发

lcq845402558 发表于 2017-2-17 12:50
感谢您的回答，W的加rel的超软赝势确实不太好找，所以我将mixing_beta调到了0.1。
scf成功job done。但 ...

也可以提高一下看看，对于W应该50 Ry还不算大。基函数尺寸过小有可能会导致S非正。
至于赝势的话你可以在这边下载到新的版本，http://www.qe-forge.org/gf/proje ... p;frs_package_id=41，进入相应的目录make ps即可，但不见得所有元素都有最新的。

lcq845402558

卡开发发 发表于 2017-2-17 13:39
也可以提高一下看看，对于W应该50 Ry还不算大。基函数尺寸过小有可能会导致S非正。
至于赝势的话你可以 ...

谢谢老师，问题已经解决。
我在system下添加了
     occupations = 'smearing',
     smearing='gaussian',
     degauss=0.01,
之前单层有算dos，所以occupations选的是tetrahedra。
虽然不知道什么原因，但是解决了。

感谢指导。