A lower bound on primordial power spectrum from halo substructure

Problema
Embora a amplitude da perturbação de curvatura primordial seja fortemente restringida em grandes escalas ($\sim 10^{-9}$) pelas observações do Fundo Cósmico de Micro-ondas (CMB) e da estrutura em grande escala, ela permanece mal restringida em escalas menores ($k \gtrsim \mathcal{O}(1) \text{ Mpc}^{-1}$). As restrições existentes sobre o poder em pequena escala, derivadas de distorções espectrais do CMB, minihalos ultra-compactos e buracos negros primordiais, fornecem primariamente limites superiores. No entanto, a observação de estruturas de pequena escala na Via Láctea, como contagens de satélites, fluxos estelares e anomalias de lente gravitacional, implica a existência de um limite inferior de flutuações primordiais necessário para originar essas estruturas. Este artigo visa derivar um limite inferior para o espectro de potência primordial ($P_\zeta$) em pequenas escalas através da análise da consistência da subestrutura de halos observada com modelos teóricos.

Metodologia
Os autores modelam o espectro de potência da perturbação de curvatura primordial, $P_\zeta$, como uma forma quase invariante de escala com um corte abrupto em $k_{cut} = 5 h \text{ Mpc}^{-1}$ e um aumento (bump) log-normal adicional em escalas menores ($k > k_{cut}$). O modelo é definido por:
$$P_\zeta = P_{CMB}^\zeta \Theta(k_{cut} - k) + P_{bump}^\zeta$$
onde $P_{bump}^\zeta$ é caracterizado pela amplitude $A$, o número de onda de pico $k_b$ e a largura $\Delta$.

Para conectar este espectro primordial à subestrutura observável, os autores empregam os seguintes passos:

1. Evolução do Halo Hospedeiro: Eles computam a evolução temporal de um halo hospedeiro do tipo Via Láctea ($M_0 = 1,3 \times 10^{12} M_\odot$ em $z=0$) usando o formalismo de Press-Schechter Estendido (EPS). Diferente de trabalhos anteriores que dependiam de realizações estocásticas, eles utilizam um método eficiente que calcula o valor esperado da massa do halo hospedeiro em cada etapa de redshift, validado contra realizações aleatórias.
2. Função de Massa de Subhalos: Usando o pacote Python SASHIMI, eles calculam a função de distribuição de massa de subhalos. Isso incorpora o formalismo EPS para a história de acreção e inclui efeitos de remoção por maré (tidal stripping) baseados em formulações analíticas.
3. Correção para o Corte: Para abordar a produção não física de halos em regiões onde o espectro de potência é suprimido (abaixo do corte), eles introduzem uma função de probabilidade de colapso, $P_{coll}(z, M)$, derivada da variância contribuída exclusivamente pelo componente do aumento (bump). Isso modifica a função de massa de subhalos para garantir que halos não se formem onde o poder primordial é efetivamente zero.
4. Comparação Observacional: As funções de massa de subhalos teóricas são comparadas contra dados observacionais de:
   • Fluxos Estelares e Lente Gravitacional: Fornecendo restrições na função de distribuição de massa de subhalos dentro de 300 kpc do centro Galáctico.
   • Contagens de Satélites: Usando dados do Dark Energy Survey e PanSTARRS1 para contar galáxias anãs esferoidais com $V_{max} > 4 \text{ km/s}$.

Principais Contribuições

• Cálculo Eficiente: Os autores demonstram que calcular a evolução da massa do halo hospedeiro via o valor esperado do formalismo EPS produz resultados consistentes com simulações estocásticas, mas reduz significativamente o custo computacional.
• Derivação do Limite Inferior: Eles derivam com sucesso um limite inferior na amplitude do espectro de potência primordial em pequenas escalas ($k \sim 10-30 \text{ Mpc}^{-1}$) ao exigir que o modelo reproduza a abundância observada de subhalos.
• Análise de Remoção por Maré: O estudo quantifica explicitamente o impacto da remoção por maré na função de massa de subhalos, mostrando que ela alarga a distribuição e melhora o ajuste aos dados observacionais, particularmente para aumentos estreitos ($\Delta \sim 0,1$).
• Restrições de Parâmetros: O artigo fornece limites de exclusão de nível de confiança de 95% para o espaço de parâmetros da amplitude do aumento ($A$) e do número de onda de pico ($k_b$) para várias larguras ($\Delta$).

Resultados

• Requisito de Amplitude: Para ser consistente com os dados de fluxos estelares e lentes gravitacionais, a amplitude do aumento log-normal ($P_\zeta$) deve ser de pelo menos $\sim 10^{-8}$ (para $\Delta = 0,1$) em números de onda $k_b \sim 20-40 h \text{ Mpc}^{-1}$. Modelos com amplitudes significativamente menores que esta falham em produzir subhalos suficientes na faixa de massa de $10^8 - 10^{10} M_\odot$.
• Contagens de Satélites: As restrições derivadas de contagens de satélites são menos severas do que as de lentes e fluxos, fornecendo um limite inferior de $P_\zeta \sim 10^{-9}$. No entanto, as contagens de satélites fornecem uma restrição sobre uma gama mais ampla de escalas.
• Dependência da Largura ($\Delta$): Para larguras pequenas ($\Delta = 0,1$), a função de massa é aguda, e a remoção por maré é crucial para ajustar os dados. À medida que $\Delta$ aumenta, a restrição sobre $k_b$ torna-se menos sensível e o espaço de parâmetros permitido se alarga.
• Regiões de Exclusão: A análise exclui regiões onde a amplitude do aumento é muito baixa para gerar a subestrutura observada. Especificamente, para $k_b \gtrsim 80 h \text{ Mpc}^{-1}$, as restrições tornam-se rigorosas porque os subhalos resultantes são pequenos demais para serem detectados como satélites com $V_{max} > 4 \text{ km/s}$.

Significância
O artigo estabelece que a existência de estruturas observadas de pequena escala na Galáxia impõe um limite inferior estrito no espectro de potência primordial em escalas menores do que aquelas sondadas pelo CMB. Este achado é significativo porque complementa as restrições superiores existentes, estreitando efetivamente o espaço de parâmetros permitido para modelos de geração de flutuações primordiais. Os resultados sugerem que mecanismos que suprimem o poder em pequena escala (ex: corrida espectral negativa) devem ser consistentes com estes limites inferiores. Além disso, o estudo valida o uso de fluxos estelares e lentes gravitacionais como sondas poderosas para restringir a física do universo primordial, especificamente o mecanismo de geração de flutuações primordiais.