为什么使用castep里的otfg 赝势不能收敛而uspp赝势则能够收敛

535219002

449  -2.45999014E+004  5.81225924E-001   2.43126999E-004   13368.18  <-- SCF
    450  -2.45998287E+004  5.79387772E-001  -1.51499200E-003   13399.16  <-- SCF
    451  -2.45998420E+004  5.78336387E-001   2.78745776E-004   13424.24  <-- SCF
    452  -2.45998411E+004  5.79172387E-001  -1.89351957E-005   13447.37  <-- SCF
    453  -2.45998057E+004  5.79066214E-001  -7.37546413E-004   13474.68  <-- SCF
    454  -2.45997846E+004  5.78306232E-001  -4.39984053E-004   13503.74  <-- SCF
    455  -2.45998050E+004  5.79645330E-001   4.25238158E-004   13530.37  <-- SCF
    456  -2.45997688E+004  5.79417688E-001  -7.53745208E-004   13555.24  <-- SCF
    457  -2.45997876E+004  5.79926536E-001   3.91114253E-004   13577.38  <-- SCF
    458  -2.45997576E+004  5.80051222E-001  -6.25157938E-004   13601.42  <-- SCF
    459  -2.45997580E+004  5.80876559E-001   8.31335582E-006   13628.24  <-- SCF
    460  -2.45997898E+004  5.81003349E-001   6.62970207E-004   13654.29  <-- SCF
    461  -2.45998042E+004  5.81016121E-001   2.99361194E-004   13680.53  <-- SCF
    462  -2.45998067E+004  5.82281506E-001   5.21298934E-005   13706.90  <-- SCF
    463  -2.45998090E+004  5.81325260E-001   4.82933375E-005   13729.96  <-- SCF
    464  -2.45998152E+004  5.82137371E-001   1.28226502E-004   13753.06  <-- SCF
    465  -2.45998276E+004  5.84487304E-001   2.59973159E-004   13777.96  <-- SCF
    466  -2.45998295E+004  5.84922712E-001   3.96550757E-005   13799.84  <-- SCF
    467  -2.45998077E+004  5.86019444E-001  -4.54735223E-004   13823.58  <-- SCF
    468  -2.45998251E+004  5.82777140E-001   3.61056775E-004   13848.23  <-- SCF
    469  -2.45998212E+004  5.79907478E-001  -7.93713494E-005   13873.97  <-- SCF
    470  -2.45998602E+004  5.79274150E-001   8.11020984E-004   13901.30  <-- SCF
    471  -2.45998593E+004  5.78824397E-001  -1.89004828E-005   13924.77  <-- SCF
    472  -2.45998826E+004  5.78480049E-001   4.86337390E-004   13947.81  <-- SCF
    473  -2.45999043E+004  5.80090801E-001   4.51433012E-004   13974.61  <-- SCF
    474  -2.45999002E+004  5.80202328E-001  -8.43025450E-005   13997.71  <-- SCF
    475  -2.45998920E+004  5.79451836E-001  -1.70619066E-004   14021.57  <-- SCF
    476  -2.45998849E+004  5.80023467E-001  -1.47755747E-004   14046.82  <-- SCF
    477  -2.45998548E+004  5.78695964E-001  -6.27340206E-004   14072.91  <-- SCF
    478  -2.45998499E+004  5.77751738E-001  -1.03353506E-004   14095.49  <-- SCF
    479  -2.45998564E+004  5.76987690E-001   1.35575327E-004   14124.24  <-- SCF
    480  -2.45998706E+004  5.77121681E-001   2.96632661E-004   14148.26  <-- SCF
    481  -2.45998735E+004  5.76618011E-001   6.00985358E-005   14170.26  <-- SCF
    482  -2.45998714E+004  5.76268474E-001  -4.34886744E-005   14195.38  <-- SCF
    483  -2.45999056E+004  5.77246707E-001   7.11689772E-004   14222.85  <-- SCF
    484  -2.45999516E+004  5.77678929E-001   9.57725060E-004   14252.15  <-- SCF
    485  -2.45999467E+004  5.77972967E-001  -1.00072633E-004   14275.66  <-- SCF
    486  -2.45999183E+004  5.76741845E-001  -5.91957290E-004   14302.16  <-- SCF
    487  -2.45999168E+004  5.74375090E-001  -3.09428816E-005   14330.09  <-- SCF
    488  -2.45999450E+004  5.73719332E-001   5.85570200E-004   14352.99  <-- SCF
    489  -2.45999450E+004  5.78141162E-001   8.37942868E-008   14381.20  <-- SCF
    490  -2.45999099E+004  5.75287002E-001  -7.29878132E-004   14407.64  <-- SCF
    491  -2.45999135E+004  5.74869076E-001   7.35495442E-005   14433.28  <-- SCF
    492  -2.45999223E+004  5.78298133E-001   1.85202666E-004   14457.84  <-- SCF
    493  -2.45999200E+004  5.78971549E-001  -4.78574797E-005   14480.88  <-- SCF
    494  -2.45999280E+004  5.79727907E-001   1.66450958E-004   14503.31  <-- SCF
    495  -2.45999348E+004  5.82217131E-001   1.39861226E-004   14529.52  <-- SCF
    496  -2.45999775E+004  5.80768664E-001   8.91251482E-004   14554.19  <-- SCF
    497  -2.45999640E+004  5.80179402E-001  -2.81120471E-004   14578.32  <-- SCF
    498  -2.45999709E+004  5.80315230E-001   1.42714548E-004   14601.33  <-- SCF
    499  -2.45999579E+004  5.81216703E-001  -2.70624015E-004   14626.65  <-- SCF
    500  -2.45999458E+004  5.81204499E-001  -2.52522168E-004   14649.96  <-- SCF
