Cosmología viscosa: cómo la viscosidad volumétrica cambia el nacimiento de agujeros negros primordiales ✛ ✎ ✦

Original: "Primordial black hole formation in bulk-viscous cosmology"
arXiv:2606.26532 · 2026-06-25 · arXiv non-exclusive · análisis propio · 5 min
Tener en cuenta la viscosidad volumétrica del fluido cósmico eleva el umbral de colapso y aumenta la masa de los agujeros negros primordiales resultantes, como se muestra en simulaciones numéricas.

Contexto

Uno de los principales candidatos para la materia oscura son los agujeros negros primordiales, hipotéticamente nacidos del colapso de regiones superdensas después de la expansión inflacionaria. Su predicción se asocia a nombres como Stephen Hawking y otros pioneros. Sin embargo, el análisis estándar supone que el fluido cósmico es ideal, mientras que en las condiciones extremas del Big Bang pueden surgir efectos disipativos, en particular, la viscosidad volumétrica. Este trabajo aplica por primera vez la simulación numérica en la relatividad general completa para averiguar cómo la viscosidad cambia el umbral de colapso y el espectro de masas de los agujeros negros primordiales, algo crucial para verificar la hipótesis de la materia oscura e interpretar los datos de microlente gravitacional.

Métodos

En el marco del formalismo tridimensional BSSN (Baumgarte–Shapiro–Shibata–Nakamura) se resolvió el sistema de ecuaciones de Einstein y las ecuaciones de hidrodinámica con presión viscosa según la teoría de Eckart. El tensor de energía-momento contenía el término Π = −ε $ρ^{1/2}$ $Θ_f$, donde ε es la fuerza adimensional de la viscosidad y $Θ_f$ es la expansión del fluido. El fondo se describió mediante las ecuaciones de Friedmann con una ecuación de estado efectiva $w_{eff}$ = w − ε√(24π). Las sobredensidades iniciales se establecieron con dos perfiles de curvatura: una función sinc (característica de un pico agudo en el espectro) y una gaussiana. La evolución se calculó en una malla radial no uniforme mediante el método explícito de Runge-Kutta de cuarto orden. La formación de un agujero negro se identificó por la aparición de una superficie atrapada, y su masa se determinó mediante la masa de Misner-Sharp. La estimación analítica del desplazamiento del umbral se basó en el argumento de Jeans sobre el equilibrio entre la compresión gravitacional y la presión del sonido.

Resultados

Se encontró que el umbral crítico µ_c aumenta con la viscosidad de forma casi lineal: para ε=0.01 y w=1/3 el incremento es de ~0.04 (de 0.60 a 0.64 para el perfil sinc, de 0.80 a 0.83 para el gaussiano). La dependencia del parámetro de la ecuación de estado mantiene un carácter monótono, similar al caso no viscoso. El ajuste de ley de potencia µ_c − µ_{c,0} = K ε^Γ da Γ ≈ 1 y K del orden de la unidad (Tabla II). Cerca del umbral, la masa del agujero negro primordial sigue el escalado $M_{PBH}$/$M_H$ = k |µ−µ_c|^γ. La viscosidad incrementa la amplitud k: por ejemplo, para el perfil sinc, k cambia de 8.30 a 10.22, mientras que el exponente γ disminuye ligeramente (de 0.37 a 0.34). Este resultado concuerda con el cuadro físico: la disipación ralentiza la compresión, requiriendo un mayor gradiente inicial de densidad para superar el colapso gravitacional, de forma análoga a la acción de la inestabilidad de Jeans en un medio viscoso.

Implicaciones

Los resultados demuestran por primera vez que incluso una pequeña viscosidad volumétrica (ε ~ 0.01) puede cambiar sistemáticamente el umbral de formación de agujeros negros primordiales y su función de masa. Esto significa que las predicciones fiables de la abundancia de agujeros negros primordiales como materia oscura deben tener en cuenta la no idealidad del plasma primordial. La viscosidad desplaza la distribución hacia masas mayores, lo que podría suavizar las tensiones con las restricciones observacionales. Además, el trabajo amplía el arsenal de la cosmología numérica, mostrando la aplicabilidad del formalismo BSSN a medios disipativos, y abre el camino hacia simulaciones más realistas de procesos gravitacionales extremos.

Desarrollo futuro

Futuras investigaciones podrían incluir la viscosidad de segundo orden (teoría de Israel-Stewart) con un tiempo de relajación finito, lo que eliminaría las acausalidades. Las simulaciones 3+1 completas con anisotropía permitirán calcular los espines y elipticidades de los agujeros negros recién nacidos. También resulta interesante el aspecto de la tensión de Hubble: la viscosidad suaviza la ecuación de estado del fondo, lo que podría reducir la discrepancia entre las mediciones locales y globales de la constante de Hubble, creando una conexión no trivial entre la microfísica del colapso y la cosmología a gran escala.

Impacto

