To evaluate the parameters of Eq. 1, knowledge of four points on the free energy surfaces is necessary. Two points are the energy of the I and J states in their respective local minima geometries. The other two points are obtained by computing the energy of the system accommodated at the equilibrium geometry of I while following the charge arrangement of J and vice versa. To compute these quantities, our calculations were performed in three steps: