Vector de Runge-Lenz

Figura 1: El vector de RL A (en rojo) para cuatro puntos (marcados como 1, 2, 3 y 4) sobre la órbita elíptica de un planeta que se mueve bajo la acción de una fuerza central que sigue la ley de la inversa del cuadrado. En centro de atracción se marca como un pequeño círculo negro a partir del cual se consideran que emanan los vectores posición (en negro). El momento angular L es perpendicular a la órbita. Los vectores coplanarios p×L y (mk/r)r se muestran en azul y en verde, respectivamente. Estas variables se definen más adelante en este artículo. El vector A es constante en dirección y magnitud.

El vector de Runge-Lenz (o vector de Laplace-Runge-Lenz) es una constante de movimiento del problema de los dos cuerpos en interacción gravitatoria mutua. La existencia de esta integral de movimiento es una de las formas más simples de probar que las trayectorias planetarias en ese caso son cónicas.

Motivación

El vector de Runge-Lenz aparece de manera natural a partir de la ecuación de movimiento. Si tomanos un sistema de referencia inercial con origen en el centro de masas y consideramos las distancias relativas de los dos cuerpos o astros respecto a él, podemos definir el vector diferencia:

$\mathbf {r} =\mathbf {r} _{1}-\mathbf {r} _{2},\qquad \qquad \mathbf {R} _{CM}={\frac {m_{1}\mathbf {r} _{1}+m_{2}\mathbf {r} _{2}}{m_{1}+m_{2}}}=0$

Donde los dos vectores que aparecen en segundo término son los vectores de posición del primer y segundo cuerpo respectivamente. En términos del vector diferencia la ecuación de movimiento puede expresarse como:

(1)${\ddot {\mathbf {r} }}=-\mu {\frac {\mathbf {r} }{r^{3}}}$

Donde μ = G(m₁+m₂). La existencia de una constante adicional además de la energía y el momento angular puede probarse muy fácilmente a partir de (1). Si se multiplica por la velocidad se obtiene:

$\mathbf {r} \times {\ddot {\mathbf {r} }}=\mathbf {0}$

Que se puede integrar sin dificultad:

$\mathbf {r} \times {\dot {\mathbf {r} }}=\mathbf {l} ={\mbox{cte.}}$

Este vector constante de hecho coincide con el momento angular por unidad de masa, y puede calcularse sin dificultad a partir de los datos iniciales. Si multiplicamos este vector constante por la aceleración nuevamente y hacemos algunas manipulaciones algebraicas tenemos:

$\mathbf {l} \times {\ddot {\mathbf {r} }}=-{\frac {\mu }{r^{3}}}\mathbf {l} \times \mathbf {r} =-{\frac {\mu }{r^{3}}}[(\mathbf {r} \times {\dot {\mathbf {r} }})\times \mathbf {r} ]=-{\frac {\mu }{r^{3}}}[r^{2}{\dot {\mathbf {r} }}-({\dot {\mathbf {r} }}\cdot \mathbf {r} )\mathbf {r} ]=-\mu \left[{\frac {\dot {\mathbf {r} }}{r}}-{\frac {\dot {r}}{r^{2}}}\mathbf {r} \right]=-\mu {\frac {d}{dt}}\left({\frac {\mathbf {r} }{r}}\right)$

De esta última igualdad reordenando términos se tiene que:

(2)${\frac {d}{dt}}\left(\mathbf {l} \times {\dot {\mathbf {r} }}+\mu {\frac {\mathbf {r} }{r}}\right)=0\qquad \Rightarrow \qquad \mathbf {A} =\left(\mathbf {l} \times {\dot {\mathbf {r} }}+\mu {\frac {\mathbf {r} }{r}}\right)={\mbox{cte.}}$

Y es este vector $\mathbf {A}$ el que se conoce como vector (Laplace-)Runge-Lenz.

Propiedades

• El vector de Runge-Lenz está contenido en el plano de la órbita.
• Además el vector de Runge-Lenz coincide con uno de los semiejes de la cónica
• El módulo del vector de Runge-Lenz de μ > 0 coincide con el valor absoluto de la excentricidad.
• Las componentes del vector de Runge-Lenz no son todas independientes ya que cumplen la relación $\mathbf {A} \cdot \mathbf {L} =0$

Forma de la órbita

De la conservación del vector Runge-Lenz, (2), se sigue que la forma de la órbita en el problema de Kepler es una cónica,​ para verlo basta multiplicar dicho vector escalarmente por el vector posición, para obtener que:

${\begin{cases}-\mathbf {r} \cdot (\mathbf {l} \times {\dot {\mathbf {r} }})=\mu r-\mathbf {A} \cdot \mathbf {r} ,\\\mathbf {l} \cdot (\mathbf {r} \times {\dot {\mathbf {r} }})=\mu r-\mathbf {A} \cdot \mathbf {r} \end{cases}}\qquad \Rightarrow \qquad \|\mathbf {l} \|^{2}=l^{2}=\mu r-\mathbf {A} \cdot \mathbf {r}$

Si la fuerza es atractiva (μ > 0) podemos reescribir la última expresión como:

(3)${\frac {l^{2}}{r}}=\mu -\mathbf {A} \cdot {\hat {\mathbf {r} }}=\mu (1+e\cos \varphi )\qquad \Rightarrow \qquad {\frac {p}{r}}=1+e\cos \varphi$

Donde e es la excentricidad de la órbita, y p = l²/μ es el semilatus rectum de la elipse, cuya ecuación en coordenadas polares viene dada por (3).

Referencias

Bibliografía

• Pérez Sánchez, Gustavo Adolfo (2015). «El Problema de Kepler. Aplicaciones del vector de Laplace-Runge-Lenz en órbitas perturbadas.». Trabajo de Fín de Grado. UNED. Consultado el 27 de junio de 2015. 

• Fernández Rañada, Antonio (2005). Fondo de Cultura Económica, ed. Dinámica Clásica (1ª edición). México, D. F. pp. 302-304. ISBN 84-206-8133-4.