Veamos primero por qué la forma cuadrada inversa es especial. Teorema de Betrand afirma que sólo dos tipos de potenciales centrales producirán órbitas estables. El potencial de oscilador armónico $V=\frac{1}{2}kr^2$ y el potencial $V=-\frac{k}{r}$ que producirá una ley de fuerza cuadrática inversa. Obviamente, la edad del universo es finita, por lo que el hecho de que las órbitas de los planetas hayan sobrevivido hasta ahora no tiene por qué implicar que siga siendo así en el futuro.
Otro argumento por el que este tipo de potencial es tan común es que, al hacer teoría cuántica de campos, el propagador (los detalles dependen de si la partícula es un bosón (gauge), un fermión o un escalar, por ahora me quedo con los escalares) tiene forma
$$\frac{1}{q^2+m^2}$$
Por lo tanto, si esta partícula donde el portador de la fuerza de su fuerza con el acoplamiento $g$ el potencial es básicamente la transformada de Fourier del propagador
$$V(r) =-g^2\frac{1}{(2\pi)^3}\int d^3k\frac{4\pi}{q^2+m^2}e^{i\vec{k}\cdot \vec{r}} = -g^2\frac{1}{r}e^{-mr}$$
Este es el famoso potencial de Yukawa. Para los portadores de fuerza sin masa, el término de amortiguación llega a 1 y la fuerza se vuelve de largo alcance con una ley de fuerza cuadrática inversa. Hasta pequeños detalles esto es análogo al caso del bosón gauge, por ejemplo, la falta de masa del fotón hace que la fuerza EM sea de largo alcance, mientras que la masividad de los bosones W,Z hace que las fuerzas débiles sean de corto alcance.
Las derivaciones anteriores utilizan las tres dimensiones espaciales. Las teorías con dimensiones extra han sugerido que las dimensiones extra grandes alterarán la ley del cuadrado inverso en algunas distancias no tan cortas (rango sub-mm). Existen resultados experimentales publicados, por ejemplo, del grupo Eöt-Wash ( http://www.npl.washington.edu/eotwash/experiments/shortRange/sr.html ) y están disponibles en arXiv.
Una posible prueba está aquí $$V(r)=-G\frac{m_1m_2}{r}(1+\alpha\exp(-r/\lambda))$$
El siguiente gráfico muestra los límites de exclusión para ambos parámetros $\alpha$ y $\lambda$