Pt(111)晶面SCF不收敛

Freeman

老师们好。
我最近要做Pt(111)晶面slab模型的优化，但是第一步的SCF就不收敛。我调了各种参数，如mixing_beta，截断能等，但都没有用。能量确实一直在下降，但是每次下降的幅度一直不变，accuracy也几乎不变。最后电子步就超过了100。我觉得单纯地增大迭代次数应该不管用。所以请大家帮忙看一下输入文件，请教一下问题出在哪里。
&CONTROL
    calculation   = "relax"
    prefix        = "ptsur"
    pseudo_dir    = "/workfs/yczhang/MyCrys/pseudo/SSSP_efficiency_pseudos"
    outdir        = "./outdir"
/
&SYSTEM
    a                         =  1.13595e+01
    c                         =  3.46375e+01
    degauss                   =  1.00000e-02
    ecutwfc                   =  20
    ecutrho                   =  200
    ibrav                     = 4
    nat                       = 48
    ntyp                      = 1
    occupations               = "smearing"
    smearing                  = "gaussian"
    nspin                     = 2
    starting_magnetization(1) = 0.3
/
&ELECTRONS
    mixing_beta               = 0.3D0
/
&IONS
/
&CELL
/
K_POINTS automatic
3 3 1 0 0 0
ATOMIC_SPECIES
Pt    195.07800  pt_pbe_v1.4.uspp.F.UPF
ATOMIC_POSITIONS {angstrom}
Pt      2.839818   1.639635   15.000116 0 0 0
Pt      0.000057   3.279172   17.318634 0 0 0
Pt      0.000000   0.000000   19.637615
Pt      1.419881   4.099039   15.000116 0 0 0
Pt      1.419994   0.819769   17.318634 0 0 0
Pt     -1.419938   2.459404   19.637615
Pt     -1.419994   4.099039   15.000116 0 0 0
Pt      4.259869   0.819769   17.318634 0 0 0
Pt      1.419938   2.459404   19.637615
Pt     -0.000057   1.639635   15.000116 0 0 0
Pt      2.839932   3.279172   17.318634 0 0 0
Pt      2.839875   0.000000   19.637615
Pt      8.519568   1.639635   15.000116 0 0 0
Pt      5.679807   3.279172   17.318634 0 0 0
Pt      5.679750   0.000000   19.637615
Pt      7.099631   4.099039   15.000116 0 0 0
Pt      7.099744   0.819769   17.318634 0 0 0
Pt      4.259813   2.459404   19.637615
Pt      4.259756   4.099039   15.000116 0 0 0
Pt      9.939619   0.819769   17.318634 0 0 0
Pt      7.099688   2.459404   19.637615
Pt      5.679693   1.639635   15.000116 0 0 0
Pt      8.519682   3.279172   17.318634 0 0 0
Pt      8.519625   0.000000   19.637615
Pt     -0.000057   6.558443   15.000116 0 0 0
Pt     -2.839818   8.197980   17.318634 0 0 0
Pt     -2.839875   4.918808   19.637615
Pt     -1.419994   9.017847   15.000116 0 0 0
Pt     -1.419881   5.738576   17.318634 0 0 0
Pt     -4.259813   7.378212   19.637615
Pt     -4.259869   9.017847   15.000116 0 0 0
Pt      1.419994   5.738576   17.318634 0 0 0
Pt     -1.419938   7.378212   19.637615
Pt     -2.839932   6.558443   15.000116 0 0 0
Pt      0.000057   8.197980   17.318634 0 0 0
Pt      0.000000   4.918808   19.637615
Pt      5.679693   6.558443   15.000116 0 0 0
Pt      2.839932   8.197980   17.318634 0 0 0
Pt      2.839875   4.918808   19.637615
Pt      4.259756   9.017847   15.000116 0 0 0
Pt      4.259869   5.738576   17.318634 0 0 0
Pt      1.419937   7.378212   19.637615
Pt      1.419881   9.017847   15.000116 0 0 0
Pt      7.099744   5.738576   17.318634 0 0 0
Pt      4.259813   7.378212   19.637615
Pt      2.839818   6.558443   15.000116 0 0 0
Pt      5.679807   8.197980   17.318634 0 0 0
Pt      5.679750   4.918808   19.637615




ptsur.pw.in (2.88 KB, 下载次数 Times of downloads: 4)

2020-11-1 22:09 上传 Uploaded

点击下载
Click to download

wzkchem5

输出文件也给一下？

Freeman

wzkchem5 发表于 2020-11-1 22:30
输出文件也给一下？

现在电脑不在身边。能量确实一直在下降，但是每次下降的幅度一直不变，accuracy也几乎不变。最后电子步就超过了100。我觉得单纯地增大迭代次数应该不管用。

卡开发发

Freeman 发表于 2020-11-1 22:43
现在电脑不在身边。能量确实一直在下降，但是每次下降的幅度一直不变，accuracy也几乎不变。最后电子步就 ...

Pt 111应该是非磁性体系，没必要ispin=2，加上磁矩初猜还给了0.3，这可能会导致难收敛。
另外，ecutwfc=20 Ry~270eV，这很可能偏低。

Freeman

卡开发发 发表于 2020-11-1 23:31
Pt 111应该是非磁性体系，没必要ispin=2，加上磁矩初猜还给了0.3，这可能会导致难收敛。
另外，ecutwfc= ...

谢谢您。去掉nspin后很快就收敛了。之前在文献上看到别人做Pt(111)时加了自旋极化，所以我也加了。看来不能盲信文献。

卡开发发

Freeman 发表于 2020-11-2 18:33
谢谢您。去掉nspin后很快就收敛了。之前在文献上看到别人做Pt(111)时加了自旋极化，所以我也加了。看来不 ...

原理上其实是影响到密度的初猜导致SCF不稳定，如果牵扯到吸附物带磁性的话，可以给Pt的磁矩猜0。

Freeman

卡开发发 发表于 2020-11-2 19:02
原理上其实是影响到密度的初猜导致SCF不稳定，如果牵扯到吸附物带磁性的话，可以给Pt的磁矩猜0。

你说得对。那篇文献就是讲Pt吸附有磁性的氧气的。

卡开发发

Freeman 发表于 2020-11-3 21:24
你说得对。那篇文献就是讲Pt吸附有磁性的氧气的。

O的自旋设置下，Pt按0猜。