计算化学公社

标题: Pt(111)晶面SCF不收敛 [打印本页]

作者
Author: Freeman 时间: 2020-11-1 22:10
标题: Pt(111)晶面SCF不收敛
本帖最后由 Freeman 于 2020-11-1 22:47 编辑

老师们好。
我最近要做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

(, 下载次数 Times of downloads: 4)

作者
Author: wzkchem5 时间: 2020-11-1 22:30
输出文件也给一下？

作者
Author: Freeman 时间: 2020-11-1 22:43

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

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

作者
Author: 卡开发发 时间: 2020-11-1 23:31
本帖最后由卡开发发于 2020-11-1 23:33 编辑

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

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

作者
Author: Freeman 时间: 2020-11-2 18:33

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

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

作者
Author: 卡开发发 时间: 2020-11-2 19:02

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

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

作者
Author: Freeman 时间: 2020-11-3 21:24

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

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

作者
Author: 卡开发发 时间: 2020-11-3 22:39

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

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

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