Problema de los n cuerpos

En física, la cuestión del problema de los n-cuerpos trata de determinar los movimientos individuales de un grupo de partículas materiales (en sus orígenes, un conjunto de objetos astronómicos) que interactúan mutuamente según las leyes de la gravitación universal de Newton.^[1]​ La resolución de este problema ha sido motivada por el deseo de predecir los movimientos del Sol, la Luna, los planetas y las estrellas visibles. En el siglo XX, el entendimiento de la dinámica de los sistemas de cúmulos globulares de estrellas se convirtió en un importante problema de n-cuerpos.^[2]​ La extensión de este problema al campo de la relatividad es considerablemente más difícil de solucionar.

El problema físico clásico puede plantearse de forma simplificada como:

El problema de los dos cuerpos ha sido completamente solucionado (como se detalla más adelante), así como el famoso problema de los tres cuerpos restringido.^[5]​

Conociendo tres posiciones de un cuerpo astronómico sobre su órbita, es posible obtener la ecuación de su movimiento. Por ejemplo, Isaac Newton (1643-1727), a partir de los datos que le facilitó el astrónomo John Flamsteed^[6]​ fue capaz de obtener mediante geometría analítica una ecuación para predecir el movimiento de un planeta y determinar sus propiedades orbitales: posición, diámetro orbital, periodo y velocidad orbital.^[7]​ Independientemente de este hecho, Newton y otros físicos pronto descubrieron en el curso de unos pocos años, que aquellas ecuaciones del movimiento no pronosticaban algunas órbitas demasiado correctamente.^[8]​ Newton comprendió que las fuerzas gravitatorias mutuas entre todos los planetas afectaban al conjunto de sus propias órbitas.

Este descubrimiento fue directamente al centro de la cuestión respecto al significado físico exacto del problema de los n-cuerpos: como Newton advirtió, no es suficiente con especificar la posición inicial y la velocidad, o tampoco tres posiciones orbitales, para determinar con certeza la órbita de un planeta: las fuerzas gravitatorias interactivas tienen que ser conocidas también.

Así llegaron el interés y las primeras reflexiones sobre el "problema de los n-cuerpos" a comienzos del siglo XVII. Estas fuerzas atractivas gravitatorias se ajustan a las leyes del movimiento de Newton y a su Ley de la Gravitación Universal, pero la complejidad de la interacción entre "n-cuerpos" hizo históricamente intratable la obtención de cualquier solución exacta. Irónicamente, esta evidencia dirigió muchos esfuerzos al hallazgo de aproximaciones incorrectas.

Inicialmente, el problema de los n-cuerpos no fue planteado correctamente porque no se incluía el efecto de las fuerzas interactivas gravitatorias. Newton no lo expresa explícitamente, pero de sus Principia se deduce que el problema de los n-cuerpos es irresoluble debido precisamente a aquellas fuerzas interactivas gravitacionales.^[9]​ En sus Principia, párrafo 21, se afirma que:^[10]​

Newton concluyó a través de su 3ª Ley que "según esta Ley, todos los cuerpos tienen que atraer cada cual a los otros." Esta última declaración, que implica la existencia de fuerzas interactivas gravitatorias, es clave.

Como se verá más adelante, el problema también se ajusta al 1º y 2º Principios (no newtonianos) de D'Alembert, y al algoritmo no lineal del problema de los n-cuerpos, el último intento de hallar una solución cerrada para el cálculo de las referidas fuerzas interactivas.

La cuestión de encontrar la solución general al problema de los n-cuerpos fue considerada muy importante y desafiante. De hecho, en el siglo XIX tardío el rey Óscar II de Suecia,^[11]​ aconsejado por Gösta Mittag-Leffler, estableció un premio para quien pudiese encontrar la solución al problema. El anuncio era bastante concreto:

En caso de que el problema no pudiera ser solucionado, cualquier otra contribución importante a la mecánica clásica sería entonces considerada para recibir un premio digno. El premio fue otorgado al matemático francés Henri Poincaré, aunque no solucionó el problema original (la primera versión de su contribución incluso contuvo un error serio). La versión finalmente impresa contenía muchas ideas importantes dirigidas al desarrollo de la teoría del caos. El problema con su planteamiento original fue finalmente solucionado por Karl F. Sundman para n = 3.

El problema de los n-cuerpos considera ${displaystyle N}$ punto de masa ${displaystyle m_{i},i=1,2,ldots ,N}$ en un sistema de referencia inercial en las tres dimensiones del espacio ${displaystyle mathbb {R} ^{3}}$ moviéndose bajo la influencia de la atracción gravitacional mutua. Cada masa ${displaystyle m_{i}}$ tiene asociado un vector de posición ${displaystyle mathbf {q} _{i}}$. La segunda ley de Newton dice que la aceleración que experimenta cada masa ${displaystyle m_{i}d^{2}mathbf {q} _{i}/dt^{2}}$ es proporcional a la suma de las fuerzas que actúan sobre la masa. La ley de gravitación universal de Newton establece que la fuerza gravitatoria que experimenta una masa ${displaystyle m_{i}}$ por el efecto de otra masa sola ${displaystyle m_{j}}$ viene dada por^[12]​

donde ${displaystyle G}$ es la constante de la gravitación universal y ${displaystyle left|mathbf {q} _{j}-mathbf {q} _{i} ight|}$ es la magnitud de la distancia entre ${displaystyle mathbf {q} _{i}}$ y ${displaystyle mathbf {q} _{j}}$ (según una métrica inducida por una norma ${displaystyle ell _{2}}$).

El sumatorio de todas las masas produce las n-ecuaciones del movimiento de cada cuerpo:

donde ${displaystyle U}$ es la energía potencial de cada uno

Definiendo el momento como ${displaystyle mathbf {p} _{i}=m_{i}dmathbf {q} _{i}/dt}$, las ecuaciones de Hamilton del movimiento para el problema de los n-cuerpos se transforman en^[13]​

Donde aparece el operador Hamiltoniano

y la energía cinética T

Las ecuaciones de Hamilton demuestran que el problema de los n-cuerpos es un sistema de ${displaystyle 6N}$ ecuaciones diferenciales de primer orden, con 6N condiciones iniciales: 3N coordenadas de las posiciones iniciales y ${displaystyle 3N}$ valores del impulso inicial.

La simetría entre las fuerzas implicadas entre cada dos masas permiten simplificar el problema,^[14]​ refiriendo los resultados del problema al centro de masas del sistema

que se desplaza con velocidad constante, por lo que ${displaystyle mathbf {C} =mathbf {L} _{0}t+mathbf {C} _{0}}$, donde ${displaystyle mathbf {L} _{0}}$ es la cantidad de movimiento y ${displaystyle mathbf {C} _{0}}$ es la posición inicial. Las constantes del movimiento ${displaystyle mathbf {L} _{0}}$ y ${displaystyle mathbf {C} _{0}}$ representan seis integrales del movimiento. El resultado de la simetría rotacional sobre el momento angular total es constante

donde ${displaystyle imes }$ es el producto cruzado. Los tres componentes del momento angular total ${displaystyle mathbf {A} }$ produce tres constantes más del movimiento. La anterior constante general del movimiento viene dada por la ley de conservación de la energía ${displaystyle H}$. Por lo tanto, cada problema de los n-cuerpos conlleva diez integrales de movimiento.

Dado que ${displaystyle T}$ y ${displaystyle U}$ son funciones homogéneas de grado 2 y −1, respectivamente, las ecuaciones de movimiento tienen un invariante escalar: si ${displaystyle mathbf {q} _{i}(t)}$ es una solución, también lo es ${displaystyle lambda ^{-2/3}mathbf {q} _{i}(lambda t)}$ para cualquier ${displaystyle lambda >0}$. ^[15]​

El momento de inercia de un sistema de n-cuerpos viene dado por

y el virial se da por ${displaystyle Q=(1/2)dI/dt}$. Entonces la fórmula de Lagrange-Jacobi establece que^[16]​

Para sistemas en equilibrio dinámico, el promedio de tiempo a largo plazo de ${displaystyle langle d^{2}I/dt^{2} angle }$ es cero. En promedio, la energía cinética total es la mitad de la energía potencial total, ${displaystyle langle T angle =langle U angle /2}$, que es un ejemplo del teorema de virial para sistemas gravitatorios.^[17]​ Si M es la masa total y R el tamaño característico del sistema (por ejemplo, el radio que contiene la mitad de la masa del sistema), entonces el tiempo crítico para que un sistema adquiera un equilibrio dinámico es ${displaystyle t_{ m {cr}}={sqrt {GM/R^{3}}}}$. ^[18]​

