Modeling Binary Lenses and Sources with the BAGLE Python Package

Resumen Técnico: Modelado de Lentes y Fuentes Binarias con el Paquete Python BAGLE

1. Planteamiento del Problema
La microlente gravitacional es una herramienta fundamental para detectar y medir la masa de objetos compactos aislados y oscuros (como agujeros negros, estrellas de neutrones y exoplanetas flotantes) que no emiten suficiente luz para ser observados directamente. Sin embargo, una limitación significativa en el análisis de eventos de microlente es la suposición frecuente de que tanto la lente como la fuente son sistemas estelares simples.

Simulaciones recientes indican que el 55% de los eventos de microlente observados involucran sistemas binarios (ya sea la lente, la fuente o ambas). La mayoría de estos eventos pueden enmascararse como eventos de lente única y fuente única, lo que lleva a interpretaciones erróneas de los parámetros físicos si no se modelan adecuadamente. Además, cuando la escala de tiempo del evento de microlente es comparable o mayor que el período orbital del sistema binario, es necesario incorporar la dinámica orbital (Kepleriana) para modelar correctamente las curvas de luz fotométricas y las trayectorias astrométricas. Aunque existen paquetes de software existentes (como VBMicrolensing, MulensModel y pyLIMA), muchos carecen de la capacidad de realizar ajustes conjuntos (joint-fitting) de datos fotométricos y astrométricos para sistemas binarios complejos, o no incluyen todas las aproximaciones de movimiento orbital necesarias.

2. Metodología
Los autores presentan la integración de modelos de sistemas binarios en el paquete de software BAGLE (Bayesian Analysis of Gravitational Lensing Events), una herramienta basada en Python diseñada para el análisis bayesiano de eventos de microlente. La metodología se centra en desarrollar formalismos matemáticos y algoritmos numéricos para tres configuraciones principales:

1. Fuente Binaria, Lente Puntual (BSPL): Dos estrellas fuente orbitando un lente único.
2. Fuente Puntual, Lente Binario (PSBL): Una fuente única siendo lenteada por dos masas (lente binaria).
3. Fuente Binaria, Lente Binaria (BSBL): Ambos sistemas son binarios.

Para cada configuración, se implementan cuatro aproximaciones de movimiento orbital para el compañero secundario:

• Estático: Separación angular fija (válido cuando el período orbital es mucho mayor que la duración del evento).
• Lineal: Movimiento relativo con velocidad constante.
• Acelerado: Movimiento con aceleración constante (aproximación de órbitas de muy largo período).
• Kepleriano Completo: Movimiento orbital elíptico o circular definido por ocho elementos orbitales (argumento del periastrón, inclinación, longitud del nodo ascendente, excentricidad, período, tiempo de periastrón y semiejes mayores proyectados).

El marco matemático utiliza constantes de Thiele-Innes para calcular las trayectorias proyectadas de los componentes primarios y secundarios alrededor del centro de masa. Para las lentes binarias, se resuelve la ecuación de lente compleja (polinomio de quinto orden) para encontrar las posiciones de las imágenes, filtrando las "imágenes fantasma" mediante un criterio de determinante jacobiano. El paquete permite el ajuste conjunto de datos fotométricos (curvas de luz) y astrométricos (desplazamiento del centroide) en un marco de referencia bariocéntrico del sistema solar (SSB).

3. Contribuciones Clave

• Implementación de Modelos Binarios en BAGLE: Se introduce por primera vez el soporte completo para geometrías BSPL, PSBL y BSBL dentro de un marco de ajuste bayesiano que integra fotometría y astrometría.
• Flexibilidad Orbital: A diferencia de otros paquetes que a menudo se limitan a órbitas circulares o requieren parámetros específicos, BAGLE ofrece un espectro completo de modelos: desde aproximaciones lineales/acceleradas (computationalesmente eficientes para órbitas largas) hasta soluciones Keplerianas completas (necesarias para órbitas cortas o excentricidades altas).
• Ajuste Conjunto (Joint-Fitting): BAGLE es uno de los pocos paquetes (junto con RTModel, aunque este último tiene limitaciones en BSBL) capaz de ajustar simultáneamente datos fotométricos y astrométricos para eventos binarios, lo cual es crucial para romper degeneraciones en los parámetros físicos.
• Herramientas de Visualización y Validación: Se incluyen funciones para generar mapas de magnificación y mapas de desplazamiento del centroide, así como la capacidad de calcular curvas críticas y causticas.
• Comparativa de Rendimiento: Se realiza una validación exhaustiva comparando BAGLE con VBMicrolensing y MulensModel. Los resultados muestran que BAGLE es extremadamente preciso (diferencias de magnitud de $10^{-11}$ a $10^{-15}$ con MulensModel) y es el más rápido para simulaciones a gran escala de eventos BSPL, aunque actualmente es más lento en la generación de modelos PSBL debido a la resolución de polinomios (un problema que se planea solucionar con JAX).

4. Resultados

• Validación de Modelos: Las comparaciones de curvas de luz y residuos para eventos simulados con paralaje muestran una concordancia excelente con otros códigos de referencia, confirmando la precisión de las ecuaciones implementadas.
• Impacto del Movimiento Orbital: Mediante el ajuste de un conjunto de datos simulado (mock dataset) que replica un evento BSPL con movimiento Kepleriano, se demostró que los modelos estáticos o lineales producen residuos correlacionados y un chi-cuadrado reducido ($\bar{\chi}^2$) alto (3.2). En contraste, los modelos con movimiento orbital circular y elíptico redujeron el $\bar{\chi}^2$ a 1.3 y 1.0 respectivamente, indicando que la inclusión de la dinámica orbital es esencial para modelar correctamente curvas de luz complejas.
• Dependencia de Causticas: Se analizó cómo la geometría de las causticas y curvas críticas varía con la relación de masas ($q$) y la separación angular. Se confirmó que la asimetría de las causticas depende directamente de $q$, mientras que para $q=1$ la simetría se mantiene independientemente de la separación, aunque las causticas se estiran a mayores distancias.
• Eficiencia Computacional: BAGLE destaca en la generación de modelos BSPL para grandes simulaciones, superando a sus competidores en velocidad cuando se instantian múltiples conjuntos de parámetros.

5. Significancia
Este trabajo representa un avance significativo en el campo de la microlente gravitacional al proporcionar una herramienta robusta y versátil para desentrañar la complejidad de los sistemas binarios.

• Preparación para Futuras Misiones: Los nuevos modelos binarios en BAGLE están diseñados específicamente para procesar los volúmenes masivos de datos que generarán el Observatorio Vera C. Rubin y el Telescopio Nancy Grace Roman.
• Caracterización de Objetos Oscuros: La capacidad de modelar con precisión la dinámica orbital de lentes binarias es crítica para medir con exactitud la masa de agujeros negros aislados y candidatos a exoplanetas flotantes, reduciendo las incertidumbres sistemáticas que surgen al asumir lentes simples.
• Futuro del Análisis: Al permitir el ajuste conjunto de fotometría y astrometría en sistemas binarios, BAGLE abre la puerta a descubrimientos más precisos sobre la población de objetos compactos oscuros en la Vía Láctea y la dinámica de sistemas estelares binarios en regiones densas como el bulbo galáctico.

En conclusión, la integración de modelos binarios en BAGLE no solo mejora la precisión de los ajustes existentes, sino que habilita nuevos tipos de análisis astrofísicos que eran previamente imposibles o poco fiables, posicionando a la herramienta como un estándar clave para la próxima generación de estudios de microlente.