Laboratorio de órbitas · Etapa 3 — La precesión secular
La elipse que gira. La elipse de Kepler es casi perfecta a corto plazo, pero los planetas se perturban unos a otros y, con los siglos, la órbita rota y se deforma lentamente. Al dejar que los elementos evolucionen, el modelo sigue a la realidad también a largo plazo. Y aparece uno de los episodios más famosos de la física: la precesión del perihelio de Mercurio.
Fundamentos teóricos
En el problema de dos cuerpos (Sol + un planeta) la elipse es fija para siempre. Pero el Sistema Solar tiene muchos cuerpos: la atracción del resto de planetas (sobre todo Júpiter) produce perturbaciones que, promediadas en el tiempo, hacen que los elementos de la órbita cambien despacio:
- El perihelio gira poco a poco (la elipse precesa):
ϖ(T) = ϖ₀ + ϖ̇·T. - El nodo se desplaza:
Ω(T) = Ω₀ + Ω̇·T. - El semieje, la excentricidad y la inclinación derivan ligeramente:
a(T),e(T),I(T).
Estas tasas seculares están medidas y tabuladas (JPL/Standish). No salen de integrar las ecuaciones de N cuerpos, sino de ajustarlas a las observaciones; pero capturan el efecto neto de las perturbaciones.
El caso de Mercurio. Su perihelio gira unos 574 segundos de arco por siglo. La mecánica de Newton, sumando el tirón de todos los planetas, explica ~531″/siglo… pero faltaban 43″/siglo. Ese pequeño resto fue un misterio durante décadas, hasta que la relatividad general de Einstein (1915) lo explicó exactamente: cerca del Sol, el espacio-tiempo curvado añade esa precesión extra. Fue una de sus primeras grandes confirmaciones.
Cómo se implementa en SusoLabs - Cosmos
Es la misma maquinaria de la etapa 2 (ecuación de Kepler, rotaciones a la eclíptica), con un cambio: en lugar de congelar los elementos en J2000, los evaluamos en el tiempo con sus tasas:
a = a₀ + ȧ·T, e = e₀ + ė·T, I = I₀ + İ·T, Ω = Ω₀ + Ω̇·T, ϖ = ϖ₀ + ϖ̇·T
Así la elipse precesa y se deforma lentamente mientras el planeta la recorre. Es exactamente la receta de Standish para posiciones aproximadas de los planetas (precisa al nivel del minuto de arco durante siglos).
Nota importante a tener en cuenta. Las tasas seculares ya incluyen la precesión observada de Mercurio (los 574″/siglo completos, relatividad incluida, porque están ajustadas a los datos). Por eso en SusoLabs - Cosmos no «vemos» el déficit de 43″/siglo como un fallo: lo contamos como la historia física que hay detrás del número. Y, como en las otras etapas, esto sigue siendo una aproximación (no integramos las ecuaciones de N cuerpos): es exacta a siglos, no a millones de años.
Qué vemos en SusoLabs - Cosmos
- Con una duración larga (décadas o un siglo) y la etapa 3, el residuo vuelve a ser pequeño y plano, mientras que con la etapa 2 (elipse congelada) crecía. Salta entre las etapas 2 y 3 con un siglo de duración para verlo.
- Elige Mercurio y observa la órbita: la elipse de la etapa 3 gira muy lentamente respecto a la de la etapa 2; su perihelio se desplaza, igual que en la realidad.
- En planetas poco perturbados o en intervalos cortos, las etapas 2 y 3 son casi indistinguibles: la precesión es un efecto lento. Por eso hace falta mirar a largo plazo para apreciarla.