On the Reparameterization Between Cartesian Position-Velocity Vectors and Orbital Elements in the Kepler Problem

Resumen Técnico: Sobre la Reparametrización entre Vectores de Posición-Velocidad Cartesianos y Elementos Orbitales en el Problema de Kepler

Enunciado del Problema
En el problema de Kepler, las órbitas ligadas se caracterizan tradicionalmente mediante elementos orbitales ($\theta$), los cuales proporcionan una descripción compacta del tamaño, la forma y la orientación. Sin embargo, los observables para binarias estelares y sistemas planetarios están a menudo más directamente relacionados con vectores de posición y velocidad cartesianos dependientes del tiempo, $(x, v)$. Si bien reparametrizar los problemas de inferencia orbital de elementos orbitales a vectores de estado cartesianos es computacionalmente ventajoso —particularmente cuando las órbitas están débilmente restringidas y la geometría posterior es no gaussiana o altamente correlacionada—, esta transformación requiere un manejo riguroso de las densidades de probabilidad a priori en marcos bayesianos. Específicamente, transformar los priores entre estos espacios exige el determinante de la matriz jacobiana $6 \times 6$, $J_{Kep} = \partial(x, v)/\partial\theta$. La literatura previa ha tratado a menudo esta transformación de forma implícita, ha utilizado aproximaciones numéricas o ha fallado en contabilizar las diferencias de prior inducidas, lo que ha llevado a distribuciones posteriores incomparables. Además, las formulaciones existentes en campos específicos, como el microlensing binario, han contenido errores conceptuales respecto a la definición de los ángulos orbitales que hacen que los jacobianos intermedios sean singulares.

Metodología
Los autores derivan expresiones analíticas, cerradas y compactas para los determinantes jacobianos de la transformación entre vectores de estado cartesianos y elementos orbitales, así como para varias variantes útiles.

1. Derivación Analítica: Los autores utilizan la naturaleza canónica de las variables de Delaunay $(l, g, h, L, G, H)$ para derivar el determinante jacobiano para el problema de Kepler completo en 3D. Al relacionar los elementos orbitales estándar $(a, e, i, \Omega, \omega, M)$ con las variables de Delaunay, obtienen una expresión simple para $\det J_{Kep}$. También derivan expresiones para transformaciones que involucran parámetros alternativos como el periodo orbital $P$ o el tiempo de paso por el periastron $\tau$, y para jacobianos reducidos donde se elimina el grado de libertad nodal (por ejemplo, fijando la coordenada $y$).
2. Corrección de la Formulación de Microlensing: El artículo revisa la formulación ampliamente utilizada para el movimiento orbital de la lente en eventos de microlensing binario por J. Skowron et al. (2011). Los autores identifican que Skowron et al. definieron inadvertidamente la longitud del nodo ascendente ($\Omega_{node}$) relativa al eje binario proyectado en el cielo en una época de referencia. Dado que el eje binario está determinado por los propios elementos orbitales, esta definición hace que el jacobiano sea formalmente singular (el determinante es idénticamente cero). Los autores proponen una derivación corregida donde la longitud del nodo ascendente se define relativa a un eje independiente de la órbita binaria (por ejemplo, la trayectoria de la fuente o el norte celeste).
3. Verificación Numérica y Muestreo: Las expresiones analíticas se verifican numéricamente utilizando el paquete jnkepler y la diferenciación automática de JAX. Para demostrar la utilidad práctica de la reparametrización, los autores realizan una comparación cuantitativa del ajuste de órbitas astrométricas utilizando un algoritmo de Cadena de Markov Monte Carlo (MCMC) basado en gradientes (Muestreador Sin Retorno, NUTS). Simulan datos astrométricos para un compañero (similar a HR 8799d) y comparan dos estrategias de muestreo: (i) muestreo directo de elementos orbitales con priores uniformes, y (ii) muestreo de vectores de estado cartesianos $(x, v)$ con un prior uniforme, aplicando el determinante jacobiano analítico para corregir la densidad del prior para los elementos orbitales. El muestreo de $(x, v)$ se realiza mediante un esquema específico sin rechazo que asegura que las propuestas permanezcan dentro de la región orbital ligada.

Contribuciones Clave

• Expresiones Jacobianas de Forma Cerrada: El artículo proporciona fórmulas analíticas explícitas y compactas para el determinante jacobiano de la transformación entre $(x, v)$ y los elementos orbitales estándar (Ecuación 5), así como variantes que utilizan el periodo $P$ o el tiempo de periastron $\tau$ (Ecuaciones 16–17), y casos reducidos de 5D y 2D (Ecuaciones 19, 23). Estas expresiones facilitan una conversión de prior transparente y eficiente en la inferencia bayesiana sin depender de la diferenciación numérica.
• Corrección de la Parametrización de Microlensing: Los autores aclaran un error conceptual en la formulación de Skowron et al. (2011) respecto a la definición de la longitud del nodo ascendente. Demuestran que, si bien la parametrización original hace que el jacobiano intermedio sea singular, la formulación corregida (utilizando un eje de referencia independiente) produce un jacobiano bien definido que numéricamente coincide con el resultado original, ya que el determinante es independiente de la longitud del nodo ascendente.
• Análisis Cuantitativo del Rendimiento: A través de una simulación controlada, los autores muestran que la reparametrización a vectores de estado cartesianos mejora significativamente la eficiencia y la robustez del muestreo MCMC basado en gradientes. Demuestran que este enfoque reduce la fracción de transiciones divergentes y aumenta el tamaño de muestra efectivo (ESS) por minuto-núcleo en un factor de 2–7 en comparación con el muestreo en elementos orbitales, incluso cuando los priores se ajustan rigurosamente.

Resultados

• El determinante jacobiano analítico derivado para el problema de Kepler completo es $\det J_{Kep} = \frac{1}{2} \mu^{3/2} a^{1/2} e \sin i$.
• En el contexto del microlensing, corregir la definición del nodo ascendente resuelve la singularidad en el cálculo del jacobiano, confirmando que el resultado numérico original de Skowron et al. fue efectivamente correcto a pesar de la derivación defectuosa.
• En el ajuste astrométrico, la reparametrización a vectores cartesianos mitiga las correlaciones complejas presentes en el espacio de elementos orbitales. Esto conduce a una geometría posterior más suave para los métodos de Monte Carlo Hamiltoniano, resultando en una convergencia más rápida y una mayor eficiencia de muestreo. El estudio confirma que las distribuciones posteriores obtenidas mediante ambas parametrizaciones son idénticas cuando se aplica la corrección jacobiana, validando el marco teórico.

Significado
El artículo establece que la reparametrización a vectores de estado cartesianos es una estrategia robusta y eficiente para la inferencia orbital, particularmente en regímenes donde las órbitas están débilmente restringidas. Al proporcionar expresiones jacobianas analíticas, los autores permiten la aplicación consistente de priores bayesianos a través de diferentes espacios de parámetros, resolviendo ambigüedades en estudios comparativos previos. El trabajo corrige un error de larga data, aunque sutil, en la literatura de microlensing y proporciona una base cuantitativa para la preferencia del muestreo de vectores de estado en algoritmos de inferencia basados en gradientes. Los autores señalan que, si bien el muestreo por rechazo puede ser más eficiente para problemas de baja dimensionalidad con posteriores simples, la reparametrización cartesiana ofrece una alternativa ampliamente aplicable para problemas de mayor dimensionalidad, como aquellos que involucran múltiples planetas o conjuntos de datos combinados, donde la geometría posterior es compleja.