Optimal Transport Reconstruction of Biased Tracers in Primordial Non-Gaussian Fields

Resumo Técnico: Reconstrução por Transporte Ótimo de Traçadores Enviesados em Campos Não-Gaussianos Primordiais

Enunciado do Problema
O Transporte Ótimo (OT) estabeleceu-se como um arcabouço eficiente para reconstruir o campo de deslocamento cósmico, mapeando posições iniciais para localizações atuais (ou vice-versa). Trabalhos anteriores demonstraram que o OT recupera com precisão os deslocamentos para traçadores enviesados sob a suposição de condições iniciais Gaussianas. No entanto, uma questão aberta crítica é: este arcabouço permanece válido quando as condições iniciais são não-Gaussianas? Este é um desafio oportuno, dado o crescente esforço para restringir a não-Gaussianidade primordial (PNG) usando levantamentos cosmológicos de próxima geração. Uma dificuldade específica surge porque as implementações de OT frequentemente assumem uma distribuição espacial inicial uniforme. Embora isso seja válido para o campo de matéria escura completo, não se aplica a subconjuntos enviesados (ex: halos). Em cenários Gaussianos, a massa não observada ("poeira") pode ser modelada efetivamente; contudo, em cenários não-Gaussianos, a própria poeira adquire vieses dependentes de escala específicos que devem ser contabilizados para garantir uma reconstrução precisa.

Metodologia
Os autores investigam a aplicabilidade da reconstrução por OT a campos não-Gaussianos primordiais "locais", onde o potencial gravitacional $\phi_{NG}$ é modelado como uma série de potências de um campo Gaussiano $\phi_G$ (truncada no segundo ordem, $f_{NL} \neq 0$). Neste modelo, o contraste de densidade $\delta_{NG}$ contém um termo proporcional ao produto de campos Gaussianos, levando a um viés dependente de escala distintivo no agrupamento de traçadores enviesados em grandes escalas (pequeno $k$).

O estudo utiliza o conjunto de simulações Quijote-PNG, consistindo em 20 caixas periódicas ($L = 1 h^{-1}$ Gpc) com $512^3$ partículas, evoluindo de $z=127$ para $z=0$ sob uma cosmologia $\Lambda$CDM plana. Os traçadores enviesados são identificados como halos com massa $> 20 m_p$.

A inovação metodológica central envolve a modelagem da "poeira" não observada (a fração de massa $1-f_b$ não presente no catálogo de halos) para satisfazer a restrição de que o campo de matéria completo deve exibir nenhum viés dependente de escala. Os autores derivam que, se os traçadores enviesados exibem um viés dependente de escala caracterizado pelos parâmetros $b_\delta$ e $b_\phi$ (onde o viés escala como $b(k) \approx b_\delta + 2f_{NL}b_\phi \alpha(k)$), a poeira deve possuir um viés dependente de escala compensatório com sinal oposto.

Duas abordagens de reconstrução são testadas:

1. OT-K: Utiliza um modelo de filtro de Wiener para estimar a distribuição da poeira com base nos traçadores observados. Isso requer supor a combinação $2f_{NL}b_\phi$ (que pode ser restringida por dados de levantamentos) e a fração de massa $f_b$.
2. OT-W: Utiliza as partículas de poeira reais da simulação (uma amostra aleatória) combinadas com os halos ponderados pela massa para realizar a reconstrução, servindo como um limite superior de desempenho.

A qualidade da reconstrução é avaliada comparando as posições reconstruídas com as posições reais dos protohalos e analisando a dependência de escala da divergência do campo de deslocamento ($\theta_{OT}$), que serve como um proxy para o campo de densidade inicial.

Principais Contribuições e Resultados

• Caracterização da Poeira: O artigo demonstra que, em modelos de PNG locais, a "poeira" não é um fundo passivo, mas exibe seu próprio viés característico dependente de escala. Especificamente, a amplitude do viés da poeira é determinada pelo requisito de que a soma dos campos de traçadores enviesados e da poeira resulte em um campo total independente de escala.
• Extensão do OT para Não-Gaussianidade: Os autores mostram que a metodologia de OT estende-se com sucesso para condições iniciais não-Gaussianas. Ao incorporar um modelo realista para a poeira que contabiliza o viés dependente de $f_{NL}$, a reconstrução dos traçadores enviesados permanece precisa.
• Precisão da Reconstrução:
  • Posições dos Traçadores: A reconstrução OT-K (usando o modelo de poeira do filtro de Wiener) recupera corretamente o viés dependente de escala dos protohalos. A reconstrução OT-W (usando a poeira real) mostra um acordo ainda melhor, indicando que o método é robusto, mas sensível à fidelidade do modelo de poeira.
  • Proxy do Campo Completo: A divergência do campo de deslocamento OT ($\theta_{OT}$) mostra quase nenhuma dependência de escala em pequeno $k$, espelhando o comportamento do verdadeiro campo inicial. O modelo OT-W estende essa independência de escala para valores de $k$ maiores em comparação ao OT-K, alcançando uma melhoria de um fator de 2 na faixa de $k$ onde $b(k) \approx 1$.
• Robustez a $f_{NL}$: O valor médio do quadrado do campo reconstruído $\langle \theta_{OT}^2 \rangle$ permanece estável através de diferentes valores de $f_{NL}$ ($+100, 0, -100$), com variações dentro do desvio padrão, sugerindo que o método não é enviesado pela presença da não-Gaussianidade em si.

Significância e Alegações
O artigo afirma que o Transporte Ótimo fornece um arcabouço útil para reconstruir estruturas cósmicas mesmo quando as condições iniciais são não-Gaussianas, desde que o modelo para a "poeça" não observada seja suficientemente realista. Os autores enfatizam que o sucesso do método reside no fato de que a assinatura distintiva dependente de escala da PNG está presente tanto nos traçadores enviesados quanto na poeira, mas com sinais opostos, permitindo que o campo completo permaneça sem viés.

Os autores afirmam modestamente que não demonstram uma melhoria direta nas restrições sobre $f_{NL}$ a partir deste trabalho isoladamente. Em vez disso, eles fornecem um arcabouço fundamental que permite que análises baseadas em reconstrução sejam aplicadas consistentemente em cosmologias de PNG. Eles sugerem que este trabalho é um passo necessário para incorporar técnicas de reconstrução em futuras análises de alta precisão de PNG, potencialmente ajudando a conectar traçadores enviesados observados mais de perto aos seus contrapartes lagrangianos primordiais e mitigando contribuições não lineares ou estocásticas sutis relevantes para levantamentos de próxima geração. O artigo conclui que o melhor modelamento da poeira (ex: via métodos iterativos ou aprendizado de máquina) poderia gerar ganhos adicionais.