Joint Estimation of Properties of the Lunar Subsurface and Galactic Foregrounds with LuSEE-Night

Resumen Técnico: Estimación Conjunta de Propiedades del Subsuelo Lunar y Foregrounds Galácticos con LuSEE-Night

1. Planteamiento del Problema

El experimento LuSEE-Night (Lunar Surface Electromagnetics Experiment), una radiotelescopio de baja frecuencia que llegará al lado oculto de la Luna en 2027, tiene como objetivo principal detectar la señal de hidrógeno neutro (21 cm) de la Edad Oscura y el Amanecer Cósmico. Sin embargo, la recuperación de esta señal intrínsecamente débil enfrenta un desafío crítico: la calibración precisa del instrumento.

• Desconocimiento del Subsuelo: Las propiedades dieléctricas del subsuelo lunar en el sitio de aterrizaje son desconocidas. Las reflexiones desde el subsuelo pueden alterar el patrón de haz principal de la antena en un 10-20%.
• Acoplamiento Antena-Suelo: El subsuelo actúa como un medio dieléctrico que puede absorber una parte significativa de la potencia de la señal del cielo y modificar la impedancia de la antena.
• Degradación de la Señal: Si no se modelan correctamente estas interacciones, el patrón de haz se complica, dificultando la separación de la señal cosmológica (21 cm) de los foregrounds galácticos (ruido de fondo de la Vía Láctea), que son más de 5 órdenes de magnitud más intensos.
• El Dilema: Existe una degeneración potencial entre los parámetros del modelo del cielo (foregrounds) y los parámetros del subsuelo lunar, lo que podría impedir la estimación precisa de ambos simultáneamente.

2. Metodología

Los autores desarrollaron un pipeline de inferencia bayesiana para estimar conjuntamente los parámetros del subsuelo lunar y los modelos de foregrounds galácticos utilizando datos simulados de LuSEE-Night.

• Simulaciones Electromagnéticas (HFSS):

  • Se utilizaron simulaciones con Ansys HFSS para modelar el patrón de haz de campo lejano de las antenas dipolo de LuSEE-Night.
  • El subsuelo se modeló como un medio estratificado de dos capas con constantes dieléctricas ($\epsilon_{r,top}$, $\epsilon_{r,bottom}$) y una profundidad de interfaz ($L$).
  • Se generaron 216 configuraciones diferentes de subsuelo (una cuadrícula $6 \times 6 \times 6$) basadas en datos inversos del sitio de aterrizaje de Apolo 17.
  • Se calcularon las fracciones de suelo (ground fraction) y la impedancia de la antena para cada configuración.

• Modelado del Foreground Galáctico:

  • Se utilizó un modelo de ley de potencia en el espacio logarítmico de frecuencia para la temperatura de brillo del cielo ($T_{sky}$), basado en el modelo ULSA (Ultralong-wavelength Sky Model).
  • Los parámetros libres incluyen el índice espectral ($\beta$), la curvatura ($\gamma$) y factores de amplitud y desplazamiento ($k_1, k_2$).

• Generación de Datos Simulados:

  • Se convolucionó el patrón de haz con los mapas del cielo para obtener la temperatura observada ($T_{obs}$).
  • Se convirtió a Densidad Espectral de Potencia (PSD) de voltaje, considerando la impedancia de la antena y el circuito receptor.
  • Se añadieron ruidos gaussianos del 5% para simular una medición realista.

• Inferencia Bayesiana y Emuladores:

  • Para evitar el costo computacional de ejecutar nuevas simulaciones HFSS en cada paso de la inferencia, se construyó un emulador (modelo sustituto) utilizando Funciones de Base Radial (RBF) para interpolar los datos simulados en el espacio de parámetros.
  • Se empleó un algoritmo de Monte Carlo Precondicionado (PMC) para maximizar la función de verosimilitud y obtener las distribuciones posteriores de los 7 parámetros (3 del subsuelo + 4 del foreground).