El estudio influirá en la interpretación de los datos de las misiones de búsqueda de materia oscura mediante microlente gravitacional (Subaru, OGLE, futuros telescopios espaciales), así como en el análisis del fondo estocástico de ondas gravitacionales, donde la función de masas de los agujeros negros primordiales es importante. Para la física de altas energías, es un ejemplo del papel significativo de la disipación en el colapso gravitacional, lo que también podría ser relevante en contextos astrofísicos.

Próximos pasos

Realizar simulaciones con un rango más amplio de parámetros ε y w, así como con no linealidades en el espectro de potencia. Desarrollar esquemas estables para la evolución de la viscosidad en una calibración arbitraria, aptos para perturbaciones negativas y grandes.

Problemas abiertos clave

El trabajo aborda directamente el problema de la naturaleza de la materia oscura: si los agujeros negros primordiales la explican por completo, la función de masas, determinada por el umbral de colapso, debe satisfacer numerosas restricciones astrofísicas. La modificación viscosa desplaza las predicciones teóricas y podría mejorar la concordancia. Además, la investigación se enmarca en el contexto más amplio del problema de las condiciones iniciales del Universo: los efectos no ideales podrían haber desempeñado un papel en la formación de estructuras, suavizando o agudizando las inhomogeneidades.

La mirada de la máquina del conocimiento

Esta sección está dirigida al autor del trabajo y a sus colegas. Los trabajos vecinos los encontró la búsqueda por significado en nuestro archivo, no por coincidencia de palabras. No evaluamos: mostramos el campo alrededor.

Qué hace el trabajo

En este trabajo se modela numéricamente la formación de agujeros negros primordiales en el universo temprano, teniendo en cuenta la viscosidad volumétrica. Se muestra que la viscosidad eleva el umbral crítico de colapso y modifica la masa de los agujeros negros formados, desplazando las relaciones de escala.

Dónde está su fuerza

Por primera vez se incluye el efecto disipativo de la viscosidad volumétrica en el marco completo de la relatividad general, lo que proporciona correcciones cuantitativas al umbral y al espectro de masas de los agujeros negros primordiales, parámetros clave para su conexión con la materia oscura.

Hacia dónde puede avanzar

Vale la pena probar

  • Verificar si el desplazamiento lineal del umbral se mantiene para perfiles con un espectro más amplio, utilizando el enfoque del trabajo sobre la presión.
  • Evaluar cómo las correcciones viscosas al espectro de masas de los agujeros negros primordiales modificarían las predicciones para la materia oscura en el modelo del trabajo sobre el hervidor cósmico.
  • Aplicar el formalismo con viscosidad al colapso de estrellas oscuras y comparar los umbrales con los obtenidos para perturbaciones primordiales.

Por qué puede importar

La inclusión de la viscosidad puede desplazar notablemente las predicciones para la población de agujeros negros primordiales, lo que es importante para la interpretación de los datos de microlente y las restricciones sobre la materia oscura.

Aunque la viscosidad se asocia a menudo con un efecto retardador, en el contexto del nacimiento de los agujeros negros desempeña un papel estabilizador: la «miel viscosa» del medio cósmico requiere un mayor exceso de densidad para colapsar. No fue hasta 2026 cuando las simulaciones numéricas permitieron cuantificar este efecto, más de medio siglo después de los primeros trabajos de Hawking y Carr.

w_{\text{eff}} = w - \epsilon \sqrt{24\pi}
La viscosidad reduce la presión, haciendo que el fluido sea más blando.
\mu_c - \mu_{c,0} = K \epsilon^\Gamma
Parametrización empírica que mostró un índice Γ ≈ 1 y una amplitud K del orden de la unidad.
M_{\text{PBH}} / M_H = k |\mu - \mu_c|^\gamma
Comportamiento universal donde la viscosidad aumenta el coeficiente k, afectando débilmente a γ.

Cifras clave

  • Elevación del umbral de colapso: Δµ_c ≈ 0.04 para ε=0.01 (w=1/3)
  • Parámetro efectivo de la ecuación de estado: w_eff ≈ 0.25 para w=1/3, ε=0.01
  • Amplitud del escalado de masa (sinc): k aumenta de 8.3 a 10.2 para ε=0.01
  • Exponente del escalado de masa (sinc): γ ≈ 0.37 (ε=0) y 0.34 (ε=0.01)
  • Índice de la ley de potencia para el umbral: Γ ≈ 1
Científicos
Stephen HawkingJacob BekensteinJames JeansFritz ZwickyKarl SchwarzschildSubrahmanyan Chandrasekhar
Conceptos
agujero negro primordial agujeros negros del universo temprano colapso crítico colisión de burbujas Entropía de Bekenstein-Hawking creación de partículas entropía de Barrow inflación de masa nucleación de burbujas cósmicas relación carga-masa formación de estructuras cósmicas ondas gravitacionales primordiales formación de núcleo de materia oscura inhomogeneidades cosmológicas cuánticas área del horizonte cuasiestrella ley de área de Hawking Límite de Chandrasekhar subhalo de materia oscura crecimiento de estructuras