Multiwfn创建cp2k输入文件时，如何定义离子？

Corrosion

目前，在模仿陈胜利老师的固液模型计算时遇到了问题，他们通过引入不同数量的 H3O+、Li+、Na+、K+和 F-，用于调节金属表面电荷密度，同时构建双电层（总体系电荷为0）。那通过Multiwfn创建cp2k输入文件时，如何定义F-离子呢？或者怎么在cp2k的输出文件中怎么修改，来定义F-？

sobereva

先在默认设置下计算（只需要确保体系净电荷和自旋多重度正确），然后考察原子电荷检验这些物质的带电状态，如果确认是对应的离子态，结果就直接可以用。如果明显不符，或者遇到SCF难收敛，尝试用&BS调节初猜波函数中的原子组态令其与期望的带电状态相符，然后再看看情况。
注意原子电荷不可能正好是整数（参考下文的相关信息），总会有电子转移，特别是纯泛函几乎都有明显离域性误差，且很多原子电荷计算方法（如Mulliken、Hirshfeld）还往往低估离子性，所以诸如F-如果发现是比如-0.65原子电荷，也不能算错。
使用Multiwfn通过LOBA方法计算氧化态
http://sobereva.com/362（http://bbs.keinsci.com/thread-5118-1-1.html）
一篇深入浅出、完整全面介绍原子电荷的综述文章已发表！
http://sobereva.com/714（http://bbs.keinsci.com/thread-46067-1-1.html）

KazusaT

恰当设定总的电荷和自旋多重度