Testing masking effectiveness using multi-line image cubes based on COSMOS2020 for [CII] line intensity mapping at $z_{[CII]} > 3.5$

Resumo Técnico: Testando a Eficácia da Máscara Usando Cubos de Imagem Multilinhas Baseados no COSMOS2020 para o Mapeamento de Intensidade da Linha [CII] em $z_{[CII]} > 3,5$

Declaração do Problema
O Mapeamento de Intensidade de Linhas (LIM) oferece uma capacidade única de sondar a distribuição tridimensional de galáxias e do meio intergaláctico durante a Época da Reionização (EoR) e as eras pós-EoR sem resolver fontes individuais. Um alvo principal para futuros experimentos de LIM, como o módulo EoR-Spec no Telescópio Submilimétrico Fred Young (FYST), é a linha de estrutura fina de carbono ionizado ([CII]) em $z \approx 3,5 - 8,2$. No entanto, as observações de LIM no regime submilimétrico são fortemente contaminadas por emissão de primeiro plano, particularmente pela transição rotacional de Monóxido de Carbono (CO) originada de galáxias de redshifts mais baixos ($z < 2$). Diferente de levantamentos tradicionais, o LIM não consegue distinguir inerentemente o sinal [CII] alvo da emissão de CO contaminante dentro dos mesmos canais de frequência. Simulações anteriores careciam de contaminação de primeiro plano realista ou não testavam totalmente a eficácia das técnicas de máscara necessárias para recuperar o espectro de potência do [CII]. Este trabalho aborda o desafio de recuperar sinais de [CII] na presença de significativa contaminação de CO e ruído instrumental, testando especificamente se as estratégias de máscara podem isolar efetivamente o sinal alvo no campo E-COSMOS.

Metodologia
Os autores estenderam o framework de simulação empírica estabelecido por Clarke et al. (2024, C24) para criar cubos de intensidade multilinhas para a faixa de frequência do EoR-Spec (205–420 GHz) e duas extensões hipotéticas de frequências mais baixas (90 e 150 GHz).

1. Construção de Dados e Amostra:

   • O estudo utilizou o catálogo de galáxias COSMOS2020 (versão 4.1.11) cobrindo uma área de 1,44 deg$^2$ (E-COSMOS).
   • Para abordar a incompletude intrínseca do catálogo, os autores geraram dois tipos de amostras: uma amostra de "limite inferior" baseada estritamente em galáxias observadas (FARMER LP) e uma amostra "extrapolada". A extrapolação envolveu a duplicação de galáxias para corresponder às estatísticas de completude de massa derivadas de CANDELS e o ajuste de funções de Schechter à função de massa estelar, adicionando assim galáxias tênues e brilhantes para representar melhor a verdadeira população do céu.
   • Propriedades globais (Massa Estelar, SFR, $L_{IR}$) foram usadas para atribuir luminosidades de linha.

2. Modelagem de Luminosidade de Linha:

   • [CII]: Três modelos foram aplicados para delimitar a incerteza: um limite inferior conservador baseado em galáxias anãs locais (De Looze et al. 2014), um modelo de meio-termo baseado em SFR (Schaerer et al. 2020) e um limite superior baseado em galáxias de alto redshift ALPINE-ALMA (C24_m3).
   • CO: Modelos para a transição $J=1-0$ foram derivados de relações $L_{IR}$ (Sargent et al. 2014; Li et al. 2016). Outras transições de CO ($J=2-1$ até $J=13-12$) foram derivadas usando templates de Distribuição de Energia de Linhas Espectrais (SLED), misturando perfis de galáxias brilhantes (NGC 253) e tênues (Via Láctea) com base na posição da galáxia na sequência principal.
   • Construção do Cubo: As luminosidades foram convertidas em intensidades e inseridas em voxels 3D. Os cubos foram convolvidos com o feixe do instrumento (Gaussiano 2D no espaço do mapa, Lorentziano 1D no espaço da frequência) e sub-gradeados por $3\times3\times3$ para facilitar a máscara.

3. Modelagem de Ruído:

   • O estudo incorporou tanto ruído branco quanto ruído correlacionado (ruído 1/f) baseado em simulações de Ferramentas de Astrofísica Escaláveis em Tempo Ordenado (TOAST) para o EoR-Spec. Os cubos de ruído foram escalonados para representar o período completo de 2000 horas de observação.

4. Técnicas de Máscara:

   • Máscara Direcionada (Targeted Masking): Utilizou o catálogo COSMOS2020 como uma máscara de primeiro plano. Galáxias brilhantes conhecidas por emitir CO foram identificadas, e voxels ao redor delas foram mascarados. Os autores testaram vários raios de máscara (0,5$\sigma$ a 2,5$\sigma$) e determinaram iterativamente uma ordem de máscara ideal baseada em um parâmetro de eficiência (remoção de CO por unidade de volume mascarado).
   • Máscara Cega (Blind Masking): Não dependeu de um catálogo externo. Em vez disso, identificou máximos locais nos cubos de intensidade total através de múltiplas bandas de frequência. Voxels foram pareados se sua separação espacial e diferença de frequência correspondiam a uma única galáxia emitindo em duas linhas diferentes (por exemplo, diferentes transições de CO ou CO e [CII]). Esta técnica foi testada com e sem as bandas hipotéticas de 90 e 150 GHz.

