求助VASP优化FeS2 (100)面上吸附构型时，吸附的IO4离子畸变

yh545930564 · 发表于 Post on 2022-8-5 23:39:50

各位老师和同学好，小弟在优化FeS2 (100)面上的S位点吸附IO4时，收敛模型总是把IO4离子分解成O + IO3。
尝试过多种初猜构型，如提高IO4的高度，调整IO4的构型等，但得到的结果都是IO4分解，I和O距离太远，不符合实际。

在这里我把IO4离子直接当成IO4分子处理了，然后竖直放到S位点上成键。
不知如何调整才能让IO4吸附在S位点上。

以下是我的INCAR/POSCAR/CONTCAR文件，非常感谢各位的帮助！

yh545930564 · 发表于 Post on 2022-8-6 08:48:09

希望有老师或者前辈能帮忙解答一下呀，网上找了大半天，不知道如何解决，谢谢各位

乐平 · 发表于 Post on 2022-8-6 10:23:51

做吸附计算不加 IVDW = 11 或者 12？

yh545930564 · 发表于 Post on 2022-8-6 10:39:26

乐平发表于 2022-8-6 10:23
做吸附计算不加 IVDW = 11 或者 12？

谢谢您的解答，大致了解了一下色散作用，我先加上运行试试，非常感谢。
我的目的是比较S位点和Fe位点的吸附能大小，Fe位点的吸附构型看起来比较合理，但现在S位点的构型奇特。

yh545930564 · 发表于 Post on 2022-8-6 14:44:39

乐平发表于 2022-8-6 10:23
做吸附计算不加 IVDW = 11 或者 12？

您好，我加了 IVDW = 12 后优化还是会断键，请问您知道是什么原因吗？谢谢

yh545930564 · 发表于 Post on 2022-8-6 21:12:41

希望各位前辈帮帮忙呀，不胜感激

QH1995 · 发表于 Post on 2022-8-7 14:45:31

试试两个O原子分别加在两个不同的S原子上？

yh545930564 · 发表于 Post on 2022-8-7 16:34:34

QH1995 发表于 2022-8-7 14:45
试试两个O原子分别加在两个不同的S原子上？

好滴，我试试，谢谢

wzkchem5 · 发表于 Post on 2022-8-7 18:07:16

这个体系实验已知I-O键不解离吗？
可能实验上本来就应该自发解离的，因为IO4-是强氧化剂，FeS2是不太弱的还原剂，可能IO4-直接把FeS2氧化了。
此外，注意体系电荷以及溶剂化问题，如果不加抗衡离子，你算的是IO4自由基，不是IO4-离子；如果不加溶剂，你算的是气相反应，但我觉得你这个实验应该做的是水相反应吧

yh545930564 · 发表于 Post on 2022-8-7 19:36:02

wzkchem5 发表于 2022-8-7 18:07
这个体系实验已知I-O键不解离吗？
可能实验上本来就应该自发解离的，因为IO4-是强氧化剂，FeS2是不太弱的 ...

   非常感谢您的解答。
   这个体系最后IO4离子会变成O自由基(羟基自由基)和一个IO3自由基。在整个计算中，我首先想得到的是吸附态的模型，然后会构建如上所述的终态模型，再计算过渡态得到反应路径和势垒。目前就是卡到这第一步了，吸附态模型一直优化不出来。最想问前辈的问题是，我如此的计算思路是否可行？
   如果加上抗衡离子，应该是Na+，这方面小弟刚学计算，可参考资料太少，不懂如何设置；也有人直接设置NELECT参数，但又说VASP带电体系计算不靠谱，想来也是对比Fe位点和S位点吸附能大小，绝对数值不用太精确，所以就忽略这个问题了，谢谢前辈指正。
   这个实验确实是水相反应，具体参数设置我先找找资料。
   最后，再次感谢前辈百忙之中解答！

yh545930564 · 发表于 Post on 2022-8-7 19:37:31

QH1995 发表于 2022-8-7 14:45
试试两个O原子分别加在两个不同的S原子上？

您好，我试了一下，最后两个O分别和不同的S原子结合，然后IO2单飞了

wzkchem5 · 发表于 Post on 2022-8-7 22:16:02

yh545930564 发表于 2022-8-7 12:36
非常感谢您的解答。
这个体系最后IO4离子会变成O自由基(羟基自由基)和一个IO3自由基。在 ...

”IO4离子会变成O自由基(羟基自由基)和一个IO3自由基“，电荷不守恒啊，负电荷去哪了？
还是那句话，电荷必须搞对，不仅反应物，生成物以及过渡态的电荷也必须对。如果该加抗衡离子不加，就等于本来你要算IO4-，结果算了IO4自由基，这不是”绝对数值不用太精确“，是整个计算都是定性错误。不能因为电子小，就不注意把电子的数目搞对，差哪怕一个电子，结果都完全是垃圾。

yh545930564 · 发表于 Post on 2022-8-7 22:29:47

本帖最后由 yh545930564 于 2022-8-7 22:32 编辑

wzkchem5 发表于 2022-8-7 22:16
”IO4离子会变成O自由基(羟基自由基)和一个IO3自由基“，电荷不守恒啊，负电荷去哪了？
还是那句话，电 ...

      谢谢您的指点。
      终态产物我再确定一下；然后电荷的话，我找下资料加个抗衡离子再试试。
      非常感谢，

。

neocc · 发表于 Post on 2022-8-8 00:31:57

本帖最后由 neocc 于 2022-8-8 09:30 编辑

wzkchem5 发表于 2022-8-7 18:07
这个体系实验已知I-O键不解离吗？
可能实验上本来就应该自发解离的，因为IO4-是强氧化剂，FeS2是不太弱的 ...

请问老师在溶剂环境中并加入抗衡离子，是不是要利用vaspsol做隐式溶剂化

wzkchem5 · 发表于 Post on 2022-8-8 00:44:50

neocc 发表于 2022-8-7 17:31
请问老师加入溶剂环境并抗衡离子，是不是要利用vaspsol做隐式溶剂化

对，并且最好再加一些显式溶剂分子

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

[VASP] 求助VASP优化FeS2 (100)面上吸附构型时，吸附的IO4离子畸变

评分 Rate

浏览过的版块