Evaporation of Primordial Black Holes with Multimodal Mass and Extended Spin Distributions: Cosmological Imprints on the Effective Number of Relativistic Species

Resumo Técnico: Evaporação de Buracos Negros Primordiais com Massa Multimodal e Distribuições de Spin Estendidas

Definição do Problema
Buracos Negros Primordiais (PBHs) são hipotetizados para se formarem via vários mecanismos no Universo primordial, cada um imprimindo uma função de massa distinta. Estudos anteriores da evaporação de PBHs e suas consequências cosmológicas basearam-se amplamente na aproximação "monocromática" (assumindo uma única massa característica) ou, no máximo, em distribuições estendidas únicas. Essa abordagem corre o risco de obscurecer as assinaturas características de mecanismos de formação específicos. Além disso, embora o impacto de distribuições de massa estendidas no número efetivo de espécies relativísticas ($N_{\text{eff}}$) tenha sido explorado, o caso de uma população genuinamente multimodal — onde múltiplos canais de formação operam simultaneamente com pesos derivados de probabilidades de colapso físico em vez de parâmetros ajustados — permanece não abordado. Adicionalmente, a interação entre distribuições de massa estendidas, o spin inicial do PBH e o spin da radiação escura emitida ($s_{\text{DR}}$) requer um tratamento sistemático para determinar a robustez das previsões observacionais.

Metodologia
Os autores introduzem o FRISHBEE, uma extensão do código público FRISBHEE, para modelar a evaporação de Hawking de populações de PBHs com distribuições de massa e spin estendidas e multimodais. A metodologia envolve:

1. Distribuições de Massa: O framework implementa quatro funções de massa distintas correspondentes a canais de formação física específicos:

   • Log-Normal (LN): Decorrente de realces localizados no espectro de potência primordial (ex: ultra-slow-roll inflation).
   • Lei de Potência (PL): Resultante de espectros invariantes de escala colapsando em uma era de equação de estado constante.
   • Colapso Crítico (CC): Caracterizado por escalonamento crítico próximo ao limiar de colapso com um corte exponencial.
   • Pré-aquecimento Métrico (MP): Gerado pela amplificação ressonante de perturbações durante as oscilações do inflaton.
   • Multimodal (MM): Uma superposição linear dos quatro canais acima.

2. Ponderação Física: Diferente de abordagens fenomenológicas anteriores, as contribuições fracionárias de densidade de energia ($\alpha_i$) de cada canal na mistura multimodal não são parâmetros livres. Elas são derivadas autoconsistentemente da probabilidade de colapso primordial, $\beta'_i$, que depende da amplitude do espectro de potência de curvatura primordial na escala característica de cada canal. Três cenários de ponderação de referência são testados: Abundância Igual (EA), Ponderada pela Massa Média (MW) e Pumentada pela Densidade Numérica (NW).

3. Distribuições de Spin: A distribuição conjunta massa-spin é fatorizada. Dois perfis de spin são considerados: um perfil Gaussiano (centrado em $a_\star = 0,7$ ou próximo ao quase-extremal $a_\star = 0,99$) e a distribuição universal impulsionada por fusões de Fishbach et al. (com pico próximo a $a_\star \approx 0,7$).

4. Equações de Evolução: O código resolve as equações acopladas de Friedmann-Boltzmann em uma formulação de espaço de fase comovente. Ele rastreia a evolução da massa ($M$) e do spin ($a_\star$) do PBH usando taxas de emissão de Hawking de Kerr (incluindo superradiância para $s_{\text{DR}} \geq 1$) e o conteúdo de partículas do Modelo Padrão. O tratamento considera explicitamente o desacoplamento de neutrinos, focando na janela pré-BBN ($M_{\text{PBH}}^{\text{in}} \in \{10^5, 10^7, 10^8\}$ g) onde a evaporação ocorre bem antes do desacoplamento de neutrinos ($T_{\text{evap}} \gg T_{\text{dec}}$), garantindo $\Delta N_{\text{eff}} > 0$.

Principais Contribuições

