Aparte de todo lo dicho por los demás me gustaría establecer el marco teórico para una solución generalizada (cualquier velocidad, cualquier masa, cualquier carga, cualquier distancia, siempre que las "bolas" no caigan en una singularidad).
Hay dos maneras de ver este problema.
Lo más fácil es elegir relatividad especial si las masas de las cargas son relativamente pequeñas, en cuyo caso podemos despreciar los efectos gravitatorios. En tal situación podemos utilizar
$$ \frac{d(m_{0}U^\mu)}{d\tau} = - e (F'')^\mu_\nu U^\nu $$
con $$ (F'')^\mu_\nu = \Lambda^\mu_\alpha(-v_{rcv}) \Lambda_\nu^\beta(-v_{rcv}) (F')^\alpha_\beta $$ y $$ (F')^\mu_\nu = \Lambda^\mu_\alpha(v_{src}) \Lambda_\nu^\beta(v_{src}) F^\alpha_\beta $$ (para el resto de ecuaciones, qué es cada una, cómo combinarlas, véase aquí )
Usamos esta ecuación para cada una de las bolas de espuma, y luego resolvemos (usando posiciones retardadas) "órbitas" a velocidades cercanas a la de la luz.
Variamos $q_1$ , $q_2$ dibujar gráficos, deducir lo que ocurre.
Por supuesto, tenemos que definir algunas condiciones de borde (superficie), que son muy importantes porque definen lo que sucede cuando las 2 bolas chocan.
¿Se dispersarán? ¿Se combinarán? ¿Se combinarán para formar un nuevo tipo de materia? ¿Se aniquilarán para crear muchas ondas EM u otro tipo de radiación?
Por ello es muy importante tener la definición adecuada de lo que hay bajo la superficie exterior de la pelota.
Simplemente suponiendo que es sólo una singularidad bajo la superficie exterior puede estar en desacuerdo cuando se llegó a comparar con los experimentos de la vida real.
La otra forma es utilizar relatividad general .
Aquí podríamos tomar dos caminos.
La más simple es suponer que una de las bolas tiene la carga y la masa mucho más pequeña que la otra : $m_1 >> m_2$ y $q_1 >> q_2$ .
Para tal caso @Void proporcionó aquí una respuesta en el marco de la métrica de Reissner-Nordström, pero intentaré responder desde una perspectiva un poco diferente la de la teoría de los puentes .
Einstein derivó la métrica en caso de simetría esférica para electricidad y gravedad combinadas de forma un poco diferente; eligió el signo del tensor de energía de tal forma que resolviendo las ecuaciones de campo obtenemos la métrica $g_{\mu\nu}$ :
$$ ds^2 = (1 - \frac{2m}{r} - \frac{q^2}{2r^2}) dt^2 - \frac{1}{1 - \frac{2m}{r} - \frac{q^2}{2r^2}}dr^2 - r^2(d\theta^2 + \sin^2{\theta} d\phi^2) $$
Por tanto, para dicha métrica, el horizonte de sucesos se definirá en $$\left(1 - \frac{2m}{r} - \frac{q^2}{2r^2} \right) = 0$$ Esto significa que incluso sin la ayuda de la masa podemos obtener un horizonte de sucesos.
Dado que utilizamos el cuadrado de la carga implica que no importa qué signo tiene la carga.
Para ello utilizando el análisis tradicional de agujeros negros llegaremos a la conclusión de que cualquier cosa que pase el horizonte de sucesos no tendrá forma de salir, por muy cerca de la velocidad de la luz que vayamos.
Por otro lado Einstein sugirió un cambio de variable que nos ayudaría a librarnos de la singularidad del horizonte de sucesos.
El primer paso sería elegir $u^2 = r^2 - \frac{q^2}{2}$ masa fijada $m = 0$ y luego aplicarlo a la métrica para obtener:
$$ ds^2 = - du^2 - \left(u^2 + \frac{q^2}{2} \right) \left(d\theta^2 + \sin^2{\theta} d\phi^2 \right) + \frac{2u^2}{2u^2+q^2} dt^2 $$
Así que como podemos ver si $u$ varía de $-\infty$ a $+\infty$ pero $r$ sólo tendrá valores positivos entre $\sqrt{\frac{q^2}{2}}$ y $+\infty$ . Nuestra bola más pequeña se moverá de una hoja de espaciotiempo a otra.
La última forma y la más complicada pero que dará respuestas para casos generales, no importa lo grandes/pequeñas, cuántas, lo rápidas/lentas que sean las "bolas".
Tenemos $p$ singularidades. Encerramos cada singularidad denotada por $s$ en una superficie cerrada.
$$ \int^s{(\Phi_{\mu k} + 2 \Lambda_{\mu k}) \cos{(x^k \cdot N)} dS } = 0$$
Asignamos a cada singularidad $s$ la posición $\overset{s}{\xi}$ con $\xi^k(x^0)$ siendo en realidad un 3-vector.
Una distancia de $s$ singularidad se define como: $$ \overset{s}{r}^2 = (x^1 - \overset{s}{\xi}^1)^2 + (x^2 - \overset{s}{\xi}^2)^2 + (x^3 - \overset{s}{\xi}^3)^2 $$
Las ecuaciones de campo generalizadas son: $$ \Phi_{\mu \nu} + 2 \Lambda_{\mu \nu} = C_{\mu \nu} $$
donde
$$ C_{mn} = - \sum^p_{s=1}{\left( \left(\frac{\overset{s}{C}_m}{\overset{s}{r}} \right)_{,n} + \left(\frac{\overset{s}{C}_n}{\overset{s}{r}} \right)_{,m} - \delta_{mn}\left(\frac{\overset{s}{C}_k}{\overset{s}{r}} \right)_{,k} \right)} $$
$$ C_{00} = - \sum^p_{s=1}{\left(\frac{\overset{s}{C}_k}{\overset{s}{r}} \right)_{,k}}$$
$$ C_{0n} = - \sum^p_{s=1}{\left(\frac{\overset{s}{C}_0}{\overset{s}{r}} \right)_{,n} + \left(\frac{\overset{s}{C}_n}{\overset{s}{r}} \right)_{,0} }$$
teniendo $$ \overset{s}{C}_{m} = \frac{1}{4\pi} \int^s{ 2 \Lambda_{mn} \cos{(x^n \cdot N)} dS }$$ $$ \overset{s}{C}_{0} - \frac{1}{3}\overset{s}{C}_{k}\overset{s}{\dot{\xi^k}} = \frac{1}{4\pi} \int^s{ 2 \Lambda_{on} \cos{(x^n \cdot N)} dS }$$
Ahora podemos aplicar un método de aproximación como el definido aquí y como se indica aquí podemos usar esto para partículas cargadas en campo electromagnético.
Ahora queremos encontrar la distancia para la cual las dos singularidades se vuelven inseparables, en el sentido de que ninguna de ellas puede salir de esa distancia.
Lo llamamos horizonte de sucesos (si es que existe). $ r_H(x^0_j) = \sqrt{(\overset{i}{\eta^k} - \overset{s}{\xi^k})(\overset{i}{\eta^k} - \overset{s}{\xi^k})} $ tal que $ r_H(x^0_c) \le r_H(x^0_j) $ para todos $ x^0_c > x^0_j$ (aquí $x^0$ es obviamente el componente temporal).
Por lo tanto, buscamos un campo que sea el máximo de $|C_{\mu\nu}|$ .
Simplificando los supuestos o utilizando distintos métodos de aproximación adecuados a cada etapa, podríamos llegar a una respuesta, o bien montar un superordenador y esperar a obtener un resultado.
Hay una nota más. Para generalizar realmente necesitamos considerar la influencia de todos los demás campos de toda la energía del universo ( $10^{53}$ Kg), que como resulta, se manifestará a distancias muy pequeñas en forma de efecto Casimir de alcance cósmico o energía/entrelazamiento del vacío, teniendo valores e influencias en el rango $ \Delta x \Delta p = h/2$ . Si estas influencias son pequeñas en comparación con los campos locales podemos confiar en nuestro resultado por tener observación experimental, podemos predecir lo mismo que los experimentos.