Cualquier análisis sobre fuerzas planetarias interactivas ha comenzado siempre históricamente con el problema de los dos cuerpos. El propósito de esta sección es mostrar la verdadera complejidad del cálculo de las fuerzas planetarias. Nótese que en esta sección también se hace referencia a otros asuntos, como la gravedad, el baricentro o las leyes de Kepler; al igual que en el epígrafe siguiente (el problema de los tres cuerpos). Estos conceptos cuentan con sus propias páginas. Sin embargo, aquí se citan exclusivamente desde la perspectiva del problema de los n-cuerpos.

El problema de los dos cuerpos ${displaystyle (N=2)}$ fue solucionado totalmente por Johann Bernoulli (1667-1748) (y no por Newton) mediante la utilización de la teoría clásica, asumiendo que una masa principal permanece fija, como se demuestra a continuación.^[19]​ Si se considera entonces el movimiento de dos cuerpos, como la pareja sol-tierra, con el sol fijo, entonces:

La ecuación que describe el movimiento de la masa ${displaystyle m_{2}}$ en relación con la masa ${displaystyle m_{1}}$ se obtiene fácilmente de las diferencias entre estas dos ecuaciones, y después de cancelar los términos comunes, se obtiene: ${displaystyle alpha +(eta /r^{3})mathbf {r} =0}$, donde

La ecuación ${displaystyle alpha +(eta /r^{3})mathbf {r} =0}$ es la ecuación diferencial fundamental del problema de los dos cuerpos que Bernoulli resolvió en 1734. Advirtió que para utilizar este enfoque, las fuerzas tienen que determinarse primero, para a continuación resolver la ecuación del movimiento. Esta ecuación diferencial tiene soluciones elípticas, parabólicas o hiperbólicas^[20]​^,.^[21]​

Es incorrecto pensar en ${displaystyle m_{1}}$ (el sol) como fijo en el espacio cuando se aplica la ley de la gravitación Universal de Newton y hacerlo así conduce a resultados erróneos. El punto fijo para dos cuerpos aislados interactuando gravitatoriamente es su baricentro mutuo. El problema de los dos cuerpos puede ser resuelto exactamente utilizando el sistema de coordenadas de Jacobi respecto al baricentro.

En la práctica, se puede calcular de forma simplificada la posición aproximada del baricentro del Sistema Solar mediante la combinación de solo las masas de Júpiter y del Sol:

El Sol se tambalea mientras gira alrededor del centro de la galaxia, arrastrando al Sistema Solar y a la Tierra junto con él. Kepler llegó a sus famosas tres ecuaciones como el mejor ajuste matemático de los movimientos aparentes de los planetas utilizando los datos de Tycho Brahe, y no del ajuste de las curvas al verdadero movimiento de los planetas alrededor del sol (ver figura). Robert Hooke y Newton eran bien conscientes de que las fuerzas asociadas a la ley de la gravitación universal de Newton no estaban en principio asociadas con órbitas elípticas.^[10]​ De hecho, la ley universal de Newton no es capaz de explicar el comportamiento de la órbita de Mercurio, el comportamiento gravitacional del cinturón de asteroides, o el de los anillos de Saturno.^[23]​ Newton afirmó (en la sección 11 de los Principia) que la razón principal, sin embargo, para no predecir las órbitas elípticas fue que su modelo matemático se limitó a una situación que apenas existía en el mundo real, es decir, a los movimientos de los cuerpos atraídos hacia un centro inmóvil. Algunos libros de texto de astronomía y física actuales no enfatizan la importancia negativa de la asunción de Newton, y que al final su modelo matemático es en efecto la realidad de la enseñanza. Debe entenderse que la solución del problema clásico de dos cuerpos anteriormente expuesta es una idealización matemática. Véase también la primera ley de Kepler.

Algunos escritores modernos han criticado el Sol fijo de Newton como el emblema de una escuela de pensamiento reduccionista (véase a continuación Ensayos sobre la historia de la Mecánica de Truesdell). Una consideración al margen: la física "Newtoniana" no incluye (entre otras cosas) el movimiento relativo y esta puede ser la raíz de la razón por la que Newton consideró "fijo" el Sol.^[24]​^[25]​

Esta sección se refiere a la importancia histórica del problema particularizado para tres cuerpos y de las simplificaciones introducidas para la posterior resolución del problema de los n-cuerpos.

No se sabe prácticamente nada acerca de posibles intentos tempranos de resolver el problema de los n-cuerpos para n igual o mayor que tres.^[26]​ Sin embargo, a partir del siglo XIX, el caso para n = 3 fue el más estudiado. Muchas tentativas anteriores de entender el problema de los tres cuerpos fueron cuantitativas, con el objetivo de encontrar soluciones explícitas para situaciones especiales.

