PFOA-与单个 H2O 在气相和 SMD(Water) 中吸附能计算是否需要 Na⁺ 中和电荷？

noob

各位老师好，我想请教一个关于 PFOA⁻ 与单个 H2O 结合能/吸附自由能计算的问题。

我的目标是比较 PFOA⁻ 离子在气相和液相环境中与单个水分子的结合能力。计算使用 ORCA，分别在气相和 SMD(Water) 隐式水模型中对 PFOA⁻、H2O 以及 PFOA⁻···H2O 复合物进行优化和频率分析，然后用类似下面的公式计算结合自由能：

ΔG_bind = G(PFOA⁻···H2O) - G(PFOA⁻) - G(H2O)

目前遇到的主要问题是：如果直接计算 PFOA⁻···H2O，体系总电荷为 -1；如果加入 Na⁺ 来中和电荷，体系总电荷为 0，但结果差异很大。Na⁺ 会明显参与羧酸根附近的相互作用，甚至在气相中把 Na⁺ 初始放到约 8 Å 远的位置后，优化过程中仍然倾向于靠近 COO⁻，并且优化收敛也比较困难。

想请教两个问题：

1. 如果目标是研究 PFOA⁻ 离子本身与单个水分子的结合能力，在气相和 SMD(Water) 中是否都应该直接按带 -1 电荷的 PFOA⁻ 体系计算，而不是加入 Na⁺ 来中和电荷？

2. 用 SMD(Water) 隐式水模型下的 ΔG_bind 与气相下的 ΔG_bind 进行比较，来讨论 PFOA⁻ 在液相和气相中与单个水分子的结合能力差异，这样的方法是否合理？还是说由于隐式水模型本身已经描述了连续溶剂环境，再加入一个显式水分子时，得到的结合自由能需要特别解释？

感谢各位老师指教！