Different effects of the Lorentz and Gaussian bump functions on the formation of primordial black holes and secondary gravitational waves

Resumo Técnico: Diferentes Efeitos das Funções de Bump de Lorentz e Gaussiana em Buracos Negros Primordiais e Ondas Gravitacionais Secundárias

Enunciado do Problema
Buracos Negros Primordiais (PBHs) são candidatos viáveis para matéria escura fria (DM) e potenciais fontes de eventos de ondas gravitacionais (GW). Sua formação é frequentemente atribuída à amplificação de perturbações escalares durante a inflação, o que pode ser alcançado introduzindo um "bump" (elevação) local ao potencial de inflação base, $V(\phi) = V_b(\phi)(1 + f(\phi))$. Embora várias funções de bump (Gaussiana, hiperbólica, Lorentz) tenham sido utilizadas na literatura, não há um princípio físico subjacente que dite a escolha específica de $f(\phi)$. Consequentemente, permanece incerto como a forma funcional específica do bump influencia o espectro de potência de curvatura resultante, a abundância de PBHs e as ondas gravitacionais secundárias induzidas (SIGWs). Este artigo aborda a falta de uma análise comparativa entre as funções de bump amplamente utilizadas, especificamente as formas Lorentz e Gaussiana, sob condições de base idênticas.

Metodologia
Os autores utilizam o modelo de inflação de Starobinsky como o potencial base $V_b(\phi)$, escolhido por sua compatibilidade com as observações do Fundo Cósmico de Micro-ondas (CMB). Eles introduzem correções de bump locais $f(\phi)$ em duas formas:

1. Bump de Lorentz: $f_L(\phi) = \frac{b}{1 + (\frac{\phi-\phi_0}{c})^2}$
2. Bump Gaussiano: $f_G(\phi) = b e^{-\frac{(\phi-\phi_0)^2}{2c^2}}$

Para isolar os efeitos da forma funcional, os autores mantêm o parâmetro de amplitude ($b$), a posição do pico ($\phi_0$) e o parâmetro de largura ($c$) idênticos para ambos os casos. Eles resolvem numericamente a equação de Mukhanov-Sasaki para computar o espectro de potência escalar $\mathcal{P}_\zeta$. Usando o formalismo de Press-Schechter, eles calculam a abundância de PBHs ($Y_{PBH}$) e o espectro de densidade de energia das ondas gravitacionais induzidas por escalares ($\Omega_{GW}$). Três conjuntos distintos de parâmetros (I, II, III) são testados para verificar a generalidade das descobertas.

Principais Contribuições e Resultados

• Aumento do Espectro de Potência: O estudo revela uma diferença nova e significativa na capacidade das duas funções de aumentar o espectro de potência. O bump de Lorentz exibe uma "cauda mais gorda" em comparação ao bump Gaussiano, que decai exponencialmente. Essa diferença estrutural cria um platô mais largo e plano no potencial de inflação, levando a um período mais estendido de ultra-slow roll (rolagem ultra-lenta). Consequentemente, o caso Lorentz produz um espectro de potência de pico ($\mathcal{P}_\zeta \sim \mathcal{O}(10^{-2})$) que é ordens de magnitude superior ao caso Gaussiano ($\mathcal{P}_\zeta \sim \mathcal{O}(10^{-7})$) sob os mesmos parâmetros.
• Abundância de Buracos Negros Primordiais: Devido ao espectro de potência significativamente maior, o bump de Lorentz gera uma abundância substancial de PBHs, enquanto o bump Gaussiano produz uma abundância negligenciável (efetivamente zero nos cenários testados).
  • Conjunto de Parâmetros I: Gera PBHs de massa estelar ($\sim 37 M_\odot$), potencialmente relevantes para eventos de fusão de sistemas binários de buracos negros detectados pelo LIGO.
  • Conjunto de Parâmetros II: Produz PBHs com massas de $\sim 10^{-5} M_\odot$, que poderiam explicar eventos de microlente de curtíssimo período observados pelo OGLE.
  • Conjunto de Parâmetros III: Resulta em PBHs de massa $\sim 10^{-12} M_\odot$ com uma abundância $Y_{PBH} \approx 1$, sugerindo que eles poderiam constituir toda a matéria escura.
• Ondas Gravitacionais Induzidas por Escalares (SIGWs): As perturbações escalares aumentadas no caso de Lorentz induzem uma densidade de energia de SIGWs significativamente maior em comparação ao caso Gaussiano. As frequências de pico desses sinais caem dentro das faixas de sensibilidade de vários detectores atuais e futuros:
  • Set I: Detectável por arrays de pulsares (EPTA, PPTA, SKA) e NANOGrav.
  • Set II: Detectável pelo SKA.
  • Set III: Detectável por interferômetros espaciais (LISA, Taiji, TianQin, DECIGO).

Significância e Alegações
O artigo afirma que a escolha da função de bump não é meramente uma conveniência matemática, mas possui consequências físicas profundas. O decaimento de lei de potência da função de Lorentz permite um aumento mais eficaz das perturbações em comparação ao decaimento exponencial da função Gaussiana. Essa distinção oferece um potencial controle observacional: experimentos futuros que detectarem intervalos específicos de massa de PBH ou sinais de SIGW poderiam distinguir entre esses dois cenários teóricos.

Os autores observam modestamente que o uso do formalismo de Press-Schechter é uma simplificação e que sua análise assume estatísticas gaussianas, negligenciando potenciais efeitos de não-gaussianidade nas caudas da função de densidade de probabilidade. Eles sugerem que trabalhos futuros incorporando a teoria de picos e não-gaussianidade, bem como investigando correções de um loop nas escalas do CMB, forneceriam uma compreensão mais precisa e abrangente dos mecanismos de formação de PBHs.