Since the definition of the ESP is inconsistent when finding equilibrium (single state) and MECoIn (averaged over crossing states) structures, PESs may not connect smoothly. However, as this study focuses on specific stationary points (with few exceptions), this issue is not of immediate concern. In cases such as molecular dynamics simulations, where surface hopping events are investigated, a weighted average solvation approach may prove useful.