• Análisis de Eigenmodos:

  • Se realizó una descomposición en valores singulares (SVD) de la matriz de covarianza de los espectros simulados para identificar si las variaciones del subsuelo y del cielo eran ortogonales (independientes).

3. Contribuciones Clave

1. Marco de Inferencia Conjunta: Se demuestra por primera vez que es posible estimar simultáneamente las propiedades dieléctricas del subsuelo lunar y los parámetros espectrales de los foregrounds galácticos sin necesidad de calibración externa previa.
2. Caracterización de Firmas Espectrales: Se identifica que las variaciones inducidas por el subsuelo son localizadas alrededor de la resonancia de la antena (aprox. 30 MHz), afectando la amplitud, posición y ancho del pico. En contraste, las variaciones de los foregrounds producen estructuras espectrales suaves y de banda ancha.
3. Validación de la Ortogonalidad: El análisis de eigenmodos confirma que las direcciones de variación dominantes para el subsuelo y el cielo son casi ortogonales, lo que permite al marco bayesiano "desenredar" ambos efectos.
4. Análisis de Precisión: Se cuantifica la capacidad del sistema para recuperar parámetros con una precisión del ~5% bajo condiciones ideales.

4. Resultados

• Recuperación de Parámetros: El algoritmo PMC convergió exitosamente a distribuciones posteriores que abarcan los valores de entrada reales ("ground truth").
  • Los parámetros del subsuelo ($\epsilon_{r,top}, \epsilon_{r,bottom}, L$) se recuperaron con desviaciones estándar muy bajas (ej. $\pm 0.0016$ para $\epsilon_{r,top}$).
  • Los parámetros del foreground ($\beta, \gamma, k_1, k_2$) también se estimaron con alta precisión, aunque $k_2$ mostró una mayor incertidumbre relativa.
• Correlaciones: Se observaron correlaciones positivas moderadas entre los parámetros del subsuelo (debido a su efecto combinado en el acoplamiento antena-suelo) y entre los parámetros espectrales del foreground, pero el espacio de parámetros global permaneció bien restringido.
• Validación del Emulador:
  • El ajuste del modelo a los datos simulados fue excelente, con residuos consistentes con el nivel de ruido asumido (5%).
  • Limitación de Interpolación: Se encontró que la cuadrícula de simulación ($6 \times 6 \times 6$) es suficiente para interpolaciones cercanas a los puntos de la cuadrícula, pero los errores de interpolación aumentan significativamente cuando los parámetros reales se alejan de la cuadrícula (especialmente cerca de la resonancia de 30 MHz). Esto indica que se necesita un muestreo más denso para una inferencia robusta en un escenario real.

5. Significado y Conclusiones

Este trabajo es fundamental para el éxito de la misión LuSEE-Night y futuras misiones de radioastronomía lunar:

• Viabilidad de la Misión: Confirma que LuSEE-Night no solo puede medir la señal del universo temprano, sino que también tiene la capacidad inherente de caracterizar su propio entorno de aterrizaje (propiedades del regolito), eliminando una fuente principal de incertidumbre sistemática.
• Robustez del Modelo: Aunque el modelo de subsuelo es simplificado (dos capas), la distinción espectral entre los efectos del suelo (localizados en la resonancia) y el cielo (suaves) sugiere que el método es robusto incluso con modelos imperfectos.
• Requisitos Futuros: Se destaca la necesidad de:
  1. Muestreos más densos y sofisticados (ej. Latin Hypercube Sampling) para mejorar la precisión del emulador.
  2. Modelos de subsuelo más complejos (multicapa) que reflejen la verdadera estructura del regolito lunar.
  3. Considerar efectos de ruido del receptor y calibración de la temperatura de brillo de la Luna, que se asumieron ideales en este estudio.

En resumen, el artículo establece un marco teórico sólido para la calibración conjunta de instrumentos lunares, transformando un problema de calibración (el subsuelo desconocido) en una oportunidad de ciencia (caracterización del subsuelo y del cielo simultáneamente).