• Framework Multimodal: A primeira implementação sistemática de uma população de PBH multimodal onde as frações da população são fisicamente derivadas de parâmetros inflacionários (amplitudes do espectro de potência) em vez de ajustadas.
• Desenvolvimento de Código: O lançamento do FRISHBEE, que incorpora as quatro funções de massa e o módulo de mistura multimodal, validado contra os resultados monocromáticos do FRISBHEE.
• Interação Spin-Massa: Uma quantificação detalhada de como a largura da distribuição de massa modula o aumento induzido pelo spin em $\Delta N_{\text{eff}}$, particularmente na presença de superradiância.

Resultos
O estudo computa $\Delta N_{\text{eff}}$ para várias distribuições, escalas de massa e configurações de spin:

• Aumento pela Distribuição de Massa: Funções de massa estendidas aumentam sistematicamente o $\Delta N_{\text{eff}}$ em relação ao benchmark monocromático na janela pré-BBN. A hierarquia de aumento é preservada através das escalas de massa: Log-Normal ($\sim 1,84\times$) > Multimodal ($\sim 1,77\times$) > Lei de Potência ($\sim 1,69\times$) > Pré-aquecimento Métrico ($\sim 1,15\times$) > Colapso Crítico ($\sim 1,03\times$). Esta hierarquia mantém-se para os três cenários de ponderação (EA, MW, NW), confirmando a robustez do sinal multimodal.
• Efeitos de Spin:
  • Para radiação escura escalar ($s_{\text{DR}}=0$) ou partículas de baixo spin, o aumento do spin inicial ($a_\star$) resulta em aumentos modestos ($R_i \lesssim 1,7$).
  • Para radiação escura de spin-2 ($s_{\text{DR}}=2$) e spins quase-extremais ($a_\star = 0,99$), a superradiância impulsiona significativamente a emissão. No entanto, o aumento é diluído em distribuições de massa largas porque os PBHs mais leves na cauda perdem o spin rapidamente antes de evaporar.
  • Inversão: Para $s_{\text{DR}}=2$ e $a_\star = 0,99$, a razão de aumento $\Delta N_{\text{eff}} / \Delta N_{\text{eff}}^{\text{mono}}$ cai abaixo de unidade para distribuições largas (Log-Normal, Lei de Potência). Neste regime específico, a aproximação monocromática superestima o $\Delta N_{\text{eff}}$ porque assume que todos os PBHs mantêm o spin quase-extremal, enquanto as populações estendidas perdem o spin mais rapidamente em média.
• Significância Observacional:
  • Para radiação escura escalar ($s_{\text{DR}}=0$) em $M_{\text{PBH}}^{\text{in}} = 10^7$ g, a previsão monocromática ($\Delta N_{\text{eff}} \approx 0,057$) está abaixo da sensibilidade projetada para o CMB-S4 ($\sigma \approx 0,06$).
  • Em contraste, as distribuições estendidas Log-Normal, Lei de Potência e Multimodal rendem $\Delta N_{\text{eff}} \approx 0,096\text{--}0,105$, atingindo o nível de detecção de $\sim 1,6\text{--}1,7\sigma$ para o CMB-S4.
  • Para radiação escura de spin-2 ($s_{\text{DR}}=2$), o sinal absoluto é suprimido por uma ordem de magnitude, tornando todas as distribuições indetectáveis nas escalas de massa atuais.

Significância
O artigo estabelece que a forma da função de massa do PBH é um driver primário da impressão cosmológica da evaporação de Hawking, não meramente um detalhe secundário. Ao ligar as frações da população multimodal diretamente à amplitude do espectro de potência primordial, o framework transforma o $\Delta N_{\text{eff}}$ de um simples limite de existência em um potencial discriminador entre cenários de formação de PBH de canal único e de múltiplos canais. Os resultados indicam que experimentos de próxima geração de CMB (CMB-S4, Simons Observatory) poderiam potencialmente distinguir entre populações de PBH monocromáticas e estendidas/multimodais, desde que a radiação escura seja escalar. O trabalho destaca que a observabilidade do sinal de Hawking é contingente tanto no canal de formação quanto no spin da radiação escura emitida, necessitando de um movimento além da aproximação monocromática para previsões fenomenológicas precisas.