Principais Resultados

• Força do Sinal CO vs. [CII]:

  • Acima de 300 GHz ($z_{[CII]} < 5$), as intensidades previstas de [CII] e CO são comparáveis (dentro de $\sim0,5$ dex), com o CO dominando em frequências mais baixas.
  • Abaixo de 300 GHz, a emissão de CO excede o [CII] em mais de 1 dex em modelos não extrapolados, tornando a recuperação via máscara extremamente difícil sem cobertura de frequência adicional.
  • Os espectros de potência de CO geralmente excedem os limites de sensibilidade atmosférica previstos para $k > 0,5$ Mpc$^{-1}$, mas modelos mais tênues caem abaixo do nível de ruído em modos de baixo $k$.

• Eficácia da Máscara Direcionada:

  • Sucesso: A máscara direcionada recuperou com sucesso os espectros de potência de [CII] acima de 300 GHz quando utilizando máscaras largas ($\ge 1,5\sigma$). Em condições ideais (catálogo de primeiro plano completo), o sinal [CII] recuperado foi distinto do CO residual.
  • Limitações: A 280 GHz, mesmo a máscara máxima falhou em recuperar o [CII] devido ao esmagador sinal de CO.
  • Sobre-mascaramento (Overmasking): O estudo identificou um risco específico de "sobre-mascaramento" (remover acidentalmente o sinal de [CII]) ao mascarar a transição de CO $J=5-4$ em 410 GHz com máscaras de galáxias brilhantes. Isso sugere que transições específicas devem ser excluídas das estratégias de máscara no campo E-COSMOS.
  • Completude do Catálogo: Se o catálogo de primeiro plano for incompleto (faltando galáxias extrapoladas tênues), a máscara direcionada torna-se significativamente menos eficaz, pois fontes de CO não mascaradas permanecem.

• Eficácia da Máscara Cega:

  • A máscara cega provou ser superior quando o catálogo de primeiro plano é incompleto, pois não depende de conhecimento prévio da localização das galáxias.
  • A técnica é mais eficaz ao utilizar as bandas hipotéticas de 90 e 150 GHz, que permitem um pareamento de voxels mais robusto entre frequências. Sem essas bandas, a máscara cega é menos eficiente, mas ainda supera a máscara direcionada com um catálogo incompleto.
  • A máscara cega é altamente sensível ao ruído; acima de frequências altas (>300 GHz), picos de ruído podem ser confundidos com emissão de linha, reduzindo a eficácia da técnica em até 0,5 dex.

• Impacto do Ruído:

  • O ruído combinado (branco e correlacionado) dificulta a recuperação de [CII] acima de 300 GHz. Alcançar uma razão Sinal-Ruído (S/N) $>5$ para o sinal de agrupamento (clustering) exigiria aproximadamente cinco vezes o tempo de observação planejado atualmente (ou melhor sensibilidade do instrumento).
  • Embora o CO possa ser recuperado abaixo de 300 GHz, a recuperação de [CII] neste regime permanece improvável sem técnicas de correlação cruzada devido à dominância do CO e do ruído.

Significância e Alegações
Os autores alegam que seu trabalho fornece uma avaliação realista e empiricamente fundamentada da viabilidade do LIM de [CII] com o instrumento EoR-Spec. O estudo demonstra que:

1. Viabilidade: Recuperar [CII] via máscara é viável para frequências acima de 300 GHz, desde que raios de máscara largos sejam usados e o catálogo de primeiro plano seja suficientemente completo.
2. Necessidade de Dados Multi-banda: A inclusão de bandas de frequências mais baixas (90/150 GHz) ou observações de cross-instrumento é crítica para o sucesso da máscara cega e para limpar efetivamente a contaminação de CO.
3. Validação Metodológica: A abordagem empírica usando COSMOS2020, combinada com extrapolação e templates SLED, produz estimativas de CO consistentes com observações de ALMA, VLA e NOEMA, validando o framework de simulação.
4. Orientação Operacional: O artigo identifica restrições operacionais específicas, como a necessidade de evitar o mascaramento da transição de CO $J=5-4$ em E-COSMOS para prevenir o sobre-mascaramento, e destaca que os níveis atuais de ruído provavelmente exigirão tempos de observação estendidos ou mitigação de ruído avançada (ex: limpeza por PCA) para alcançar detecções robustas.

Os autores concluem que, embora a mitigação de contaminantes via máscara direcionada e cega seja uma estratégia viável para o EoR-Spec, o sucesso do experimento depende fortemente da completude dos catálogos de primeiro plano, da disponibilidade de dados de múltiplas frequências e do gerenciamento do ruído instrumental. Eles posicionam este trabalho como uma fundação para futuros estudos de correlação cruzada e técnicas de desconfusão mais avançadas.