La solución de Moulton es más fácil de visualizar (y definitivamente más fácil de resolver) si se considera el cuerpo más masivo (por ejemplo, el Sol) como "estacionario" en el espacio, y el cuerpo menos masivo (por ejemplo, Júpiter) en órbita alrededor de él, con los puntos de equilibrio (puntos de Lagrange) manteniendo a 60 grados por delante y por detrás. El cuerpo menos masivo pŕacticamente en una órbita fija (aunque en realidad, ninguno de los cuerpos son verdaderamente fijos, dependiendo sus órbitas en realidad del baricentro del sistema completo). Para relaciones entre masas suficientemente pequeñas respecto a las dos principales, estos puntos de equilibrio triangular son estables, tales que partículas (casi) sin masa orbitarán en estos puntos establemente alrededor de la masa principal (Sol). Los cinco puntos de equilibrio del problema circular se conocen como los puntos de Lagrange. (Véase la figura siguiente):

En el problema de los tres cuerpos restringido, el modelo matemático de la figura anterior (Ref. Moulton), muestra los puntos de Lagrange L₄ y L₅ donde se ubican los asteroides troyanos; el punto m₁ es ocupado por el Sol; y en el punto m₂ se sitúa Júpiter. L₂ es un punto dentro del cinturón de asteroides. En este modelo, el sistema completo Sol-Júpiter gira sobre su baricentro. La solución del problema de los tres cuerpos restringido predijo la ubicación de los asteroides troyanos antes de que fueran divisados por primera vez. Los círculos y bucles cerrados h son el eco de los flujos electromagnéticos emitidos desde el Sol y Júpiter. Se ha conjeturado que los dos puntos h₁ son sumideros de gravedad (donde las fuerzas gravitatorias son cero), motivo por el que permanecen atrapados los asteroides Troyanos. Se desconoce la masa total de este conjunto de asteroides.

El problema restringido de los tres cuerpos asume que la masa de uno de los tres cuerpos es despreciable. Para un análisis del caso en el que el cuerpo de masa insignificante es un satélite del cuerpo de menor masa, véase esfera de Hill; para sistemas binarios, consúltese el artículo dedicado al lóbulo de Roche. También se conocen simulaciones específicas para el resultado del problema de los tres cuerpos en movimiento, sin ningún signo evidente de trayectorias repetitivas.

El problema restringido (circular y elíptico) fue trabajado extensivamente por muchos famosos matemáticos y físicos, en particular por Poincaré a finales del siglo XIX. El trabajo de Poincaré sobre el problema restringido de los tres cuerpos supuso la fundación de la teoría del caos determinístico. Para el problema restringido existen cinco puntos de equilibrio. Tres son colineales con las masas (en el marco rotatorio) y son inestables. Los dos restantes se encuentran en el tercer vértice de dos triángulos equiláteros con los dos cuerpos situados en uno de sus lados.

El problema planetario es un problema de n-cuerpos en el caso de que una de las masas es mucho mayor que todas las demás. Un ejemplo prototípico de un problema planetario es el sistema Sol-Júpiter-Saturno, donde la masa del sol es aproximadamente 1000 veces mayor que las masas de Júpiter o Saturno.^[15]​ Una solución aproximada al problema debe descomponerse en ${displaystyle N-1}$ parejas de problemas de Kepler "estrella-planeta", tratando las interacciones entre los planetas como perturbaciones. La teoría de perturbaciones funciona bien cuando no hay resonancias orbitales en el sistema, es decir, cuando ninguno de los cocientes de las frecuencias estables de Kepler es un número racional. Las resonancias aparecen como pequeños denominadores en la fórmula desarrollada.

La existencia de resonancias y de denominadores pequeños condujeron a la importante cuestión de la estabilidad en el problema planetario: "Los planetas, en órbitas casi circulares alrededor de una estrella, ¿permanecen en órbitas estables o limitadas en el tiempo?"^[15]​^[29]​ En 1963, Vladimir Arnold demostró con la teoría KAM la existencia de un tipo de estabilidad del problema planetario: existe un conjunto de medida positiva de órbitas cuasiperiódicas en el caso del problema planetario restringido al plano. En la teoría KAM,^[29]​ las órbitas planetarias caóticas cuasiperiódicas quedan confinadas a regiones de volumen tórico. El resultado de Arnold fue ampliado a un teorema más general por Féjoz y Herman en 2004.^[30]​

