请问你是怎么解决结构跑散掉的问题？我也有这样的困扰！诚心请教一下抱拳！

wzkchem5 发表于 2023-2-10 14:50
最简单的方法是不要用VASP，改用算孤立体系的程序。
一定要用VASP的话，可能只能比较磺酸根和碘离子在一 ...

OK收到，感谢老师的仔细回复！

YZH_jesse 发表于 2023-2-10 03:06
非常感谢老师回答！关于加抗衡离子计算binding energy我还有一个小问题：

比如当前体系是：1个磺酸根 ...

最简单的方法是不要用VASP，改用算孤立体系的程序。
一定要用VASP的话，可能只能比较磺酸根和碘离子在一起时的能量和离得很远的能量，而不能计算磺酸根和碘离子完全分开、各自的能量。然后靠把周期性盒子取得足够大来让误差尽量小

wzkchem5 发表于 2023-2-10 05:39
1. 那样算的是小分子的自由基，不是离子。举例：你算水溶液里的乙酸根离子，但通过改变电子数使得体系是 ...

非常感谢老师回答！关于加抗衡离子计算binding energy我还有一个小问题：

比如当前体系是：1个磺酸根小分子+1个碘离子+2个抗衡离子（比如钠离子）+水分子显式溶剂，计算磺酸根小分子和碘离子的结合能公式为：

E = E_tot - E_so3_solution - E_I_solution + E_solution

下标solution代表水分子加抗衡离子。那么在计算等式右边后三项的时候，它们各自的体系又不可避免偏离电中性（从左到右分别为+1/+1/+2），此时可能又没办法继续引入抗衡离子，这个时候要怎么处理呢？

YZH_jesse 发表于 2023-2-9 13:37
ezkchem5老师您好！我想问一下：

1. 在带电小分子溶液体系中，直接在NELECT中修改价电子数量（不添加 ...

1. 那样算的是小分子的自由基，不是离子。举例：你算水溶液里的乙酸根离子，但通过改变电子数使得体系是中性，那你算的就是乙酸自由基，跑个AIMD就会发现脱羧了，根本不能反映乙酸根离子的行为。
2. 能代表在你算的这个结合构型里排斥的能量。但是要确定你算的这个结合构型是有代表性的，随手摆一个构型然后优化是不行的

wzkchem5 发表于 2022-12-4 19:48
抗衡离子不会导致悬键。悬键是指的共价键没有饱和，因为有单电子导致收敛困难。抗衡离子和待研究体系之间 ...

ezkchem5老师您好！我想问一下：

1. 在带电小分子溶液体系中，直接在NELECT中修改价电子数量（不添加抗衡离子）得到的结果可靠吗？（我看很多二维材料相关问题的回复中指出这种做法得到的结果无意义）

2. 我计算带有磺酸根的小分子和碘离子在显式溶剂中的binding energy结果为正值，这个正值是否有意义？它能表示两者之间排斥的能量吗？

感谢老师！

wzkchem5 发表于 2022-12-4 19:48
抗衡离子不会导致悬键。悬键是指的共价键没有饱和，因为有单电子导致收敛困难。抗衡离子和待研究体系之间 ...

学习到了！感谢wzkchem5大佬的细心解答！

wang5360848 发表于 2022-12-4 11:12
我明白了，请问溶剂化模型除了更接近实际情况以外，会有助于此情况的收敛吗？因为加了抗衡离子感觉会有悬 ...

抗衡离子不会导致悬键。悬键是指的共价键没有饱和，因为有单电子导致收敛困难。抗衡离子和待研究体系之间是离子键，或者假如电离得比较开的话可能只是长程库伦相互作用而已，不会导致悬键。
如果是只有溶剂里稳定、气相里不稳定的结构，加溶剂倾向于有利于收敛，但不一定。不过即便溶剂模型让收敛变差了，也必须加，否则结果不可信。
吸附能肯定要加溶剂模型，除非你实验做的是气相里的吸附。把吸附反应方程式左右两边溶剂分子配平就行了

wzkchem5 发表于 2022-12-4 16:58
实验上的抗衡离子是什么，计算就用什么；实验已知是强电解质就把抗衡离子放在离阴离子远一点的地方，否则 ...

我明白了，请问溶剂化模型除了更接近实际情况以外，会有助于此情况的收敛吗？因为加了抗衡离子感觉会有悬键不容易收敛。
另外计算吸附能需要加了溶剂化模型吗？那是不是减去H2O分子的能量就可以呢？

wang5360848 发表于 2022-12-4 09:13
wzkchem5大佬您好！请问加抗衡离子有什么要求吗？比如我处理S2O8^(2-)离子，是在其中两个O上分别连一个H ...

实验上的抗衡离子是什么，计算就用什么；实验已知是强电解质就把抗衡离子放在离阴离子远一点的地方，否则就紧挨着阴离子。
此外注意必须做显式溶剂化，必须做构型搜索，否则结果是不能发表的

wzkchem5 发表于 2022-11-19 01:12
计算时加了抗衡阴离子了吗？如果没加的话，这个碘分子会被处理成自由基阴离子，因此自然会键长严重偏长。

wzkchem5大佬您好！请问加抗衡离子有什么要求吗？比如我处理S2O8^(2-)离子，是在其中两个O上分别连一个H，还是在盒子的任意位置散放2个H原子呢？

为什么不用Gaussian？

非常感谢老师回答！

我目前没有加抗衡阴离子，体系整体上是带一个正电荷的。我马上引入抗衡阴离子试试看

计算时加了抗衡阴离子了吗？如果没加的话，这个碘分子会被处理成自由基阴离子，因此自然会键长严重偏长。