It is worth noting that in our protocol of computing the Poisson potential of liquid phase (Fig. ), we did not sample ionic configurations of any liquid-air interface. Therefore, the Poisson potential that we obtained does not contain the surface dipole contribution and there is no double counting in Eq. 6.