La configuración central ${displaystyle mathbf {q} _{1}(0),ldots ,mathbf {q} _{N}(0)}$ es una configuración inicial tal que si las partículas se liberan con velocidad cero, todas se contraen hacia el centro de masas del sistema ${displaystyle mathbf {C} }$. ^[29]​ Tal movimiento se denomina una homotecia. La configuración central también pueden dar lugar a movimientos homográficos en los que todas las masas se mueven a lo largo de trayectorias Keplerianas (elípticas, circulares, parabólicas o hiperbólicas), con todas las trayectorias con la misma excentricidad ${displaystyle e}$. Para trayectorias elípticas, ${displaystyle e=1}$ corresponde a un movimiento de homotecia y ${displaystyle e=0}$ a un movimiento de equilibrio relativo, en el que la configuración sigue siendo una isometría de la configuración inicial, como si se tratase de la configuración de un cuerpo rígido. Configuraciones de ^[31]​ El caso de la configuración central también ha desempeñado un papel importante en la comprensión de la topología de las variedades invariantes creadas al fijar las primeras integrales del sistema de ecuaciones.

Aquellas situaciones en las que todas las masas se mueven cada una en su misma curva sin colisiones, se llaman coreografías. Una coreografía para ${displaystyle N=3}$ fue descubierta por Lagrange en 1772,^[32]​ en la que tres cuerpos están situados en los vértices de un triángulo equilátero que se mantiene en rotación. Otra coreografía con forma de lemniscata (en forma de "ocho") para ${displaystyle N=3}$ fue hallada numéricamente por C. Moore en 1993 y generalizada y probada por A. Chenciner y R. Montgomery en el año 2000. Desde entonces, se han encontrado muchas otras coreografías para ${displaystyle Ngeq 3}$.

Para cada solución del problema, no solo aplicando una isometría o un cambio de tiempo, sino también invirtiendo el sentido del flujo del tiempo (a diferencia de los casos de fricción, donde no se conserva la energía), se pueden deducir otras soluciones.

En la literatura física sobre el problema de los n-cuerpos (para ${displaystyle n}$ ≥ 3), a veces se hace referencia a la imposibilidad de resolver el problema de los ${displaystyle n}$ cuerpos (empleando el método anterior). Sin embargo, debe tenerse cuidado al hablar de la imposibilidad de una solución, ya que esto solo se refiere al método de integrales primeras (es una afirmación con ciertos paralelismos a los teoremas de Abel y de Galois sobre la imposibilidad de resolver ecuaciones algebraicas de quinto grado o superior por medio de fórmulas que solo utilicen raíces).

Una manera clásica de resolver el problema de los n-cuerpos es mediante las series de Taylor, que es una implementación de la solución de una ecuación diferencial mediante una serie de potencias.

Se comienza por definir el sistema de ecuaciones diferenciales:

Como x_i (t = t₀) y dx_i(t)/dt_t=t0 se dan como condiciones iniciales, cada d²x_i(t)/dt² es conocido. Se diferencia la d²x_i(t)dt² resultando d³x_i(t)/dt³ que en t₀ también es conocido. La serie de Taylor se construye iterativamente.

Con el fin de generalizar el resultado de Sundman del caso n > 3 (o n = 3 y c = 0) hay que superar dos obstáculos:

Por último, el resultado de Sundman fue generalizado para el caso de n > 3 cuerpos por Q. Wang en la década de 1990. Puesto que la estructura de las singularidades es más complicada, Wang tuvo que abandonar completamente las cuestiones relacionadas con estas singularidades. El punto central de su enfoque es transformar, de manera apropiada, las ecuaciones en un sistema nuevo, de forma que el intervalo de existencia de las soluciones de este nuevo sistema esté comprendido en el intervalo ${displaystyle [0,infty )}$.

Puede haber dos tipos de singularidades en el problema de los n-cuerpos:

Este supuesto se denomina conjetura de Painlevé (con singularidades sin colisiones). Su existencia ha sido conjeturada para n > 3 por Painlevé (véase conjetura de Painlevé). Ejemplos de este comportamiento para n = 5 se han construido por Xia^[33]​ y un modelo heurístico para n = 4 por Gerver.^[34]​ Donald G. Saari ha demostrado que para 4 o menos cuerpos, el conjunto de datos iniciales que da lugar a estas singularidades tiene medida de Lebesgue cero.^[35]​

Si bien existen soluciones analíticas para el clásico problema de dos cuerpos (es decir, no relativista) y para las configuraciones con ${displaystyle N>2}$, en general los problemas de n-cuerpos deben ser resueltos o simulados utilizando métodos numéricos.^[18]​

