计算化学公社

标题: 为什么使用castep里的otfg 赝势不能收敛而uspp赝势则能够收敛 [打印本页]

作者
Author: 535219002 时间: 2021-7-22 21:45
标题: 为什么使用castep里的otfg 赝势不能收敛而uspp赝势则能够收敛
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不是应该更精确吗

作者
Author: 卡开发发 时间: 2021-7-23 09:59
精确和容不容易SCF收敛没啥直接关系。参数合理（如动能截断OTFG高于Vanderbilt老的超软赝势）的情况下，试着调整density mixing的选项。

作者
Author: 535219002 时间: 2021-7-26 09:00

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

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

作者
Author: 卡开发发 时间: 2021-7-26 09:12

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

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

作者
Author: 535219002 时间: 2021-7-26 15:03

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

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

作者
Author: 卡开发发 时间: 2021-7-26 19:41

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

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

作者
Author: 535219002 时间: 2021-7-26 20:32

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

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

作者
Author: 卡开发发 时间: 2021-7-26 21:59

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

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

作者
Author: 535219002 时间: 2021-7-27 09:35
标题: hou
本帖最后由 535219002 于 2021-7-27 09:36 编辑

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

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

作者
Author: 卡开发发 时间: 2021-7-27 09:44

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

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

作者
Author: 535219002 时间: 2021-7-27 10:04

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

(, 下载次数 Times of downloads: 60) (, 下载次数 Times of downloads: 50) 我是这样设置的参数，您看看有什么问题吗

作者
Author: 卡开发发 时间: 2021-7-27 10:46

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

没啥问题。

作者
Author: 535219002 时间: 2021-7-27 10:47

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

那问题是出在哪里呢

作者
Author: 卡开发发 时间: 2021-7-27 13:43

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

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

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

作者
Author: 535219002 时间: 2021-7-28 17:30

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

(, 下载次数 Times of downloads: 52) [attach]37471[/attach]卡老师您看参数设置有问题嘛

作者
Author: 卡开发发 时间: 2021-7-28 21:51

535219002 发表于 2021-7-28 17:30
卡老师您看参数设置有问题嘛

应该没啥问题，混合方式和参数只能适当调整做尝试。

作者
Author: 535219002 时间: 2021-7-29 17:35
卡老师是这样的，我的体系在扩展晶胞之前是可以收敛的，但是扩展成2x2x1的超晶胞之后就不能收敛了，这是为什么呢。

作者
Author: cccwl 时间: 2023-4-14 18:54
楼主，你之前的SCF收敛问题按倒上方截图操作后解决了吗？我也遇到了类似问题，我是fecr合金

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