qzupc 发表于 2024-8-19 10:29 理论上直接磁矩设置为0,然后精度fine应该没问题,如果SCF发生一些震荡可以稍微调整一下density mixing相关的参数。 |
卡开发发 发表于 2024-8-17 10:52 谢谢卡老师 最后收敛成功了 解决办法是 设置这个体系磁矩为零,没有随意初猜磁矩;另外先降低scf tolerance 从coarse先开始优化 然后提高精度到fine; 没有试过一开始设置磁矩为零然后直接精度fine计算 |
本帖最后由 卡开发发 于 2024-8-17 10:55 编辑 2楼提到的问题,你可以把磁矩初猜都设置成0看下。另外如果实在比较难收敛可以试一下老的超软赝势,OTFG有时候是很难收敛。 |
尝试把empty bands提高应该可以 |
WilliamH 发表于 2024-8-8 14:33 你好 我目前先降低精度到coarse ,几何优化成功了 ,我看输出结果确实和你说的一样 整个体系没有磁性,现在我想的是在这个优化好的结构上提高精度到medium cutoff截断能我就还是保持450ev 感谢大佬的指点 附上我输出结果 Atomic Populations (Mulliken) ----------------------------- Species Ion Spin s p d f Total Charge(e) Spin(hbar/2) ========================================================================================= O 1 up: 0.94 2.41 0.00 0.00 3.35 -0.71 -0.00 1 dn: 0.94 2.41 0.00 0.00 3.35 O 2 up: 0.93 2.41 0.00 0.00 3.34 -0.68 -0.00 2 dn: 0.93 2.41 0.00 0.00 3.34 O 3 up: 0.93 2.40 0.00 0.00 3.33 -0.66 0.00 3 dn: 0.93 2.40 0.00 0.00 3.33 O 4 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 4 dn: 0.94 2.41 0.00 0.00 3.35 O 5 up: 0.93 2.41 0.00 0.00 3.34 -0.68 0.00 5 dn: 0.93 2.41 0.00 0.00 3.34 O 6 up: 0.93 2.40 0.00 0.00 3.33 -0.66 0.00 6 dn: 0.93 2.40 0.00 0.00 3.33 O 7 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 7 dn: 0.94 2.41 0.00 0.00 3.35 O 8 up: 0.93 2.41 0.00 0.00 3.34 -0.68 -0.00 8 dn: 0.93 2.41 0.00 0.00 3.34 O 9 up: 0.93 2.40 0.00 0.00 3.33 -0.66 0.00 9 dn: 0.93 2.40 0.00 0.00 3.33 O 10 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 10 dn: 0.94 2.41 0.00 0.00 3.35 O 11 up: 0.93 2.41 0.00 0.00 3.34 -0.68 -0.00 11 dn: 0.93 2.41 0.00 0.00 3.34 O 12 up: 0.93 2.40 0.00 0.00 3.33 -0.66 0.00 12 dn: 0.93 2.40 0.00 0.00 3.33 O 13 up: 0.93 2.40 0.00 0.00 3.33 -0.66 0.00 13 dn: 0.93 2.40 0.00 0.00 3.33 O 14 up: 0.93 2.41 0.00 0.00 3.34 -0.68 -0.00 14 dn: 0.93 2.41 0.00 0.00 3.34 O 15 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 15 dn: 0.94 2.41 0.00 0.00 3.35 O 16 up: 0.93 2.40 0.00 0.00 3.33 -0.66 0.00 16 dn: 0.93 2.40 0.00 0.00 3.33 O 17 up: 0.93 2.41 0.00 0.00 3.34 -0.68 0.00 17 dn: 0.93 2.41 0.00 0.00 3.34 O 18 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 18 dn: 0.94 2.41 0.00 0.00 3.35 O 19 up: 0.93 2.40 0.00 0.00 3.33 -0.66 -0.00 19 dn: 0.93 2.40 0.00 0.00 3.33 O 20 up: 0.93 2.41 0.00 0.00 3.34 -0.68 -0.00 20 dn: 0.93 2.41 0.00 0.00 3.34 O 21 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 21 dn: 0.94 2.41 0.00 0.00 3.35 O 22 up: 0.93 2.40 0.00 0.00 3.33 -0.66 -0.00 22 dn: 0.93 2.40 0.00 0.00 3.33 O 23 up: 0.93 2.41 0.00 0.00 3.34 -0.68 0.00 23 dn: 0.93 2.41 0.00 0.00 3.34 O 24 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 24 dn: 0.94 2.41 0.00 0.00 3.35 O 25 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 25 dn: 0.94 2.41 0.00 0.00 3.35 O 26 up: 0.93 2.41 0.00 0.00 3.34 -0.68 -0.00 26 dn: 0.93 2.41 0.00 0.00 3.34 O 27 up: 0.93 2.40 0.00 0.00 3.33 -0.66 0.00 27 dn: 0.93 2.40 0.00 0.00 3.33 O 28 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 28 dn: 0.94 2.41 0.00 0.00 3.35 O 29 up: 0.93 2.41 0.00 0.00 3.34 -0.68 -0.00 29 dn: 0.93 2.41 0.00 0.00 3.34 O 30 up: 0.93 2.40 0.00 0.00 3.33 -0.66 0.00 30 dn: 0.93 2.40 0.00 0.00 3.33 O 31 up: 0.94 2.41 0.00 0.00 3.35 -0.71 -0.00 31 dn: 0.94 2.41 0.00 0.00 3.35 O 32 up: 0.93 2.41 0.00 0.00 3.34 -0.68 -0.00 32 dn: 0.93 2.41 0.00 0.00 3.34 O 33 up: 0.93 2.40 0.00 0.00 3.33 -0.66 0.00 33 dn: 0.93 2.40 0.00 0.00 3.33 O 34 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 34 dn: 0.94 2.41 0.00 0.00 3.35 O 35 up: 0.93 2.41 0.00 0.00 3.34 -0.68 -0.00 35 dn: 0.93 2.41 0.00 0.00 3.34 O 36 up: 0.93 2.40 0.00 0.00 3.33 -0.66 0.00 36 dn: 0.93 2.40 0.00 0.00 3.33 O 37 up: 0.93 2.40 0.00 0.00 3.33 -0.66 0.00 37 dn: 0.93 2.40 0.00 0.00 3.33 O 38 up: 0.93 2.41 0.00 0.00 3.34 -0.68 -0.00 38 dn: 0.93 2.41 0.00 0.00 3.34 O 39 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 39 dn: 0.94 2.41 0.00 0.00 3.35 O 40 up: 0.93 2.40 0.00 0.00 3.33 -0.66 -0.00 40 dn: 0.93 2.40 0.00 0.00 3.33 O 41 up: 0.93 2.41 0.00 0.00 3.34 -0.68 -0.00 41 dn: 0.93 2.41 0.00 0.00 3.34 O 42 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 42 dn: 0.94 2.41 0.00 0.00 3.35 O 43 up: 0.93 2.40 0.00 0.00 3.33 -0.66 0.00 43 dn: 0.93 2.40 0.00 0.00 3.33 O 44 up: 0.93 2.41 0.00 0.00 3.34 -0.68 -0.00 44 dn: 0.93 2.41 0.00 0.00 3.34 O 45 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 45 dn: 0.94 2.41 0.00 0.00 3.35 O 46 up: 0.93 2.40 0.00 0.00 3.33 -0.66 0.00 46 dn: 0.93 2.40 0.00 0.00 3.33 O 47 up: 0.93 2.41 0.00 0.00 3.34 -0.68 -0.00 47 dn: 0.93 2.41 0.00 0.00 3.34 O 48 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 48 dn: 0.94 2.41 0.00 0.00 3.35 O 49 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 49 dn: 0.94 2.41 0.00 0.00 3.35 O 50 up: 0.93 2.41 0.00 0.00 3.34 -0.68 -0.00 50 dn: 0.93 2.41 0.00 0.00 3.34 O 51 up: 0.93 2.40 0.00 0.00 3.33 -0.66 0.00 51 dn: 0.93 2.40 0.00 0.00 3.33 O 52 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 52 dn: 0.94 2.41 0.00 0.00 3.35 O 53 up: 0.93 2.41 0.00 0.00 3.34 -0.68 -0.00 53 dn: 0.93 2.41 0.00 0.00 3.34 O 54 up: 0.93 2.40 0.00 0.00 3.33 -0.66 0.00 54 dn: 0.93 2.40 0.00 0.00 3.33 O 55 up: 0.94 2.41 0.00 0.00 3.35 -0.71 -0.00 55 dn: 0.94 2.41 0.00 0.00 3.35 O 56 up: 0.93 2.41 0.00 0.00 3.34 -0.68 -0.00 56 dn: 0.93 2.41 0.00 0.00 3.34 O 57 up: 0.93 2.40 0.00 0.00 3.33 -0.66 0.00 57 dn: 0.93 2.40 0.00 0.00 3.33 O 58 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 58 dn: 0.94 2.41 0.00 0.00 3.35 O 59 up: 0.93 2.41 0.00 0.00 3.34 -0.68 -0.00 59 dn: 0.93 2.41 0.00 0.00 3.34 O 60 up: 0.93 2.40 0.00 0.00 3.33 -0.66 0.00 60 dn: 0.93 2.40 0.00 0.00 3.33 O 61 up: 0.93 2.40 0.00 0.00 3.33 -0.66 0.00 61 dn: 0.93 2.40 0.00 0.00 3.33 O 62 up: 0.93 2.41 0.00 0.00 3.34 -0.68 -0.00 62 dn: 0.93 2.41 0.00 0.00 3.34 O 63 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 63 dn: 0.94 2.41 0.00 0.00 3.35 O 64 up: 0.93 2.40 0.00 0.00 3.33 -0.66 -0.00 64 dn: 0.93 2.40 0.00 0.00 3.33 O 65 up: 0.93 2.41 0.00 0.00 3.34 -0.68 -0.00 65 dn: 0.93 2.41 0.00 0.00 3.34 O 66 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 66 dn: 0.94 2.41 0.00 0.00 3.35 O 67 up: 0.93 2.40 0.00 0.00 3.33 -0.66 0.00 67 dn: 0.93 2.40 0.00 0.00 3.33 O 68 up: 0.93 2.41 0.00 0.00 3.34 -0.68 -0.00 68 dn: 0.93 2.41 0.00 0.00 3.34 O 69 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 69 dn: 0.94 2.41 0.00 0.00 3.35 O 70 up: 0.93 2.40 0.00 0.00 3.33 -0.66 -0.00 70 dn: 0.93 2.40 0.00 0.00 3.33 O 71 up: 0.93 2.41 0.00 0.00 3.34 -0.68 -0.00 71 dn: 0.93 2.41 0.00 0.00 3.34 O 72 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 72 dn: 0.94 2.41 0.00 0.00 3.35 O 73 up: 0.94 2.41 0.00 0.00 3.35 -0.71 -0.00 73 dn: 0.94 2.41 0.00 0.00 3.35 O 74 up: 0.93 2.41 0.00 0.00 3.34 -0.68 -0.00 74 dn: 0.93 2.41 0.00 0.00 3.34 O 75 up: 0.93 2.40 0.00 0.00 3.33 -0.66 0.00 75 dn: 0.93 2.40 0.00 0.00 3.33 O 76 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 76 dn: 0.94 2.41 0.00 0.00 3.35 O 77 up: 0.93 2.41 0.00 0.00 3.34 -0.68 0.00 77 dn: 0.93 2.41 0.00 0.00 3.34 O 78 up: 0.93 2.40 0.00 0.00 3.33 -0.66 0.00 78 dn: 0.93 2.40 0.00 0.00 3.33 O 79 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 79 dn: 0.94 2.41 0.00 0.00 3.35 O 80 up: 0.93 2.41 0.00 0.00 3.34 -0.68 -0.00 80 dn: 0.93 2.41 0.00 0.00 3.34 O 81 up: 0.93 2.40 0.00 0.00 3.33 -0.66 0.00 81 dn: 0.93 2.40 0.00 0.00 3.33 O 82 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 82 dn: 0.94 2.41 0.00 0.00 3.35 O 83 up: 0.93 2.41 0.00 0.00 3.34 -0.68 -0.00 83 dn: 0.93 2.41 0.00 0.00 3.34 O 84 up: 0.93 2.40 0.00 0.00 3.33 -0.66 0.00 84 dn: 0.93 2.40 0.00 0.00 3.33 O 85 up: 0.93 2.40 0.00 0.00 3.33 -0.66 0.00 85 dn: 0.93 2.40 0.00 0.00 3.33 O 86 up: 0.93 2.41 0.00 0.00 3.34 -0.68 -0.00 86 dn: 0.93 2.41 0.00 0.00 3.34 O 87 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 87 dn: 0.94 2.41 0.00 0.00 3.35 O 88 up: 0.93 2.40 0.00 0.00 3.33 -0.66 0.00 88 dn: 0.93 2.40 0.00 0.00 3.33 O 89 up: 0.93 2.41 0.00 0.00 3.34 -0.68 0.00 89 dn: 0.93 2.41 0.00 0.00 3.34 O 90 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 90 dn: 0.94 2.41 0.00 0.00 3.35 O 91 up: 0.93 2.40 0.00 0.00 3.33 -0.66 -0.00 91 dn: 0.93 2.40 0.00 0.00 3.33 O 92 up: 0.93 2.41 0.00 0.00 3.34 -0.68 -0.00 92 dn: 0.93 2.41 0.00 0.00 3.34 O 93 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 93 dn: 0.94 2.41 0.00 0.00 3.35 O 94 up: 0.93 2.40 0.00 0.00 3.33 -0.66 -0.00 94 dn: 0.93 2.40 0.00 0.00 3.33 O 95 up: 0.93 2.41 0.00 0.00 3.34 -0.68 0.00 95 dn: 0.93 2.41 0.00 0.00 3.34 O 96 up: 0.94 2.41 0.00 0.00 3.35 -0.71 0.00 96 dn: 0.94 2.41 0.00 0.00 3.35 Ce 1 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 1 dn: 1.12 3.04 0.76 0.40 5.32 Ce 2 up: 1.11 3.05 0.76 0.39 5.32 1.36 0.00 2 dn: 1.11 3.05 0.76 0.39 5.32 Ce 3 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 3 dn: 1.12 3.04 0.76 0.40 5.32 Ce 4 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 4 dn: 1.12 3.04 0.76 0.40 5.32 Ce 5 up: 1.11 3.05 0.76 0.39 5.32 1.36 0.00 5 dn: 1.11 3.05 0.76 0.39 5.32 Ce 6 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 6 dn: 1.12 3.04 0.76 0.40 5.32 Ce 7 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 7 dn: 1.12 3.04 0.76 0.40 5.32 Ce 8 up: 1.11 3.05 0.76 0.39 5.32 1.36 0.00 8 dn: 1.11 3.05 0.76 0.39 5.32 Ce 9 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 9 dn: 1.12 3.04 0.76 0.40 5.32 Ce 10 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 10 dn: 1.12 3.04 0.76 0.40 5.32 Ce 11 up: 1.11 3.05 0.76 0.39 5.32 1.36 0.00 11 dn: 1.11 3.05 0.76 0.39 5.32 Ce 12 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 12 dn: 1.12 3.04 0.76 0.40 5.32 Ce 13 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 13 dn: 1.12 3.04 0.76 0.40 5.32 Ce 14 up: 1.11 3.05 0.76 0.39 5.32 1.36 0.00 14 dn: 1.11 3.05 0.76 0.39 5.32 Ce 15 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 15 dn: 1.12 3.04 0.76 0.40 5.32 Ce 16 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 16 dn: 1.12 3.04 0.76 0.40 5.32 Ce 17 up: 1.11 3.05 0.76 0.39 5.32 1.36 -0.00 17 dn: 1.11 3.05 0.76 0.39 5.32 Ce 18 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 18 dn: 1.12 3.04 0.76 0.40 5.32 Ce 19 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 19 dn: 1.12 3.04 0.76 0.40 5.32 Ce 20 up: 1.11 3.05 0.76 0.39 5.32 1.36 0.00 20 dn: 1.11 3.05 0.76 0.39 5.32 Ce 21 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 21 dn: 1.12 3.04 0.76 0.40 5.32 Ce 22 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 22 dn: 1.12 3.04 0.76 0.40 5.32 Ce 23 up: 1.11 3.05 0.76 0.39 5.32 1.36 0.00 23 dn: 1.11 3.05 0.76 0.39 5.32 Ce 24 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 24 dn: 1.12 3.04 0.76 0.40 5.32 Ce 25 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 25 dn: 1.12 3.04 0.76 0.40 5.32 Ce 26 up: 1.11 3.05 0.76 0.39 5.32 1.36 0.00 26 dn: 1.11 3.05 0.76 0.39 5.32 Ce 27 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 27 dn: 1.12 3.04 0.76 0.40 5.32 Ce 28 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 28 dn: 1.12 3.04 0.76 0.40 5.32 Ce 29 up: 1.11 3.05 0.76 0.39 5.32 1.36 0.00 29 dn: 1.11 3.05 0.76 0.39 5.32 Ce 30 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 30 dn: 1.12 3.04 0.76 0.40 5.32 Ce 31 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 31 dn: 1.12 3.04 0.76 0.40 5.32 Ce 32 up: 1.11 3.05 0.76 0.39 5.32 1.36 0.00 32 dn: 1.11 3.05 0.76 0.39 5.32 Ce 33 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 33 dn: 1.12 3.04 0.76 0.40 5.32 Ce 34 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 34 dn: 1.12 3.04 0.76 0.40 5.32 Ce 35 up: 1.11 3.05 0.76 0.39 5.32 1.36 0.00 35 dn: 1.11 3.05 0.76 0.39 5.32 Ce 36 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 36 dn: 1.12 3.04 0.76 0.40 5.32 Ce 37 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 37 dn: 1.12 3.04 0.76 0.40 5.32 Ce 38 up: 1.11 3.05 0.76 0.39 5.32 1.36 0.00 38 dn: 1.11 3.05 0.76 0.39 5.32 Ce 39 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 39 dn: 1.12 3.04 0.76 0.40 5.32 Ce 40 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 40 dn: 1.12 3.04 0.76 0.40 5.32 Ce 41 up: 1.11 3.05 0.76 0.39 5.32 1.36 0.00 41 dn: 1.11 3.05 0.76 0.39 5.32 Ce 42 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 42 dn: 1.12 3.04 0.76 0.40 5.32 Ce 43 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 43 dn: 1.12 3.04 0.76 0.40 5.32 Ce 44 up: 1.11 3.05 0.76 0.39 5.32 1.36 0.00 44 dn: 1.11 3.05 0.76 0.39 5.32 Ce 45 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 45 dn: 1.12 3.04 0.76 0.40 5.32 Ce 46 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 46 dn: 1.12 3.04 0.76 0.40 5.32 Ce 47 up: 1.11 3.05 0.76 0.39 5.32 1.36 0.00 47 dn: 1.11 3.05 0.76 0.39 5.32 Ce 48 up: 1.12 3.04 0.76 0.40 5.32 1.37 -0.00 48 dn: 1.12 3.04 0.76 0.40 5.32 ======================================================================================= |
就目前的信息,很明显,CeO2中Ce4+不存在单电子,整体是非磁性的(可以开启自旋极化,但也不要随意初猜磁矩);450eV的动能截断对于这个体系和赝势,我想大概率是偏低的。 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2024-11-24 04:38 , Processed in 0.180522 second(s), 26 queries , Gzip On.