Para un pequeño número de cuerpos, un problema de n-cuerpos puede ser resuelto utilizando métodos directos, también llamados métodos de partícula a partícula. Estos métodos integran numéricamente las ecuaciones diferenciales del movimiento. La integración numérica de este problema puede ser un desafío por varias razones. En primer lugar, el potencial gravitacional presenta la singularidad de que crece hasta el infinito cuando la distancia entre dos partículas tiende a cero. El potencial gravitatorio puede ablandarse para eliminar esta singularidad en distancias pequeñas:^[18]​

En segundo lugar, en general para N > 2, el problema de N cuerpos es caótico,^[36]​ lo que significa que incluso pequeños errores en la integración pueden crecer exponencialmente con el tiempo. En tercer lugar, en una simulación sobre grandes extensiones de tiempo de un modelo (por ejemplo, millones de años), los errores numéricos se acumulan con la integración cuando el tiempo aumenta.

Hay una serie de técnicas para reducir los errores en la integración numérica.^[18]​ Para ello, se utilizan sistemas de coordenadas locales en el tratamiento a diferentes escalas en algunos problemas, por ejemplo, un sistema de coordenadas de la Luna respecto a la Tierra en el contexto de una simulación del sistema solar. Los métodos variacionales y la teoría de las perturbaciones pueden producir trayectorias analíticas aproximadas en las que la integración numérica puede ser corregida. El uso de un integrador simpléctico asegura que la simulación obedece a las ecuaciones de Hamilton con un alto grado de precisión y en particular que la energía se conserva.

Los métodos directos mediante integración numérica requieren cómputos del orden de ${displaystyle N^{2}/2}$ operaciones para evaluar la energía potencial de todos los pares de partículas, y por lo tanto tienen una complejidad de cálculo de orden ${displaystyle O(N^{2})}$. Para las simulaciones con muchas partículas, el factor ${displaystyle O(N^{2})}$ requiere cálculos a gran escala con tiempos de ejecución especialmente lentos.^[18]​

En este sentido, se han desarrollado una serie de métodos aproximados que reducen la complejidad de lls cálculos en comparación con lls métodos directos:^[18]​

En sistemas astrofísicos con fuertes campos gravitacionales, como los cercanos del horizonte de eventos de un agujero negro, las simulaciones de n-cuerpos deben tomar en cuenta la relatividad general; tales simulaciones son el dominio de relativitidad numérica. Sistemas numéricos que simulan las ecuaciones de campo de Einstein son muy prometedores^[18]​ y un formalismo parametrizado post-Newtoniano (PPN), como el de las ecuaciones de Einstein–Infeld–Hoffmann, se utiliza si es posible. Actualmente, el problema de los dos cuerpos en la relatividad general solo es analíticamente resoluble para el caso de Kepler, en el que una masa se supone que es mucho más grande que la otra.^[37]​

La mayoría del trabajo realizado sobre el problema de los n-cuerpos se ha centrado en el campo gravitatorio. Pero existen otros sistemas para los que la matemática de los n-cuerpos y las técnicas de simulación se han probado útiles.

En problemas de electrostática de gran escala, como en la simulación de proteínas y uniones celulares en biología estructural, el potencial eléctrico tiene la misma forma que el potencial gravitatorio, excepto en que las cargas pueden ser indistintamente positivas o negativas, produciendo tanto fuerzas repulsivas como atractivas.^[38]​ Los Fast Coulomb solvers son unas contrapartidas electrostáticas a los métodos de simuladores rápidos multipolo. A menudo son utilizados con condiciones de frontera periódica en simulación de regiones, y se utilizan técnicas de sumatorios de Ewald para acelerar los cálculos.^[39]​

En estadística y aprendizaje automático, algunos modelos tienen funciones de pérdida de una forma similar a como se comporta el potencial gravitatoria: la suma de un núcleo de funciones sobre todos los pares de objetos, donde el núcleo de la función depende de la distancia parametrizada entre los distintos objetos.^[40]​

Ejemplos de problemas que encajan en esta tipología son la búsqueda de vecinos próximos en aprendizaje de captación de datos, estimación de densidad de núcleos, y optimizaciones núcleos de cálculo. Se han desarrollado alternativas para reducir la complejidad del tiempo de cálculo de ${displaystyle O(N^{2})}$ a ${displaystyle O(N)}$, tales como los algoritmos de doble árbol, que tienen aplicación para el problema gravitatorio de los n-cuerpos.

---