Question 1

¿Qué es el problema de los N-cuerpos?

Accepted Answer

El problema gravitacional de los N-cuerpos plantea lo siguiente: dadas N masas puntuales con posiciones y velocidades iniciales conocidas, calcular sus posiciones futuras bajo la gravedad newtoniana mutua F = Gm₁m₂/r². Para N = 2 la solución es una sección cónica (órbita kepleriana); para N ≥ 3 no existe una solución cerrada general y se requiere integración numérica. Cada paso temporal cuesta ingenuamente O(n²) evaluaciones de fuerza, lo que hace inviables las simulaciones grandes —modelar millones de estrellas— sin un algoritmo de aproximación.

Question 2

¿Cómo funciona el quadtree de Barnes-Hut?

Accepted Answer

En cada paso temporal el algoritmo inserta todos los cuerpos en un quadtree —una partición recursiva del espacio 2D en cuatro cuadrados iguales—. Para calcular la fuerza sobre el cuerpo i, el árbol se recorre desde la raíz: si la celda de un nodo está "suficientemente lejos" (s/d < θ, donde s es el ancho de la celda y d la distancia a i), toda la celda se aproxima por su masa total y su centro de masa. En caso contrario, el nodo se abre y se prueban sus hijos. Esto reduce el número típico de interacciones de O(n) por cuerpo a O(log n).

Question 3

¿Qué controla el ángulo de apertura θ?

Accepted Answer

θ (theta) es el parámetro clave de precisión/velocidad. Con θ = 0 cada celda se abre siempre sin importar la distancia, reproduciendo la suma directa exacta O(n²). Con θ = 0,5 (un valor habitual en producción) el algoritmo es aproximadamente un 95 % tan preciso como la suma directa, mientras corre entre 30 y 100 veces más rápido para n = 1.000. Con θ = 1,5 la velocidad es máxima, pero la estructura de galaxias distantes se difumina notablemente. La simulación muestra la razón de aceleración en vivo para que puedas observar este compromiso directamente.

Question 4

¿Qué es el suavizado gravitacional (ε)?

Accepted Answer

El suavizado gravitacional reemplaza la ley de fuerza 1/r² por 1/(r² + ε²), lo que evita que la fuerza diverja a infinito cuando dos cuerpos se acercan mucho. Sin suavizado, los encuentros cercanos requerirían un paso temporal infinitesimalmente pequeño para integrarse con precisión, lo que haría que la simulación colapsara o produjera velocidades de "tirachinas" sin sentido. La longitud de suavizado ε se fija habitualmente como una fracción de la separación media entre partículas; en simulaciones cosmológicas reales ronda 1 kpc para partículas de materia oscura.

Question 5

¿Por qué el algoritmo de Barnes-Hut es O(n log n)?

Accepted Answer

El quadtree tiene O(n) nodos y profundidad O(log n). Para cada uno de los n cuerpos, el recorrido del árbol visita en promedio O(log n) nodos antes de que suficientes celdas cumplan el criterio θ. Multiplicando se obtiene O(n log n) evaluaciones de fuerza por paso. Construir el árbol también cuesta O(n log n). En la práctica los factores constantes son pequeños, y para n = 10.000 la aceleración respecto a O(n²) ronda entre 100 y 1.000 veces.

Question 6

¿Qué ocurre durante los preajustes de fusión de dos galaxias?

Accepted Answer

Cuando dos galaxias se atraviesan mutuamente, las fuerzas de marea arrancan estrellas de sus extremos y las lanzan formando largas "colas de marea" y puentes que se extienden por millones de años luz. El proceso está gobernado casi por completo por la gravedad; los discos estelares se atraviesan en gran medida sin obstáculos porque las estrellas están muy separadas entre sí. Las galaxias Vía Láctea y Andrómeda están en curso de colisión y se espera que se fusionen en unos 4.500 millones de años, produciendo una estructura de colas de marea similar a la visible en este preajuste.

Question 7

¿Cómo se compara la aceleración con el método ingenuo O(n²) en la práctica?

Accepted Answer

La simulación muestra el número de evaluaciones de Barnes-Hut junto a las n(n-1)/2 comparaciones ingenuas necesarias para la suma directa. Para n = 500 y θ = 0,7, Barnes-Hut típicamente evalúa entre 3.000 y 5.000 pares de fuerza por cuerpo por paso en lugar de 499, dando una aceleración de 50 a 100 veces. A medida que n crece la ventaja se acumula: con n = 2.000 el método ingenuo necesita casi 2 millones de evaluaciones de pares por paso, mientras que Barnes-Hut necesita solo unas 20.000.

Question 8

¿Qué son la energía cinética y potencial, y por qué fluctúan?

Accepted Answer

La energía cinética KE = Σ ½mᵢvᵢ² y la energía potencial gravitacional PE = -Σᵢ≠ⱼ Gmᵢmⱼ/rᵢⱼ (sumada sobre todos los pares). En una simulación verdaderamente conservativa, la energía total E = KE + PE debería ser constante, pero la integración numérica introduce errores. El integrador de salto de rana (Verlet) usado aquí es simpléctico, lo que significa que conserva un hamiltoniano modificado y mantiene muy pequeña la deriva de energía a largo plazo. Las fluctuaciones que ves son en parte intercambios físicos genuinos (por ejemplo, un cuerpo en una órbita excéntrica intercambia KE y PE) y en parte error de truncamiento numérico por un dt finito.

Question 9

¿Qué es el destacado de "sensibilidad a las condiciones iniciales"?

Accepted Answer

Al marcar "Resaltar el descenso θ para el cuerpo elegido", al hacer clic en un cuerpo se muestra qué celdas del quadtree aproxima frente a cuáles recorre. La misma celda del quadtree puede tratarse de forma distinta para cuerpos en posiciones diferentes, lo que ilustra que Barnes-Hut es un algoritmo adaptativo. Para observar el caos en las trayectorias reales de N-cuerpos, usa la comparación de dos trayectorias en las simulaciones relacionadas: incluso una perturbación de Δ = 10⁻⁵ en las condiciones iniciales provoca trayectorias que divergen exponencialmente, demostrando que el problema de los N-cuerpos es en general caótico.

Question 10

¿Se usa el algoritmo de Barnes-Hut en códigos reales de astrofísica?

Accepted Answer

Sí — la familia de códigos GADGET (Springel 2001, 2005) usa un árbol jerárquico muy similar a Barnes-Hut y se ha empleado para simular la Millennium Simulation (10.000 millones de partículas) y las corridas de formación de galaxias de IllustrisTNG. Los códigos modernos suelen combinar Barnes-Hut con métodos multipolares rápidos (FMM), técnicas de partícula-malla (PM) o paralelismo en GPU para simular hasta 10¹² partículas. El algoritmo también se usa fuera de la astrofísica, en dinámica molecular, simulación de fluidos y renderizado de sistemas de partículas densos.

Question 11

¿Qué significa "centro de masa" en el contexto del quadtree?

Accepted Answer

Cada nodo interno del quadtree almacena la masa total M = Σmᵢ y la posición promedio ponderada por masa r_cm = Σmᵢrᵢ / M de todos los cuerpos dentro de esa celda. Cuando se aproxima una celda distante, la fuerza gravitacional sobre el cuerpo i se calcula como si toda la masa de esa celda estuviera concentrada en r_cm. Esta aproximación de un solo polo es precisa a primer orden en s/d; se pueden incluir expansiones multipolares de orden superior (cuadrupolo, octupolo) para mayor precisión a costa de más almacenamiento y cómputo.

🌌 N-cuerpos de Barnes–Hut