PHANTOM: A MATLAB and Octave Toolbox Connecting Linear Field Statistics to Dark Matter Halo Observables

Resumo Técnico: PHANTOM – Um Toolbox de MATLAB e Octave para Análise de Halos de Matéria Escura

Definição do Problema
A estrutura interna e a abundância de halos de matéria escura (DM) são críticas para ligar as condições iniciais cosmológicas a observáveis de galáxias e aglomerados. Embora o toolkit em Python colossus tenha padronizado muitas dessas cálculos (Diemer 2018), uma parte significativa dos modelos de instrumentos astronômicos, pipelines de levantamento e códigos de ajuste depende do MATLAB. Atualmente, não existe um toolbox público e validado de MATLAB ou Octave que reproduza o conjunto de modelos do colossus ou que estenda estes para cenários alternativos de matéria escura, como Matéria Escura Morna (WDM) e Matéria Escura Difusa (FDM)/Matéria Escura de Campo Escalar (SFDM). Consequentemente, pesquisadores que utilizam fluxos de trabalho baseados em MATLAB enfrentam o fardo de portar manualmente os outputs do colossus ou implementar funções de ajuste fragmentadas da literatura, correndo riscos de inconsistências e erros de validação.

Metodologia e Arquitetura
O artigo apresenta o phantom (Profile and Halo Analysis for Numerous Theoretical dark Matter Observables), um toolbox de MATLAB e um pacote Octave projetados para fazer a ponte entre estatísticas de campos lineares e observáveis de halos. A arquitetura do software é construída em torno de uma única estrutura de cosmologia que armazena o background de expansão, o crescimento e os handles do espectro de potência linear. Este objeto é construído uma única vez e propagado através do grafo de chamadas, garantindo consistência em todos os cálculos subsequentes.

O pacote é organizado em cinco camadas lógicas:

1. Cosmologia: Implementa um framework padrão de Friedmann–Lemaître–Robertson–Walker (FLRW), suportando geometrias planas e não planas, vários modelos de energia escura (incluindo parametrização CPL) e componentes relativísticos. Ele computa taxas de expansão, tempos de lookback e distâncias via integração numérica direta.
2. Estatísticas de Campo: Computa o espectro de potência de matéria linear, a variância de matéria e funções de correlação. Suporta múltiplos funções de transferência (Eisenstein & Hu, CAMB, axionCAMB) e filtros de suavização (top-hat, Gaussian, sharp-k, smooth-k, vsmk).
3. Estatísticas de Halo: Calcula a função de massa de halo (HMF) e o viés (bias) linear de halo usando um sistema de despacho que aceita chaves de string para vários modelos (ex: Press–Schechter, Sheth–Tormen, Tinker, e prescrições específicas de WDM/FDM).
4. Estrutura de Halo: Determina relações de concentração–massa e constrói perfis de densidade normalizados (NFW, Hernquist, Einasto, Diemer–Kravtsov e perfis compostos de sol de FDM/solitão).
5. Observáveis: Deriva quantidades físicas como massa contida, velocidade circular, densidade superficial projetada e convergência de lente de densidade a partir dos perfis de densidade.

Uma característica fundamental de design é o uso de funções de despacho que permitem que modelos sejam trocados ou estendidos sem modificar os scripts voltados ao usuário. O código prioriza a transparência numérica; por exemplo, realiza integração numérica direta para variância e busca de raízes para relações de concentração–massa em vez de depender apenas de tabelas de interpolação pré-computadas, garantindo a autossuficiência em diferentes contextos de invocação.

Principais Contribuições

• Implementação Nativa: Fornece a primeira implementação nativa validada de MATLAB/Octave do conjunto de modelos do colossus, cobrindo cenários CDM, WDM e FDM.
• Física Estendida: Incorpora física específica para modelos de matéria escura alternativa, incluindo:
  • WDM: Funções de transferência de supressão dependentes de escala e as correspondentes modificações de HMF/viés (ex: Schneider et al. 2012, Lovell et al. 2014).
  • FDM: Implementação do núcleo de solitão de Schrödinger–Poisson, perfis de densidade compostos de solitão-NFW e relações de concentração–massa específicas (ex: Schive et al. 2016, Laroche et al. 2022).
• Validação: Todas as rotinas principais são validadas contra colossus, hmf e halomod. O artigo demonstra concordância de sub-percentual para modelos compartilhados em espectros de potência, variâncias, funções de correlação e funções de massa de halo e perfis de densidade.
• Flexibilidade: Suporta grades arbitrárias definidas pelo usuário para estatísticas de campo e observáveis de halo, facilitando a integração em pipelines existentes de modelagem de instrumentos e previsões de levantamentos.

Resultos
O artigo valida o phantom através de comparações extensas:

• Espectros de Potência & Variância: Comparações com o colossus para CDM (Eisenstein & Hu) mostram resíduos fracionários dentro de $\pm 7 \times 10^{-5}$ (0,007%) para modelos de zero-bárion e $\sim 0,03\%$ para modelos de bárions completos em escalas de BAO. Os cálculos de variância concordam dentro de 0,2% em uma ampla gama de escalas.
• Funções de Correlação: A concordância com o colossus é mantida ao longo de quatro décadas de amplitude, com desvios limitados a $\sim 1,5\%$ próximo ao cruzamento de zero da função de correlação, consistente com limitações numéricas conhecidas.
• Funções de Massa de Halo & Viés: Cálculos de HMF e viés para CDM concordam com colossus e hmf em nível de sub-percentual. Modelos não-frios (WDM/FDM) mostram a supressão esperada em baixas massas consistente com suas respectivas funções de transferência.
• Perfis de Densidade: Perfis NFW e Hernquist concordam com o colossus em melhor de 0,1%. Perfis Einasto mostem um offset de $\sim 0,5\%$ atribuído a convenções de definição de altura de pico. Perfis compostos de FDM reproduzem com sucesso núcleos de solitão e transição para halos externos NFW, correspondendo a dados de simulação.
• Desempenho: Benchmarks de tempo de execução em um único núcleo mostram que o phantom é comparável ou mais rápido que o colossus para espectros de potência, embora ligeiramente mais lento para cálculos de variância e concentração devido à escolha de integração direta sobre buscas em tabelas. Os autores observam que este trade-off prioriza a transparência numérica.

Significância e Alegações
O artigo alega que o phantom preenche uma lacuna crítica no ecossistema de software para astrônomos que trabalham com MATLAB e Octave. Ao fornecer um ponto de entrada único e validado para análise de halos através de CDM, WDM e FDM, ele elimina a necessidade de portabilidade manual de dados ou implementações customizadas não verificadas. Os autores enfatizam que a ferramenta foi desenhada para integrar-se perfeitamente em pipelines de análise existentes para modelagem de instrumentos, previsões de levantamentos e pós-processamento de simulações.

A significância do trabalho reside na sua capacidade de estender o alcance da modelagem avançada de halos (incluindo teorias alternativas de matéria escura) para uma base de usuários que anteriormente carecia de ferramentas nativas. O artigo enquadra modestamente a ferramenta como um "ponto de entrada validado" que permite aos usuários explorar diferentes observáveis e cenários de matéria escura dentro de um framework consistente, em vez de alegar resolver novos problemas físicos por si só. Extensões futuras mencionadas incluem um módulo completo de modelo de halo para espectros de potência não lineares, modelos para matéria escura auto-interagente (SIDM) e backends numéricos aprimorados.