------------------------------------------------------------------------ <-- SCF

*Warning* max. SCF cycles performed but system has not reached the groundstate.

Current total energy, E     =  -24599.93427521     eV
Current free energy (E-TS)  =  -24599.94577734     eV
(energies not corrected for finite basis set)

NB est. 0K energy (E-0.5TS)      =  -24599.94002628     eV

  ********************************************************
  Warning: electronic minimisation did not converge during
           finite basis set correction.
  ********************************************************
Checkpoint file cannot be written.
Error calculate_finite_basis : Convergence failed when doing finite basis set correction.
Current trace stack:
calculate_finite_basis_corr
check_elec_ground_state
castep
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 3
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 5
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 6
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 9
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 18
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 1
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 2
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 4
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 7
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 8
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 10
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 11
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 12
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 13
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 14
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 15
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 16
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 17
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 19

Last known process information:
===============================
Name:        castepexe.exe
Umask:        0002
State:        R (running)
Tgid:        371090
Ngid:        371090
Pid:        371090
PPid:        371086
TracerPid:        0
Uid:        1041        1041        1041        1041
Gid:        1041        1041        1041        1041
FDSize:        64
Groups:        1003 1041
VmPeak:         1330504 kB
VmSize:         1293168 kB
VmLck:               0 kB
VmPin:               0 kB
VmHWM:          938224 kB
VmRSS:          900680 kB
RssAnon:          861808 kB
RssFile:            4904 kB
RssShmem:           33968 kB
VmData:          881924 kB
VmStk:             624 kB
VmExe:          150104 kB
VmLib:           28324 kB
VmPTE:            2280 kB
VmSwap:               0 kB
Threads:        1
SigQ:        0/1541074
SigPnd:        0000000000000000
ShdPnd:        0000000000000000
SigBlk:        0000000000000000
SigIgn:        0000000000100000
SigCgt:        00000001800046ae
CapInh:        0000000000000000
CapPrm:        0000000000000000
CapEff:        0000000000000000
CapBnd:        0000001fffffffff
CapAmb:        0000000000000000
NoNewPrivs:        0
Seccomp:        0
Speculation_Store_Bypass:        thread vulnerable
Cpus_allowed:        00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000001,00000000
Cpus_allowed_list:        32
Mems_allowed:        00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,00000000,0000000f
Mems_allowed_list:        0-3
voluntary_ctxt_switches:        9703
nonvoluntary_ctxt_switches:        14909

