La fuerza del análisis del problema se realiza gracias a @freecharly. Ahora trabajo usando otros dos métodos: trabajo virtual y energía del campo. Supongamos cargo $q_1$ es fijo y carga en $q_2$ se mueve a lo largo de la cuerda por un pequeño desplazamiento virtual $\delta r$$q_1$. El Coulomb campo debido a la carga neta $q_1/\epsilon_r$ repele $q_2$ y atrae la envolvente de carga $-q_2(1-1/\epsilon_r)$ circundante $q_2$. Al $q_2$ se mueve por la distancia a la $\delta r$, el límite de carga no se mueva con $q_2$. Por definición, la envolvente de carga no se puede mover. Lo que en realidad sucede es que la envolvente de carga en la posición original de $q_2$ despolariza a la neutralidad, mientras que algunos de los nuevos envolvente de carga de la misma cantidad vuelve a aparecer en la nueva posición de $q_2$. Por lo que no se realiza un trabajo de la envolvente de carga $-q_2(1-1/\epsilon_r)$ porque no la envolvente de carga se movió en realidad la distancia $\delta r$. Por lo tanto,
$$F\delta r=\frac{q_1q_2}{4\pi\epsilon_0\epsilon_rr^2}\delta r,$$
que significa que la primera fórmula es correcta. Observe que el visionario "desplazamiento" de la envolvente de los cargos de $-q_2(1-1/\epsilon_r)$ no es un verdadero desplazamiento. Por lo tanto no se trabaja.
La energía del campo método no distingue libre de cargos y obligado cargos o realizar un seguimiento de cómo se mueven las cargas. Utiliza la energía de los condensadores $\,W=\frac{1}{2}CU^2\,$ $\,C=\epsilon_0\epsilon_rS/d\,$ $\,E=U/d\,$ obtener
$$W=\frac{1}{2}\epsilon_0\epsilon_rE^2Sd.$$
Por lo tanto, la densidad de energía de un $E$-campo en un medio dieléctrico es mayor que el mismo $E$-campo en el vacío por un factor de $\epsilon_r$ debido a la polarización de los medios. Cuando la distancia entre las cargas netas $q_1/\epsilon_r$ $q_2/\epsilon_r$ aumenta por $\delta r$, si estos eran los cargos en el vacío, el $E$-campo de energía se reduce a una cantidad igual al trabajo realizado por la fuerza en la segunda fórmula veces $\delta r$. Pero ahora el $E$-campo de los cargos netos $q_1/\epsilon_r$ $q_2/\epsilon_r$ están en el medio. Por lo tanto, el trabajo realizado es $\epsilon_r$ veces mayor. En la fórmula tenemos
$$F\delta r=\frac{(q_1/\epsilon_r)(q_2/\epsilon_r)}{4\pi\epsilon_0r^2}\delta r\times\epsilon_r=\frac{q_1q_2}{4\pi\epsilon_0\epsilon_rr^2}\delta r.$$
$E$-campos en medios dieléctricos contienen más energía de la que el mismo $E$-campos en el vacío por $\epsilon_r$ veces.