It is known that interaction energy = $-\vec{p}.\vec{E}$ where $\vec{p}$ is dipole moment and $\vec{E}$ is the electric field.
I have to calculate the interaction energy of a system whose dipole moment and potential are the only available data.
My assumption is since the negative gradient of potential is electric field, can we calculate the dipolar component of the potential as an alternative to electric field, i.e.
$$\begin{align} \phi_{dx}&=\phi_x(x-x_c) \\ \phi_{dy}&=\phi_y(y-y_c) \\ \phi_{dz}&=\phi_z(z-z_c) \end{align}$$ where $x$, $y$, $z$ are the positions of atoms and $x_c$, $y_c$, $z_c$ is the geometrical centre, and then calculate the dot product with the dipole moment $(p_x,p_y,p_z)$.
Does it sounds meaningful?? Kindly provide some suggestions in solving my problem.