Uma estrutura de simulação para testar as capacidades de imagem de matrizes de rádio em larga escala na superfície lunar é apresentada. Os principais componentes de ruído são discutidos, e um pipeline de software é acompanhado com detalhes sobre como personalizá-lo para novos usos científicos.
Nos últimos anos, houve um interesse renovado em retornar à Lua por razões científicas e exploratórias por natureza. A Lua fornece o campo de treinamento perfeito para construir bases em grande escala que se pode aplicar a outros planetas como Marte. A existência de uma zona de rádio tranquila no lado lunar tem prometido para estudos iniciais do universo e pesquisas de exoplanetas, enquanto o lado próximo fornece uma base estável que pode ser usada para observar emissões de baixa frequência da magnetosfera da Terra que podem ajudar a medir sua resposta ao clima espacial que está chegando. A construção de uma matriz de rádio em larga escala proporcionaria grandes retornos científicos, bem como agiria como um teste da capacidade da humanidade de construir estruturas em outros planetas. Este trabalho se concentra em simular a resposta de pequenas e grandes matrizes de rádio na Lua, consistindo de centenas ou milhares de antenas. A resposta do array depende da estrutura da emissão, juntamente com a configuração e sensibilidade do array. Um conjunto de locais são selecionados para os receptores de rádio simulados, usando modelos de elevação digital do instrumento Lunar Orbiter Laser Altimeter no Lunar Reconnaissance Orbiter para caracterizar a elevação dos locais receptores. Um código personalizado de aplicativos de software de astronomia é descrito e usado para processar os dados dos receptores simulados, alinhando os quadros de coordenadas lunares e do céu usando SPICE para garantir que as projeções adequadas sejam usadas para imagens. Esta estrutura de simulação é útil para iterar o design de matriz para imagem de qualquer alvo científico em um pequeno campo de visão. Este quadro não suporta atualmente todas as imagens do céu.
O campo da radioastronomia começou em 1932 com a detecção acidental de emissão de rádio galáctica por Karl G. Jansky1 a 20 MHz, em uma faixa agora comumente chamada de rádio de baixa frequência. Desde então, a radioastronomia tem crescido rapidamente, alcançando observações ópticas de maior frequência que vêm acontecendo há séculos. Outro avanço foi a utilização da interferometria de rádio, onde grupos de antenas separadas por grandes distâncias são usados para criar uma abertura sintética, proporcionando uma maneira de ampliar a sensibilidade e resolução das observações de rádio2,3. Isso pode ser intuitivamente pensado como uma extensão da fórmula de resolução regular para observações ópticas:
Para um prato de observação de medidores de tamanho D, e um comprimento de onda de observação de λ metros, ΦHPBW é o tamanho angular em radians da Half Power Beam Width (HPBW), definindo a resolução no céu. Este processo de sintetizar uma fração de um prato grande completo com apenas pontos espalhados em uma área principalmente vazia também é chamado de síntese de abertura. No reino da interferometria de rádio, a resolução de uma matriz é determinada pela distância mais distante entre quaisquer dois receptores na matriz, e esta distância é usada como D na Equação 1.
A matemática por trás da interferometria foi bem documentada em textos clássicos como Interferometria de Thompson e Síntese em RadioAstronomia3. A visão básica pode ser comunicada informalmente como “(para matrizes planares observando um pequeno campo de visão) a correlação cruzada de sinais entre 2 receptores (uma visibilidade) produzirá informações sobre um coeficiente 2D Fourier do padrão de brilho do céu.”. O modo Fourier é amostrado depende da separação dos receptores (a linha de base),normalizada pelo comprimento de onda de observação. Receptores mais distantes (no sistema de coordenadas UVW padrão orientado para o alvo de imagem) amostram características de frequência espacial mais altas, produzindo detalhes de resolução mais elevados em escalas menores. Por outro lado, receptores que estão próximos na mesma amostra de quadro UVW têm frequências espaciais mais baixas, fornecendo informações de estruturas de maior escala em menor resolução.
Para as frequências de rádio mais baixas, elétrons livres na ionosfera da Terra impedem ondas de rádio abaixo de 10 MHz de viajar do espaço para o solo, e vice-versa. Este chamado “corte ionosférico” há muito impede observações terrestres do céu para esta faixa de frequência. A resposta óbvia para essa limitação é colocar receptores de rádio no espaço onde eles podem registrar dados livres da influência da atmosfera da Terra e elétrons livres em sua ionosfera. Isso já foi feito antes com antenas únicas em naves espaciais como Wind4 e STEREO5,que revelaram muitos processos astrofísicos que produzem emissões nesta faixa de rádio de baixa frequência. Isso inclui emissões das interações dos elétrons com a magnetosfera da Terra, aceleração eletrônica de erupções solares, e da própria galáxia. Observações de antena única podem medir a densidade total de fluxo de tais eventos, mas não podem identificar de onde vem a emissão. Para localizar essa emissão de baixa frequência e fazer imagens nesse regime de frequência pela primeira vez, muitas antenas terão que ser enviadas ao espaço e ter seus dados combinados para fazer uma abertura sintética.
Fazer isso abriria uma nova janela através da qual a humanidade pode observar o universo, permitindo uma série de medições científicas que requerem imagens do céu nessas frequências mais baixas. A Lua é um local possível para uma abertura sintética no espaço, e vem com prós e contras quando comparado com matrizes em órbita voadoras livres. O lado lunar tem uma zona de silêncio de rádio única que bloqueia toda a interferência usual proveniente de sinais feitos pelo homem, enquanto o lado próximo fornece um lugar estático para matrizes de observação da Terra, e se construído no ponto lunar sub-Terra, a Terra estará sempre no zênite do céu. Com uma matriz estática, é mais fácil obter linhas de base curtas para medir as emissões em larga escala, pois elas não correm o risco de colidir, ao contrário das matrizes de voo livre. As desvantagens de uma matriz lunar são principalmente dificuldades em custo e energia. Uma matriz em grande escala na Lua exigiria uma quantidade substancial de infraestrutura e dinheiro, enquanto matrizes menores em órbita exigiriam muito menos recursos. Há também a questão do poder; a maioria dos lugares na Lua são expostos à luz solar suficiente para a geração de energia solar para 1/3 de cada dia lunar. Sobreviver às grandes oscilações na temperatura do dia para a noite lunar também é uma preocupação de engenharia. Deixando de lado essas dificuldades, ainda há o problema de garantir que o projeto de matriz proposto seja adequado para seus alvos científicos especificados. A resposta de qualquer matriz depende da estrutura da emissão que está sendo observada juntamente com a configuração e sensibilidade da matriz.
Várias matrizes conceituais para ir na superfície lunar foram elaboradas ao longo das décadas. Os primeiros desenhos não foram os mais detalhados, mas ainda reconheceram os avanços científicos que poderiam ser alcançados por tais matrizes6,7,8,9,10. Mais arrays também foram apresentados nos últimos anos, alguns dos quais, como FARSIDE11, DEX12, e DALI13 buscam medir os cochos de absorção do sinal de hidrogênio neutro de 21 cm na faixa de 10-40 MHz para sondar a chamada “Idade das Trevas” e restringir modelos cosmológicos do universo primitivo. Outros como o ROLSS14 chamam o rastreamento de ondas de rádio solar brilhante tipo II para a heliosfera para identificar o local da aceleração de partículas energéticas solares dentro de ejeções de massa coronal como seu caso científico convincente. Matrizes de escala menor também foram descritas como o interferômetro RIF15de 2 elementos, que usaria um único lander e um rover em movimento para provar muitas linhas de base à medida que se move para fora do módulo de pouso. O RIF concentra-se na capacidade de fazer um mapa do céu dessas baixas frequências pela primeira vez, e calcula a cobertura UV e o feixe sintetizado para observações integradas.
Matrizes de rádio baseadas no espaço também poderiam permitir imagens de baixa frequência de galáxias de rádio distantes para determinar campos magnéticos e medições astrométricas16. Imagens de baixa frequência desses corpos forneceriam uma imagem mais completa da física que rege esses sistemas, em particular produzindo dados de emissão síncrotron para a extremidade inferior da distribuição de energia eletrônica. Há também uma gama de várias emissões magnetosféricas que ocorrem nessas baixas frequências, fornecendo tanto assinaturas globais (emissões síncrotrons constantes) quanto locais (rajadas, radiação quilométrica auroral) da dinâmica eletrônica que não são detectáveis a partir do solo17. As emissões mais brilhantes registradas desses tipos vieram da Terra e de Júpiter, pois estes são os planetas mais próximos com fortes magnetosferas. No entanto, matrizes com sensibilidade e resolução suficientes poderiam observar a emissão magnetosférica de outros planetas externos, ou mesmo planetas extrasolares18. Este tema, em particular, foi chamado como uma área de interesse no recente workshop Planetry Sciences Vision 2050.
Este trabalho se concentra em simular a resposta de matrizes de rádio na Lua que consistem em qualquer lugar, desde apenas algumas antenas, até centenas ou milhares de antenas. Esta estrutura de simulação é útil para iterar o design de matriz para imagem de qualquer alvo científico em um pequeno campo de visão (alguns graus quadrados), mas atualmente não suporta todas as imagens do céu. Estimativas precisas dos mapas de brilho previstos, juntamente com perfis de ruído realistas devem ser usadas para garantir que um determinado tamanho/configuração da matriz seja suficiente para observar o alvo a um determinado nível de ruído ou resolução. A geometria da matriz também deve ser conhecida em alto grau para que as linhas de base sejam calculadas com precisão para permitir a correta imagem dos dados. Atualmente, os melhores mapas da superfície lunar são Modelos de Elevação Digital (DEMs) do Lunar Reconnaissance Orbiter (LRO’s)19 Lunar Orbiter Laser Altimeter (LOLA)20. O pipeline de simulação aceita coordenadas de latitude de longitude para cada receptor e interpola a elevação nesses pontos dos DEMs existentes para calcular a posição 3D completa.
A partir dessas coordenadas, as linhas de base são computadas e inseridas em um arquivo Common Astronomy Software Applications (CASA)21 Measurement Set (MS). O formato MS pode ser usado com muitos algoritmos de análise e imagem existentes, e contém informações sobre a configuração do array, dados de visibilidade e alinhamento com o céu. No entanto, muitas dessas rotinas de software são codificadas para trabalhar com matrizes que giram com a superfície da Terra, e não funcionam para matrizes orbitais ou lunares. Para contornar isso, este pipeline calcula manualmente as linhas de base e as visibilidades para um determinado array e alvo de imagem, e insere os dados no formato MS. A biblioteca SPICE22 é usada para alinhar corretamente os sistemas de coordenadas lunares e céus e rastrear os movimentos da Lua, Terra e Sol.
A estrutura de simulação descrita aqui segue Hegedus et al.17, e o software é arquivado pela biblioteca da Universidade de Michigan no arquivo Deep Blue23, armazenado em https://deepblue.lib.umich.edu/data/concern/data_sets/bg257f178?locale=en. Quaisquer patches ou atualizações para este software arquivado podem ser encontrados em https://github.com/alexhege/LunarSynchrotronArray. A seção a seguir descreverá os requisitos deste software e percorrerá o processo de formação de uma matriz, definindo os níveis de ruído apropriados, alimentando a matriz de uma imagem simulada da verdade da emissão alvo e simulando as reconstruções sonoras e barulhentas da matriz da emissão usando um script CASA.
Cada etapa do pipeline de simulação é necessária e se alimenta na próxima, pegando uma configuração de matriz na superfície lunar, alinhando o quadro de referência corretamente para orientar a matriz para a área alvo no céu, calculando os dados de visibilidade, adicionando os níveis de ruído apropriados e executando algoritmos de imagem nos dados resultantes.
Para cada etapa, personalizações podem ser feitas. Na etapa 2, a configuração do array definido pelo usuário pode ser qualquer lista de longitudes e latitudes. Isso então alimenta o script SPICE na etapa 3, onde pode-se escolher a hora exata das medições planejadas, bem como onde no céu a matriz deve ser focada. Na etapa 3, pode-se especificar a emissão simulada de verdade que a matriz está tentando imaginar fornecendo um arquivo CASA .truth adequado. Em seguida, na etapa 4, pode-se alterar o nível esperado de ruído dependendo da frequência de observação e dos recursos de hardware esperados. Este conjunto de códigos constitui uma estrutura de simulação flexível que pode ser usada para iterar o design de matriz para qualquer número de usos, dependendo da ciência alvo. Esses códigos podem ser todos executados em um laptop ou estação de trabalho médio, embora o tempo de computação aumente com o número de antenas. As partes mais lentas do processo estão prevendo a visibilidade, seguida por imagens. Para pequenas matrizes, todo o processo pode ser feito em minutos, enquanto para matrizes maiores de algumas centenas ou milhares de receptores, horas ou dias podem ser necessários.
Alguns dos próximos passos que poderiam ser dados com este pipeline para aumentar seu realismo incluem a adição de um sistema de remoção em primeiro plano dependente do canal. Isso requer a construção de um modelo de céu global, dominado em baixas frequências por emissão síncrotron galáctica e algumas fontes brilhantes como Cas A, rastreando qual parte do céu é visível para os receptores, e convolving esse padrão de brilho com o feixe primário, com o centro de fase da matriz alinhado ao alvo de imagem. Por tempos de integração mais longos, rastrear o movimento aparente do céu também é um problema. Outra melhoria que pode ser adicionada é um sistema de sinalização de interferência de evento/radiofrequência transitória (RFI) que pode remover canais sinalizados de imagens normais e enviá-los para um pipeline especializado que imagens e caracteriza os dados sinalizados. Este pipeline de eventos transitórios poderia então usar algoritmos especiais como o uvmodelfit que pode aproveitar a alta relação sinal-ruído desses eventos para caracterizá-los melhor do que a resolução normal da matriz27.
Há também efeitos adicionais que precisam ser levados em conta para uma calibração completa da matriz, uma das quais é o acoplamento mútuo. Como discutido em Ellingson28, isso pode levar a uma diminuição da sensibilidade em matrizes se eles tiverem receptores que estão dentro de alguns comprimentos de onda um do outro. Isso é visto em uma diminuição da sensibilidade para a matriz, ou equivalentemente, um aumento no SEFD. Isso é especialmente verdade para feixes superiores a 10 graus de distância do zênite. O exemplo de matriz neste trabalho tem como alvo a Terra, que está sempre perto do zênite por design, de modo que o acoplamento mútuo não deve afetar esse alvo de imagem em particular, mas estudos do SEFD sobre toda a gama de ângulos de elevação e frequências precisarão ser feitos no comissionamento de qualquer matriz real para desbloquear todo o seu potencial. Outra deficiência deste pipeline de simulação de matriz está nos mapas imperfeitos da superfície lunar usados. OS DEMs de medições LRO LOLA têm, na melhor das hipóteses, uma resolução de 60×60 metros/pixel nos mapas de 512 pixels/grau. Pode-se interpolar esses dados para matrizes simuladas, mas para matrizes reais será necessário um período de comissionamento/calibração onde fontes com uma posição conhecida serão usadas para determinar as separações relativas entre toda a antena e alta precisão. As possíveis fontes de calibração incluem Cas A, emissão periódica de baixa frequência de Júpiter ou Terra, ou potencialmente o Gateway Lunar29.
Há também a resposta da superfície lunar a considerar. Há uma camada de topsoil lunar chamada regolito que age como um dielétrico lossy que pode refletir a emissão de entrada com alguma eficiência, acima da rocha lunar que também pode refletir a emissão de entrada com alguma melhor eficiência30,31. Esta resposta depende da temperatura ambiente e da frequência de entrada, bem como da composição química do regolito. Estudos30,31 descobriram que em temperaturas mais baixas abaixo de 100 K, o regolito é quase transparente para a emissão de rádio, e a reflexão ocorre no nível da base com um coeficiente de reflexão de cerca de 0,5-0,6. Em temperaturas mais altas de 150-200 K, o regolito pode absorver emissões e refletir a radiação de entrada na superfície com um coeficiente de reflexão em torno de 0,2-0,3. Em temperaturas acima de 200 K, constata-se que as propriedades dielétricas do regolito são diminuídas, e a variação da reflexão pode ser ignorada. Esses efeitos podem reduzir a área efetiva do array, reduzindo a sensibilidade e exigindo tempos de integração mais longos. Este efeito pode ser modelado com pacotes de software de simulação eletromagnética, como NEC4.232, dado modelos de relativa permissão/constante dielétrica em função da profundidade lunar. Isso irá produzir o SEFD de um receptor para uma determinada frequência, que pode ser dado ao pipeline de simulação de array para calcular o ruído correto para adicionar ao sinal simulado. Adicionar uma grade de aterramento entre o receptor e a superfície lunar pode ajudar a diminuir o efeito das ondas refletidas, mas adiciona seu próprio conjunto de complicações na forma de implantação.
Muitos dos detalhes hipotéticos ou confusos em torno da implementação de um receptor de rádio na superfície lunar finalmente se solidificarão em realidade com o financiamento recente de projetos de antenas de baixa frequência, como observações de ondas de rádio na superfície lunar da bainha de fotoEletrron (ROLSES) e o Experimento de Eletromagnética da Superfície Lunar (LuSEE)33. A LuSEE foi recentemente financiada pela NASA como parte do programa Commercial Lunar Payload Services. Ambas as suítes de antena consistirão principalmente de peças de voo para instrumentos passados como STEREO/WAVES ou PSP FIELDS e estão previstas para uma entrega em 2021. As medições desses receptores finalmente solidificarão o nível de ruído quase isentérmico da baia fotoeletrrona da poeira ionizada na superfície lunar e como ela muda ao longo de um dia lunar. Essas medidas também caracterizarão o nível de reflexão e absorção da superfície lunar, e quantificarão como ele altera o SEFD do receptor. Eles também fornecerão estatísticas sobre o número de eventos transitórios ou RFI que são recebidos na superfície lunar. Essas missões abrirão caminho para conjuntos de antenas que finalmente serão capazes de fazer uma infinidade de novas observações científicas, como emissão de baixa frequência de rajadas de rádio solar, galáxias distantes e magnetosferas planetárias. O pipeline de simulação descrito neste trabalho fornece uma maneira flexível de iterar o projeto dessas futuras matrizes para uma variedade de alvos científicos.
The authors have nothing to disclose.
Graças às equipes do Lunar Reconnaissance Orbiter (LRO) e do Lunar Orbiter Laser Altimeter (LOLA) para fornecer os Mapas de Elevação Digital Lunar. Este trabalho foi apoiado diretamente pelo acordo cooperativo do Instituto Virtual de Pesquisa de Exploração do Sistema Solar da NASA número 80ARC017M0006, como parte da equipe da Network for Exploration and Space Science (NESS).
No physical materials are needed, this is a purely computational work. |