金属表面吸附并且体系加了5个水一个水合氢离子显式溶剂化，需要加NELECT参数吗

cdhhh

请教各位前辈，我现在在算一个CO2RR反应的自由能，对于固液界面处，这个体系加了5个水分子和1个水合氢离子，请问这个多出的电荷怎么处理？刚开始我结构优化和自由能计算的时候加了NELECT这个参数，就价电子数-1，但是算出来的自由能台阶图走势不对，反应方向发生改变，然后我删掉NELECT这个参数后台阶图走势基本正确但是产物选择性跟文献所报道的又相悖。请教各位老师，对于我这个体系需要在INCAR中还输入哪些参数考虑这个带电体系，如果我换成另一个体系，比如显示溶剂化然后表面又有阳离子如K+，该如何处理呢？

reid

vasp不能算带电体系。文献上所谓的吸附各种离子，其实都是原子

cdhhh

reid 发表于 2025-8-5 23:22
vasp不能算带电体系。文献上所谓的吸附各种离子，其实都是原子

感谢您的回复，那我该如何处理这个带电体系呢，是直接按中性体系算吗？不管带不带电荷

reid

cdhhh 发表于 2025-8-7 17:49
感谢您的回复，那我该如何处理这个带电体系呢，是直接按中性体系算吗？不管带不带电荷

我曾经用cp2k计算金属吸附金属离子，发现吸附能-10几eV，文献用vasp算的都是几个eV。我就改成用cp2k算吸附金属原子，发现也是几个eV了。