计算化学公社

标题: 求助计算WSe2加SOC后，自洽计算总能变得特别大 [打印本页]

作者
Author: lcq845402558 时间: 2017-2-16 20:24
标题: 求助计算WSe2加SOC后，自洽计算总能变得特别大
本帖最后由 lcq845402558 于 2017-2-16 20:46 编辑

在下刚开始接触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埃。
希望各位老师给予指导！！

谢谢！

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

作者
Author: lcq845402558 时间: 2017-2-17 12:50

卡开发发发表于 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。

希望给予指导。

作者
Author: 卡开发发 时间: 2017-2-17 13:39
本帖最后由卡开发发于 2017-2-17 13:43 编辑

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即可，但不见得所有元素都有最新的。

作者
Author: lcq845402558 时间: 2017-2-17 16:13

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

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

感谢指导。

欢迎光临计算化学公社 (http://ccc.keinsci.com/)