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

Resumen Técnico: Diferentes Efectos de las Funciones de Bulto de Lorentz y Gaussianas en Agujeros Negros Primordiales y Ondas Gravitacionales Secundarias

Planteamiento del Problema
Los Agujeros Negros Primordiales (PBH, por sus siglas en inglés) son candidatos viables para la materia oscura fría (DM) y fuentes potenciales de eventos de ondas gravitacionales (GW). Su formación se atribuye a menudo a la amplificación de las perturbaciones escalares durante la inflación, lo cual puede lograrse introduciendo un "bulto" local al potencial de inflación base, $V(\phi) = V_b(\phi)(1 + f(\phi))$. Aunque en la literatura se han utilizado diversas funciones de bulto (Gaussianas, hiperbólicas, Lorentz), no existe un principio físico subyacente que dicte la elección específica de $f(\phi)$. Por consiguiente, no está claro cómo la forma funcional específica del bulto influye en el espectro de potencia de curvatura resultante, la abundancia de PBH y las ondas gravitacionales secundarias inducidas (SIGW). Este artículo aborda la falta de un análisis comparativo entre las funciones de bulto ampliamente utilizadas, específicamente las formas Lorentz y Gaussiana, bajo idénticas condiciones base.

Metodología
Los autores emplean el modelo de inflación de Starobinsky como el potencial base $V_b(\phi)$, elegido por su compatibilidad con las observaciones del Fondo Cósmico de Microondas (CMB). Introducen correcciones de bulto local $f(\phi)$ en dos formas:

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

Para aislar los efectos de la forma funcional, los autores mantienen idénticos el parámetro de amplitud ($b$), la posición del pico ($\phi_0$) y el parámetro de ancho ($c$) para ambos casos. Resuelven numéricamente la ecuación de Mukhanov-Sasaki para computar el espectro de potencia escalar $\mathcal{P}_\zeta$. Utilizando el formalismo de Press-Schechter, calculan la abundancia de PBH ($Y_{PBH}$) y el espectro de densidad de energía de las ondas gravitacionales inducidas por escalares ($\Omega_{GW}$). Se prueban tres conjuntos de parámetros distintos (I, II, III) para verificar la generalidad de los hallazgos.

Contribuciones Clave y Resultados

• Mejora del Espectro de Potencia: El estudio revela una diferencia novedosa y significativa en la capacidad de las dos funciones para mejorar el espectro de potencia. El bulto de Lorentz exhibe una "cola más gruesa" en comparación con el bulto Gaussiano, que decae exponencialmente. Esta diferencia estructural crea una meseta más amplia y plana en el potencial de inflación, lo que conduce a un periodo más extendido de rodadura ultra-lenta (ultra-slow roll). En consecuencia, el caso de Lorentz produce un espectro de potencia pico ($\mathcal{P}_\zeta \sim \mathcal{O}(10^{-2})$) que es órdenes de magnitud superior al caso Gaussiano ($\mathcal{P}_\zeta \sim \mathcal{O}(10^{-7})$) bajo los mismos parámetros.
• Abundancia de Agujeros Negros Primordiales: Debido al espectro de potencia significativamente mayor, el bulto de Lorentz genera una abundancia sustancial de PBH, mientras que el bulto Gaussiano produce una abundancia insignificante (efectivamente cero en los escenarios probados).
  • Conjunto de Parámetros I: Genera PBH de masa estelar ($\sim 37 M_\odot$), potencialmente relevantes para los eventos de fusión de agujeros negros binarios de LIGO.
  • Conjunto de Parámetros II: Produce PBH con masas de $\sim 10^{-5} M_\odot$, que podrían explicar los eventos de microlente de escala de tiempo ultracorta de OGLE.
  • Conjunto de Parámetros III: Produce PBH con masas de $\sim 10^{-12} M_\odot$ con una abundancia $Y_{PBH} \approx 1$, lo que sugiere que podrían constituir toda la materia oscura.
• Ondas Gravitacionales Inducidas por Escalares (SIGWs): Las perturbaciones escalares mejoradas en el caso de Lorentz inducen una densidad de energía de SIGWs significativamente mayor en comparación con el caso Gaussiano. Las frecuencias pico de estas señales caen dentro de los rangos de sensibilidad de varios detectores actuales y futuros:
  • Set I: Detectable por arreglos de temporización de púlsares (EPTA, PPTA, SKA) y NANOGrav.
  • Set II: Detectable por SKA.
  • Set III: Detectable por interferómetros espaciales (LISA, Taiji, TianQin, DECIGO).

Significancia y Reivindicaciones
El artículo afirma que la elección de la función de bulto no es meramente una conveniencia matemática, sino que tiene profundas consecuencias físicas. El decaimiento de ley de potencia de la función de Lorentz permite una mejora más efectiva de las perturbaciones en comparación con el decaimiento exponencial de la función Gaussiana. Esta distinción ofrece un posible control observacional: los experimentos futuros que detecten rangos específicos de masa de PBH o señales de SIGW podrían distinguir entre estos dos escenarios teóricos.

Los autores señalan modestamente que su uso del formalismo de Press-Schechter es una simplificación, y que su análisis asume estadísticas gaussianas, omitiendo posibles efectos de no-gaussianidad en las colas de la función de densidad de probabilidad. Sugieren que trabajos futuros que incorporen la teoría de picos y la no-gaussianidad, así como la investigación de correcciones de un lazo (one-loop) en las escalas del CMB, proporcionarían una comprensión más precisa y completa de los mecanismos de formación de PBH.