Se presenta un marco de simulación para probar las capacidades de imagen de matrices de radio a gran escala en la superficie lunar. Se discuten los principales componentes de ruido, y una canalización de software se camina a través de detalles sobre cómo personalizarlo para usos científicos novedosos.
En los últimos años ha habido un renovado interés en regresar a la Luna por razones de naturaleza científica y exploratoria. La Luna proporciona el campo de entrenamiento perfecto para construir bases a gran escala que uno puede aplicar a otros planetas como Marte. La existencia de una zona radioeléctuda en el lado lejano lunar tiene promesas de primeros estudios del universo y búsquedas de exoplanetas, mientras que el lado cercano proporciona una base estable que puede ser utilizada para observar emisiones de baja frecuencia de la magnetosfera de la Tierra que pueden ayudar a medir su respuesta al clima espacial entrante. La construcción de una matriz de radio a gran escala proporcionaría grandes retornos científicos, además de actuar como una prueba de la capacidad de la humanidad para construir estructuras en otros planetas. Este trabajo se centra en simular la respuesta de matrices de radio pequeñas a grandes escalas en la Luna que constan de cientos o miles de antenas. La respuesta de la matriz depende de la estructura de la emisión junto con la configuración y sensibilidad de la matriz. Se selecciona un conjunto de ubicaciones para los receptores de radio simulados, utilizando modelos de elevación digital del instrumento Altímetro láser orbitador lunar en lunar Reconnaissance Orbiter para caracterizar la elevación de las ubicaciones del receptor. Un código personalizado de aplicaciones de software de astronomía común se describe y se utiliza para procesar los datos de los receptores simulados, alineando los marcos de coordenadas lunares y del cielo usando SPICE para garantizar que las proyecciones adecuadas se utilizan para la toma de imágenes. Este marco de simulación es útil para iterar el diseño de matrices para crear imágenes de cualquier objetivo científico determinado en un pequeño campo de visión. Este marco no admite actualmente todas las imágenes del cielo.
El campo de la radioastronomía comenzó en 1932 con la detección accidental de emisiones de radio galácticas por Karl G. Jansky1 a 20 MHz, en un rango ahora comúnmente llamado la radio de baja frecuencia. Desde entonces, la radioastronomía ha crecido rápidamente, poniéndose al día con observaciones ópticas de mayor frecuencia que han estado sucediendo durante siglos más. Otro avance fue la utilización de la interferometría radioeléctmica, donde se utilizan grupos de antenas separadas por grandes distancias para crear una abertura sintética, proporcionando una manera de ampliar la sensibilidad y resolución de las observaciones radioeléctosas2,3. Esto puede considerarse intuitivamente como una extensión de la fórmula de resolución regular para observaciones ópticas:
Para un plato de observación de metros de tamaño D, y una longitud de onda de observación de λ metros, ḤHPBW es el tamaño angular en radianes de la Anchura del haz de media potencia (HPBW), definiendo la resolución en el cielo. Este proceso de sintetizar una fracción de un plato completo grande con sólo puntos dispersos a través de un área en su mayoría vacía también se llama síntesis de apertura. En el reino de la interferometría radioeléctrica, la resolución de una matriz viene determinada por la distancia más lejana entre los dos receptores de la matriz, y esta distancia se utiliza como D en la ecuación 1.
Las matemáticas detrás de la interferometría han sido bien documentadas en textos clásicos como Interferometría y Síntesis de Thompson en Radio Astronomía3. La información básica se puede comunicar informalmente como “(para matrices planas que observan un pequeño campo de visión) la correlación cruzada de señales entre cualquier 2 receptores (una visibilidad)producirá información sobre un coeficiente 2D Fourier del patrón de brillo del cielo.” Lo que se muestrea en el modo Fourier depende de la separación de los receptores (la línea base),normalizada por la longitud de onda de observación. Los receptores que están más separados (en el sistema de coordenadas UVW estándar orientados hacia el objetivo de imagen) muestran entidades de mayor frecuencia espacial, lo que produce detalles de mayor resolución a escalas más pequeñas. Por el contrario, los receptores que están unidos en la misma muestra de trama UVW menores frecuencias espaciales, dando información de estructuras a mayor escala a una resolución más baja.
Para las frecuencias de radio más bajas, los electrones libres en la ionosfera de la Tierra evitan que las ondas de radio inferiores a 10 MHz viajen del espacio al suelo, yviceversa. Este llamado “corte ionosférico” ha impedido durante mucho tiempo las observaciones terrestres del cielo para este rango de frecuencia. La respuesta obvia a esta limitación es poner receptores de radio en el espacio donde puedan registrar datos libres de la influencia de la atmósfera terrestre y electrones libres en su ionosfera. Esto se ha hecho antes con antenas individuales en naves espaciales como Wind4 y STEREO5,que han revelado muchos procesos astrofísicos que producen emisiones en este rango de radio de baja frecuencia. Esto incluye las emisiones de las interacciones de los electrones con la magnetosfera de la Tierra, la aceleración de electrones de las erupciones solares y de la propia galaxia. Las observaciones de antena única pueden medir la densidad total de flujo de tales eventos, pero no pueden precisar de dónde proviene la emisión. Con el fin de localizar esta emisión de baja frecuencia y hacer imágenes en este régimen de frecuencia por primera vez, muchas antenas tendrán que ser enviadas al espacio y tener sus datos combinados para hacer una apertura sintética.
Hacer esto abriría una nueva ventana a través de la cual la humanidad puede observar el universo, permitiendo una serie de mediciones científicas que requieren imágenes del cielo en estas frecuencias más bajas. La Luna es un sitio posible para una abertura sintética en el espacio, y viene con pros y contras en comparación con las matrices de órbita de vuelo libre. El lado lejano lunar tiene una zona radioeléctosa única que bloquea toda la interferencia habitual procedente de señales artificiales, mientras que el lado cercano proporciona un lugar estático para la Tierra observando matrices, y si se construye en el punto lunar de la subtierra, la Tierra siempre estará en el cenit del cielo. Con una matriz estática, es más fácil obtener líneas base cortas para medir las emisiones a gran escala, ya que no están en peligro de colisionar, a diferencia de las matrices de vuelo libre. Los inconvenientes de una matriz lunar son principalmente dificultades en el costo y el poder. Una matriz a gran escala en la Luna requeriría una cantidad sustancial de infraestructura y dinero, mientras que las matrices en órbita más pequeñas requerirían muchos menos recursos. También está el tema del poder; la mayoría de los lugares de la Luna están expuestos a suficiente luz solar para la generación de energía solar durante 1/3 de cada día lunar. Sobrevivir a los grandes columpios de temperatura desde el día lunar hasta la noche también es una preocupación de ingeniería. Dejando a un lado estas dificultades, todavía existe el problema de asegurarse de que el diseño de matriz propuesto es adecuado para sus objetivos científicos especificados. La respuesta de cualquier matriz dada depende de la estructura de la emisión que se observa junto con la configuración y sensibilidad de la matriz.
A lo largo de las décadas se han elaborado varias matrices conceptuales para ir a la superficie lunar. Los primeros diseños no fueron los más detallados, pero aún así reconocieron los avances científicos que podrían lograr dichos conjuntos6,7,8,9,10. En los últimos años también se han presentado más matrices, algunas de las cuales, como FARSIDE11,DEX12y DALI13, buscan medir los valles de absorción de la señal de hidrógeno neutro de 21 cm en el rango de 10-40 MHz para sondear las llamadas “Edades Oscuras” y restringir los modelos cosmológicos del universo temprano. Otros como ROLSS14 llaman a rastrear ráfagas de radio solares brillantes tipo II en la heliosfera para identificar el sitio de aceleración de partículas energéticas solares dentro de eyecciones de masa coronal como su convincente caso científico. También se han descrito matrices a menor escala como el interferómetro de 2 elementos RIF15,que utilizaría un único lander y un rover en movimiento para muestrear muchas líneas base a medida que se mueve hacia afuera desde el lander. RIF se centra en la capacidad de hacer un mapa del cielo de estas frecuencias bajas por primera vez, y calcula la cobertura UV y el haz sintetizado para observaciones integradas.
Las matrices de radio basadas en el espacio también podrían permitir imágenes de baja frecuencia de galaxias de radio distantes para determinar campos magnéticos y mediciones astrométricas16. Las imágenes de baja frecuencia de estos cuerpos proporcionarían una imagen más completa de la física que rige estos sistemas, en particular produciendo datos de emisiones de sincrotrón para el extremo inferior de la distribución de energía de electrones. También hay una gama de varias emisiones magnetosféricas que se producen a estas frecuencias bajas, proporcionando tanto firmas globales (emisión constante de sincrotrón) como locales (ráfagas, radiación kilométrica auroral) de dinámicas de electrones que no son detectables desde el suelo17. Las emisiones registradas más brillantes de este tipo han venido de la Tierra y Júpiter, ya que estos son los planetas más cercanos con fuertes magnetosferas. Sin embargo, las matrices con suficiente sensibilidad y resolución podrían observar emisiones magnetosféricas de otros planetas exteriores, o incluso planetas extrasolares18. Este tema en particular fue calificado como un área de interés en el reciente taller Planetary Sciences Vision 2050.
Este trabajo se centra en simular la respuesta de las matrices de radio en la Luna que consisten en cualquier lugar, desde sólo unas pocas antenas, hasta cientos o miles de antenas. Este marco de simulación es útil para iterar el diseño de matrices para crear imágenes de cualquier objetivo científico determinado en un pequeño campo de visión (unos pocos grados cuadrados), pero actualmente no admite todas las imágenes del cielo. Las estimaciones precisas de los mapas de brillo predichos junto con perfiles de ruido realistas deben utilizarse para garantizar que un tamaño/configuración de matriz determinado sea suficiente para observar el objetivo a un cierto nivel de ruido o resolución. La geometría de la matriz también debe conocerse en un alto grado para que las líneas base se calcule con precisión para habilitar la correcta toma de imágenes de los datos. Actualmente, los mejores mapas de la superficie lunar son los modelos de elevación digital (DEM) del Orbitador de Reconocimiento Lunar (LRO)19 Lunar Orbiter Altimeter (LOLA)20. La canalización de simulación acepta coordenadas de latitud longitud para cada receptor e interpola la elevación en estos puntos desde los DEMs existentes para calcular la posición 3D completa.
A partir de estas coordenadas, las líneas base se calculan e insertan en un archivo common astronomy software applications (CASA)21 Measurement Set (MS). El formato ms se puede utilizar con muchos algoritmos de análisis e imágenes existentes, y contiene información sobre la configuración de la matriz, los datos de visibilidad y la alineación con el cielo. Sin embargo, muchas de estas rutinas de software están codificadas arduamente para trabajar con matrices que giran con la superficie de la Tierra, y no funcionan para orbitar o matrices lunares. Para evitar esto, esta canalización calcula manualmente las líneas base y las visibilidades de una matriz determinada y un destino de imágenes, e inserta los datos en el formato ms. La biblioteca SPICE22 se utiliza para alinear correctamente los sistemas de coordenadas lunares y del cielo y rastrear los movimientos de la Luna, la Tierra y el Sol.
El marco de simulación descrito aquí sigue a Hegedus et al.17, y el software es archivado por la biblioteca de la Universidad de Michigan en el archivo Deep Blue23,almacenado en https://deepblue.lib.umich.edu/data/concern/data_sets/bg257f178?locale=en. Cualquier parche o actualización de este software archivado se puede encontrar en https://github.com/alexhege/LunarSynchrotronArray. La siguiente sección describirá los requisitos de este software y recorrerá el proceso de formación de una matriz, estableciendo los niveles de ruido adecuados, alimentando a la matriz con una imagen de verdad simulada de la emisión objetivo y simulando las reconstrucciones ruidosas y ruidosas de la matriz de la emisión utilizando un script CASA.
Cada paso de la canalización de simulación es necesario y se alimenta en el siguiente, tomando una configuración de matriz en la superficie lunar, alineando el marco de referencia correctamente para orientar la matriz al área de destino en el cielo, calculando los datos de visibilidad, agregando los niveles de ruido adecuados y ejecutando algoritmos de imágenes en los datos resultantes.
Para cada paso, se pueden realizar personalizaciones. En el paso 2, la configuración de matriz definida por el usuario puede ser cualquier lista de longitudes y latitudes. Esto luego se alimenta en el script SPICE en el paso 3, donde uno puede elegir la hora exacta de las mediciones planificadas, así como donde en el cielo la matriz debe ser enfocada. En el paso 3, se puede especificar la emisión de verdad simulada que la matriz está intentando crear imágenes proporcionando un archivo CASA .truth adecuado. A continuación, en el paso 4 uno puede cambiar el nivel esperado de ruido dependiendo de la frecuencia de observación y las capacidades de hardware esperadas. Este conjunto de códigos constituye un marco de simulación flexible que se puede utilizar para iterar el diseño de matrices para cualquier número de usos, dependiendo de la ciencia de destino. Todos estos códigos se pueden ejecutar en un portátil o estación de trabajo promedio, aunque el tiempo de cálculo aumenta con el número de antenas. Las partes más lentas del proceso están prediciendo las visibilidades, seguidas de imágenes. Para matrices pequeñas, todo el proceso se puede hacer en minutos, mientras que para matrices más grandes de unos pocos cientos o miles de receptores, pueden ser necesarias horas o días.
Algunos pasos siguientes que podrían tomarse con esta canalización para aumentar su realismo incluyen la adición de un sistema de eliminación en primer plano dependiente del canal. Esto requiere la construcción de un modelo de cielo global, dominado a bajas frecuencias por la emisión de sincrotrón galáctico y algunas fuentes brillantes como cas A, el seguimiento de qué parte del cielo es visible para los receptores, y la convolvencia de ese patrón de brillo con el haz primario, con el centro de fase de la matriz alineado hacia el objetivo de imagen. Durante tiempos de integración más largos, el seguimiento del movimiento aparente del cielo también es un problema. Otra mejora que podría agregarse es un sistema transitorio de marcado de interferencia de eventos/radiofrecuencia (RFI) que puede eliminar los canales marcados de las imágenes normales y enviarlos a una canalización especializada que imágenes y caracteriza los datos marcados. Esta canalización de eventos transitorios podría entonces utilizar algoritmos especiales como uvmodelfit que pueden aprovechar la alta relación señal-ruido de estos eventos para caracterizarlos mejor que la resolución normal de la matriz27.
También hay efectos adicionales que deben tenerse en cuenta para una calibración completa de la matriz, uno de los cuales es el acoplamiento mutuo. Como se explica en Ellingson28,esto puede conducir a una disminución de la sensibilidad en las matrices si tienen receptores que están dentro de unas pocas longitudes de onda entre sí. Esto se ve en una disminución de la sensibilidad para la matriz, o equivalentemente, un aumento en el SEFD. Esto es especialmente cierto para las vigas de más de 10 grados de distancia del cenit. La matriz de ejemplo de este trabajo se dirige a la Tierra, que siempre está cerca del cenit por diseño, por lo que el acoplamiento mutuo no debe afectar a este objetivo de imagen en particular, pero los estudios del SEFD sobre toda la gama de ángulos y frecuencias de elevación tendrán que hacerse en la puesta en marcha de cualquier matriz real para desbloquear todo su potencial. Otra deficiencia de esta tubería de simulación de matriz se encuentra en los mapas de superficie lunares imperfectos utilizados. Los DEM de las mediciones LRO LOLA tienen en el mejor de los casos una resolución de 60×60 metros/píxel en los mapas de 512 píxeles/grado. Uno puede interpolar estos datos para matrices simuladas, pero para matrices reales tendrá que haber un período de puesta en marcha / calibración donde se utilizarán fuentes con una posición conocida para determinar las separaciones relativas entre toda la antena a alta precisión. Las posibles fuentes de calibración incluyen la Cas A, la emisión periódica de baja frecuencia de Júpiter o la Tierra, o potencialmente la Puerta Lunar29.
También está la respuesta de la superficie lunar a considerar. Hay una capa de topsoil lunar llamada regolito que actúa como un dieléctrico con pérdida que puede reflejar la emisión entrante con cierta eficiencia, por encima de la roca lunar que también puede reflejar la emisión entrante con una mejor eficiencia30,31. Esta respuesta depende de la temperatura ambiente y la frecuencia entrante, así como de la composición química del regolito. Los estudios30,31 han encontrado que a temperaturas más bajas por debajo de 100 K, el regolito es casi transparente para la emisión de radio, y la reflexión se produce a nivel básico con un coeficiente de reflexión de alrededor de 0.5-0.6. A temperaturas más altas 150-200 K, el regolito puede absorber la emisión y reflejar la radiación entrante en la superficie con un coeficiente de reflexión de alrededor de 0.2-0.3. A temperaturas superiores a 200 K, se encuentra que las propiedades dieléctricas del regolito disminuyen, y la variación de la reflexión puede ser ignorada. Estos efectos pueden reducir el área efectiva de la matriz, reduciendo la sensibilidad y requiriendo tiempos de integración más largos. Este efecto se puede modelar con paquetes de software de simulación electromagnética como NEC4.232 dados modelos de permisividad relativa / constante dieléctrica en función de la profundidad lunar. Esto generará el SEFD de un receptor para una frecuencia determinada, que se puede dar a la canalización de simulación de matriz para calcular el ruido correcto para agregar a la señal simulada. Agregar una rejilla de puesta a tierra entre el receptor y la superficie lunar puede ayudar a disminuir el efecto de las ondas reflejadas, pero añade su propio conjunto de complicaciones en forma de despliegue.
Muchos de los detalles hipotéticos o difusos en torno a la implementación de un receptor de radio en la superficie lunar finalmente se consolidarán en realidad con la financiación reciente de proyectos de antenas de baja frecuencia como observaciones de ondas de radio en la superficie lunar de la vaina fotoelectron (ROLSES) y el Experimento electromagnético de superficie lunar (LuSEE)33. LuSEE fue financiado recientemente por la NASA como parte del programa Commercial Lunar Payload Services. Ambas suites de antena consistirán principalmente en repuestos de vuelo para instrumentos pasados como STEREO/WAVES o PSP FIELDS y están previstas para una entrega en 2021. Las mediciones de estos receptores finalmente consolidarán el nivel de ruido cuastérmico de la vaina fotoelectronal del polvo ionizado en la superficie lunar y cómo cambia en el transcurso de un día lunar. Estas mediciones también caracterizarán el nivel de reflexión y absorción de la superficie lunar, y cuantificarán cómo cambia el SEFD del receptor. También proporcionarán estadísticas sobre el número de eventos transitorios o RFI que se reciben en la superficie lunar. Estas misiones allanarán el camino para conjuntos de antenas que finalmente podrán hacer una multitud de observaciones científicas novedosas como la emisión de baja frecuencia de ráfagas de radio solar, galaxias lejanas y magnetosferas planetarias. La canalización de simulación descrita en este trabajo proporciona una manera flexible de iterar el diseño de estos futuros arreglos de discos para una variedad de objetivos científicos.
The authors have nothing to disclose.
Gracias a los equipos Lunar Reconnaissance Orbiter (LRO) y Lunar Orbiter Laser Altimeter (LOLA) por proporcionar los mapas de elevación digital lunar. Este trabajo fue apoyado directamente por el acuerdo de cooperación del Instituto Virtual de Investigación del Sistema Solar de la NASA número 80ARC017M0006, como parte del equipo de la Red de Exploración y Ciencia Espacial (NESS).
No physical materials are needed, this is a purely computational work. |