优化收敛问题

lujinhui396 · 发表于 Post on 2014-11-24 21:16:18

卡开发发发表于 2014-11-24 17:28
这里我没明白你的意思，你计算的结果体相到底铁磁性稳定还是反铁磁稳定？表面呢？

加U计算时，体相是铁磁稳定。所以表面也应该是稳定的

卡开发发 · 发表于 Post on 2014-11-24 21:28:29

lujinhui396 发表于 2014-11-24 21:16
加U计算时，体相是铁磁稳定。所以表面也应该是稳定的

恩，那就意味着还是磁性的问题，我也觉得磁序应该不会变化这么剧烈。结构优化是否收敛？问题是否解决？

lujinhui396 · 发表于 Post on 2014-11-25 10:49:10

卡开发发发表于 2014-11-24 21:28
恩，那就意味着还是磁性的问题，我也觉得磁序应该不会变化这么剧烈。结构优化是否收敛？问题是否解决？

结构优化不收敛。
我优化的步骤是：（1）收敛判据：EDIFFG=-0.2.EDIFF=0.001
                        (2)得到CONTCAR后，再变成POSCAR,收敛判据：EDIFFG=-0.05.EDIFF=0.0001
                        (3)再得到CONTCAR后，再变成POSCAR，收敛判据：EDIFFG=-0.02.EDIFF=0.0001
遇到的问题：在进行第二步后还能正常的收敛。
         在进行第三步时，出现下面的错误：
ZBRENT: fatal internal in brackting
      system-shutdown; contact gK immediately
此外，优化的结果和前面第二步优化的结构相差特别大。优化步数很长。
我采用的一直是IBRION=2.

卡开发发 · 发表于 Post on 2014-11-25 14:01:21

lujinhui396 发表于 2014-11-25 10:49
结构优化不收敛。
我优化的步骤是：（1）收敛判据：EDIFFG=-0.2.EDIFF=0.001
...

如果ibrion=2的hua体相有变化就有点问题了，如果体系不是很大的话把poscar传来我看看。

lujinhui396 · 发表于 Post on 2014-11-25 14:33:49

这个是第一步优化的结果

lujinhui396 · 发表于 Post on 2014-11-25 14:34:38

这个是第二步优化的结果

lujinhui396 · 发表于 Post on 2014-11-25 14:35:07

这个是第三步优化的结果

卡开发发 · 发表于 Post on 2014-11-25 23:48:03

lujinhui396 发表于 2014-11-25 14:33
这个是第一步优化的结果

看样子应该是初始结构不好，明天后天我去学校的时候帮你跑一下看看了。

lujinhui396 · 发表于 Post on 2014-11-26 15:12:31

卡开发发发表于 2014-11-25 23:48
看样子应该是初始结构不好，明天后天我去学校的时候帮你跑一下看看了。

好的，谢谢

卡开发发 · 发表于 Post on 2014-11-30 13:24:01

lujinhui396 发表于 2014-11-26 15:12
好的，谢谢

问一下您做的是Co3O4哪个表面？是完美的表面吗，还是有缺陷的？

lujinhui396 · 发表于 Post on 2014-11-30 19:55:06

卡开发发发表于 2014-11-30 13:24
问一下您做的是Co3O4哪个表面？是完美的表面吗，还是有缺陷的？

10 0面，因为它有极性，我就人为的去掉了个二价钴，满足极性补偿。

卡开发发 · 发表于 Post on 2014-11-30 20:06:13

lujinhui396 发表于 2014-11-30 19:55
10 0面，因为它有极性，我就人为的去掉了个二价钴，满足极性补偿。

我看了一下结果计量比是对的，也就是没有缺原子，但似乎是从上表面移动了一个Co到了下表面。

卡开发发 · 发表于 Post on 2014-11-30 20:36:35

lujinhui396 发表于 2014-11-30 19:55
10 0面，因为它有极性，我就人为的去掉了个二价钴，满足极性补偿。

您这个表面上还负载了H2，这样即便上下表面极性抵消好像意义也不是很大啊？如果就按照无缺陷的情况不知是否会有所改善？

lujinhui396 · 发表于 Post on 2014-11-30 20:44:50

卡开发发发表于 2014-11-30 20:36
您这个表面上还负载了H2，这样即便上下表面极性抵消好像意义也不是很大啊？如果就按照无缺陷的情况不知是 ...

之所以采用这个去掉二价钴的模型，是因为，有关于四氧化三铁的100 表面通过低能电子衍射，发现表面是有铁缺陷。
另外，对于没有去掉二价钴的情形，我也考虑了。就是对于完美的表面，我采用默认的收敛判据和采用EDIFFG=-0.02，优化得到的结果不一样，所以，我才考虑是不是因为出现的这个ZBRENT警告的原因。

lujinhui396 · 发表于 Post on 2014-11-30 20:50:49

卡开发发发表于 2014-11-30 20:06
我看了一下结果计量比是对的，也就是没有缺原子，但似乎是从上表面移动了一个Co到了下表面。

也可以这么理解，我实际的操作是，上下表面都去掉了一个二价钴。Fe3O4和MgAl2O4（与四氧化三钴结构相同，尖晶石，四氧化三铁为特殊的反尖晶石）在实验上都观察到表面有Fe和Mg缺陷，但是我没有查到关于四氧化三钴这个表面是否有钴缺陷的实验文献。另外，模拟四氧化三钴的文献，对100面的很少，有两个课题组，一个采用了去掉二价钴的模型，一个采用了完美的表面。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[VASP] 优化收敛问题

浏览过的版块