这是使用otfg的计算结果，为什么otfg不好收敛呢，按道理说otfg不是应该更精确吗

卡开发发

精确和容不容易SCF收敛没啥直接关系。参数合理（如动能截断OTFG高于Vanderbilt老的超软赝势）的情况下，试着调整density mixing的选项。

535219002

卡开发发 发表于 2021-7-23 09:59
精确和容不容易SCF收敛没啥直接关系。参数合理（如动能截断OTFG高于Vanderbilt老的超软赝势）的情况下，试 ...

老师，我用的otfg的截断能和超软的数值大小一样大

卡开发发

535219002 发表于 2021-7-26 09:00
老师，我用的otfg的截断能和超软的数值大小一样大

otfg按道理会略高一点。难收敛的情况建议对迭代方面的参数做调整。

535219002

卡开发发 发表于 2021-7-26 09:12
otfg按道理会略高一点。难收敛的情况建议对迭代方面的参数做调整。

我已经改了density mixing里的charge和spin，把他们都改小了。

卡开发发

535219002 发表于 2021-7-26 15:03
我已经改了density mixing里的charge和spin，把他们都改小了。

先看看你做的是什么体系比较好。

535219002

卡开发发 发表于 2021-7-26 19:41
先看看你做的是什么体系比较好。

我做的体系是含有过渡金属氧化物的，用了+u进行计算。

卡开发发

535219002 发表于 2021-7-26 20:32
我做的体系是含有过渡金属氧化物的，用了+u进行计算。

那么体系是否有磁性？如果有是否经过初始猜测？

535219002

卡开发发 发表于 2021-7-26 21:59
那么体系是否有磁性？如果有是否经过初始猜测？

卡老师，我是直接用的文献里的u值，先没有加自旋极化，后面才加的自旋极化，加了自旋极化后不好收敛了。
初始猜测应该怎么做呢

卡开发发

535219002 发表于 2021-7-27 09:35
卡老师，我是直接用的文献里的u值，先没有加自旋极化，后面才加的自旋极化，加了自旋极化后不好收敛了。
...

Modify->Electronic Configuration里面有原子的Formal Spin可以修改，然后计算选项卡勾选Use Formal Spin。

535219002

卡开发发 发表于 2021-7-27 09:44
Modify->Electronic Configuration里面有原子的Formal Spin可以修改，然后计算选项卡勾选Use Formal Spin ...

我是这样设置的参数，您看看有什么问题吗

卡开发发

535219002 发表于 2021-7-27 10:04
我是这样设置的参数，您看看有什么问题吗

没啥问题。

535219002

卡开发发 发表于 2021-7-27 10:46
没啥问题。

那问题是出在哪里呢

卡开发发

535219002 发表于 2021-7-27 10:47
那问题是出在哪里呢

如果其他精度方面的设置以及体系相关的参数没问题，那么只要调整迭代方面参数就行了。例如Denstiy Mixing和DIIS相关的参数可以再调整下，以及迭代的算法也可以做调整，比如Pulay可以换Broyden等。

另外，可能你的体系有磁性，建议勾选metal或者不要使用fix occ并给一定的展宽这样会比较好。

535219002

卡开发发 发表于 2021-7-27 13:43
如果其他精度方面的设置以及体系相关的参数没问题，那么只要调整迭代方面参数就行了。例如Denstiy Mixing ...

卡老师您看参数设置有问题嘛