Atractor Multi-Cuerpo

Caos determinista: el futuro depende del presente con precisión infinita

Una partícula atraída por N cuerpos fijos dispuestos en círculo. Las ecuaciones son perfectamente deterministas: dado el estado inicial, el futuro está completamente determinado. Sin embargo, la trayectoria es impredecible en la práctica: una perturbación microscópica en la posición inicial produce una trayectoria completamente diferente. Este es el caos determinista.

El Sistema Físico

N atractores (centros de fuerza) están fijos en los vértices de un polígono regular. Una partícula libre se mueve bajo la acción de la suma de sus fuerzas más una fricción viscosa:

F_i = K / r_i^n · r̂_i
Fuerza de cada atractor sobre la partícula. r̂ es el vector unitario hacia el atractor.
mẍ = Σ F_i − γ ẋ
Ecuación de movimiento total. γ es el coeficiente de fricción.

La suma vectorial de N fuerzas que compiten crea un campo de potencial complejo con múltiples mínimos (cuencas de atracción). La trayectoria de la partícula depende críticamente de por cuál cuenca "cae" en cada paso.

El Exponente n: Geometría del Potencial

El exponente n controla cómo decae la fuerza con la distancia y cambia radicalmente la topología del campo:

n = 2 (gravitacional)

F ∝ 1/r². Las cuencas de atracción son suaves y las fronteras entre ellas son relativamente regulares. El atractor más cercano casi siempre "gana".

n = 3 (umbral del caos)

Las fronteras entre cuencas se vuelven fractales. Cerca de las líneas de separación, una pequeña variación en posición puede enviar la partícula a cualquier atractor.

n = 3.5 (caos profundo)

Las fronteras fractales se densifican. La predicción a largo plazo se vuelve imposible incluso con alta precisión numérica.

n = 4 (ultra local)

La fuerza decae muy rápidamente. El atractor más cercano domina casi siempre, pero las cuencas forman patrones de alta frecuencia fractal en las fronteras.

Los N Atractores

La simulación permite elegir entre 3 y 8 atractores. La simetría del polígono determina el número de cuencas y la complejidad del diagrama:

3
triángulo — cuencas simples
4
cuadrado — simetría doblada
5
pentágono — default, caos moderado
6
hexágono — patrón tipo copo
7
heptágono — sin simetría cuadrática
8
octágono — máxima complejidad

Sensibilidad a Condiciones Iniciales

El experimento de las dos partículas es la demostración más directa del caos. Actívalo con el checkbox "Dos partículas". Las dos partículas parten de posiciones casi idénticas (diferencia de 10⁻⁶ unidades). Durante unos ciclos siguen juntas; luego sus trayectorias divergen exponencialmente hasta ser completamente independientes.

El indicador de divergencia en el panel mide la distancia entre ambas partículas. En un sistema no caótico, esta distancia crecería linealmente o permanecería acotada. En el atractor multi-cuerpo, crece exponencialmente:

|δx(t)| ≈ |δx₀| · e^(λt)
λ es el exponente de Lyapunov. Si λ > 0, el sistema es caótico.

El Campo de Potencial

El fondo de la simulación muestra el campo de potencial V(x,y) como un mapa de color. Las zonas oscuras son pozos de potencial (donde están los atractores); las zonas claras son crestas de potencial (fronteras entre cuencas).

La partícula se mueve cuesta abajo en este paisaje de potencial, pero como tiene inercia (velocidad inicial), puede "saltar" de una cuenca a otra antes de asentarse. Con suficiente fricción, eventualmente cae en el atractor más cercano; con poca fricción, puede explorar todo el campo durante mucho tiempo.

Controles de la Simulación

Experimentos Guiados

Experimento 1 — El umbral del caos en n

  1. Selecciona 5 atractores, fricción media (0.003), velocidad inicial 5.
  2. Ajusta el exponente a n = 2. Observa la trayectoria durante 30 segundos. Anota en cuál atractor termina.
  3. Cambia a n = 3. Resetea. Repite varias veces con las mismas condiciones. ¿Termina siempre en el mismo atractor?
  4. Con n = 4, el comportamiento es más errático cerca de las fronteras. Observa cómo los colores del campo de potencial se vuelven más "nervioso" en las fronteras.
  5. El umbral del caos para este tipo de sistema ocurre alrededor de n = 2.5–3.

Experimento 2 — Divergencia de Lyapunov

  1. Activa "Dos partículas". Selecciona n = 3, 5 atractores, fricción baja (0.001).
  2. Observa el indicador de divergencia. Durante los primeros ciclos, las dos trayectorias son indistinguibles.
  3. En el momento en que la partícula pasa cerca de una frontera entre cuencas, la divergencia se dispara.
  4. A partir de ahí, las dos partículas van a atractores distintos: el sistema "eligió" destinos diferentes a partir de condiciones microscópicamente distintas.
  5. Resetea y repite: cada vez, el momento y el destino de la bifurcación puede ser diferente.

Experimento 3 — Fricción y predictibilidad

  1. Selecciona n = 3.5, 6 atractores.
  2. Con fricción alta (~0.01): la partícula cae rápidamente en el atractor más cercano. El sistema es predecible.
  3. Con fricción cero (~0.0): la partícula orbita caóticamente durante largo tiempo sin asentarse. La trayectoria llena el espacio.
  4. Encuentra el valor de fricción donde el comportamiento pasa de predecible a impredecible. Este es el umbral práctico de caos para este sistema.

Fractalidad de las Cuencas

La figura que forman las fronteras entre cuencas de atracción es un fractal: tiene la misma estructura a todas las escalas. Si pudieras hacer zoom infinito en la frontera entre dos cuencas, siempre encontrarías más estructura detallada, nunca una línea suave.

La frontera fractal es la firma del caos. Exactamente en la frontera entre cuencas, la más pequeña perturbación decide a cuál atractor va la partícula. Este conjunto de puntos "indecidibles" tiene dimensión fractal mayor que 1 pero menor que 2. Cuanto mayor sea esa dimensión, más caótico es el sistema.

Conexiones

🕰
Péndulo Simple
Oscilador sin caos
🌀
Chaos Lab
Lorenz, Mandelbrot y más caos
💨
Gas Ideal
Caos estadístico: muchas partículas
🪐
Órbitas de Kepler
N-cuerpos gravitacionales reales